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?
1 answer
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 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.
Log in to answer this question.
-
how to find the Allele Number of specific variant from gnomAD
written by Sara
-
How to find all the region where one read is mapped in a bam file
written by prasoongama •
-
how to find the amplified region from aln.bam file
written by trisha •
-
Extract the mapped contig sequences from SAM/BAM file
written by BDK_compbio
-
Calling allele-specific methylation from bisulfite sequencing reads
written by Jautis
-
For bam and mutation signature I need coverage and number of observations
written by kulvait
-
Why Iontorrent PGM generates high percent of duplicates?
written by murali
-
How to over come "TypeError: argument of type 'int' is not iterable" in python?
written by murali
-
Get the positive and negative strands of reads in bam file for particular nucleotide position?
written by murali
-
Extracting Reads That Are Properly Paired From Bam File
written by Abhi
How are you calling the variants? I mean samtools, GATK, varscan etc.