Colleagues,

I am having trouble running cufflinks with an annotation file on the public Galaxy.  The assembled transcripts gtf file has all the FPKM at 0 although the gene expression and transcript expression tab files have values for FPKM.

I have seen the SEQanswers threads about the compatibility of tophat bam files relative to chromosomes labeled as 1,2,3… versus Chr1, Chr2, Chr3…   I am using the iGenomes bovine UMD3.1 genome and annotation file (chromosomes are  1,2,3) from the history.  I altered the gtf file to Chr1, Chr2, Chr3… but it did not help.

Another potential discrepency/conflict is that the genome and gtf file have the bosTau6 database attribute from when I uploaded them. However I am running them from the history (bosTau6 is not an option for tophat).  I do not seem to be able to remove the attribute.

Am I missing something else?

Here is the command line

Info: cufflinks v1.0.3
cufflinks -q --no-update-check -I 50000 -F 0.050000 -j 0.050000 -p 8 -G /galaxy/main_database/files/003/142/dataset_3142240.dat -N -b ref.fa

Here is the details page

Tool: Cufflinks
Name: Cufflinks on data 69, data 4, and data 5: assembled transcripts
Created: Nov 09, 2011
Filesize: 44.6 Mb
Dbkey: bosTau6
Format: gtf
Tool Version: 

Input Parameter  Value
SAM or BAM file of aligned RNA-Seq reads  4: Tophat for Illumina on data 4 and data 69: accepted_hits
Max Intron Length  50000
Min Isoform Fraction  0.05
Pre MRNA Fraction  0.05
Perform quartile normalization  Yes
Conditional (reference_annotation)  1
Reference Annotation  5: iGen_UMD3_1_genes.gtf
Conditional (bias_correction)  0
Conditional (seq_source)  1
Using reference file  69: UMD31_iGen_1-29X.fa
Conditional (singlePaired)  0

Cordially,

Chris