Question: Tophat Error Unable to find file in temporary folder that is clearly there
gravatar for baunruh
2.1 years ago by
baunruh0 wrote:

So I am running an older version of tophat in order to validate some work that was done, unfortunately I am receiving an error that I cannot figure out. Most of the errors similar to this one online were an issue with RAM and running on a cluster. I am running on a cluster but I have determined it is not a RAM problem or a problem with the underlying cluster (I ran it on several different clusters that I have access to and received the same error). It appears that there is a problem calling back to the temporary file, but the file is clearly present in the temporary file!!

Code Ran:

tophat --library-type fr-secondstrand --color --quals -p 4 -G TandemUTR.gtf -o TandemUTR_GTF_out mm9 TandemUTR.csfasta TandemUTR.qual

Using Versions:

Error Message Below

[Mon Feb 12 16:49:13 2018] Beginning TopHat run (v1.3.3)
[Mon Feb 12 16:49:13 2018] Preparing output location /home/baunruh/ALEwork/Vollmer/CT0_Vollmer_gtf/
[Mon Feb 12 16:49:13 2018] Checking for Bowtie index files
[Mon Feb 12 16:49:13 2018] Checking for reference FASTA file
[Mon Feb 12 16:49:13 2018] Checking for Bowtie
    Bowtie version:
[Mon Feb 12 16:49:13 2018] Checking for Samtools
    Samtools Version: 0.1.19
[Mon Feb 12 16:49:13 2018] Generating SAM header for /home/baunruh/ALEwork/ZhangKerry/mm9/mm9
[Mon Feb 12 16:49:13 2018] Preparing reads
    format:      fasta
[Mon Feb 12 16:49:13 2018] Reading known junctions from GTF file
    Left  reads: min. length=50, count=32461692
[Mon Feb 12 16:56:37 2018] Mapping left_kept_reads against mm9 with Bowtie 

gzip: stdout: Broken pipe
[Mon Feb 12 16:56:37 2018] Processing bowtie hits
Traceback (most recent call last):
  File "/home/baunruh/tophat1.3.3/tophat", line 2604, in <module>
  File "/home/baunruh/tophat1.3.3/tophat", line 2563, in main
  File "/home/baunruh/tophat1.3.3/tophat", line 2218, in spliced_alignment
  File "/home/baunruh/tophat1.3.3/tophat", line 1820, in split_reads
    zreads = ZReader(reads_filename, params.system_params, False)
  File "/home/baunruh/tophat1.3.3/tophat", line 1190, in __init__
IOError: [Errno 2] No such file or directory: '/home/baunruh/ALEwork/Vollmer/CT0_Vollmer_gtf/tmp/left_kept_reads_missing.fq'
ADD COMMENTlink modified 2.1 years ago by genomax80k • written 2.1 years ago by baunruh0

There's no reason to use tophat anymore. Please switch to a more recent and supported aligner, such as STAR, hisat2, or BBMap.

ADD REPLYlink written 2.1 years ago by Devon Ryan94k

I know, its I am trying to repeat some steps to make sure we don't have errors.

ADD REPLYlink written 2.1 years ago by baunruh0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 782 users visited in the last hour