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