VCF format: what does AD 1,0 denote?

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

• 6,549 views • link

updated 8.9 years ago by Pierre Lindenbaum

• written 8.9 years ago by jtwalker •

1 answer

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)

• 0 views • link

written 8.9 years ago by Pierre Lindenbaum

More posts like this

Pooled-heterozigosity (SECTION2)
written 4.7 years ago by reza

Following the [question][1] I asked before and unfortunately did not receive any help, I have another question that I hope will benefit from your help. …
sed (or similar) command to remove part of the FORMAT from a vcf file?
written 6.4 years ago by Jautis

Hello, I'm looking for a quick way to remove PGT, PID, and PS from the FORMAT of a vcf file output by GATK. Currently, some …
Tools to filter the heterozygous SNPs with biased allele depth (AD) values in vcf file.
written 6.6 years ago by rimgubaev

Hello everyone! There is the following problem in my vcf file produced by GATK: a large proportion of heterozygous sites has a bias in allele …
Tassel5v2 vcf file with 'N' in reference and alternate allele
written 6.8 years ago by Tm

Hello All, Can anyone please help me understand why VCF file obtained using tassel5v2 pipeline "DiscooverySNPcallerV2" is giving 'N' in reference as well as alternate …
Manually calculate Minor Allele Frequency in this VCF
written 7.1 years ago by cookersjs •

I have a VCF with the following lines: ##fileformat=VCFv4.2 ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype"> ##FORMAT=<ID=AD,Number=G,Type=Integer,Description="Allelic Depths of REF and ALT(s) in the order listed"> ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth"> #CHROM POS ID …
Allele Depth (AD) / Allele Balance (AB) Filtering in GATK 4
written 8.2 years ago by maegsul

Hi, I am trying to filter my GATK 4.0.3 - HaplotypeCaller generated multi-sample VCF for allele depth (AD) annotation at sample genotype-level (so available in …
GT (genotype information) in haploid organism has a value of 2 .
written 8.3 years ago by shubhra.bhattacharya

Dear all, Although this question may seem a bit trivial but can someone please explain what the GT field (0,1,2) imply in case of haploid …
Access array in bcftools
written 8.6 years ago by dariober

Hi- I'm filtering a VCF file with `bcftools view` according the a FORMAT field with two items. How can I access the *n-th* item in …
one variant's AD tag in VCF file shows 0,0
written 9.2 years ago by sxzhuxu •

I find a variant in my vcf file. This variant's DPis 232, however the AD is '0,0'. I can't understand why REF and ALT allele …
How to identify Homozygous or Heterozygous variants from the annotation file ?
written 11.2 years ago by deepue

<p>Hi,</p> <p>I am new to the NGS analysis and wondering how one can identify whether a variant is a Homozygous or Heterozygous from the annotated …