Hi,

I'd like to calculate RNA-Seq meta read coverage over different subsets of genes similar to the genebody coverage function in RSeQC (http://rseqc.sourceforge.net/#genebody-coverage-py). I saw a number of packages that will do this for genomic coordinates, but I'm looking for something that will take into account just exons when scaling the region. Are there any R packages that will do this?

Thanks