How to get per sample read depth from VCF with no DP field
2
0
Entering edit mode
4.8 years ago
Rubal ▴ 340

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


vcf coverage depth • 2.5k views
1
Entering edit mode
4.8 years ago
chen ★ 2.1k

AD = Allele Depth, is not enough to calculate depth

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

0
Entering edit mode

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.

0
Entering edit mode
4.8 years ago
chen ★ 2.1k

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

1
Entering edit mode

This belongs under the comment from @Rubal above.

0
Entering edit mode