Dear bioconductor community,

is there an R tool that generates a coverage plot around gene bodies like the one created by the RSeQC tool? Preferably working on region information stored in GenomicRanges objects. I just found a similar question here link which is 6 years ago.

Some input on this is appreciated! Thank you!

I always use a pipeline in which functions of the RSeQC tool are called (incl. the gene coverage plot).

Having said this, in the past I came across an R package called genecovr (link) that based on the available info might do what you are looking for ("The genecovr package provides functionality for evaluating gene body coverage and other metrics in genome assemblies. The central measurement unit is a gene or transcript model"). Also see the vignette.

Also, I know that the R package accompanying the QoRTs tool is used to make visualizations of QC metrics, including gene coverage plots. Again, I don't have any hands-on experience with it, so I cannot comment on its compatibility with your output files...

<<forgot to add>> When browsing Bioconductor packages I also noticed the package CoverageView. Although primarily designed for the visualization of ChIp-seq-like experiments, "additionally, CoverageView can also be used in RNA-seq experiments in order to compare how the coverage is distributed along the transcriptome."

The calculateCoverageProfileList function from the RCAS package can generate summary coverage values across genomic features (5' UTR, cds, 3' UTR, introns, etc.) in an RNA aware manner.