Hi guys,
there are a lot of log files and stuff in the tophat output directory...but is there any way to get easy statistics such as the number of unique mapped reads...number of repeat mapped reads, unmappablem?! And where can i find some infroamtions about all these log files??
Gives "bowtie.left/rightkeptreads.fixmap.log" the information i am looking for??
In the manual there is nothing about this!
Cheers, Steve!
•
written
by
Stevelor
1 answer
The simplest way is to use the picard tools or samtools libraries to inform you about the resulting alignments. In particular, you can use samtools flagstat to give you some mapping numbers. If you want to look at duplicates, you'll need to sort and them mark duplicates. You might also look at the picard tool CollectRnaSeqMetrics for a tool specific to RNA-seq.
Log in to answer this question.
-
valid UMIs in STARsolo output?
written by hamarillo •
-
RNA-seq expression profiling
written by umn_bist •
-
Get # of unique mapped reads from Tophat output
written by colin.kern
-
Variation In The Unique Reads Stats
written by alok.helix
-
Tophat Mapping Rate With Rna-Seq.
written by C Shao •
-
Does The New Tophat2 - Bowtie2 Pipeline Really Map 100% Of The Rna-Seq Reads ? Is This Real ?
written by biorepine
-
Basic Question : #Reads Used By Tophat
written by Abhi
-
Counting Repeat And Unique Reads Of Tophat Output
written by Stevelor
-
Human Transcriptome Mapping
written by Stevelor
-
Tophat Question: No Accepted_Hits.Bam Generated
written by Rnaer •
