SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
what wrong with boost? feng Bioinformatics 3 10-16-2011 09:50 AM
something wrong when using consed dingkai0564 Bioinformatics 0 08-25-2011 12:05 AM
BWA weird thing CNVboy Bioinformatics 1 06-16-2011 07:34 PM
PubMed: Massively parallel sequencing: the next big thing in genetic medicine. Newsbot! Literature Watch 0 08-15-2009 05:05 AM
Spam Increase: I guess it's a good thing? ECO Site Announcements 0 06-27-2009 05:48 PM

Reply
 
Thread Tools
Old 08-11-2011, 12:51 AM   #1
chenyao
Member
 
Location: Beijing

Join Date: Jul 2011
Posts: 74
Default Any thing wrong with my cufflink process?

Recently, I run cufflink in galaxy on the internet. I want to compare two samples, However, I found no transcript or gene passed the significant level, even many of them have large FPKM in one sample and 0 FPKM in another sample.

Any thoughts?

Below is my cufflink process:

I have four samples belong to two group. the test have three samples, and the control has one sample.

First, using accept_hit.bam from tophat, I run cufflink without annotation on each sample.

Then, for the four "gtf" files from four samples, I run cuffcompare to combine these transcript and compare to the annotation genome. However, at this step, I found the transcript accuracy is very low.
See one example:
Missed exons: 10673/11776 ( 90.6%)
Wrong exons: 1254/2007 ( 62.5%)
Missed introns: 8529/8637 ( 98.7%)
Wrong introns: 2/5 ( 40.0%)
Missed loci: 0/504 ( 0.0%)
Wrong loci: 1248/2002 ( 62.3%)

at last, I run cufdiff between this two group sample.

Thank you.
chenyao is offline   Reply With Quote
Old 08-12-2011, 05:07 AM   #2
chenyao
Member
 
Location: Beijing

Join Date: Jul 2011
Posts: 74
Default

Any thoughts?
chenyao 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 08:43 AM.


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