Question: Rookie question-running Bowtie from Terminal
0
gravatar for benkirk86
2.1 years ago by
benkirk860
benkirk860 wrote:

Hi, I am running into trouble on the first line of the Bowtie tutorial. This question has been asked before but I haven't found an answer to it. I am absolutely in the correct directory. Can anybody help with how to fix this? Thanks!

When i enter >>

./bowtie e_coli reads/e_coli_1000.fq

this gets returned:

Traceback (most recent call last): File "./bowtie", line 91, in <module> main() File "./bowtie", line 88, in main os.execv(bin_spec, bowtie_args) OSError: [Errno 2] No such file or directory

bowtie rna-seq osx • 953 views
ADD COMMENTlink written 2.1 years ago by benkirk860

The path (relative) that you are providing does not appear to be correct (e_coli reads/e_coli_1000.fq). You could provide the full path /path_to/e_coli reads/e_coli_1000.fq and see if that helps (change path_to part to actual value on your machine).

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by genomax67k

From webpage tutorial: "The Bowtie source and binary packages come with a pre-built index of the E. coli genome, and a set of 1,000 35-bp reads simulated from that genome. To use Bowtie to align those reads, issue the following command. If you get an error message "command not found", try adding a ./ before the bowtie.

bowtie e_coli reads/e_coli_1000.fq

The first argument to bowtie is the basename of the index for the genome to be searched. The second argument is the name of a FASTQ file containing the reads."

Yeah, my pwd is ~/bowtie-1.2 where the folder is located.

ADD REPLYlink written 2.1 years ago by benkirk860

Do you see a listing when you do ls -lh e_coli reads/e_coli_1000.fq from that directory?

ADD REPLYlink written 2.1 years ago by genomax67k

No it is not listed

ls -lh e_coli reads/e_coli_1000.fq
ls: e_coli: No such file or directory
-rw-r--r--@ 1 benkirk86  staff    78K Dec 19 22:29 reads/e_coli_1000.fq

is what it returns.

I've already changed my bash_profile. I am lost on what to do.

ADD REPLYlink modified 2.1 years ago by genomax67k • written 2.1 years ago by benkirk860

What you are trying to do is run bowtie with pre-made index for e_coli genome. Those index files are in a separate directory called indexes. So your command needs to be

./bowtie indexes/e_coli reads/e_coli_1000.fq

You may want to spend a couple of hours at this site to familiarize yourself with basics of unix. It will make things much clearer.

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by genomax67k

Thanks for the link! And your help. The argument is tripping me up, I will get it figured out

ADD REPLYlink written 2.1 years ago by benkirk860
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: 1707 users visited in the last hour