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

Thank in advance!

vcf coverage depth • 2.5k views
ADD COMMENT
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

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

ADD REPLY
0
Entering edit mode
4.8 years ago
chen ★ 2.1k

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

ADD COMMENT
1
Entering edit mode

This belongs under the comment from @Rubal above.

ADD REPLY
0
Entering edit mode

I would like to discuss more about the allele frequency spectrum thread at (https://www.biostars.org/p/101218/). Could you share your email address. It would help a lot. Thanks

ADD REPLY

Login before adding your answer.

Traffic: 1095 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