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
problem running in CASAVA for RNA-seq data with multiple indexes clariet Bioinformatics 0 06-06-2011 01:30 PM
problem with tophat commnd directing bowtie indexes repinementer Bioinformatics 2 06-16-2010 08:05 AM
Tophat --segment-length Pepe Bioinformatics 0 05-31-2010 10:33 PM
seed length in TopHat? joseph Bioinformatics 5 09-23-2009 09:23 PM

Thread Tools
Old 09-09-2009, 11:06 AM   #1
Location: Californica

Join Date: Sep 2009
Posts: 19
Default TopHat questions (indexes and seq length)

Hi all,

I just downloaded and installed tophat with bowtie.

Using the test_data, it worked fine and found some junctions.

I went ahead and ran Chris Burge's cerebellum data (12M reads, 36bp) on h_sapiens_37_asm (hg19). It ran about 2.5 hours and to my surprise there is no junctions.

Should I just use h_sapiens (hg18) instead of h_sapiens_37_asm?
What are the difference between assembly and contig indexes?

I know that the tophat chops the seqeunce into small pieces and tries to map the small ones. Is there a way we can specify the length of small pieces? I think, in my run, it did not chop the initial sequence (36bp).

Thank you,
statsteam is offline   Reply With Quote
Old 09-09-2009, 06:03 PM   #2
Cole Trapnell
Senior Member
Location: Boston, MA

Join Date: Nov 2008
Posts: 212

This is unfortunately a known bug in TopHat v1.0.10, caused by the fact that the h_sapiens_37_asm index contains spaces in the index record names, and TopHat parses the Bowtie output and the adjoining fasta files for the index inconsistently. The workaround is to index hg19 yourself, simplifying the fasta record names to chr1, chrX, etc. This will be fixed in an upcoming release. Sorry for the inconvenience.
Cole Trapnell is offline   Reply With Quote

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:30 PM.

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