![]() |
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
bowtie index problem (bowtie-build and then bowtie-inspect) | tgenahmet | Bioinformatics | 4 | 09-10-2013 12:51 PM |
SpliceMap error: cannot find chromosomes in bowtie index | jgw | RNA Sequencing | 2 | 02-15-2013 07:37 AM |
SpliceMap error: cannot find chromosomes in bowtie index (crossposted) | jgw | Bioinformatics | 0 | 02-03-2012 06:24 AM |
Bowtie Index Files | zun | Bioinformatics | 0 | 10-21-2010 07:55 PM |
problem with tophat commnd directing bowtie indexes | repinementer | Bioinformatics | 2 | 06-16-2010 09:05 AM |
![]() |
|
Thread Tools |
![]() |
#1 | |
ng
Location: LA Join Date: May 2008
Posts: 78
|
![]()
Hello all,
I am using TopHat for the first time and it does not find my bowtie indexes. I have an old version of Bowtie installed; all the indexes that I need are there (I used it before, no problem). When I run TopHat, it gives the following error: Quote:
Any help? Thanks! |
|
![]() |
![]() |
![]() |
#2 |
Member
Location: US Join Date: Oct 2010
Posts: 14
|
![]()
The problem is that you need to specify in the command the name of the bowtie index in addition to the folder.
So, something like this instead: Code:
tophat -r 20 /bowtie-0.12.5/indexes/genome_file_name /data_with_reads/sequence.txt |
![]() |
![]() |
![]() |
#3 | |
ng
Location: LA Join Date: May 2008
Posts: 78
|
![]() Quote:
hg18.1.ebwt hg18.2.ebwt hg18.3.ebwt hg18.4.ebwt hg18.rev.2.ebwt hg18.rev.1.ebwt |
|
![]() |
![]() |
![]() |
#4 |
Member
Location: US Join Date: Oct 2010
Posts: 14
|
![]()
You actually won't specify a file, but rather the base file name. So, in this case your command would be:
Code:
tophat -r 20 /bowtie-0.12.5/indexes/hg18 /data_with_reads/read_1.fastq /data_with_reads/read_2.fastq |
![]() |
![]() |
![]() |
#5 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,480
|
![]()
Just "hg18", without the quotes. From the TopHat manual, you just specify the "basename", which I recall being hg18, or mm9, or whatever you're using. You can then set the PATH in the BOWTIE_INDEXES environment variable. If you're using bash:
Code:
export BOWTIE_INDEXES=/bowtie-0.12.5/indexes Code:
tophat -r 20 hg18 /path/to/reads.fastq |
![]() |
![]() |
![]() |
#6 |
ng
Location: LA Join Date: May 2008
Posts: 78
|
![]()
NM_010117, dpryan,
Thank you, this helped! |
![]() |
![]() |
![]() |
#7 |
Junior Member
Location: Los Angels Join Date: Aug 2012
Posts: 1
|
![]()
Thank you, NM 010117. Your answer to rebrendi helped me.
|
![]() |
![]() |
![]() |
#8 | |
Member
Location: Boston Join Date: Jun 2012
Posts: 33
|
![]() Quote:
Is there a command line option to specify the bowtie index path? Our indexes are stored in a different location from bowtie and have different names from the ususal. Using the methods you listed tophat --fusion... still cannot find its indexes --no matter what path and filename I list. Joe White |
|
![]() |
![]() |
![]() |
#9 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,480
|
![]()
You can just specify the path along with the base, so "/I'm/some/random/path/some_base" will work.
|
![]() |
![]() |
![]() |
#10 |
Senior Member
Location: East Coast USA Join Date: Feb 2008
Posts: 7,080
|
![]()
*All* index files with that specific base name have to be in the path (directory) that Devon mentions above.
|
![]() |
![]() |
![]() |
#11 |
Member
Location: Brazil Join Date: May 2016
Posts: 14
|
![]()
tophat --color --quals fr-secondstrand Bowtie-directory/BOWTIE/ebwt_base -o NGR1_NGR2_1_thout genome NFN_run6_Transcriptoma_29102013_bcSample1_F3_NGR1.csfasta NFN_run6_Transcriptoma_29102013_bcSample1_F3_NGR2.csfasta NFN_run6_Transcriptoma_29102013_bcSample1_F3_QV_NGR1.qual NFN_run6_Transcriptoma_29102013_bcSample1_F3_QV_NGR2.qual
Beginning TopHat run (v2.1.0) Checking for Bowtie Bowtie version: 1.1.1.0 [2016-06-22 13:23:45] Checking for Bowtie index files (genome).. Error: Could not find Bowtie index files (fr-secondstrand.*.ebwt) the bowtie don't find the file index! |
![]() |
![]() |
![]() |
#12 |
Devon Ryan
Location: Freiburg, Germany Join Date: Jul 2011
Posts: 3,480
|
![]()
Look at the actual error message. Sure, it says it can't find them, but look at the names of what you're telling it your indexes are called. Hint: you can't just randomly put "fr-secondstrand" anywhere in the command.
|
![]() |
![]() |
![]() |
Thread Tools | |
|
|