SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
visualizing scripture output BED chrisbala Bioinformatics 6 12-06-2013 11:05 AM
cufflinks funny output, scripture comparison rcorbett Bioinformatics 11 12-25-2010 09:55 AM
scripture problem zwzhu Bioinformatics 0 08-27-2010 04:46 AM
scripture trouble memo Bioinformatics 1 08-13-2010 05:20 AM
Scripture output wenhuang Bioinformatics 1 07-28-2010 03:36 AM

Reply
 
Thread Tools
Old 10-21-2011, 01:29 PM   #1
nimmi
Member
 
Location: Bethesda

Join Date: Jul 2010
Posts: 15
Default scripture output

We did RNA-seq on 5 cases and 5 controls. Following the documentation given on scripture website I got some results for Chr1. Below are few lines from the output. There are a total of 64K rows for Chr1, of which ~62K have FWER<0.05 and last column (Poisson pvalue) is significant for all the rows. Why are so many transcripts significant? Did anyone else observe this trend in their data?

chrom chromStart chromEnd name score strand thickStart thickEnd itemRgb blockCount blockSizes blockStarts Family wise error rate (FWER) at the observed score Read enrichment (average reads per base over transcript / genome average) Total reads aligned transcript Average reads per base over transcript RPKM Genome average reads per base Transcript length Poisson nominal p-value for the observed count given the genome-wide average
chr1 569913 570307 chr1:569913-570307 40848.5882 + 569913 570307 0,0,0 6 1,2,1,7,4,53, 0,3,7,9,17,341, 1.11E-16 10181.53598 109494 1610.205882 40848.5882 0.158149604 68 0
chr1 16266622 16302620 chr1:16266622-16302620 68.37972491 - 16266622 16302620 0,0,0 15 2125,111,210,131,117,88,301,143,613,192,141,189,203,87,143, 0,2203,2401,2936,3471,3680,4172,4569,4809,5521,6026,6807,8167,32911,35855, 1.11E-16 17.0436889 12922 2.695452649 68.37972491 0.158149604 4794 0
chr1 16266622 16302620 chr1:16266622-16302620 68.36914145 - 16266622 16302620 0,0,0 15 2125,111,210,131,122,83,301,143,613,192,141,189,203,87,143, 0,2203,2401,2936,3471,3685,4172,4569,4809,5521,6026,6807,8167,32911,35855, 1.11E-16 17.04105096 12920 2.695035461 68.36914145 0.158149604 4794 0

thanks in advance,
Nirmala
nimmi 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 12:43 AM.


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