Question: Variant callers reporting different read depth on same alignment
0
gravatar for Lauren
18 months ago by
Lauren20
Lauren20 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 18 months ago by Pierre Lindenbaum113k • written 18 months ago by Lauren20
1

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

ADD REPLYlink written 18 months ago by Devon Ryan85k

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

ADD REPLYlink written 18 months ago by WouterDeCoster32k
2
gravatar for Pierre Lindenbaum
18 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum113k 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 18 months ago by Pierre Lindenbaum113k
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: 1786 users visited in the last hour