How To Plot Rarefaction/Saturation Curve Based On Bam File Generated By Tophat?
2
1
Entering edit mode
8.3 years ago
huangqp1984 ▴ 10

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 • 4.8k views
ADD COMMENT
0
Entering edit mode
8.3 years ago

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 COMMENT
0
Entering edit mode
3.2 years ago
achamess ▴ 60

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 COMMENT

Login before adding your answer.

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