Error RNA-seq data, TOPHAT
1
0
Entering edit mode
7.5 years ago
David_emir ▴ 430

Hi All,

I tried aligning the RNA-seq reads to the genome using tophat, used the following command

[root@BIO-DT-415 RNA_SEQ]# tophat -p 8 -G genes.gtf -o C1_R1_thout genome GSM794483_C1_R1_1.fq GSM794483_C1_R1_2.fq


It's throwing error:

/usr/local/bin/tophat -p 8 -G genes.gtf -o C1_R1_thout genome GSM794483_C1_R1_1.fq GSM794483_C1_R1_2.fq
/usr/local/bin/gtf_juncs genes.gtf  > C1_R1_thout/tmp/genes.juncs
2014-11-26 16:25:31] Reporting output tracks
[FAILED]



I am not able to decode this, Please let me know if you have any solution for this.

Thanks a lot

-Ateeq Khaliq

RNA-Seq • 2.7k views
1
Entering edit mode

Try running that command manually (you'll need to be in the same directory) and then see what the actual error message is. Often this means you've run out of memory, though tophat isn't without its own bugs.

0
Entering edit mode

Thanks Devon , Thanks for coming to my rescue again !!! I have increased the RAM size to 16 gigs, now its running fine . Hoping to finish this protocol by this week!!! I am sure I will be back again to this forum.. Thank you so much ...its great help for me !!!

3
Entering edit mode
7.5 years ago
Manvendra Singh ★ 2.2k

how much memory are you providing on 8 threads?

1
Entering edit mode

Thanks Manu, you pointed it right, I have increased the ram size, now its working fine. Thanks a lot !!!

0
Entering edit mode

Welcome !

Good!