Sometime while doing Cufflinks for transcriptome assembly it pauses at certain point of time. Does anyone have faced it, please help me
•
written
by
qudrat
2 answers
For how long have you waited while it has been paused? It also pauses for me at a certain point during each run, but it always comes back provided that I have enough memory.
What specific command did you run/execute?
How much RAM have you got?
3 CPU cores is very little. In the past, I used as many as 32 CPU cores (human genome) and it still took time.
For a bacterial genome 1,000th the size of a human genome, it can still take 2-3 hours for Cufflinks to finish.
This is most probably because of highly expressed transcripts (like rRNA or those coming from mtDNA). To quickly check it you can summarize reads by featurecounts and see which genes have very abnormally high depth. Then you can mask that genes using –mask-file option in cufflinks http://cole-trapnell-lab.github.io/cufflinks/cufflinks/.
Overexpressed genes are not necessarily pseudogenes. I my case out of 100mln reads 30mln reads were mapped to mitochondrial gene and cufflinks was getting stuck on that one. But after masking it the problem was solved.
Log in to answer this question.
-
How to remove the isoform in the trasncript.fa file from salmon output?
written by mathavanbioinfo
-
de novo assembly of raw illumina reads of RNA seq experiment (24gb data)
written by dileepbiochem77 •
-
SpliceDB and SAM file
written by qudrat
-
cuffcompare and cuffmerge
written by qudrat
-
Cufflinks analysis of Tophat BAM file
written by qudrat
-
RNA-Seq Data analysis FPKM values mostly showing 0
written by simang5c •
-
Cufflink on multiple bam files
written by 1769mkc
-
StringTie: To use a reference gene-set or not?
written by i.sudbery
-
Gene ontology level
written by Sudhir Jadhao
-
Memory requirement for the transcriptome data analysis from ion PGM
written by iamtuttu5 •
