Seqanswers Leaderboard Ad

**Suimye** · 05-17-2012, 08:42 PM

I also got this message when Im running TopHat v2.0.0 (Bowtie version:2.0.0.5).
Do you have any solution for this error?

Suimye

**jdanderson** · 05-17-2012, 08:51 PM

Sorry, I actually can't recall how I resolved this error, although I wasn't using Tophat v2.0 when I ran into it.

Maybe Cole, or another Tophat programmer could be of more help. It also might help if you posted more info about the run you were attempting to perform and what preceeded the error message.

Cheers,
Johnathon

**Suimye** · 05-17-2012, 09:54 PM

Hi, Johnathon

Thank you for your rapidly replying and suggestion.
I am just trying by using the older version of Tophat and bowtie.

Sincerely,
Suimye

**sameet** · 06-11-2012, 11:31 AM

The error persists!

I am using the latest tophat version, but the error persists. I am also trying with the oldest version available in our system to see if that works.

**sameet** · 06-11-2012, 12:55 PM

[Mon Jun 11 15:35:24 2012] Beginning TopHat run (v1.4.0)
-----------------------------------------------
[Mon Jun 11 15:35:24 2012] Preparing output location wt_Deli_thout/
[Mon Jun 11 15:35:24 2012] Checking for Bowtie index files
[Mon Jun 11 15:35:24 2012] Checking for reference FASTA file
[Mon Jun 11 15:35:24 2012] Checking for Bowtie
Bowtie version: 0.12.5.0
[Mon Jun 11 15:35:24 2012] Checking for Samtools
Samtools Version: 0.1.18
[Mon Jun 11 15:35:24 2012] Generating SAM header for /fdb/bowtie/indexes/genome
format: fastq
quality scale: phred33 (default)
[Mon Jun 11 15:35:24 2012] Preparing reads
left reads: min. length=45, count=9213949
right reads: min. length=42, count=7830358
[Mon Jun 11 15:41:45 2012] Mapping left_kept_reads against genome with Bowtie
[Mon Jun 11 16:06:18 2012] Processing bowtie hits
[Mon Jun 11 16:12:37 2012] Mapping left_kept_reads_seg1 against genome with Bowtie (1/2)
[Mon Jun 11 16:13:57 2012] Mapping left_kept_reads_seg2 against genome with Bowtie (2/2)
[Mon Jun 11 16:16:01 2012] Mapping right_kept_reads against genome with Bowtie
[Mon Jun 11 16:39:09 2012] Processing bowtie hits
[Mon Jun 11 16:44:42 2012] Mapping right_kept_reads_seg1 against genome with Bowtie (1/1)
[Mon Jun 11 16:53:43 2012] Searching for junctions via segment mapping
[FAILED]
Error: segment-based junction search failed with err =-11
Loading left segment hits...

I still get the same error with even older release of tophat.

**dchen** · 06-12-2012, 05:02 AM

I received the same error message, with Tophat 2.0.0 + Bowtie2.0.0.6,
but could fix it by correcting the Bowtie genome index, in my case hg19.
While the downloaded index contained the random chromosomes, the hg19.fa fasta file did not. After re-building the index without the random chromosomes,
the error message disappeared.

Hoping that this was helpful,
Doris

**sameet** · 06-12-2012, 05:20 AM

Thanks Doris, I will try to rebuild the index and see if that helps. I will update the post accordingly.

**hmortens** · 08-23-2012, 08:03 AM

How did you re-build the index without the random chromosomes????

**bruce01** · 01-09-2013, 09:59 AM

By random chromosomes do you mean scaffolds (eg chrG*)? Do you just remove all from ref.fasta and build your index from this? I thought I had them removed.

Does anyone know specifically what is causing the error?

**apredeus** · 08-08-2013, 07:53 AM

I also get the error with the same ID, with a very established tophat 1.4.1.1

These are shorter reads I'm working with (25 bp), so I wonder it that is the problem?

I dug in the code a bit but couldn't find what exactly "-11" refers to. I wish there'd be a more verbose mode to tophat that would give more information on errors.

Anyhow, if anybody knows more about this, I'd love to know it.

PS yeah, and segment_juncs.log has no error at all (and it's the last log written). It looks like the program just stopped.

**apredeus** · 08-08-2013, 08:14 AM

and the command I use is the following:

tophat -p 8 -G $DATA/mm9.mrna.10.31.gtf --prefilter-multihits -o $WDIR/tophat_$TAG --segment-length 12 --max-multihits 15 $BOWTIE/mm9 $FQFILE

**alabadorf** · 02-04-2014, 09:33 AM

I encountered this error too and it was indeed related to a mismatch between the index and fasta file. Instead of rebuilding the index without random chromosomes, I used bowtie2-inspect to create a new fasta file from the index. The error went away.

The exit status code -11 corresponds to a segmentation fault suffered by the segment_juncs program, where I suspect it tries to find a fasta sequence but can't. The source of the error is actually in the SeqAn distribution, which they have bundled with the tophat release.

**1520191** · 04-23-2014, 03:49 AM

Originally posted by alabadorf View Post

I encountered this error too and it was indeed related to a mismatch between the index and fasta file. Instead of rebuilding the index without random chromosomes, I used bowtie2-inspect to create a new fasta file from the index. The error went away.

The exit status code -11 corresponds to a segmentation fault suffered by the segment_juncs program, where I suspect it tries to find a fasta sequence but can't. The source of the error is actually in the SeqAn distribution, which they have bundled with the tophat release.

Hi,I would like to know how you use bowtie2-inspect create a new fasta file from the index?whats your command?
my reference file is :musref
my index file are : musref.1.bt2 musref.2.bt2 musref.3.bt2 musref.4.bt2 musref.rev.1.bt2 musref.rev.2.bt2
what command can help me create a .fa file?
Thank you a lot!

**mastal** · 04-23-2014, 04:01 AM

see the Bowtie 2 Manual,

the section titled "The bowtie2-inspect index inspector"

Bowtie 2: Manual

http://bowtie-bio.sourceforge.net/bowtie2/manual.shtml#the-bowtie2-inspect-index-inspector

here is a quote from the manual:

"bowtie2-inspect extracts information from a Bowtie index about what kind of index it is and what reference sequences were used to build it. When run without any options, the tool will output a FASTA file containing the sequences of the original references (with all non-A/C/G/T characters converted to Ns). It can also be used to extract just the reference sequence names using the -n/--names option or a more verbose summary using the -s/--summary option."

Topics	Statistics	Last Post
Evaluating Genome Sequencing for ECMO Patients in the NICU by seqadmin Started by seqadmin, 12-17-2024, 10:28 AM	0 responses 26 views 0 likes	Last Post by seqadmin 12-17-2024, 10:28 AM
New Genetic Toolkit Refines Studies on Gene Function and Disease by seqadmin Started by seqadmin, 12-13-2024, 08:24 AM	0 responses 42 views 0 likes	Last Post by seqadmin 12-13-2024, 08:24 AM
Study Links Brain Mechanism to Emotional Responses in Animals and Humans by seqadmin Started by seqadmin, 12-12-2024, 07:41 AM	0 responses 28 views 0 likes	Last Post by seqadmin 12-12-2024, 07:41 AM
Study Identifies Ribosomal RNA Fingerprints as Early Cancer Biomarkers by seqadmin Started by seqadmin, 12-11-2024, 07:45 AM	0 responses 42 views 0 likes	Last Post by seqadmin 12-11-2024, 07:45 AM

Seqanswers Leaderboard Ad

Announcement

Tophat: segment-based junction search failed with err =-11

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Comment

Latest Articles

ad_right_rmr

News