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
Do you have a genome file? What's it look like? Bowtie should have a creator script for indexing a genome file.
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.
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.
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.