Question: How To Plot Rarefaction/Saturation Curve Based On Bam File Generated By Tophat?
1
gravatar for huangqp1984
5.8 years ago by
huangqp198410
huangqp198410 wrote:

I want to determine whether the sequencing depth is enough to cover most of the transcripts in the sample. SO, is there some tools or scripts to plot the saturation curve with the X axis being unique reads(randomly extracted in million from the BAM file generated by Tophat) mapped on genome and the Y axis being transcripts or genes detected with given uniquely aligned reads? ...thanks!

tophat • 3.5k views
ADD COMMENTlink modified 8 months ago by achamess40 • written 5.8 years ago by huangqp198410
0
gravatar for Mikael Huss
5.8 years ago by
Mikael Huss4.6k
Stockholm
Mikael Huss4.6k wrote:

This is probably not exactly what you are looking for, but almost :-) RSeQC has a few functions to plot saturation curves from BAM files given a gene annotation. There is one for RPKM saturation (how the RPKM estimates change as a function of sequencing depth) and one for junction saturation.

ADD COMMENTlink written 5.8 years ago by Mikael Huss4.6k
0
gravatar for achamess
8 months ago by
achamess40
United States
achamess40 wrote:

Are there any readily available tools or scripts to perform this analysis? RSeqQC doesn't do it for me.

I'd like to produce a graph like this: https://cofactorgenomics.com/1-graph-will-give-you-new-perspective-sequencing-experiment/

ADD COMMENTlink modified 8 months ago • written 8 months ago by achamess40
Please log in to add an answer.

Help
Access

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