Why different read depth on IGV and in VCF for a variant
1
0
Entering edit mode
2.3 years ago
pragnapcu ▴ 10

For a variant, the read depth in the VCF file (hardfiltered) is 10 whereas, on IGV, the total count is 79. What might be the reason?

depth GATK read haplotypecaller IGV • 1.4k views
ADD COMMENT
1
Entering edit mode
2.3 years ago
LChart 4.6k

Both the variant caller and IGV have default thresholds for mapping quality; could it be that the variant caller was run with a higher mapping quality threshold, and that many of the reads viewed in IGV fall below that threshold? (Also, depending on the caller, the DP field may reflect the filtered or raw depth).

ADD COMMENT
0
Entering edit mode

Yes. It makes sense but, how can I validate the variant in this case (as the depth is low)?

ADD REPLY
2
Entering edit mode

You'll have to re-sequence that locus. Or look at replicate samples that are expected to have the same variants.

ADD REPLY

Login before adding your answer.

Traffic: 2531 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6