Question: From cuffmerge to cuffdiff: What input for cuffdiff?
gravatar for YaGalbi
4.1 years ago by
Biocomputing, MRC Harwell Institute, Oxford, UK
YaGalbi1.5k wrote:

I want to compare 2 conditions: controls (4 biological replicates) v treated (3 biological replicates). I have run cufflinks successfully on all 7 BAM files independently,

I then ran cuffmerge twice to produce 2 transctipt.gtf files:

1) merged_control_transcrpts.gtf

2) merged_treated_transcripts.gtf

I was expecting that the two transctipts.gtf files would then be used as input (with BAM files) for cuffdiff. But, cuffdiff only allows 1 transctipts.gtf. Where have I gone wrong? Am I supposed to merge everything? Can cuffdiff tell the treated group from the control group?

Thanks in advance. Kenneth

ADD COMMENTlink modified 4.1 years ago by onspotproductions140 • written 4.1 years ago by YaGalbi1.5k
gravatar for onspotproductions
4.1 years ago by
United States
onspotproductions140 wrote:

You need to merge all your cufflinks assemblies into a single gtf file using cuffmerge. Both conditions will be contained in a single file. You can then use this within cuffdiff and your original BAM files to test for your conditions.

ADD COMMENTlink written 4.1 years ago by onspotproductions140

You sir....are a legend.

ADD REPLYlink written 4.1 years ago by YaGalbi1.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1401 users visited in the last hour