Visualizing average read coverage across features with varying length
1
1
Entering edit mode
6.4 years ago
rxzlmn ▴ 10

Hi,

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?

RNA-Seq R • 2.1k views
ADD COMMENT
1
Entering edit mode

For "base tools", are we including just writing your own functions? This shouldn't be too terrible to do yourself (it's just a matter of smoothing things so interpolation makes sense). Have a look at the approx() function.

ADD REPLY
0
Entering edit mode

I think that's exactly what I was looking for! Many thanks

ADD REPLY
2
Entering edit mode
6.4 years ago

Have a look at theĀ  deeptools suite in particular computeMatrix.

something like:

computeMatrix scale-regions -R exons.bed -S coverage.bw -out exons.mat.gz --outFileNameData profile.tab -B 100 -A 100

will put in profile.tab the raw numbers you can use in R.

ADD COMMENT
0
Entering edit mode

Thanks, I'll try that!

ADD REPLY

Login before adding your answer.

Traffic: 2412 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6