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
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...Yesterday, 05:41 AM
Comment