SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat 2.0.8b installation error canhu RNA Sequencing 2 05-24-2013 03:09 AM
Tophat 2.0.6 installation error sudheer.svss Bioinformatics 0 02-18-2013 10:50 PM
tophat installation error bob-loblaw Bioinformatics 2 12-17-2012 01:59 PM
Installation error of tophat-2.0.0 fortunezhang General 0 04-18-2012 06:02 AM
TopHat installation vinay052003 Bioinformatics 5 11-22-2011 01:54 PM

Reply
 
Thread Tools
Old 06-20-2014, 04:05 AM   #1
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default Error when testing Tophat installation

I'm just starting out with RNA-seq, and I'm trying to get the installation of Tophat to work. I'm following the "Getting started" instructions on the site, and am currently testing the installation. When I run

Code:
tophat -r 20 test_ref reads_1.fq reads_2.fq
... this is what I get:

Code:
[2014-06-20 13:01:22] Beginning TopHat run (v2.0.11)
-----------------------------------------------
[2014-06-20 13:01:22] Checking for Bowtie
		  Bowtie version:	 2.2.3.0
[2014-06-20 13:01:22] Checking for Samtools
		Samtools version:	 0.1.19.0
[2014-06-20 13:01:22] Checking for Bowtie index files (genome)..
	Found both Bowtie1 and Bowtie2 indexes.
[2014-06-20 13:01:22] Checking for reference FASTA file
[2014-06-20 13:01:22] Generating SAM header for test_ref
[2014-06-20 13:01:22] Preparing reads
	 left reads: min. length=75, max. length=75, 100 kept reads (0 discarded)
	right reads: min. length=75, max. length=75, 100 kept reads (0 discarded)
[2014-06-20 13:01:22] Mapping left_kept_reads to genome test_ref with Bowtie2 
	[FAILED]
Error running:
/Users/erikfasterius/bin/bam2fastx --all ./tophat_out/tmp/left_kept_reads.bam|/Users/erikfasterius/bin/bowtie2 -k 20 -D 15 -R 2 -N 0 -L 20 -i S,1,1.25 --gbar 4 --mp 6,2 --np 1 --rdg 5,3 --rfg 5,3 --score-min C,-14,0 -p 1 --sam-no-hd -x test_ref -|/Users/erikfasterius/bin/fix_map_ordering --bowtie2-min-score 15 --read-mismatches 2 --read-gap-length 2 --read-edit-dist 2 --read-realign-edit-dist 3 --index-outfile ./tophat_out/tmp/left_kept_reads.mapped.bam.index --sam-header ./tophat_out/tmp/test_ref_genome.bwt.samheader.sam - ./tophat_out/tmp/left_kept_reads.mapped.bam ./tophat_out/tmp/left_kept_reads_unmapped.bam
... and I don't really understand the error. It doesn't say what type of error it is, just says [FAILED] and give the long string of whatever caused the error, which I don't understand.

Any help would be greatly appreciated!
Erik
ErikFas is offline   Reply With Quote
Old 06-20-2014, 05:54 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

There should be a "logs" directory in the folder you specified for tophat output. Start looking at the tophat.log and run.log to see if there is some indication in there as to why the run is failing.
GenoMax is offline   Reply With Quote
Old 06-20-2014, 05:59 AM   #3
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Ok, in the tophat.log I can see this as different from the output in Terminal:

Code:
[sam_read1] missing header? Abort!
Full tophat.log:

Code:
[2014-06-20 14:57:28] Beginning TopHat run (v2.0.11)
-----------------------------------------------
[2014-06-20 14:57:28] Checking for Bowtie
		  Bowtie version:	 2.2.3.0
[2014-06-20 14:57:28] Checking for Samtools
		Samtools version:	 0.1.19.0
[2014-06-20 14:57:28] Checking for Bowtie index files (genome)..
[2014-06-20 14:57:28] Checking for reference FASTA file
[2014-06-20 14:57:28] Generating SAM header for test_ref
[2014-06-20 14:57:28] Preparing reads
	 left reads: min. length=75, max. length=75, 100 kept reads (0 discarded)
	right reads: min. length=75, max. length=75, 100 kept reads (0 discarded)
