Question: Variant callers reporting different read depth on same alignment
0
gravatar for Lauren
22 months ago by
Lauren30
Lauren30 wrote:

I generated two VCFs, one with freebayes and one with GATK on the same BAM file.

At one position, freebayes called an MNV. At the same position, GATK called an SNV. What confuses me is that the depths are different in the VCF file. for example,

freebayes called:

chr 9 35906602 C A,T 0/1 6,7,6:19

and GATK called:

chr 9 35906602 C A 0/1 8,7:21

This might be obvious, but I'm new to all this. How can two callers report different total read depths like this? And how would one caller report a higher reference than the other?

ADD COMMENTlink modified 22 months ago by Pierre Lindenbaum117k • written 22 months ago by Lauren30
1

It's quite likely that they have slightly different defaults for filtering.

ADD REPLYlink written 22 months ago by Devon Ryan88k

And how does it look like in e.g. igv?

ADD REPLYlink written 22 months ago by WouterDeCoster37k
2
gravatar for Pierre Lindenbaum
22 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum117k wrote:

What confuses me is that the depths are different in the VCF file. f

different tools use different treshold; for example GATK will not consider

ADD COMMENTlink written 22 months ago by Pierre Lindenbaum117k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 831 users visited in the last hour