tophat Error running running 'prep_reads'

victoryhe · 10-12-2011, 10:18 AM

Hi, Naive here, i'm getting an error from tophat with mac os x 10.6

-----
tophat -o /tophatoutput/ --solexa1.3-quals /bowtie-0-1.12.7/indexes/mm9 a.fastq

-----

[Wed Oct 12 13:31:31 2011] Beginning TopHat run (v1.3.2)
-----------------------------------------------
[Wed Oct 12 13:31:31 2011] Preparing output location /tophatoutput/
[Wed Oct 12 13:31:31 2011] Checking for Bowtie index files
[Wed Oct 12 13:31:31 2011] Checking for reference FASTA file
Warning: Could not find FASTA file /bowtie-0-1.12.7/indexes/mm9.fa
[Wed Oct 12 13:31:31 2011] Reconstituting reference FASTA file from Bowtie index
Executing: /bowtie-0-1.12.7/indexes/mm9 > /tmp/mm9.fa
[Wed Oct 12 13:34:42 2011] Checking for Bowtie
Bowtie version: 0.12.7.0
[Wed Oct 12 13:34:42 2011] Checking for Samtools
Samtools Version: 0.1.18
[Wed Oct 12 13:34:42 2011] Generating SAM header for /bowtie-0-1.12.7/indexes/mm9
[Wed Oct 12 13:35:36 2011] Preparing reads
format: fastq
quality scale: phred64 (reads generated with GA pipeline version >= 1.3)
[FAILED]
Error running 'prep_reads'
terminate called after throwing an instance of 'int'

-------

i have never seen such an error before, something wrong on command line?

cedance · 10-13-2011, 05:09 PM

quality scale: phred64 (reads generated with GA pipeline version >= 1.3)
FAILED

Here lies the problem. Are you sure --solexa1.3-quals is right? Could you paste here a complete fastq read (header + read + quality)?

victoryhe · 10-17-2011, 04:53 AM

thanks! without the solexa option, it works.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
prep_reads error when running Tophat	bzhang	Bioinformatics	35	10-02-2015 01:49 AM
Tophat running error 1?	day_sleeperz	RNA Sequencing	2	10-04-2011 07:19 AM
Error when running TopHat	clk	RNA Sequencing	3	09-29-2011 12:10 PM
Tophat Running Error	Wei-HD	Bioinformatics	12	02-11-2010 02:42 PM
TopHat running error	pfranchini	Bioinformatics	2	08-10-2009 06:46 AM

10-12-2011, 10:18 AM	#1
victoryhe Junior Member Location: boston Join Date: Oct 2011 Posts: 7	tophat Error running running 'prep_reads' Hi, Naive here, i'm getting an error from tophat with mac os x 10.6 ----- tophat -o /tophatoutput/ --solexa1.3-quals /bowtie-0-1.12.7/indexes/mm9 a.fastq ----- [Wed Oct 12 13:31:31 2011] Beginning TopHat run (v1.3.2) ----------------------------------------------- [Wed Oct 12 13:31:31 2011] Preparing output location /tophatoutput/ [Wed Oct 12 13:31:31 2011] Checking for Bowtie index files [Wed Oct 12 13:31:31 2011] Checking for reference FASTA file Warning: Could not find FASTA file /bowtie-0-1.12.7/indexes/mm9.fa [Wed Oct 12 13:31:31 2011] Reconstituting reference FASTA file from Bowtie index Executing: /bowtie-0-1.12.7/indexes/mm9 > /tmp/mm9.fa [Wed Oct 12 13:34:42 2011] Checking for Bowtie Bowtie version: 0.12.7.0 [Wed Oct 12 13:34:42 2011] Checking for Samtools Samtools Version: 0.1.18 [Wed Oct 12 13:34:42 2011] Generating SAM header for /bowtie-0-1.12.7/indexes/mm9 [Wed Oct 12 13:35:36 2011] Preparing reads format: fastq quality scale: phred64 (reads generated with GA pipeline version >= 1.3) [FAILED] Error running 'prep_reads' terminate called after throwing an instance of 'int' ------- i have never seen such an error before, something wrong on command line? Last edited by victoryhe; 10-12-2011 at 10:22 AM.

10-13-2011, 05:09 PM	#2
cedance Senior Member Location: Germany Join Date: Feb 2011 Posts: 108	quality scale: phred64 (reads generated with GA pipeline version >= 1.3) FAILED Here lies the problem. Are you sure --solexa1.3-quals is right? Could you paste here a complete fastq read (header + read + quality)?

10-17-2011, 04:53 AM	#3
victoryhe Junior Member Location: boston Join Date: Oct 2011 Posts: 7	thanks! without the solexa option, it works.