VarDictJava not reporting EGFR variant, perl version does
See original GitHub issueL861Q_chr7_55259523-55259525.zip
See difference between:
vardict.pl -G /path/to/hg19.fa -f 0.001 -N abc -c 1 -S 2 -E 3 -g 4 -b L861Q.bam -R chr7:55259523-55259525
and
vardict-java -G /path/to/hg19.fa -f 0.001 -N abc -c 1 -S 2 -E 3 -g 4 -b L861Q.bam -R chr7:55259523-55259525
Issue Analytics
- State:
- Created 6 years ago
- Comments:10 (10 by maintainers)
Top Results From Across the Web
AstraZeneca-NGS/VarDictJava: VarDict Java port
The original Perl VarDict is a sensitive variant caller for both single and paired sample variant calling from BAM files. VarDict implements several...
Read more >VarDict: a novel and versatile variant caller for next ... - NCBI
MuTect does not report quality and depth was used in its place. There are total 21 913 of synthetic somatic SNV and InDel...
Read more >VarDictJava
VarDictJava is a variant discovery program written in Java and Perl. It is a partial Java port of VarDict variant caller.
Read more >(PDF) VarDict: A novel and versatile variant caller for next- ...
PDF | Accurate variant calling in next generation sequencing (NGS) is critical to understand cancer genomes better. Here we present VarDict, ...
Read more >(PDF) VarDict: a novel and versatile variant caller for next- ...
Accurate variant calling in next generation sequencing (NGS) is critical to ... MuTect does not report quality and VarDict while the y-axis shows...
Read more >Top Related Medium Post
No results found
Top Related StackOverflow Question
No results found
Troubleshoot Live Code
Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start FreeTop Related Reddit Thread
No results found
Top Related Hackernoon Post
No results found
Top Related Tweet
No results found
Top Related Dev.to Post
No results found
Top Related Hashnode Post
No results found
Top GitHub Comments
Miika, thanks for spotting this. Polina – thanks for this detective work and identifying the underlying issue. Is it possible to have a random assignment of CIGAR trimming to the forward or reverse strands for overlaps so we don’t introduce apparent bias? There would still be some decrease in frequency relative to consensus building, but it would resolve the most apparent change.
(while I’m writing, looks like Miika had a similar idea with R1/R2 based “randomization”)
@mjafin, thanks, I close the issue.