Cuffmerge not working with reference genome

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Cuffmerge not working with reference genome

Mark Lindsay
Dear All

when I run Cuffmerge using the latest GenCode v19 GTF as a reference genome….this fails to run…tried this using multiple approaches and datasets.

Of note…...Cuffmerge works fine when the Gencode GTF is omitted as a reference genome. CuffCompare also works fine when the same Gencode GTF is used as a reference genome.

Presumably there is something wrong with the Cuffmerge set-up relating to the reference genome? Or am I doing something wrong.

Best wishes

Mark




___________________________________________________________
The Galaxy User List is being replaced by the Galaxy Biostar
User Support Forum at https://biostar.usegalaxy.org/

Posts to this list will be disabled in May 2014.  In the
meantime, you are encouraged to post all new questions to
Galaxy Biostar.

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:

  http://lists.bx.psu.edu/

To search Galaxy mailing lists use the unified search at:

  http://galaxyproject.org/search/mailinglists/
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Cuffmerge not working with reference genome

Jen Hillman-Jackson
Hello Mark,

Do you mean the GenCode reference annotation gtf dataset? This is not a reference genome. If you want to use a custom reference genome, this would be provided as a fasta file, using these instructions: https://wiki.galaxyproject.org/Support#Custom_reference_genome

My guess is that you have done the following (since a gtf file would not be accepted as a custom reference genome unless you changed the datatype to be "fasta"):

-> provided the GenCode gtf file as a "reference annotation" dataset and for mapping used a built-in "reference genome" or mapped somewhere else. When there are issues, it is most often a reference genome chromosome identifier mismatch problem. This wiki section explains this issue in more detail:
https://wiki.galaxyproject.org/Support#Reference_genomes.

Sometimes this issue can be corrected by altering the identifiers in the gtf file to match those in the reference genome (what your fastq data was mapped to, convert BAM->SAM to see the identifiers if necessary). This FAQ shows one method:
https://wiki.galaxyproject.org/Support#Tools_on_the_Main_server:_RNA-seq
-> https://usegalaxy.org/u/jeremy/p/transcriptome-analysis-faq#faq5

Hopefully the resources here help to resolve the issue!

Best,

Jen
Galaxy team

ps. This did not post to the mailing list because the "to" was not to just the mailing list. Please post new questions that way, or much better (since this mailing list will be retired very soon), please activate your account at Galaxy Biostar and post your question there (this has replaced the galaxy-user mailing list). Here is how: https://wiki.galaxyproject.org/Support#Biostar

On 5/20/14 12:55 AM, Mark Lindsay wrote:
Dear All

when I run Cuffmerge using the latest GenCode v19 GTF as a reference genome….this fails to run…tried this using multiple approaches and datasets.

Of note…...Cuffmerge works fine when the Gencode GTF is omitted as a reference genome. CuffCompare also works fine when the same Gencode GTF is used as a reference genome.

Presumably there is something wrong with the Cuffmerge set-up relating to the reference genome? Or am I doing something wrong.

Best wishes

Mark





-- 
Jennifer Hillman-Jackson
http://galaxyproject.org

___________________________________________________________
The Galaxy User List is being replaced by the Galaxy Biostar
User Support Forum at https://biostar.usegalaxy.org/

Posts to this list will be disabled in May 2014.  In the
meantime, you are encouraged to post all new questions to
Galaxy Biostar.

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:

  http://lists.bx.psu.edu/

To search Galaxy mailing lists use the unified search at:

  http://galaxyproject.org/search/mailinglists/
Jennifer Hillman-Jackson
http://galaxyproject.org
Loading...