Cufflinks takes too long time on my PC
0
0
Entering edit mode
5.1 years ago
takoyaki ▴ 120

I'm now running expression analysis from RNA-Seq sequence data.

I want to calculate FPKM value and make heat map by using R and I've already finish QC & mapping. So, I used cufflinks to calculate, and wanted to continue to next analysis like cuffmerge and cuffdiff.

However, cufflinks don't run smoothly ( it proceeds only 10% in an hour ).

QC was run by FastQC and fastq_quality_trimmer, and mapping was run by Hisat2. Also, Hisat2's result was sam, so I used samtools to convert sam to bam.

・This is commands that I actually ran.

"cufflinks -p 2 [annotation's gtf file] [bam file] -o [output folder]"

・This is the outputs displayed on screen

BAM record error: found spliced alignment without XS attribute

Processing Locus 6:160042173-160042246 [* ] 31%BAM record error: found spliced alignment without XS attribute

This is my computer's environment

・MacBook Pro, mac OS Mojave

・memory 8GB : CPU 2.9GHz, Intel Core i7, 2 cores

・annotation's gtf file 102MB , mapping results bam file 3.7 GB

Does anyone know why cufflinks takes so much time or is it usual in situation like my PC ?

Are there any mistake points or bad combination in my procedure ?

If my command is wrong, please tell me which is wrong or how should I write ? Also, do you think cufflinks is running normally considering about the outputs displayed on screen.

Thanks

RNA-Seq next-gen • 1.5k views
ADD COMMENT
2
Entering edit mode

A pipeline with low computational requirements is Salmon-tximport-DESeq2

ADD REPLY
0
Entering edit mode

Thank you so much ! I'll try it

ADD REPLY

Login before adding your answer.

Traffic: 1071 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