We use the tophat-cufflinks-cuffmerge-cuffdiff pipeline for RNA Seq analysis. I am wondering if there is a way to have "cufflinks" - the command not the package - just do the transcript assembly since we never use the fpkm outputs from it anyway?
-
-
Just run Cufflinks alone. Cufflinks provides a transcript.gtf file which consists of all possible transcripts. You could use gffread then and provide your reference genome and this gtf to obtain a transcript fasta file. For your case, you just wouldn't continue with any of the other modules they have for downstream analysis. Is this what you were asking for? -
@cedance: Thanks for your response. My post perhaps was not clear. After the cufflinks run, we use cuffmerge to generate merged gtf and then run cuffdiff on it on each of our bam files. But, I see that cufflinks also spends time doing abundance estimates using the GTF file provided for guidance assembly (RABT) - I dont think these abundance results are used by any downstream application. I was wondering if there was a way to avoid this.Comment
-
Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.
We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing...
Comment