error after running tophat
1
2
Entering edit mode
7.8 years ago
Vasu ▴ 720

I had run the tophat script and it went through everything but in the last after "Reporting output tracks" it gave me the following error. Please help me. Thanks in Advance.

Error: [Errno 2] No such file or directory: 'SRR493360/thout/tmp/accepted_hits0_sorted.bam'
Found 46454 junctions from happy spliced reads
RNA-Seq tophat • 1.7k views
ADD COMMENT
2
Entering edit mode
7.8 years ago

Did you run out of space on the drive? That's the most common cause of this.

ADD COMMENT
0
Entering edit mode

Is it ? I'm not sure about that. Anyways thanks for the reply. I'l check once.

ADD REPLY
0
Entering edit mode

I had checked now. Still it has some 3G left. Is that free space enough for that?

ADD REPLY
1
Entering edit mode

It depends on how big your fastq files are. Tophat has to have room for ~3x the output, since it creates an unsorted BAM file and then sorts that (requiring creation of tmp files approximately the same size as the resulting sorted BAM files). So unless your fastq files are <300 megs or so, I wouldn't expect you to have enough room.

ADD REPLY
0
Entering edit mode

ok. Thankyou for the reply.

ADD REPLY
0
Entering edit mode

Subset 1000 reads and see if it runs fine. This helps to check if the issue is with hard drive space.

ADD REPLY
0
Entering edit mode

hw I need to run for 1000 reads ? Can you please help me?

ADD REPLY
2
Entering edit mode
zcat R1.fastq.gz | head -4000 > R1_subset.fastq
zcat R2.fastq.gz | head -4000 > R2_subset.fastq
ADD REPLY
0
Entering edit mode

Hi, could you show me how to do top hat on cluster, I have same problem and with huge file, so even I run "head 8000000 > 1_subset.fastq. But the remaining reads are more than 2 million, so I tried to run tail 8000000 > 2_subset.fastq. How should I deal with the middle part? And how to output the results, separately or combine together?

Thanks

ADD REPLY
0
Entering edit mode

You are right! The error is due to less space in the drive. Today I did it on cluster, it worked! Thank you

ADD REPLY
0
Entering edit mode

Hi, could you show me how to do top hat on cluster, I have same problem and with huge file, so even I run head 8000000 > 1_subset.fastq. But the remaining reads are more than 2 million, so I tried to run tail 8000000 > 2_subset.fastq. How should I deal with the middle part? And how to output the results, separately or combine together?

Thanks

ADD REPLY

Login before adding your answer.

Traffic: 1655 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6