[2019-04-23 20:55:59] Beginning TopHat run (v2.0.8) ----------------------------------------------- [2019-04-23 20:55:59] Checking for Bowtie Bowtie version: 0.12.9.0 [2019-04-23 20:55:59] Checking for Samtools Samtools version: 0.1.19.0 [2019-04-23 20:55:59] Checking for Bowtie index files [2019-04-23 20:55:59] Checking for Bowtie index files [2019-04-23 20:55:59] Checking for reference FASTA file [2019-04-23 20:55:59] Generating SAM header for /oak/stanford/groups/akundaje/marinovg/genomes/dm3/bowtie-indexes/dm3 format: fastq quality scale: phred33 (default) [2019-04-23 20:56:00] Reading known junctions from GTF file [2019-04-23 20:56:01] Preparing reads left reads: min. length=36, max. length=36, 4077336 kept reads (1373 discarded) [2019-04-23 20:56:36] Using pre-built transcriptome index.. [2019-04-23 20:56:37] Mapping left_kept_reads to transcriptome FlyBase with Bowtie [2019-04-23 20:57:35] 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 [2019-04-23 20:57:35] Mapping left_kept_reads.m2g_um to genome dm3 with Bowtie [2019-04-23 20:58:42] Searching for junctions via segment mapping [2019-04-23 20:59:21] Retrieving sequences for splices [2019-04-23 20:59:27] Indexing splices [2019-04-23 20:59:33] 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-Rep1-Celniker-ENCFF988NQG-TopHat-2.0.8/tmp/left_kept_reads.m2g_um_unmapped.bam [2019-04-23 20:59:37] Joining segment hits [2019-04-23 20:59:52] Reporting output tracks ----------------------------------------------- [2019-04-23 21:00:39] Run complete: 00:04:40 elapsed