I am trying to figure out a way to plot read densities over features of varying length. I have a set of ~100 genomic features that are between ~3kb to 7kb in length, and I want to plot per-base coverage or raw read coverage across those elements on the same x scale, i.e. from 0 to 100% of feature length.
The only tool so far I found is ngsplot, but it is kind of limited as it is stand-alone and does not have many options.
Is there a way to do that with existing packages or base tools in R?