Tophat: Se Reads

Hi All, Can anybody help me how to run tophat for SE data?

I ran the following command but I got an warning and error

tophat --bowtie1 -p 4 -G mouse.gtf -o a_thout mouse/mouse_index/mouse a.fastq

Warning: you have only one segment per read.
If the read length is greater than or equal to 45bp,
we strongly recommend that you decrease --segment-length to about half the read length because TopHat will work better with multiple segments

Mapping left_kept_reads.m2g_um against mouse with Bowtie
[FAILED]
Error: segment-based junction search failed with err =127

Thank you

rna-seq

• 3,161 views • link

updated 13.7 years ago by Pavel Senin

• written 13.7 years ago by RandManP •

1 answer

(1) I would try to run without -p 4, i.e. using a single thread (2) I would configure the segment length so it could infer splice sites/junctions

• 0 views • link

written 13.7 years ago by Pavel Senin

More posts like this

tophat fusion error
written 7.4 years ago by dimitrischat

Hi all. Trying to run tophat fusion with 2 fastq files from rna-seq. I am using the example command from here: https://ccb.jhu.edu/software/tophat/fusion_tutorial.shtml . tophat -o …
Error while mapping raw reads to genome using tophat2
written 9.6 years ago by Bioinfonext

Deal All, I getting error while mapping raw reads to genome using tophat2: Please suggest how to solve it. [2017-01-15 05:09:52] Resuming TopHat pipeline with …
Am I using Tophat2 correctly? Segment-based junction search error keeps coming up.
written 9.9 years ago by aswartz85 •

Although I'm very inexperienced with bioinformatics, what I"m trying to do is very straightforward. I want to align my Miseq mRNAseq reads to the mouse …
Tophat2 "error running tophat_reports"
written 10.3 years ago by Bibi •

I am using Tophat2 (v2.1.1) to find fusions on RNA-seq files. To test the --fusion-search function, I launch the example in the manual : tophat …
tophat failed: Warning: junction database is empty!
written 11.4 years ago by biolab

Hi everyone I embedded tophat command in a perl script. When running the script, I found a tophat failure error. I searched similar problem posts, …
How to run tophat without gff file for bacterial genome
written 11.6 years ago by bassanio •

I would like to find the RNA seq differential expression against bacterial genomes with out gff /gtf files. I used the following commands: ``` tophat …
Tophat2 : Error: Segment-Based Junction Search Failed With Err =-11
written 13.9 years ago by Nicolas Rosewick

Hi, I tried yesterday tophat 2.0.4 with fusion-search on several on my samples and I have an error in a few of them at the …
Tophat : Read Length Warning
written 13.9 years ago by narges •

Hello, I am running tophat over reads with length equal to 36 to 38. And while running it the below warning appears when checking the …
Tophat_Reports "Failed" Error
written 14.1 years ago by vinay kumar •

Hi, Has anyone get tophat2 to work for the fusion search? I am using tophat-2.0.4 and it terminates "Reporting output tracks" step. The error says: …
Tophat Stalled Error = -9
written 14.7 years ago by Craig Anderson •

Heya, I've just tried out tophat on all of my QC'd RNAseq data using the the transcriptome as the index, but unfortunately its become snagged! …