SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Cuffdiff splicing.diff Greg Bioinformatics 1 11-04-2013 04:04 PM
More diff expressed genes than transcripts with Cuffdiff drosoform RNA Sequencing 2 01-18-2013 02:28 PM
Diff Exp Post Cuffdiff studentofscience Bioinformatics 0 12-03-2012 01:06 AM
Cuffdiff NOTEST for splicing.diff Greg Bioinformatics 5 11-14-2011 05:22 AM
multiple FPKM problem for single gene in gene_exp.diff after running cuffdiff ngs RNA Sequencing 4 03-30-2011 01:55 PM

Reply
 
Thread Tools
Old 03-18-2013, 01:00 AM   #1
leifive
Member
 
Location: Nanjing, China

Join Date: Mar 2013
Posts: 10
Default Problem about splicing.diff from cuffdiff

Hi,
This problem has puzzled me one week long!
I followed the workflow:Tophat - Cufflinks - Cuffmerge - Cuffdiff,and everything is OK except the splicing.diff file generated by cuffdiff. After running Cuffdiff, it print "Performed 0 splicing tests", and splicing.diff file is filled up with LOWDATA and NOTEST without even one OK. I used merged.gtf as input, which has tss_id and p_id. And I also tried to use combined.gtf from Cuffcompare as input, but the result was the same. Can anyone help me?

Here is my command:
1).cuffmerge -o ./merged_asm/ -g genes.gtf -s genome.fa -p 8 assemblies.txt
2).cuffdiff -o ./diff_asm/ -b ./genome.fa -p 8 -L L1,L2,L3 -u ./merged_asm/merged.gtf ./Sample_out/accepted_hits.bam ./4h_out/accepted_hits.bam ./S48h12_out/accepted_hits.bam

Then is some lines of my merged.gtf
Code:
PHYCAscaffold_1	Cufflinks	exon	41144	42265	.	+	.	gene_id "XLOC_000001"; transcript_id "TCONS_00000001"; exon_number "1"; oId "transcriptId502542"; nearest_ref "proteinId502542"; class_code "o"; tss_id "TSS1";
PHYCAscaffold_1	Cufflinks	exon	41195	41689	.	+	.	gene_id "XLOC_000001"; transcript_id "TCONS_00000002"; exon_number "1"; oId "proteinId502542"; contained_in "TCONS_00000001"; nearest_ref "proteinId502542"; class_code "="; tss_id "TSS1"; p_id "P1";
PHYCAscaffold_1	Cufflinks	exon	42789	45015	.	+	.	gene_id "XLOC_000002"; transcript_id "TCONS_00000003"; exon_number "1"; oId "transcriptId502544"; nearest_ref "proteinId502544"; class_code "="; tss_id "TSS2"; p_id "P2";
leifive 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 06:18 AM.


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