Hi,
I have processed a WGS sample including alignment (bwa-mem2), variant calling (GATK HaplotypeCaller) and annotation (ANNOVAR). In the annotated file, a variant fitting the phenotype was identified. However, on visualizing the bam in IGV, this variant was not there.
What could be the possible reasons for this?
Will the newly aligned intermediate bam file generated by HaplotypeCaller contains this variant?
Please let me know.
Thank you
1 answer
Depends on what you did while processing your bam file.
What was the the DP,PL,GQ, and read counts for each allele of the variant?
Did you remove duplicates, run base quality score recalibration and local realignment around indels in your bam file? All these can remove good and bad variants from your bam file.
Log in to answer this question.
-
Seeking Advice on Multisample Variant Calling Pipeline with BWA-MEME and GATK DRAGEN Mode
written by George •
-
IGV is not showing the low frequency variants
written by k.manojkumar2991 •
-
Active regions determination in HaplotypeCaller
written by Nawar •
-
multiple processing with single .bam file gonna make problems??
written by hs960201 •
-
gatk HaplotyperCaller bam reference
written by Peter Chung
-
No mutations in BAM (IGV) but a mutation in final VCF?
written by Clark_BioMorgan •
-
GATK; Variant calling RNA-Seq best practices: bam processing steps
written by user230613
-
Gene specific region visualization in Alignment (.bam) or Variant calling file (.vcf) in IGV
written by ravi.uhdnis
-
Gatk - Haplotypecaller Is So Slow, What Is Faster And As Good?
written by newDNASeqer
-
How To Prepare Gatk .Vcf Files Output For Annovar ?
written by newDNASeqer
