[2018-05-07 05:38:45] Beginning TopHat run (v2.0.8) ----------------------------------------------- [2018-05-07 05:38:45] Checking for Bowtie Bowtie version: 0.12.9.0 [2018-05-07 05:38:45] Checking for Samtools Samtools version: 0.1.19.0 [2018-05-07 05:38:45] Checking for Bowtie index files [2018-05-07 05:38:45] Checking for Bowtie index files [2018-05-07 05:38:45] Checking for reference FASTA file [2018-05-07 05:38:45] Generating SAM header for /oak/stanford/groups/akundaje/marinovg/genomes/dm6/bowtie-indexes/dm6 format: fastq quality scale: phred33 (default) [2018-05-07 05:38:46] Reading known junctions from GTF file [2018-05-07 05:38:49] Preparing reads left reads: min. length=36, max. length=36, 8003194 kept reads (1747 discarded) [2018-05-07 05:39:51] Using pre-built transcriptome index.. [2018-05-07 05:39:54] Mapping left_kept_reads to transcriptome Drosophila_melanogaster.BDGP6.90 with Bowtie [2018-05-07 05:41:57] Resuming TopHat pipeline with unmapped reads Warning: you have only one segment per read. If the read length is greater than or equal to 45bp, we strongly recommend that you decrease --segment-length to about half the read length because TopHat will work better with multiple segments [2018-05-07 05:41:57] Mapping left_kept_reads.m2g_um to genome dm6 with Bowtie [2018-05-07 05:44:07] Searching for junctions via segment mapping [2018-05-07 05:45:10] Retrieving sequences for splices [2018-05-07 05:45:16] Indexing splices [2018-05-07 05:45:25] Mapping left_kept_reads.m2g_um_unmapped to genome segment_juncs with Bowtie (1/1) open: No such file or directory Error: bam2fastx failed to open BAM file S2-CAGE-Rep2-Celniker-ENCFF173DNO-TopHat-2.0.8/tmp/left_kept_reads.m2g_um_unmapped.bam [2018-05-07 05:45:29] Joining segment hits [2018-05-07 05:45:53] Reporting output tracks ----------------------------------------------- [2018-05-07 05:47:13] Run complete: 00:08:28 elapsed