How do I setup the pre-built indexed reference for Bowtie2
0
0
Entering edit mode
3.7 years ago
field654 ▴ 30

Hi,

This is really basic but how to I set up the prebuilt reference?

http://bowtie-bio.sourceforge.net/bowtie2/index.shtml

I got the reference from the above link. The reference I tried is named "grch38_1kgmaj_bt2.zip"

After I unzipped I got 6 files, with the same base name and suffixes of .1 .2 .3 .4 .rev.1 .rev.2

That's exactly like the described, built reference. So I copied them to a new folder /BOWTIE2-INDEXES/grch38_1kgmaj_bt2

But when I tried either the align with the reference or inspect, the bowtie2 program won't recognize the folder as a reference.

Like I used the command within the dir to where I copied the 6 files:

bowtie2-inspect grch38_1kgmaj

I simply got the error

Could not locate a Bowtie index corresponding to basename "grch_38_1kgmaj"
Error: Encountered internal Bowtie 2 exception (#1)
Command: /usr/bin/bowtie2-inspect-s --wrapper basic-0 grch_38_1kgmaj

It seems that I have to tell the program somehow that I had the reference here.

Thank you.

Field

Problem solved. It didn't work for me either because I mistyped the file name or I didn't type the base name after the directory

assembly • 2.5k views
ADD COMMENT
0
Entering edit mode

Do you have a genome file? What's it look like? Bowtie should have a creator script for indexing a genome file.

ADD REPLY
0
Entering edit mode

Yes it does. It's bowtie2-build "genome" "export basename for reference".

I'm really new to bioinformatics and bowtie2. I should definitely have tried that.

I started with the pre-built reference though.

Thank you for you kind help.

ADD REPLY
0
Entering edit mode

You will need to add basename for the index after the directory name. If you downloaded pre-built indexes from Bowtie site they should be fine to use as is.

ADD REPLY
0
Entering edit mode

As a matter of fact, when I entered the folder of the index and inspected it worked. I assume the previous time I failed because I mistyped the file name where I missed a "_" or so.

Thank you so much.

ADD REPLY

Login before adding your answer.

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