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?