|
gwajnberg@lbbc02:~/tophat-1.3.1.Linux_x86_64$ ./tophat -o /storage/gwajnberg/outputtophat/ /storage/gwajnberg/bowtie_splitted/hg19/hg19 /storage/gwajnberg/SRP000614/SRX002932/SRR013983.fastq |
|
[Tue Jul 5 12:59:34 2011] Beginning TopHat run (v1.3.1) |
|
[Tue Jul 5 12:59:34 2011] Preparing output location /storage/gwajnberg/outputtophat// |
|
[Tue Jul 5 12:59:34 2011] Checking for Bowtie index files |
|
[Tue Jul 5 12:59:34 2011] Checking for reference FASTA file |
|
[Tue Jul 5 12:59:34 2011] Checking for Bowtie |
|
Bowtie version: 0.12.7.0 |
|
[Tue Jul 5 12:59:34 2011] Checking for Samtools |
|
Samtools Version: 0.1.16 |
|
[Tue Jul 5 12:59:34 2011] Generating SAM header for /storage/gwajnberg/bowtie_splitted/hg19/hg19 |
|
[Tue Jul 5 12:59:58 2011] Preparing reads |
|
format: fastq |
|
quality scale: phred33 (default) |
|
Left reads: min. length=25, count=121868 |
|
[Tue Jul 5 13:00:06 2011] Mapping left_kept_reads against hg19 with Bowtie |
|
[Tue Jul 5 13:00:38 2011] Processing bowtie hits |
|
[Tue Jul 5 13:02:24 2011] Mapping left_kept_reads_seg1 against hg19 with Bowtie (1/33) |
|
[Tue Jul 5 13:04:09 2011] Mapping left_kept_reads_seg2 against hg19 with Bowtie (2/33) |
|
[Tue Jul 5 13:06:00 2011] Mapping left_kept_reads_seg3 against hg19 with Bowtie (3/33) |
|
[Tue Jul 5 13:07:48 2011] Mapping left_kept_reads_seg4 against hg19 with Bowtie (4/33) |
|
[Tue Jul 5 13:09:35 2011] Mapping left_kept_reads_seg5 against hg19 with Bowtie (5/33) |
|
[Tue Jul 5 13:11:20 2011] Mapping left_kept_reads_seg6 against hg19 with Bowtie (6/33) |
|
[Tue Jul 5 13:13:00 2011] Mapping left_kept_reads_seg7 against hg19 with Bowtie (7/33) |
|
[Tue Jul 5 13:14:36 2011] Mapping left_kept_reads_seg8 against hg19 with Bowtie (8/33) |
|
[Tue Jul 5 13:16:04 2011] Mapping left_kept_reads_seg9 against hg19 with Bowtie (9/33) |
|
[Tue Jul 5 13:17:21 2011] Mapping left_kept_reads_seg10 against hg19 with Bowtie (10/33) |
|
[Tue Jul 5 13:18:16 2011] Mapping left_kept_reads_seg11 against hg19 with Bowtie (11/33) |
|
[Tue Jul 5 13:18:47 2011] Mapping left_kept_reads_seg12 against hg19 with Bowtie (12/33) |
|
[Tue Jul 5 13:18:56 2011] Mapping left_kept_reads_seg13 against hg19 with Bowtie (13/33) |
|
[Tue Jul 5 13:19:00 2011] Mapping left_kept_reads_seg14 against hg19 with Bowtie (14/33) |
|
[Tue Jul 5 13:19:03 2011] Mapping left_kept_reads_seg15 against hg19 with Bowtie (15/33) |
|
[Tue Jul 5 13:19:06 2011] Mapping left_kept_reads_seg16 against hg19 with Bowtie (16/33) |
|
[Tue Jul 5 13:19:08 2011] Mapping left_kept_reads_seg17 against hg19 with Bowtie (17/33) |
|
[Tue Jul 5 13:19:11 2011] Mapping left_kept_reads_seg18 against hg19 with Bowtie (18/33) |
|
[Tue Jul 5 13:19:14 2011] Mapping left_kept_reads_seg19 against hg19 with Bowtie (19/33) |
|
[Tue Jul 5 13:19:16 2011] Mapping left_kept_reads_seg20 against hg19 with Bowtie (20/33) |
|
[Tue Jul 5 13:19:19 2011] Mapping left_kept_reads_seg21 against hg19 with Bowtie (21/33) |
|
[Tue Jul 5 13:19:21 2011] Mapping left_kept_reads_seg22 against hg19 with Bowtie (22/33) |
|
[Tue Jul 5 13:19:24 2011] Mapping left_kept_reads_seg23 against hg19 with Bowtie (23/33) |
|
[Tue Jul 5 13:19:27 2011] Mapping left_kept_reads_seg24 against hg19 with Bowtie (24/33) |
|
[Tue Jul 5 13:19:29 2011] Mapping left_kept_reads_seg25 against hg19 with Bowtie (25/33) |
|
[Tue Jul 5 13:19:32 2011] Mapping left_kept_reads_seg26 against hg19 with Bowtie (26/33) |
|
[Tue Jul 5 13:19:34 2011] Mapping left_kept_reads_seg27 against hg19 with Bowtie (27/33) |
|
[Tue Jul 5 13:19:37 2011] Mapping left_kept_reads_seg28 against hg19 with Bowtie (28/33) |
|
[Tue Jul 5 13:19:39 2011] Mapping left_kept_reads_seg29 against hg19 with Bowtie (29/33) |
|
[Tue Jul 5 13:19:42 2011] Mapping left_kept_reads_seg30 against hg19 with Bowtie (30/33) |
|
[Tue Jul 5 13:19:45 2011] Mapping left_kept_reads_seg31 against hg19 with Bowtie (31/33) |
|
[Tue Jul 5 13:19:47 2011] Mapping left_kept_reads_seg32 against hg19 with Bowtie (32/33) |
|
[Tue Jul 5 13:19:50 2011] Mapping left_kept_reads_seg33 against hg19 with Bowtie (33/33) |
|
[Tue Jul 5 13:19:52 2011] Searching for junctions via segment mapping |
|
[Tue Jul 5 13:21:31 2011] Retrieving sequences for splices |
|
[Tue Jul 5 13:25:29 2011] Indexing splices |
|
[Tue Jul 5 13:25:30 2011] Mapping left_kept_reads_seg1 against segment_juncs with Bowtie (1/33) |
|
[Tue Jul 5 13:25:34 2011] Mapping left_kept_reads_seg2 against segment_juncs with Bowtie (2/33) |
|
[Tue Jul 5 13:25:39 2011] Mapping left_kept_reads_seg3 against segment_juncs with Bowtie (3/33) |
|
[Tue Jul 5 13:25:43 2011] Mapping left_kept_reads_seg4 against segment_juncs with Bowtie (4/33) |
|
[Tue Jul 5 13:25:48 2011] Mapping left_kept_reads_seg5 against segment_juncs with Bowtie (5/33) |
|
[Tue Jul 5 13:25:52 2011] Mapping left_kept_reads_seg6 against segment_juncs with Bowtie (6/33) |
|
[Tue Jul 5 13:25:57 2011] Mapping left_kept_reads_seg7 against segment_juncs with Bowtie (7/33) |
|
[Tue Jul 5 13:26:01 2011] Mapping left_kept_reads_seg8 against segment_juncs with Bowtie (8/33) |
|
[Tue Jul 5 13:26:05 2011] Mapping left_kept_reads_seg9 against segment_juncs with Bowtie (9/33) |
|
[Tue Jul 5 13:26:08 2011] Mapping left_kept_reads_seg10 against segment_juncs with Bowtie (10/33) |
|
[Tue Jul 5 13:26:12 2011] Mapping left_kept_reads_seg11 against segment_juncs with Bowtie (11/33) |
|
[Tue Jul 5 13:26:13 2011] Mapping left_kept_reads_seg12 against segment_juncs with Bowtie (12/33) |
|
[Tue Jul 5 13:26:14 2011] Mapping left_kept_reads_seg13 against segment_juncs with Bowtie (13/33) |
|
[Tue Jul 5 13:26:15 2011] Mapping left_kept_reads_seg14 against segment_juncs with Bowtie (14/33) |
|
[Tue Jul 5 13:26:16 2011] Mapping left_kept_reads_seg15 against segment_juncs with Bowtie (15/33) |
|
[Tue Jul 5 13:26:16 2011] Mapping left_kept_reads_seg16 against segment_juncs with Bowtie (16/33) |
|
[Tue Jul 5 13:26:17 2011] Mapping left_kept_reads_seg17 against segment_juncs with Bowtie (17/33) |
|
[Tue Jul 5 13:26:17 2011] Mapping left_kept_reads_seg18 against segment_juncs with Bowtie (18/33) |
|
[Tue Jul 5 13:26:18 2011] Mapping left_kept_reads_seg19 against segment_juncs with Bowtie (19/33) |
|
[Tue Jul 5 13:26:19 2011] Mapping left_kept_reads_seg20 against segment_juncs with Bowtie (20/33) |
|
[Tue Jul 5 13:26:19 2011] Mapping left_kept_reads_seg21 against segment_juncs with Bowtie (21/33) |
|
[Tue Jul 5 13:26:20 2011] Mapping left_kept_reads_seg22 against segment_juncs with Bowtie (22/33) |
|
[Tue Jul 5 13:26:21 2011] Mapping left_kept_reads_seg23 against segment_juncs with Bowtie (23/33) |
|
[Tue Jul 5 13:26:21 2011] Mapping left_kept_reads_seg24 against segment_juncs with Bowtie (24/33) |
|
[Tue Jul 5 13:26:22 2011] Mapping left_kept_reads_seg25 against segment_juncs with Bowtie (25/33) |
|
[Tue Jul 5 13:26:22 2011] Mapping left_kept_reads_seg26 against segment_juncs with Bowtie (26/33) |
|
[Tue Jul 5 13:26:23 2011] Mapping left_kept_reads_seg27 against segment_juncs with Bowtie (27/33) |
|
[Tue Jul 5 13:26:24 2011] Mapping left_kept_reads_seg28 against segment_juncs with Bowtie (28/33) |
|
[Tue Jul 5 13:26:24 2011] Mapping left_kept_reads_seg29 against segment_juncs with Bowtie (29/33) |
|
[Tue Jul 5 13:26:25 2011] Mapping left_kept_reads_seg30 against segment_juncs with Bowtie (30/33) |
|
[Tue Jul 5 13:26:25 2011] Mapping left_kept_reads_seg31 against segment_juncs with Bowtie (31/33) |
|
[Tue Jul 5 13:26:26 2011] Mapping left_kept_reads_seg32 against segment_juncs with Bowtie (32/33) |
|
[Tue Jul 5 13:26:27 2011] Mapping left_kept_reads_seg33 against segment_juncs with Bowtie (33/33) |
|
[Tue Jul 5 13:26:27 2011] Joining segment hits |
|
sort: write failed: standard output: Broken pipe |
|
sort: write error |
|
[Tue Jul 5 13:29:01 2011] Reporting output tracks |
|
[FAILED] |
|
Error: [Errno 2] No such file or directory |
|
Found 0 junctions from happy spliced reads |
A quick guess: you forgot to provide the path to your reference fasta file.
Haven't run tophat (yet), but it looks like the "joining segment hits" step failed. Broken Pipe means one command was expecting to eat the output of another command via a Unix pipe, and one of the two commands failed. Then tophat tried to write some output at "Reporting output tracks", but that probably failed because of the earlier problem.
Well, there isn´t any option to specify where is the fasta file....the command is :
somewhere I read that the reference file ( fa or fasta) must be in the same directory of bowtie indexes... and that´s what I did!
I built the indexes from this specified index! Somwhere I also read that if there isn´t a ref file, the programa builds a reference file from the recognized indexes files. I´ll try some othere things later, s