Question: Bowtie2 index file problem in RNA-Seq mapping
0
gravatar for rekren
20 months ago by
rekren20
rekren20 wrote:

Hello colleagues,

I have .fastq files obtained by RNA-Seq and I'm trying to map them into customly merged and generated miRNA reference which is taken from miRBase.

I am taking these error.

Could not locate a Bowtie index corresponding to basename "db"
Error: Encountered internal Bowtie 2 exception (#1)
Command: /usr/bin/bowtie2-align-s --wrapper basic-0 -q -x db --passthrough -U /home/graug/Desktop/mirna/2_R.mirna.cleaned.fastq 
(ERR): bowtie2-align exited with value 1

To generate index files, I've typed the following code.

bowtie-build mirna_db.fa db

To map RNA-Seq .fastq files into custom miRNA database of mine, following code was used.

bowtie2 -x db -U -q  /mirna/2_R.mirna.cleaned.fastq -S /mirna/2_R.mirna.sam --no-unal

I've checked the db index withbowtie-inspect command and index seems pretty sanity. Couldn't figure out what am I missing ? Can you help me to find it, please.

index rna-seq bowtie2 • 1.4k views
ADD COMMENTlink modified 20 months ago • written 20 months ago by rekren20
1

Use bowtie2-build instead of bowtie-build.

ADD REPLYlink written 20 months ago by arup1.8k

I wasn't aware that bowtie2 had different index format. Thanks a lot.

On the other hand, I noticed that the order of command has to be like this;

bowtie2 -x db -S /mirna/2_R.mirna.sam --no-unal -U -q  /mirna/2_R.mirna.cleaned.fastq

It gave another error on the previous order of my command which says;

Warning: Output file '/home/graug/Desktop/miRNA/2_R1.cleaned.fastq' was specified without -S.  This will not work in future Bowtie 2 versions.  Please use -S instead.
Warning: Could not open read file "-q" for reading; skipping...
Error: No input read files were valid
(ERR): bowtie2-align exited with value 1
ADD REPLYlink modified 20 months ago • written 20 months ago by rekren20
3
gravatar for genomax
20 months ago by
genomax73k
United States
genomax73k wrote:

Did you build the indexes with bowtie-build (which is v.1) instead of bowtie2-build?

Bowtie 2’s .bt2 index format is different from Bowtie 1’s .ebwt format, and they are not compatible with each other.

You should use bowtie (v.1.x) to align miRNA data though this post claims otherwise: A: Mapping miRNA to the Human Genome with Bowtie2.

In any case make sure you use the correct indexes with correct aligner version.

ADD COMMENTlink modified 20 months ago • written 20 months ago by genomax73k
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: 1182 users visited in the last hour