View Single Post
Old 10-25-2010, 03:50 AM   #2
Enrico Palazzo
Junior Member
 
Location: Hamburg, Germany

Join Date: Jul 2010
Posts: 9
Default same problem

Hi,

did you find a reason for this error? I get the same error running

tophat -p 8 --output-dir output/26_1 genome/Homo
_sapiens.hg19/hg19 readfiles/s_100617_3_1.fastq

******************************************************
Error log:
[Mon Oct 25 10:44:10 2010] Beginning TopHat run (v1.0.13)
-----------------------------------------------
[Mon Oct 25 10:44:10 2010] Preparing output location output/26_1/
[Mon Oct 25 10:44:10 2010] Checking for Bowtie index files
[Mon Oct 25 10:44:10 2010] Checking for reference FASTA file
[Mon Oct 25 10:44:10 2010] Checking for Bowtie
Bowtie version: 0.12.1.0
[Mon Oct 25 10:44:10 2010] Checking reads
seed length: 51bp
format: fastq
quality scale: phred33 (default)
[Mon Oct 25 10:54:50 2010] Mapping reads against hg19 with Bowtie
[Mon Oct 25 11:30:57 2010] Joining segment hits
Splitting reads into 2 segments
[Mon Oct 25 11:42:02 2010] Mapping reads against hg19 with Bowtie
[Mon Oct 25 11:47:40 2010] Mapping reads against hg19 with Bowtie
[Mon Oct 25 11:52:33 2010] Searching for junctions via segment mapping
[Mon Oct 25 12:01:35 2010] Retrieving sequences for splices
[Mon Oct 25 12:04:01 2010] Indexing splices
[Mon Oct 25 12:23:04 2010] Mapping reads against segment_juncs with Bowtie
[Mon Oct 25 12:26:55 2010] Mapping reads against segment_juncs with Bowtie
[Mon Oct 25 12:30:42 2010] Joining segment hits
Traceback (most recent call last):
File "tophat", line 1635, in <module>
sys.exit(main())
File "tophat", line 1595, in main
user_supplied_juncs)
File "tophat", line 1502, in spliced_alignment
stdout=open(merged_map,"w"))
File "/usr/lib64/python2.6/subprocess.py", line 444, in call
return Popen(*popenargs, **kwargs).wait()
File "/usr/lib64/python2.6/subprocess.py", line 595, in __init__
errread, errwrite)
File "/usr/lib64/python2.6/subprocess.py", line 1106, in _execute_child
raise child_exception
OSError: [Errno 2] No such file or directory
******************************************************

Last entry in run.log was:

sort -k 1,1n --temporary-directory=output/26_1/tmp/ output/26_1/tmp/file
vsmUly output/26_1/left_kept_reads.fq.candidate_hits.sam > output/26_1/tmp/file5jC6td

file5jC6td have been created in tmp/ but it's empty.

Any suggestions?

Thanks
Enrico Palazzo is offline   Reply With Quote