How to get Allele Coverage from bam file or vcf file?
1
0
Entering edit mode
6.8 years ago
murali ▴ 100

I am working with targeted sequencing, I want to find the allele coverage?

What is the allele coverage, and how can I extract it from bam or vcf file?

(If I understood correct, for particular allele how many number of reads mapped) Is that correct?

sequencing allele coverage • 3.1k views
ADD COMMENT
0
Entering edit mode

How are you calling  the variants? I mean samtools, GATK, varscan etc.

 

ADD REPLY
3
Entering edit mode
6.8 years ago
Ram 36k

The term you're looking for is Allele Depth. It can be found in the INFO field for every sample. Google Allele Depth + <target_operation>  (calculate, filter, etc) and you will see tools to get your target operation done.

This (4 year old) thread [ What Is Ad (Allelic Depth) In 1000Genomes Vcf? ]suggests a possible catch to this though - Sum(AD) need not be equal to DP. Not sure if this applies today, and to all variant callers.

Search the forum for Allele Depth - there's plenty of stuff in here.

 

ADD COMMENT

Login before adding your answer.

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