Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
cufflinks - having difficulty reading reference GTF tguo Bioinformatics 31 07-07-2019 06:03 AM
cufflinks without reference stopped at chrM arrchi Bioinformatics 5 01-28-2016 01:08 AM
Cufflinks cufflinks v1.0.3 - segmentation fault bias correction chrNT annotations adrian Bioinformatics 0 06-08-2011 01:28 PM
cufflinks --reference-seq qusestion sumfry RNA Sequencing 1 04-18-2011 11:39 AM
transferring annotations from reference genome to the resequenced genome mike.t Bioinformatics 1 09-17-2010 05:35 AM

Thread Tools
Old 01-18-2012, 06:59 PM   #1
Location: Cleveland, OH

Join Date: Aug 2011
Posts: 11
Default Question about Cufflinks and reference annotations


I ran cufflinks with the -r option for performing reference annotation based transcript assembly. I used a reference GTF file downloaded from the UCSC genome browser (refFlat table). Here is an example of my output in the transcripts.gtf file:

chr1 Cufflinks transcript 1236828 1249909 1000 - . gene_id "CUFF.65"; transcript_id "CPSF3L"; FPKM "78.3979674437"; frac "0.481787"; conf_l
o "76.995976"; conf_hi "79.799959"; cov "1213.450041"; full_read_support "yes";
chr1 Cufflinks exon 1236828 1237167 1000 - . gene_id "CUFF.65"; transcript_id "CPSF3L"; exon_number "1"; FPKM "78.3979674437"; frac "0.481787
"; conf_lo "76.995976"; conf_hi "79.799959"; cov "1213.450041";
chr1 Cufflinks exon 1237261 1237390 1000 - . gene_id "CUFF.65"; transcript_id "CPSF3L"; exon_number "2"; FPKM "78.3979674437"; frac "0.481787
"; conf_lo "76.995976"; conf_hi "79.799959"; cov "1213.450041";
chr1 Cufflinks exon 1237469 1237611 1000 - . gene_id "CUFF.65"; transcript_id "CPSF3L"; exon_number "3"; FPKM "78.3979674437"; frac "0.481787
"; conf_lo "76.995976"; conf_hi "79.799959"; cov "1213.450041";

Cufflinks identifies the transcript as belonging to the CPSF3L gene. However, the gene id is listed as "CUFF.65".

Shouldn't the gene ID also be "CPSF3L"?

This is actually happening for the VAST majority of assembled transcripts; the program assigns the gene IDs "CUFF.x", but the transcript IDs are true gene identifiers.

Is this normal?
mhayes is offline   Reply With Quote
Old 03-23-2015, 06:40 PM   #2
Junior Member
Location: china

Join Date: Oct 2013
Posts: 5

I have the same question with you! I am confused...
If you figured it out, please let me know, thank you.
hubery_Bio 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 07:09 PM.

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