Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
cuffcompare or cuffmerge for cuffdiff upper Bioinformatics 6 04-07-2015 04:06 PM
Questions about Cuffmerge and Cuffdiff aquleaf Bioinformatics 1 08-08-2012 10:41 AM
Total RNAseq and cufflinks/cuffmerge/cuffdiff/cummeRbund pipeline JWC Bioinformatics 8 08-01-2012 06:41 AM
cufflinks XLOC identifiers golharam Bioinformatics 0 05-31-2012 11:42 AM
Thoughts on running cufflinks with cuffmerge output? hcosta Bioinformatics 0 03-23-2012 07:43 PM

Thread Tools
Old 08-31-2012, 04:56 AM   #1
Location: uk

Join Date: Apr 2010
Posts: 43
Default XLOC identifiers from cufflinks/cuffmerge/cuffdiff

Is there a way to stop cufflinks replacing Ensembl gene IDs from the reference gtf with its own XLOC identifiers? Or is it possible to convert all the XLOCs back to Ensembl gene IDs after running cuffmerge/cuffdiff?
AdamB is offline   Reply With Quote
Old 09-01-2012, 10:58 AM   #2
Location: Sweden

Join Date: Jan 2012
Posts: 45

Use -G option i.e GTF file for tophat and cufflinks.
Krish_143 is offline   Reply With Quote
Old 09-03-2012, 09:21 AM   #3
Location: uk

Join Date: Apr 2010
Posts: 43

I should have mentioned - this is when I use the -g or -G option.
AdamB is offline   Reply With Quote
Old 03-05-2015, 04:27 AM   #4
Senior Member
Location: UK

Join Date: Feb 2014
Posts: 206

I'd like to go back to .gtf file of reference genome and parse the Ensembl ID. But it seems that not clever. Does anyone have better or easier idea?
super0925 is offline   Reply With Quote
Old 01-24-2016, 08:12 PM   #5
Location: Auckland

Join Date: Sep 2011
Posts: 72

That's not what the -G option does. It means "only quantitate against genes supplied in the GTF file". While -g means "assemble new genes as well, if there is read support, and quantitate them too".

I believe that all genes get an XLOC number in the gene_id field, but this is for cufflinks to internally be able to distinguish genes. However any genes that have a name supplied in an associated GTF/GFF file will also have a name (from the gtf file) in the gene_name field. You can see this if/when you interrogate your data in cummeRbund /R.
danwiththeplan is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 06:40 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO