Question: bowtie problem: cannot locate bowtie index
0
gravatar for biolab
4.9 years ago by
biolab1.1k
biolab1.1k wrote:

Hi everyone,

I actually have three minor questions about bowtie.  I wrote a run.sh file as follows:

bowtie2-build rrna.fa rRNA
bowtie2 --very-sensitive  -x rRNA  -f reads.fa --un a.fa
bowtie2-build genome.fa GENO
bowtie2 --very-sensitive  -x GENO  -f a.fa  --un b.fa

However, when running bash run.sh, I was encountered a problem as follows.  This is my first question.

Could not locate a Bowtie index corresponding to basename "GENO"
Error: Encountered internal Bowtie 2 exception (#1)

My second question is related to bowtie --very-sensitive option.  Does this option map more reads on the genome?

My third question is about bowtie vs tophat. My understanding is: bowtie can not detect exon splicing, so abandoned  some mapping reads; in contrast, Tophat overcomes this problem and enables more reads  mapping on genome.  Is my understanding right?

THANK YOU very much!

bowtie • 6.5k views
ADD COMMENTlink modified 20 months ago by Biostar ♦♦ 20 • written 4.9 years ago by biolab1.1k
1
gravatar for Devon Ryan
4.9 years ago by
Devon Ryan93k
Freiburg, Germany
Devon Ryan93k wrote:
  1. Presumably something was mistyped in the bowtie2-build command when you tried to build the GENO index. Just look and see what files are in that directory and how big they are. In practice, it's simplest to just run the bowtie2-build command yourself and ensure that it doesn't throw an error message.
  2. Not necessarily. I've had it map fewer with that option, though the difference is quite small. The --very-sensitive and --very-sensitive-local presets mostly affect reads that can align with a high edit distance. However, since the seed size changes with these, you may occasionally exhaust the number of seed extensions before finding the actual best alignment (this is true with any setting, but is more likely to occur as the seed size is decreased). Having said that, you'll typically get a few more reads mapped with these settings, though it'll take longer to do so.
  3. Correct, bowtie2 isn't intended to map mRNA to the genome and will miss splice junctions. If you have that sort of data and sufficient RAM, then use STAR instead. It's MUCH faster than tophat2 and produces equally good results.
     
ADD COMMENTlink written 4.9 years ago by Devon Ryan93k

Hi Devon,

Thank you very much for your detailed explanations, all of which are helpful for my work.

ADD REPLYlink written 4.9 years ago by biolab1.1k

I'm accepting this answer so Biostars bot doesn't bump it again.

ADD REPLYlink written 20 months ago by RamRS25k
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: 772 users visited in the last hour