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!
•
written
by
huangqp1984 •
2 answers
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.
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/
Log in to answer this question.
-
Count The Coverage Rate Of A Bacterial Genome'S Bam File
written by lanyuhao1991 •
-
How To Get The Correct Alignment Statistics After Using Bwa Mem With -M Option
written by skm770 •
-
Variation In The Unique Reads Stats
written by alok.helix
-
How To Find Unique Reads?
written by alok.helix
-
Split A Bam File Into Several Files Containing All The Alignments For X Number Of Reads.
written by gaelgarcia05
-
Shrimp Uniquely Aligned Reads
written by vj
-
New Tophat Adds Multimappers As Entirely New Reads
written by Allpowerde
-
Running Tophat --Fusion-Search On A Filtered Bam That Was Generated By Tophat
written by Malachi Griffith
-
Sample Sam File
written by Steffi
-
Reads Per Contig From Sam File?
written by Newvin •
