Question: VCF format: what does AD 1,0 denote?
gravatar for jtwalker
9 weeks ago by
jtwalker10 wrote:

I'm looking at one of my VCF files and I noticed that allele depth (AD) isn't just one number, but rather two values seperated by a comma

e.g.: GT:AD:DP:GQ:PL 0/0:1,0:1:3:0,3,27

The GATK doc doesn't really seem all that helpful here, as it doesn't explicitly describe what the two numbers are used for.

Can anyone explain? Thanks!

vcf • 163 views
ADD COMMENTlink modified 9 weeks ago by Pierre Lindenbaum100k • written 9 weeks ago by jtwalker10
gravatar for Pierre Lindenbaum
9 weeks ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum100k wrote:

the definition of the attribute is defined in the header

 grep '=AD,'  test_vcf01.vcf
##FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic depths for the ref and alt alleles in the order listed">

see also the VCF spec

GT:AD:DP:GQ:PL 0/0:1,0:1:3:0,3,27

so we have a number for each allele (REF and all ALT alleles). Here one short read carries the REF allele and none carry the ALT. DP (depth) is 1 (1+0==1)

ADD COMMENTlink written 9 weeks ago by Pierre Lindenbaum100k
Please log in to add an answer.


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