<p>I have a bam file of whole exome human data. I now want to check coverage of a list of genes in my data. Basically, I want to get output like, for example:</p>
<pre><code>Gene    Percentage of coding region covered
A1BG               100 %
A1BG-AS1           100%
A1CF                99.94%
NEK8                100%
FANCI              100%
A2ML1             100%
</code></pre>
<p>I want to find out at which depth how much percentage of given genes are covered. for example AT 30X NEK8 gene covered 100%.
Interested in the canonical transcript and if canonical transcript is not present then we can consider other transcript in that case.</p>

CC BY-SA 4.0
2023-04-26T11:10:09.007Z
Please [edit] your question and add more detail. How do you define "coverage"? Is one read enough? Do you want a minimum coverage? What transcripts of the gene? All of them? Only the canonical? What if there is no canonical? How do you define a gene? Are you also interested in UTRs? How about non-exonic regulatory regions which will probably not be included in WES data?
terdon
CC BY-SA 4.0
2023-04-26T11:41:50.130Z
Please clarify your specific problem or provide additional details to highlight exactly what you need. As it's currently written, it's hard to tell exactly what you're asking.
Community
CC BY-SA 4.0
2023-04-28T07:35:01.943Z
I have edited the question and provide the information hope now its clear
Ani
CC BY-SA 4.0
2023-04-28T10:37:14.140Z
Thanks. So, what data do you have? How can we know the locations of these genes? How do we know what transcript to use? Do you have a bed file with the regions of interest? Also, again, how do we define coverage? Do you want the average coverage over the length of the gene? The maximum? The minimum? Do you want a report by exon? By the entire gene length? What is your minimum coverage threshold to consider something as "covered"?
terdon
CC BY-SA 4.0
2023-04-28T12:53:35.377Z
I have a bam file of human whole exome data and also a bed file of all regions of interest.  Yes i want the average coverage of gene at threshold 30x and report by gene.
Ani
<gene><coverage><exome>
How can I calculate the coverage of a gene in exome sequencing / calculate coverage per gene?
# How can I calculate the coverage of a gene in exome sequencing / calculate coverage per gene?

(Asked by: username_0 on 2023-04-26)

I have a bam file of whole exome human data. I now want to check coverage of a list of genes in my data. Basically, I want to get output like, for example:
Gene    Percentage of coding region covered
A1BG               100 %
A1BG-AS1           100%
A1CF                99.94%
NEK8                100%
FANCI              100%
A2ML1             100%

I want to find out at which depth how much percentage of given genes are covered. for example AT 30X NEK8 gene covered 100%.
Interested in the canonical transcript and if canonical transcript is not present then we can consider other transcript in that case.

**Tags:** gene, coverage, exome
> **Comments:**
> **username_1**: Please [edit] your question and add more detail. How do you define "coverage"? Is one read enough? Do you want a minimum coverage? What transcripts of the gene? All of them? Only the canonical? What if there is no canonical? How do you define a gene? Are you also interested in UTRs? How about non-exonic regulatory regions which will probably not be included in WES data?  
> **username_2**: Please clarify your specific problem or provide additional details to highlight exactly what you need. As it's currently written, it's hard to tell exactly what you're asking.  
> **username_0**: I have edited the question and provide the information hope now its clear  
> **username_1**: Thanks. So, what data do you have? How can we know the locations of these genes? How do we know what transcript to use? Do you have a bed file with the regions of interest? Also, again, how do we define coverage? Do you want the average coverage over the length of the gene? The maximum? The minimum? Do you want a report by exon? By the entire gene length? What is your minimum coverage threshold to consider something as "covered"?  
> **username_0**: I have a bam file of human whole exome data and also a bed file of all regions of interest.  Yes i want the average coverage of gene at threshold 30x and report by gene.  

---

## 0 answers

