BCFtools Allelic Depth format nowhere explained?
1
0
Entering edit mode
3.1 years ago
chrisgr ▴ 20

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 • 1.5k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode
3.1 years ago

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

ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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...

ADD REPLY
0
Entering edit mode

Oh I see!

Thank you!

ADD REPLY

Login before adding your answer.

Traffic: 2371 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6