Question: VCF format: what does AD 1,0 denote?
0
gravatar for jtwalker
8 months ago by
jtwalker20
jtwalker20 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 • 425 views
ADD COMMENTlink modified 8 months ago by Pierre Lindenbaum107k • written 8 months ago by jtwalker20
1
gravatar for Pierre Lindenbaum
8 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum107k 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 https://samtools.github.io/hts-specs/VCFv4.2.pdf

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 8 months ago by Pierre Lindenbaum107k
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: 734 users visited in the last hour