This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BCFtools Allelic Depth format nowhere explained?

Hi, I had a question about Allelic depth (AD) of BCFtools.

It has this format e.g. AD=262,18,0 - What number shows the depth of the REF and what is the ALT? and what is the '0'? I found 1 form where someone said 262 is the REF and 18 is the ALT. Is that correct?

I would also like to know why this format is explained nowhere in the manuals.. It seems to me like something everyone would be interested in.

Thank you in advance,

-Chris

allelic samtools bcftools depth ad
• 2,553 views • link
written by chrisgr • chrisgr

I noticed that it follows the same format as the ref and alt collumn of the vcf so, 262 indecates the REF, the 18 the most occuring alt, and the 0 the depth of the 2nd most occuring alt

• 0 views • link
written by chrisgr • chrisgr

1 answer

I'm sure AD is described in your VCF header.....

AD=262,18,0

is a tri-allelic variant. 262 reads are REF, 18 are ALT n°1, 0 are ALT n°2

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thanks for the reply!

But no, not in the VCF header. It says only this:

FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic depths">

INFO=<ID=AD,Number=R,Type=Integer,Description="Total allelic depths">

Thank you for the explanation!

-Chris

• 0 views • link
written by chrisgr • chrisgr

FORMAT=<ID=AD,Number=R,Type=Integer,Description="Allelic depths"> . So all the information you need is here. See the VCF spec for the meaning of Number=R...

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this