Question: Show RNASeq coverage distribution over transcript length
gravatar for Federico Giorgi
19 months ago by
Columbia University
Federico Giorgi580 wrote:

Hi all!

I have RNASeq data aligned over a transcriptome in BAM format. Provided the gene annotation (GTF, GFF3), do you know of any tool that can show global transcript coverage with respect to the transcript length?

I would like to quantify any uneven transcriptome-wide transcript coverage, specifically between the 5' and 3' ends of transcripts.

The plot should look like this (I made it up with random data, but you get the idea):

Relative Global Transcript Coverage Pseudo Plot

ADD COMMENTlink modified 19 months ago by Carlo Yague4.8k • written 19 months ago by Federico Giorgi580
gravatar for Carlo Yague
19 months ago by
Carlo Yague4.8k
Carlo Yague4.8k wrote:

You can use deeptools plotProfile along with deeptools computeMatrix and the scale-regions option. This can give you something like this (TSS= 5' end, TES= 3' end). enter image description here

ADD COMMENTlink written 19 months ago by Carlo Yague4.8k

Nice but what's on the y-axis? Is it relative per-transcript read coverage? Why is it negative? It is a nice tool, but what I'm really looking for is a transcriptome-wide assessment of uneven transcript coverage, as a quality assessment of 5' degradation. I guess I will just leave this question open for a few days then I will write a package for it.

ADD REPLYlink written 19 months ago by Federico Giorgi580
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1685 users visited in the last hour