[2014-06-20 14:57:28] Mapping left_kept_reads to genome test_ref with Bowtie2 
[sam_read1] missing header? Abort!
	[FAILED]
Error running:
/Users/erikfasterius/bin/bam2fastx --all ./tophat_out/tmp/left_kept_reads.bam|/Users/erikfasterius/bin/bowtie2 -k 20 -D 15 -R 2 -N 0 -L 20 -i S,1,1.25 --gbar 4 --mp 6,2 --np 1 --rdg 5,3 --rfg 5,3 --score-min C,-14,0 -p 1 --sam-no-hd -x test_ref -|/Users/erikfasterius/bin/fix_map_ordering --bowtie2-min-score 15 --read-mismatches 2 --read-gap-length 2 --read-edit-dist 2 --read-realign-edit-dist 3 --index-outfile ./tophat_out/tmp/left_kept_reads.mapped.bam.index --sam-header ./tophat_out/tmp/test_ref_genome.bwt.samheader.sam - ./tophat_out/tmp/left_kept_reads.mapped.bam ./tophat_out/tmp/left_kept_reads_unmapped.bam
The run.log looks like this:

Code:
/Users/erikfasterius/bin/tophat -r 20 test_ref reads_1.fq reads_2.fq
#>prep_reads:
/Users/erikfasterius/bin/prep_reads --min-anchor 8 --splice-mismatches 0 --min-report-intron 50 --max-report-intron 500000 --min-isoform-fraction 0.15 --output-dir ./tophat_out/ --max-multihits 20 --max-seg-multihits 40 --segment-length 25 --segment-mismatches 2 --min-closure-exon 100 --min-closure-intron 50 --max-closure-intron 5000 --min-coverage-intron 50 --max-coverage-intron 20000 --min-segment-intron 50 --max-segment-intron 500000 --read-mismatches 2 --read-gap-length 2 --read-edit-dist 2 --read-realign-edit-dist 3 --max-insertion-length 3 --max-deletion-length 3 -z gzip --inner-dist-mean 20 --inner-dist-std-dev 20 --no-closure-search --no-coverage-search --no-microexon-search --aux-outfile=./tophat_out/prep_reads.info --index-outfile=./tophat_out/tmp/%side%_kept_reads.bam.index --sam-header=./tophat_out/tmp/test_ref_genome.bwt.samheader.sam --outfile=./tophat_out/tmp/%side%_kept_reads.bam reads_1.fq reads_2.fq
#>map_start:
/Users/erikfasterius/bin/bam2fastx --all ./tophat_out/tmp/left_kept_reads.bam|/Users/erikfasterius/bin/bowtie2 -k 20 -D 15 -R 2 -N 0 -L 20 -i S,1,1.25 --gbar 4 --mp 6,2 --np 1 --rdg 5,3 --rfg 5,3 --score-min C,-14,0 -p 1 --sam-no-hd -x test_ref -|/Users/erikfasterius/bin/fix_map_ordering --bowtie2-min-score 15 --read-mismatches 2 --read-gap-length 2 --read-edit-dist 2 --read-realign-edit-dist 3 --index-outfile ./tophat_out/tmp/left_kept_reads.mapped.bam.index --sam-header ./tophat_out/tmp/test_ref_genome.bwt.samheader.sam - ./tophat_out/tmp/left_kept_reads.mapped.bam ./tophat_out/tmp/left_kept_reads_unmapped.bam
ErikFas is offline   Reply With Quote
Old 06-20-2014, 06:34 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

Does your account have write permissions to the directory where your test data is?

Can you run bam2fastx and see if that is properly installed (i.e. do you see help printed to screen)?
GenoMax is offline   Reply With Quote
Old 06-20-2014, 06:37 AM   #5
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Well, being on an admin account I most certainly hope so! =P

This is what I get when I run bam2fastx, which I assume means it's working correctly:

Code:
Usage: bam2fastx [--fasta|-a|--fastq|-q] [--color] [-Q] [--sam|-s|-t]
   [-M|--mapped-only|-A|--all] [-o <outfile>] [-P|--paired] [-N] <in.bam>
   
