SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat Error: Could not find Bowtie index files /bowtie-0.12.5/indexes/. rebrendi Bioinformatics 11 06-22-2016 09:55 AM
bowtie build does not create .rev index files plasticdeath Bioinformatics 2 10-01-2014 11:24 AM
Bowtie can't read index files anecsulea Bioinformatics 9 02-20-2013 09:31 PM
bowtie index files jay2008 Bioinformatics 0 08-20-2011 06:17 PM
Multiple index files in Bowtie? spb2003 Bioinformatics 3 10-13-2009 04:22 PM

Reply
 
Thread Tools
Old 10-21-2010, 06:55 PM   #1
zun
Member
 
Location: Japan

Join Date: Oct 2010
Posts: 26
Unhappy Bowtie Index Files

Hi, everyone

I am a little confused about bowtie index files.


I am treating Arabidopsis and Oryza RNA-seq data.
First, to get index file for Arabidopsis reference genome, I run this command.


bowtie-build -f TAIR9_chr_all.fas TAIR9_chr_all


It output these index files,

TAIR9_chr_all.1.ebwt
TAIR9_chr_all.2.ebwt
TAIR9_chr_all.3.ebwt
TAIR9_chr_all.4.ebwt
TAIR9_chr_all.rev.1.ebwt
TAIR9_chr_all.rev.2.ebwt


Then, I run tophat command and the result was good.

BUT, for Oryza there are still index files which were made by my collegue.
These index files are below,

IRGSPb5.1.ebwt
IRGSPb5.2.ebwt
IRGSPb5.3.ebwt
IRGSPb5.4.ebwt
IRGSPb5.fa.masked.index.amb
IRGSPb5.fa.masked.index.ann
IRGSPb5.fa.masked.index.bwt
IRGSPb5.fa.masked.index.fmv
:
:
IRGSPb5.fa.masked.index.sa
IRGSPb5.fa.masked.index.sai
IRGSPb5.rev.1.ebwt
IRGSPb5.rev.2.ebwt
# I wonder what *.index.* file is...


when running tophat with these index files, I got such warning .

Warning: junction database is empty!


At first, I doubt RNA-seq data which may be low quality.
I calcurated the both of Arabidpsis and Oryza QV , however,
their QV were almost same.


I moved *.1~4.ebwt and *.rev.1~2.ebwt files from orysa index directory
to another directory.
then, I tried Tophat again pointing new index directory.
This time I got junction database from Oryza.


Are *.index.* files not neccessary and should I not include them in the index directory?
So, what are these *.index.*?

It is earlier to ask my colleague about this,but he is now in Brazil for conference ..


I am appreciated your help!

zun
zun is offline   Reply With Quote
Reply

Tags
bowtie, index, junction, tophat rna-seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:29 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO