Entering edit mode
5.8 years ago
k.kathirvel93
▴
300
Hi Everyone,
I am working in RNA-seq data analysis, I tried Cufflinks for Transcript Quantification from STAR mapped BAM file, but its not showing any progress, can anyone find me solution for this ? am using version 2.2.1 .
Thanks in advance.
Here is the Command line :
cufflinks -p 16 --library-type fr-unstranded -o STAR_Output/
-G /home/bioinformatics/Documents/Kathir/Glaucoma/RNA/Reference/ensemble_GRCh38.77/Homo_sapiens.GRCh38.77.gtf
-u /home/bioinformatics/Documents/Kathir/Glaucoma/RNA/STAR/STAR_Output--RNAAligned.sortedByCoord.out.bam
How did you reach that conclusion? It is not possible to give you any guidance with just that. Are you able to see the process run (using
top)? Is there an output file being produced. Is there a space missing in this last bit of the command line somewhere/home/bioinformatics/Documents/Kathir/Glaucoma/RNA/STAR/STAR_Output--RNAAligned.sortedByCoord.out.bamor perhapscufflinksis confused by--you have in that name (if it is legit).I have been waiting for more than 3 hours, No outputs created in the folder -o STAR_Output/ . and
topcommand Doesn't showing cufflinks is running. I have been used cufflinks three months back successfully.Are you running this directly on the command line (without a job scheduler)? Does
ps -ef | grep cuffshow a process that looks like cufflinks?Are you running this directly on the command line (without a job scheduler)?
Yes.
and
ps -ef | grebcuff showing this.I have a feeling that the process may be a zombie, if nothing is being written to output dir. You should kill it and start over. I am still not sure about the strange
--inclusion in the file name. I would suggest replacing that with an_.You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon followed by DESEq2 or edgeR.
Hi WouterDeCoster, Are you suggesting me to use StringtTie instead of Cufflinks for Quantification?.
Yes, see also here: How does Cufflinks compare to StringTie?