
<<Problem>> The cufflinks and cuffdiff results are not consistent with each other. This is killing me. Does it make sense? << Obseration >> We have 3 control samples and 3 treated sample. For many genes, their FPKM in cufflinks and cuffdiff are far from consistent. In cufflinks result, for a gene’s FPKM are control group (sample: 1,2,3):0, 0, 4.8 treated group(sample: 1,2,3):0, 0, 6.0 In cuffdiff, the estimated FPKM are control group: 12.6 treated group:2.0 <<Method>> Use ucsc gene annotation gtf file, mm9, downloaded from UCSC table database Use cufflinks on each individual sample. Cufflinks: galaxy mirror at cistrome, minimal count:10, no quantile normalization, use gtf as reference, no background correction Use cufflinks on treated groups (3 biological replicates) and control groups (3 biological replicates) Cuffdiff: galaxy mirror at cistrome, minimal count:10, no normalization, use gtf as reference, no background correction <<Additional Comments>> Cufflinks returns 55350 transcripts, while cuffdiff return 55418 transcripts, even though they use the same gene annotation gtf file. For the 6 cufflinks results (corresponding to 6 samples), the transcript ids are all the same, but the order are not, <<Question>> Does it make sense? Or did I do anything wrong?