SEQanswers

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

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

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

Hi,

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
hubery_Bio
Junior Member
 
Location: china

Join Date: Oct 2013
Posts: 5
Default

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
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 07:09 PM.


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