SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
extract gene/transcript sequence from ttophat/cufflink/cuffdiff output mehtaaditya RNA Sequencing 2 02-24-2014 04:06 AM
Cuffdiff 2 progress reporting barkasn Bioinformatics 2 01-10-2014 01:14 AM
DNA amplification differences between samples & between regions (WGA & single locus) colinveal Sample Prep / Library Generation 0 11-09-2012 04:01 AM
using Cuffdiff without gene and transcript discory narges RNA Sequencing 7 09-07-2012 04:14 AM
Existence Genetics - analysis & reporting of close to 1,500 phenotypes and growing DoctorDNA Vendor Forum 3 11-17-2011 07:03 PM

Reply
 
Thread Tools
Old 11-25-2015, 08:40 AM   #1
Marta123
Junior Member
 
Location: Barcelona

Join Date: Nov 2015
Posts: 2
Default Cuffdiff reporting same locus several times &/or more than one gene per transcript

Hi all,

I am doing RNA-Seq of pig (Sus scrofa 10.2) and I got some weird results in cuffdiff output and I wonder if someone did ever find them. (I am using the lastest versions of TopHat and Cuffflinks as well as the latest genome and gtf files from ensembl).

1. More than one gene_id reported in one transcript

2. More than one transcript with exactly the same locus

e.g Some of the output that I obtained in genes.fpkm_tracking look like this:

XLOC_000006 ENSSSCG00000004008,ENSSSCG00000025112,ENSSSCG00000030155 TSS10,TSS11,TSS12,TSS13,TSS6,TSS7,TSS8,TSS9 1:935011-1369891
XLOC_000007 ENSSSCG00000021794 TSS14 1:935011-1369891

XLOC_000026 ENSSSCG00000026205,PARK2 TSS45,TSS46,TSS47,TSS48,TSS49 1:7195765-8262315
XLOC_000027 - TSS50 1:7195765-8262315
XLOC_000028 - TSS51 1:7195765-8262315

I obtained up to 246 transcripts with > 1 gene annotated per transcript.
Is it a program bug? Running edgeR parallely did not report this output... So I don't think they are fusion genes...

Any idea?
Marta123 is offline   Reply With Quote
Reply

Tags
cuffmerge, duplicate transcripts, rna-seq

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 11:17 AM.


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