Tophat Error: Bam2Fastx Failed To Open Bam

I ran Tophat and near the end of the run it gave the following error:

[2012-09-28 15:51:58] Mapping left_kept_reads.m2g_um to genome mm9 with Bowtie 
[2012-09-28 16:08:04] Searching for junctions via segment mapping
[2012-09-28 16:12:27] Retrieving sequences for splices
[2012-09-28 16:14:34] Indexing splices
[2012-09-28 16:15:02] 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 /mydata/left_kept_reads.m2g_um_unmapped.bam
[2012-09-28 16:15:08] Processing bowtie hits
[2012-09-28 16:26:43] Reporting output tracks

It then warned that it couldn't remove all the temporary files. Any idea what might be happening? I have no idea what the "open: No such file or directory" error means. It must not be finding some files but I cannot tell which.

tophat rna-seq bowtie mapping

• 5,991 views • link

updated 10.7 years ago by Biostar

• written 13.9 years ago by user

2 answers

Really sorry for hijacking this thread, but I am facing the exact same error. Tried using relative or fullpath for the output destination but without any change. This error is persistant. The only relevant log I think is this:

bowtie.left_kept_reads.m2g_um_unmapped.log

Warning: Could not find any reads in "-"

reads processed: 0

reads with at least one reported alignment: 0 (0.00%)

reads that failed to align: 0 (0.00%)

No alignments

Any ideas?

Thank you very much in advance

• 224 views • link

written 13.2 years ago by IV

You have the answer in your post already. The very next line after "open: No such file or directory", tells you exactly which file is missing

Error: bam2fastx failed to open BAM file /mydata/leftkeptreads.m2gumunmapped.bam

Your file

leftkeptreads.m2gumunmapped.bam

is missing.

• 1 views • link

written 13.9 years ago by VS

OK, but since this file is generated by Tophat and not by me, how can I detect why it is missing and how to fix it?

• 0 views • link

written 13.9 years ago by user

Tophat is trying to find this missing file in /mydata folder. Can you give the fullpath of the folder from where you are running Tophat? Are you on Unix or windows machine?

• 0 views • link

written 13.9 years ago by VS

I am facing something similar. I will add it as a comment below.

• 0 views • link

written 13.3 years ago by IV

More posts like this

Speed up Tophat2 running
written 8.7 years ago by yuabrahamliu •

Hi, all. I'm starting to use tophat2 to map dUTP strand-specific, paired-end RNA-seq data to mm9 genome. Each sample pair includes 2 samples about 500M …
tophat2 "Broken pipe error"
written 8.9 years ago by ahmad mousavi

Dear all, I have downloaded the most recent versions of Samtools, Tophat2 and Bowtie2 on a Ubuntu 16.04 server. I have got below error, How …
Tophat runs through without errors but no read mapped
written 10.1 years ago by jrxu.bioinf •

Hello, I am a new user of tophat. The version in use is v2.1.1 (tophat2). I used default parameters except for --no-coverage-search to save time. …
tophat2 Error: segment-based junction search failed with err =-9
written 10.2 years ago by marongiu.luigi

Dear all, I am aligning HISeq RNA-Seq data to the human reference genome. I downloaded the Homo_sapiens.GRCh38.dna.toplevel.fa.gz and the Homo_sapiens.GRCh38.84.gtf.gz reference files and I built …
Low final mapping rate by TopHat in spite of high initial Bowtie mapping
written 11.0 years ago by drtamermansour •

Hi, I am running TopHat2/2.0.14 using this code: ``` tophat \ --transcriptome-index "${transcriptome_index}" \ --num-threads 4 \ --output-dir "${output_dir}" \ --coverage-search \ --microexon-search \ --prefilter-multihits …
samtools: writing to standard output failed: Broken pipe
written 11.9 years ago by mcff23 •

Hi everyone! I'm getting this error from tophat: samtools: writing to standard output failed: Broken pipe Log file: ``` [2014-09-15 16:31:23] Beginning TopHat run (v2.0.12) …
Tophat-Fusion Error 'Long_Spanning_Reads'
written 13.5 years ago by samsara

I have run <a href="http://tophat.cbcb.umd.edu/fusion_index.html">Tophat fusion</a> with paired-end RNA-Seq data, but my run was unsuccessful. The error says <code>running 'long_spanning_reads'</code>. I found a post in …
Tophat Error When "Reporting Output Tracks"
written 13.7 years ago by Pierre Lindenbaum

I'm learning tophat . here is the Makefile I'm using: <pre><code>include tools.mk REF=./chr22 REF.bowtie=$(foreach S,1.bt2 2.bt2 3.bt2 4.bt2 rev.1.bt2 rev.2.bt2, $(addsuffix .$S,$(REF)) ) GTF=Homo_sapiens.GRCh37.69.gtf export …
Tophat2 : Error: Segment-Based Junction Search Failed With Err =-11
written 13.9 years ago by Nicolas Rosewick

Hi, I tried yesterday tophat 2.0.4 with fusion-search on several on my samples and I have an error in a few of them at the …
Tophat Error : Oserror: [Errno 2] No Such File Or Directory
written 14.3 years ago by Sukhi Singh

Hi, I am running tophat and getting a lot of different errors, I could understand some of them, but this one, I don't know the …