From cuffmerge to cuffdiff: What input for cuffdiff?
1
0
Entering edit mode
5.3 years ago
YaGalbi ★ 1.5k

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

rna-seq cuffmerge cuffdiff treated controls • 1.9k views
ADD COMMENT
3
Entering edit mode
5.3 years ago

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 COMMENT
0
Entering edit mode

You sir....are a legend.

ADD REPLY

Login before adding your answer.

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