Attributes required for cuffdiff annotation file
0
0
Entering edit mode
9.7 years ago
bdeonovic ▴ 210

I want to use cufflinks to do differential expression analysis (8 samples, 4 cases, 4 controls). I have a reference annotation file (refFlat.gtf). In the manual for cuffdiff it states that cuffdiff requires certain attributes (i.e. tss_id, p_id) and that those can be added to the annotation file using the command

cuffcompare -s /path/to/genome_seqs.fa -CG -r annotation.gtf annotation.gtf

When I run this command I get a lot less lines in my annotation file:

[bdeonovic@neon-login-0-2 stoltz]$ wc -l cuffcmp.combined.gtf
35768 cuffcmp.combined.gtf

[bdeonovic@neon-login-0-2 stoltz]$ wc -l refFlat.gtf
78642 refFlat.gtf

Is this necessary? Is there a way to append these cufflink annotations without merging together my gene annotations?

cufflinks cuffdiff RNA-Seq • 2.2k views
ADD COMMENT

Login before adding your answer.

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