SEQanswers

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 03:06 PM
Questions about Cuffmerge and Cuffdiff aquleaf Bioinformatics 1 08-08-2012 09:41 AM
Total RNAseq and cufflinks/cuffmerge/cuffdiff/cummeRbund pipeline JWC Bioinformatics 8 08-01-2012 05:41 AM
cufflinks XLOC identifiers golharam Bioinformatics 0 05-31-2012 10:42 AM
Thoughts on running cufflinks with cuffmerge output? hcosta Bioinformatics 0 03-23-2012 06:43 PM

Reply
 
Thread Tools
Old 08-31-2012, 03:56 AM   #1
AdamB
Member
 
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, 09:58 AM   #2
Krish_143
Member
 
Location: Sweden

Join Date: Jan 2012
Posts: 45
Default

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

Join Date: Apr 2010
Posts: 43
Default

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

Join Date: Feb 2014
Posts: 206
Default

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, 07:12 PM   #5
danwiththeplan
Member
 
Location: Auckland

Join Date: Sep 2011
Posts: 72
Default

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
Reply

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 03:02 AM.


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