What is wrong with my bowtie2 index?

Hi all. I am trying to map mouse rRNA reads and so I made a bowtie2 index for the mus musculus rRNA genome. I downloaded all rRNA sequences from GenBank and made a FASTA file.

Here is my index code:

bowtie2-build -f RNA_musmusculus.fasta rRNA

The files for index were made (small index because it is just rRNAs):

  4.1M Jun 11 09:49 rRNA.1.bt2
 7.7K Jun 11 09:49 rRNA.2.bt2
  746 Jun 11 09:49 rRNA.3.bt2
 7.7K Jun 11 09:49 rRNA.4.bt2
 4.1M Jun 11 09:49 rRNA.rev.1.bt2
 7.7K Jun 11 09:49 rRNA.rev.2.bt2

And my mapping code is here:

index=Bowtie2Index/rRNA

bowtie2 -t -x $index --un 6234-S10_norRNA.fq.gz -1 $dir'6234-S10_paired_1.fq.gz' -2 $dir'6234-S10_paired_2.fq.gz' -S 6234-S10_mapped.sam

The directories are correct but I am getting this error message:

Bowtie2 2.3.5.1
(ERR): "Bowtie2Index/rRNA" does not exist or is not a Bowtie 2 index
Exiting now ...

Does anyone have any ideas? Thanks!

bowtie2 rrna

• 2,305 views • link

written 6.2 years ago by ksh4718 •

1 answer

The path to it is most likely messed up. Is the script in a directory that contains a folder Bowtie2Index in which the index is located?

• 0 views • link

written 6.2 years ago by ATpoint

Ah, that helped. Thanks!

• 0 views • link

written 6.2 years ago by ksh4718 •

More posts like this

no input BAM records found in ./tophat_out_SRR5481725/tmp/right_kept_reads.m2g_um.candidates15.bam …
written 7.7 years ago by ALI Hkz •

Hi here is my commands and unfortunately i give this error massage in tophat2! tophat2 -p 16 --keep-fasta-order -o ./tophat_out_${NAME}\ --read-mismatches 4\ --read-gap-length 4\ --read-edit-dist …
Tophat runs through without errors but no read mapped
written 10.1 years ago by jrxu.bioinf •

Hello, I am a new user of tophat. The version in use is v2.1.1 (tophat2). I used default parameters except for --no-coverage-search to save time. …
Low final mapping rate by TopHat in spite of high initial Bowtie mapping
written 11.0 years ago by drtamermansour •

Hi, I am running TopHat2/2.0.14 using this code: ``` tophat \ --transcriptome-index "${transcriptome_index}" \ --num-threads 4 \ --output-dir "${output_dir}" \ --coverage-search \ --microexon-search \ --prefilter-multihits …
Tophat2 output error
written 12.2 years ago by K.Nijbroek •

I am trying to run a basic Tophat2 command but it goes wrong somewhere. Perhaps anyone has experience with this error message? I'm running Tophat2 …
Tophat / Cufflinks Not Using Gene_Id Or Gene_Name From Gtf File
written 13.7 years ago by Owen S.

I am trying to get tophat / cufflinks to work with Ensembl annotations. Although this looks like a long post, my problem is quite …
Why Does Tophat Fail On Me?
written 13.9 years ago by learnwell •

Hi BioStar, I am fairly new to RNA-seq. I am trying to use Tophat for RNA-seq read alignment, and am trying to save the Bowtie …