This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get per sample read depth from VCF with no DP field

Hello Everyone,

I have a VCF file that has no DP info only AD (allelic depth). Does anyone know of a way that I can get per sample mean read depth using only the AD filter?

I could write a custom script that sums the two AD values for each site, but I'd be more comfortable if there was a standard tool available (and it would save a lot of time).

example for two individuals:

GT:AD   0|1:6,5 0|0:12,0

Thank in advance!

vcf coverage depth
• 4,442 views • link
updated by chen chen • written by Rubal Rubal

2 answers

AD = Allele Depth, is not enough to calculate depth

If you don't have DP, you need RD = Reference Depth

• 1 views • link
written by chen chen

Thanks for the reply that's good to know. Why isn't Allele Depth sufficient, I assumed if you add the counts from both alleles you would get the read depth.

• 1 views • link
written by Rubal Rubal

Allele Depth is the depth for ALT, and REF is not included

• 1 views • link
written by chen chen

This belongs under the comment from @Rubal above.

• 1 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this