This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How many reads cover the exon per gene from a BAM? The gene (or exon) is totally covered by BAM file?

Hi! I want to calculate a percent of coverage (or the read count witch cover it) of the exons that code for a gene. How can I do? Which is the pipeline? Is it possible to get this information starting from the BAM file?

I'll like to get a result that seems like this:

chr | start | end | reads | exons | gene_name | percent_of_coverage

I would like to use python tools

Thanks a lot

Andrea

coverage gene bam
• 2,671 views • link
updated by Devon Ryan Devon Ryan • written by Andrea Spinelli • Andrea Spinelli

2 answers

Hello Andrea,

if you don't want to reinvent the wheel have a look at picard CollectHsMetrics with the PER_TARGET_COVERAGE parameter. The output looks pretty much as you like to have it.

fin swimmer

• 1 views • link
written by finswimmer finswimmer

Thank you Pierre, but I don't know which is the bed file for matching.

• 1 views • link
written by Andrea Spinelli • Andrea Spinelli

Log in to answer this question.

More posts like this