problem with cuffdiff
Hy guys, I am performing some RNA sequencing. I am using a gtf file as a reference annotation for the cufflink - cuffmerge steps. Then when I do the cuffdiff I find something strange. Where there is the gene name, the corresponding values (all of them!) are 0, where there is no gene name there are values of expression, fold change etc. Does anyone have an idea of what is going on? I haven't find any similar question in the Galaxy archive. Also, I tried to visualize my outputs (TopHat Cufflink Cuffmerge) using Trackster but it gives me an error. Is anyone experiencing the same problem? Thanks, Giuseppe
Hello Giuseppe, If you haven't resolved this yet, then I would suggested double checking that the chromosome identifiers in the reference annotation file are /exactly the same/ as in the reference genome. You can run Picard's tools to get summary stats on Tophat's BAM outout files to see the reference genome idenifiers, or use SAMTools' BAM->SAM to examine. In general, it is simpler to change the reference annotation's identifiers to match those in reference genome. As the file is small, a simple find/replace in a text editor is a common solution. Hope this helps, Jen Galaxy team On 9/10/13 3:22 PM, Ianiri, Giuseppe wrote:
Hy guys, I am performing some RNA sequencing. I am using a gtf file as a reference annotation for the cufflink - cuffmerge steps. Then when I do the cuffdiff I find something strange. Where there is the gene name, the corresponding values (all of them!) are 0, where there is no gene name there are values of expression, fold change etc. Does anyone have an idea of what is going on? I haven't find any similar question in the Galaxy archive. Also, I tried to visualize my outputs (TopHat Cufflink Cuffmerge) using Trackster but it gives me an error. Is anyone experiencing the same problem? Thanks,
Giuseppe
___________________________________________________________ The Galaxy User list should be used for the discussion of Galaxy analysis and other features on the public server at usegalaxy.org. Please keep all replies on the list by using "reply all" in your mail client. For discussion of local Galaxy instances and the Galaxy source code, please use the Galaxy Development list:
http://lists.bx.psu.edu/listinfo/galaxy-dev
To manage your subscriptions to this and other Galaxy lists, please use the interface at:
To search Galaxy mailing lists use the unified search at:
-- Jennifer Hillman-Jackson http://galaxyproject.org
participants (2)
-
Ianiri, Giuseppe
-
Jennifer Jackson