Tophat error: Bowtie not found on this system
0
0
Entering edit mode
5.5 years ago
rm16 • 0

I am trying to align RNA-seq reads to my reference but I keep getting the error:

[2016-06-04 13:58:45] Beginning TopHat run (v2.0.9) ----------------------------------------------- [2016-06-04 13:58:45] Checking for Bowtie Bowtie 2 not found, checking for older version.. Error: Bowtie not found on this system.

I have a bowtie2 2.2.5 and tophat 2.1.1 installed on my machine, both are located in the same directory, both are added to my $PATH. I can independently run a bowtie2 script to align reads to the reference. The problem only arises when I try to run a tophat script.

Pseudocode for my tophat script: tophat /path/to/bowtie2-2.2.5/indexes/genomeIndexPrefix /path/to/tophat-2.1.1/RNAseqReads.fastq

Any idea what could be going wrong?

tophat bowtie RNA-Seq alignment • 3.8k views
ADD COMMENT
0
Entering edit mode

What happens if you type bowtie2?

ADD REPLY
0
Entering edit mode

This is what is printed when I type in "bowtie2" with no other arguments (I edited out some of the stuff in the middle):

No index, query, or output file specified! Bowtie 2 version 2.2.5 by Ben Langmead (langmea@cs.jhu.edu, www.cs.jhu.edu/~langmea) Usage:
bowtie2 [options]* -x <bt2-idx> {-1 <m1> -2 <m2> | -U <r>} [-S <sam>]

<bt2-idx> Index filename prefix (minus trailing .X.bt2). NOTE: Bowtie 1 and Bowtie 2 indexes are not compatible. <m1> Files with #1 mates, paired with files in <m2>. Could be gzip'ed (extension: .gz) or bzip2'ed (extension: .bz2). ... ... ... (ERR): bowtie2-align exited with value 1

ADD REPLY
0
Entering edit mode

Not sure if this explains your error, but I noticed you give the full path for the index, while the manual states:

<genome_index_base> The basename of the genome index to be searched. The basename is the name of any of the index files up to but not including the first period. Bowtie first looks in the current directory for the index files, then looks in the indexes subdirectory under the directory where the currently-running bowtie executable is located, then looks in the directory specified in the BOWTIE_INDEXES (or BOWTIE2_INDEXES) environment variable.

So you should only give the basename and make sure that Tophat can find it, e.g. with the environmental variable.

ADD REPLY
0
Entering edit mode

I just gave that a shot, but I am still getting the same error. I also ensured that the path to my bowtie2 index directory is in my PATH variable

ADD REPLY
0
Entering edit mode

It should be in the BOWTIE_INDEXES variable. Not PATH.

ADD REPLY
0
Entering edit mode

Ok, just fixed that. Still getting the error though.

ADD REPLY
0
Entering edit mode

Didn't expect that would solve it, the error message is quite specific. Can you elaborate on your "tophat script"? You don't call it simply from the command line?

ADD REPLY
0
Entering edit mode

Sure. I was unclear: I do call it from command line. Read file is located within tophat directory.

tophat genomeIndexPrefix RNAseqReads.fastq

ADD REPLY
0
Entering edit mode

I just spotted that in your first post you have this:

Beginning TopHat run (v2.0.9)

That is not the version you describe you downloaded. I suspect you have multiple installations... What do you get if you type

which tophat

It will tell you in which location the tophat executable is.

ADD REPLY
0
Entering edit mode

Presumably your script is mucking with the PATH somehow. If you execute tophat directly then do you still get an error?

For what it's worth, specifying a direct path to the basename of the index is the normal practice (ignore the mention of environment variables in the tophat manual).

ADD REPLY
0
Entering edit mode

Good to know, saves me the trouble of specifying that environmental variable ^_^

ADD REPLY

Login before adding your answer.

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