Error when testing Tophat installation

[ErikFas]

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).

[GenoMax]

Did you compile the software from source or are you using the pre-compiled binaries? What flavor of linux is this?

[ErikFas]

New copy of test_data, still the same results, and the logs look the same.

[GenoMax]

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.

[ErikFas]

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:

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

[GenoMax]

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)?

[ErikFas]

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

[sam_read1] missing header? Abort!

Full tophat.log:

[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:

/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

[GenoMax]

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.

[ErikFas]

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

tophat -r 20 test_ref reads_1.fq reads_2.fq

... this is what I get:

[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