Note: By default, reads flagged as not passing quality controls are
   discarded; the -Q option can be used to ignore the QC flag.
   
Use the -N option if the /1 and /2 suffixes should be appended to
   read names according to the SAM flags
   
Use the -O option to ignore the OQ tag, if present, when writing quality values
ErikFas is offline   Reply With Quote
Old 06-20-2014, 06:44 AM   #6
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

That looks good.

Can you download the test data again and try (it is small) with that new copy? I just tried it and did not have any problems with my copy of TopHat 2.0.11.
GenoMax is offline   Reply With Quote
Old 06-20-2014, 06:50 AM   #7
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

New copy of test_data, still the same results, and the logs look the same.
ErikFas is offline   Reply With Quote
Old 06-20-2014, 06:51 AM   #8
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

Did you compile the software from source or are you using the pre-compiled binaries? What flavor of linux is this?
GenoMax is offline   Reply With Quote
Old 06-20-2014, 06:53 AM   #9
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Both Tophat and Bowtie were precompiled, but not Samtools or Boost - I followed the instructions here for them (as well as for the rest of the testing, of course).
ErikFas is offline   Reply With Quote
Old 06-20-2014, 06:57 AM   #10
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

Let us look at some more log files. Can you check prep_reads.log, bowtie.*.log to see if there are any errors there?
GenoMax is offline   Reply With Quote
Old 06-20-2014, 07:00 AM   #11
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

bowtie.left_kept_reads.log:

Code:
100 reads; of these:
  100 (100.00%) were unpaired; of these:
    59 (59.00%) aligned 0 times
    41 (41.00%) aligned exactly 1 time
    0 (0.00%) aligned >1 times
41.00% overall alignment rate
prep_reads.log:

Code:
prep_reads v2.0.11 (4203)
---------------------------
0 out of 100 reads have been filtered out
0 out of 100 read mates have been filtered out
Looks error-free to me.
ErikFas is offline   Reply With Quote
Old 06-20-2014, 07:30 AM   #12
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

True. Not sure what to say next.

Do you also have all these other bowtie related log files?

bowtie_build.log
bowtie.left_kept_reads.log
bowtie.left_kept_reads_seg1.log
bowtie.left_kept_reads_seg2.log
bowtie.left_kept_reads_seg3.log
bowtie.right_kept_reads.log
bowtie.right_kept_reads_seg1.log
bowtie.right_kept_reads_seg2.log
bowtie.right_kept_reads_seg3.log
GenoMax is offline   Reply With Quote
Old 06-20-2014, 07:40 AM   #13
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Nope, only got the four we've already been through... So, something related to that, mayhap?
ErikFas is offline   Reply With Quote
Old 06-20-2014, 07:56 AM   #14
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

Definitely.

If you do not have the bowtie-build.log then that means the indexes are missing. Do you have bowtie2-build in your path?
GenoMax is offline   Reply With Quote
Old 06-20-2014, 07:59 AM   #15
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Yeah, a bunch of them (bowtie2-build-l/s/l-debug/s-debug).
ErikFas is offline   Reply With Quote
Old 06-20-2014, 08:22 AM   #16
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

But no plain bowtie2-build (the -l is for 64-bit indexes but they are not supported by TopHat as yet)?
GenoMax is offline   Reply With Quote
Old 06-20-2014, 08:56 AM   #17
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

No, I have a plain bowtie2-build as well.
ErikFas is offline   Reply With Quote
Old 06-20-2014, 09:05 AM   #18
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,990
Default

I am stumped.

Give a real sample a try to see if that works.
GenoMax is offline   Reply With Quote
Old 06-20-2014, 09:15 AM   #19
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Damn =P I'll do that and see if it works. Thanks for all your help, regardless!
ErikFas is offline   Reply With Quote
Old 06-21-2014, 12:01 AM   #20
ErikFas
Member
 
Location: Sweden

Join Date: Jun 2014
Posts: 86
Default

Downgrading to Bowtie 2.2.2 did the trick! Must have been a bug or something in 2.2.3...
ErikFas is offline   Reply With Quote
Reply

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 08:45 AM.


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