Show RNASeq coverage distribution over transcript length
1
0
Entering edit mode
2.8 years ago

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

RNA-Seq BAM sequencing Transcriptomics • 936 views
ADD COMMENT
1
Entering edit mode
2.8 years ago

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 COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 2588 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