Question: Tophat error: Bowtie not found on this system
0
gravatar for rm16
3.5 years ago by
rm160
rm160 wrote:

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?

rna-seq bowtie alignment tophat • 2.6k views
ADD COMMENTlink written 3.5 years ago by rm160

What happens if you type bowtie2?

ADD REPLYlink written 3.5 years ago by Devon Ryan92k

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 REPLYlink modified 3.5 years ago • written 3.5 years ago by rm160

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 REPLYlink written 3.5 years ago by WouterDeCoster42k

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 REPLYlink written 3.5 years ago by rm160

It should be in the BOWTIE_INDEXES variable. Not PATH.

ADD REPLYlink written 3.5 years ago by WouterDeCoster42k

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

ADD REPLYlink written 3.5 years ago by rm160

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 REPLYlink written 3.5 years ago by WouterDeCoster42k

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

tophat genomeIndexPrefix RNAseqReads.fastq

ADD REPLYlink written 3.5 years ago by rm160

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 REPLYlink written 3.5 years ago by WouterDeCoster42k

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 REPLYlink written 3.5 years ago by Devon Ryan92k

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

ADD REPLYlink written 3.5 years ago by WouterDeCoster42k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 976 users visited in the last hour