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

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.

-Chris

Allelic Samtools BCFtools Depth AD • 349 views
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

1
Entering edit mode
11 weeks ago

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

0
Entering edit mode

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

Thank you for the explanation!

-Chris

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

0
Entering edit mode

Oh I see!

Thank you!