changeset 0:5f784a95ce13 draft

planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/fisher_test commit 7d954797362ea5de843d812cb823ac1dd869e1df
author artbio
date Mon, 28 Aug 2017 13:25:56 -0400
parents
children 3a67bee2ef63
files fisher_test.xml test-data/counts.tab test-data/fisher.tab
diffstat 3 files changed, 1074 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/fisher_test.xml	Mon Aug 28 13:25:56 2017 -0400
@@ -0,0 +1,74 @@
+<tool id="fishertest" name="Fisher's exact test" version="1.0.0">
+	<description>on two gene hit lists</description>
+        <requirements>
+<!--            <requirement type="package" version="3.1.2">R</requirement> -->
+            <requirement type="package" version="2.4.2=r3.3.1_0">bioconductor-qvalue</requirement>
+        </requirements>
+        <command><![CDATA[
+            Rscript '$fisher_test' "\${GALAXY_SLOTS:-1}"
+    ]]></command>
+  <configfiles>
+    <configfile name="fisher_test">
+    <![CDATA[
+      ## Setup R error handling to go to stderr
+      options( show.error.messages=F,  error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
+      options(warn=-1)
+      suppressMessages(library(qvalue))
+      library(parallel)
+      args = commandArgs(trailingOnly = TRUE)
+      slots = as.numeric(args[1])
+      countsTable = read.delim("${input}", header=TRUE, check.names=FALSE, stringsAsFactor = TRUE)
+      depth1 = sum(countsTable[,2])
+      depth2 = sum(countsTable[,3])
+      float_table=data.frame(countsTable[,2], countsTable[,3])
+
+      calc_pvalue <- function(row, depth1, depth2, ... ){
+        thearray = array( c(row, (depth1 - row[1]), (depth2 - row[2])), dim=c(2,2))
+        current_test = fisher.test( thearray )
+        return(current_test\$p.value)
+      }
+
+      cl <- makePSOCKcluster(slots)
+      clusterExport(cl=cl, varlist=c("calc_pvalue", "depth1", "depth2"))
+      ptm <- proc.time()
+      p_val = parApply(cl, float_table, 1, function(x) calc_pvalue(x, depth1, depth2))
+      stopCluster(cl)
+      proc.time() - ptm
+      p_val[p_val>1]=1
+      p = qvalue(p_val)
+      finalTable = cbind(countsTable, data.frame(p\$pvalues), data.frame(p\$qvalues))
+      write.table ( finalTable, file = "${output}", row.names=FALSE, col.names=TRUE, quote= FALSE, dec = ".", sep = "\t", eol = "\n")
+    ]]>
+    </configfile>
+  </configfiles>
+        <inputs>
+                <param name="input" type="data" format="tabular" label="gene hit lists, 2 samples"/>
+        </inputs>
+        <outputs>
+                <data name="output" format="tabular" label="Fisher test p-values" />
+        </outputs>
+<tests>
+    <test>
+        <param name="input" value="counts.tab" ftype="tabular"/>
+        <output name="output" file="fisher.tab" ftype="tabular"/>
+    </test>
+</tests>
+<help>
+
+**What it does**
+
+Runs Fisher's exact test for testing the null of independence of rows and columns in a contingency table of two columns.
+
+p.pvalues: the chance of getting this data if it is independent between columns (false negative); the p-value.
+
+q.qvalues: FDR (Faslse Detection Rate) adjusted p-values; a q-value of 0.05 implies that 5% of significant tests will result in false positives.
+
+Be aware that this test does not take into account the biological noise that would be visible if replicates were available.
+
+
+</help>
+<citations>
+    <citation type="doi">10.1111/1467-9868.00346</citation>
+</citations>
+    
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/counts.tab	Mon Aug 28 13:25:56 2017 -0400
@@ -0,0 +1,500 @@
+Geneid	2: Oocyte 1: accepted_hits	4: Pronuclei 1: accepted_hits
+NM_001195025	0	0
+NM_028778	0	0
+NM_001290393	0	0
+NM_008922	3133	2288
+NM_175370	4	7
+NM_001290390	0	0
+NM_175642	977	1013
+NM_001290392	0	0
+NM_027671	0	0
+NM_178884	0	0
+NM_199021	332	1436
+NM_198680	0	0
+NM_011465	0	0
+NM_130890	0	0
+NM_001113393	0	0
+NM_001113392	0	0
+NM_001113391	19	88
+NM_031162	0	0
+NR_103716	0	0
+NM_001160371	0	0
+NM_177834	1	2
+NM_145381	606	367
+NM_177173	0	0
+NM_001160369	0	0
+NM_001160370	2	9
+NM_001289497	0	0
+NM_027974	0	6
+NM_133235	74	0
+NM_028096	77	2
+NM_201365	1	0
+NM_001126046	0	0
+NM_001081198	31	144
+NM_001160386	890	1226
+NM_009418	456	774
+NM_001114663	369	353
+NM_001252070	346	463
+NR_045151	0	0
+NM_001081050	93	35
+NM_007642	0	0
+NM_025728	0	0
+NM_029160	0	0
+NM_175210	290	346
+NM_010154	31	23
+NM_175291	0	0
+NM_010472	265	1666
+NM_023314	0	20
+NM_027456	0	0
+NM_030184	5	11
+NM_009183	0	0
+NM_011066	40	12
+NM_145519	32	140
+NM_173760	8045	6907
+NM_172851	0	0
+NM_019933	125	54
+NM_176957	0	0
+NM_001081756	0	0
+NM_181750	1448	1651
+NM_019777	0	0
+NM_172484	0	0
+NM_028320	46	70
+NM_021605	790	949
+NM_173437	1290	1131
+NM_001287179	2333	2559
+NM_175443	0	0
+NM_007495	0	0
+NM_022312	0	0
+NM_001081290	3128	1875
+NM_013862	79	236
+NM_001205204	0	0
+NM_001159525	0	0
+NM_023041	1	2
+NM_022563	464	427
+NM_011880	0	0
+NM_026626	31	0
+NM_001199003	0	0
+NM_027188	248	359
+NM_001243792	0	1
+NM_133684	2924	2392
+NM_001011874	88	28
+NM_001159750	0	0
+NM_001159751	3	0
+NM_001195662	0	1
+NM_001177795	0	0
+NM_133826	2265	2743
+NM_011541	0	0
+NM_021374	0	0
+NM_001290372	0	0
+NM_001244693	0	0
+NM_001244692	0	0
+NR_045189	0	0
+NR_045188	0	0
+NM_173868	0	0
+NR_015470	62	217
+NM_133220	0	0
+NM_145692	2	6
+NM_026493	37	10
+NR_110971	0	0
+NM_177547	0	0
+NM_029525	0	0
+NM_001033636	11	11
+NM_001077695	0	0
+NM_001302702	0	0
+NM_001081209	0	2
+NM_172841	76	564
+NM_172294	0	0
+NM_001198565	0	0
+NM_001198566	0	0
+NM_008678	2	1
+NM_010164	12	10
+NM_001252192	0	0
+NM_001098528	62	55
+NM_025303	4	6
+NM_133832	7337	4876
+NM_001033288	15	8
+NM_001111272	0	0
+NM_026456	1919	2756
+NM_020604	35	27
+NR_130314	7	1
+NM_153154	8	18
+NM_153179	30	130
+NM_001077353	0	0
+NM_001288617	0	0
+NM_010356	0	0
+NM_023872	0	0
+NM_001160139	0	0
+NM_183018	0	0
+NM_001012623	0	0
+NM_053270	0	0
+NM_001012624	0	0
+NM_001012625	0	0
+NR_110972	0	0
+NM_028534	22	21
+NM_001290683	0	0
+NM_007733	111	280
+NM_007740	0	0
+NM_026719	1136	514
+NM_001290691	0	0
+NM_001033764	0	0
+NM_001013750	0	0
+NM_177235	43	42
+NM_133833	0	0
+NM_134448	0	0
+NM_001276764	7	28
+NM_183019	9	2
+NM_001277076	12	68
+NM_001039551	0	0
+NM_053186	0	0
+NM_198006	817	1015
+NM_001289766	0	0
+NM_009539	0	0
+NM_018872	3195	2254
+NM_001289765	0	0
+NM_001290721	0	0
+NM_001290814	15	332
+NM_019570	822	545
+NM_027111	0	0
+NM_207228	0	0
+NM_001290720	2	0
+NR_110976	31	36
+NM_008719	579	349
+NM_018775	9019	10143
+NM_145142	90	9
+NM_001123382	0	0
+NM_001252202	0	0
+NM_001252200	0	0
+NM_001294171	0	0
+NM_001252201	1	0
+NM_008696	0	0
+NM_008362	0	0
+NM_010743	0	0
+NM_001025602	0	0
+NM_001033289	0	2
+NR_027826	0	13
+NM_010879	0	0
+NM_001013025	1593	2602
+NM_023645	1984	999
+NM_028450	191	89
+NM_007737	0	0
+NM_028599	2428	2519
+NM_019790	35	105
+NM_010863	0	0
+NM_001290982	0	0
+NM_001161817	0	0
+NM_001142647	0	0
+NM_153556	495	503
+NM_146108	0	0
+NR_024590	0	0
+NM_030025	25	53
+NM_172655	0	0
+NM_001001883	78	147
+NM_001081433	13714	20555
+NM_031179	28668	20953
+NM_139146	319	229
+NM_199007	130	41
+NM_001177867	0	0
+NM_001164566	0	0
+NM_025824	1651	1530
+NM_023617	0	0
+NM_144882	0	0
+NM_172513	4659	7514
+NM_001013771	83	56
+NM_001033373	558	310
+NM_028717	0	0
+NM_172406	9	346
+NM_001159948	0	0
+NM_175200	38	14
+NM_001045513	458	714
+NM_173444	444	408
+NM_139150	0	0
+NM_007561	3461	3400
+NM_030013	81	461
+NM_001037725	7010	6779
+NM_001285463	0	0
+NM_001077404	0	0
+NM_001077405	0	0
+NM_001077406	0	0
+NM_001077407	0	0
+NM_010939	0	0
+NM_001077403	0	0
+NR_045915	18	26
+NM_001285937	0	0
+NM_001285936	0	0
+NM_001081658	0	0
+NR_045914	0	2
+NM_011780	0	2
+NM_001177600	0	0
+NM_009952	0	0
+NM_001097644	8358	3205
+NM_001039493	1595	1557
+NM_001037726	0	0
+NM_133828	0	0
+NM_033563	4320	2839
+NM_011086	1063	819
+NM_139270	0	0
+NM_175510	105	116
+NM_021285	0	0
+NM_008632	0	0
+NM_001039934	0	0
+NM_177645	3	6
+NM_001122738	0	0
+NM_001113387	0	0
+NM_001080809	242	444
+NM_011770	185	205
+NM_177164	137	61
+NM_001276408	0	0
+NM_001276412	0	0
+NM_001276411	0	0
+NM_001276410	0	0
+NM_001276409	0	0
+NM_010233	0	0
+NM_001276413	0	0
+NM_001005423	109	135
+NM_008342	0	2
+NM_001045533	627	848
+NM_176980	0	7
+NM_001170641	0	0
+NM_027884	0	0
+NR_045311	0	0
+NM_001034060	0	0
+NM_001289895	0	0
+NM_029342	121	618
+NM_009518	0	0
+NM_001039495	0	0
+NM_001159884	0	0
+NM_001159885	0	0
+NM_001159883	0	0
+NM_178055	0	0
+NM_020266	0	0
+NM_001014974	258	657
+NM_176972	4705	5534
+NR_037976	0	0
+NM_007936	5937	4482
+NM_008781	0	0
+NM_026713	0	0
+NM_001159520	0	0
+NM_001004173	149	46
+NM_001285927	0	0
+NM_009255	23375	13582
+NM_172849	41	206
+NM_029777	0	0
+NM_007735	0	0
+NM_007734	145	233
+NM_172430	0	17
+NM_133975	22051	37255
+NM_152915	98	4
+NM_001003948	2771	2921
+NM_001081746	17	11
+NM_013673	30	70
+NM_001144993	0	0
+NM_001144992	0	0
+NM_001172157	0	0
+NM_001081082	0	0
+NM_008801	252	1115
+NM_153530	0	0
+NM_201644	0	0
+NM_001110192	0	0
+NM_001205391	0	0
+NM_019867	0	2
+NM_001205392	0	0
+NM_201410	0	0
+NM_029846	0	0
+NM_001110212	0	0
+NM_001111314	0	23
+NM_146112	0	0
+NM_010566	0	0
+NR_122078	0	0
+NM_201642	0	0
+NM_201641	0	0
+NM_001110193	0	0
+NM_177646	1676	2127
+NM_145079	0	0
+NM_029269	0	0
+NM_139152	0	0
+NM_029122	7144	9137
+NM_178119	2	11
+NM_001037136	0	0
+NM_001111312	0	0
+NM_001033408	7	337
+NM_001243009	0	0
+NM_001243008	0	0
+NM_023046	0	0
+NM_001039126	0	0
+NM_207225	338	255
+NM_007855	0	0
+NM_001271884	0	0
+NM_207031	0	1
+NM_172463	33	16
+NM_008798	0	0
+NM_174874	198	623
+NR_040467	0	0
+NM_028942	0	0
+NR_120500	0	64
+NM_027584	0	0
+NM_001164233	0	0
+NM_013626	5	14
+NM_011800	468	279
+NM_173187	0	0
+NM_009399	0	0
+NM_013784	1048	1091
+NM_178779	0	0
+NM_029349	0	0
+NM_001160368	0	0
+NM_009741	2	37
+NM_133821	96	34
+NM_198028	0	0
+NM_001160307	0	0
+NR_102307	0	0
+NR_102306	0	0
+NM_172853	0	31
+NM_001081386	196	135
+NM_001077425	0	25
+NM_029709	0	0
+NM_001293300	0	1
+NM_001293301	0	0
+NM_145506	51	74
+NM_022327	0	0
+NM_001081125	561	386
+NM_001113416	1958	2652
+NM_001172074	3	0
+NM_133186	0	0
+NM_001085409	0	0
+NM_028439	0	0
+NM_178082	0	0
+NM_029955	0	0
+NM_133748	0	0
+NM_001271532	0	0
+NM_001271531	0	0
+NM_001025156	0	0
+NM_027677	0	0
+NM_145100	0	0
+NM_028135	0	47
+NM_001033041	0	0
+NM_001285945	0	0
+NM_145128	0	41
+NM_027678	0	0
+NM_177445	331	743
+NM_026390	7243	6281
+NM_172485	329	620
+NM_026976	0	0
+NM_009581	8	104
+NM_007799	0	0
+NM_175382	3678	1498
+NM_172510	17	40
+NM_001114662	96	89
+NM_175294	0	0
+NM_001145804	0	0
+NM_001081011	107	390
+NM_001160318	0	0
+NM_001099276	501	655
+NM_001160317	0	0
+NM_001160316	0	0
+NM_182716	0	0
+NM_172516	4014	2948
+NM_144530	0	0
+NM_001276767	0	0
+NM_001145300	0	0
+NM_001081307	54	88
+NM_053106	0	0
+NM_001033409	60	74
+NM_009307	32	96
+NM_001291015	0	0
+NM_177642	0	0
+NM_030676	0	0
+NM_001159769	0	0
+NM_011210	0	0
+NM_001111316	0	0
+NM_001268286	0	0
+NM_181347	0	0
+NM_001166501	0	0
+NM_133239	208	490
+NR_040455	0	0
+NM_001081027	0	0
+NR_040604	0	0
+NM_020025	4	0
+NM_013835	246	188
+NM_145991	821	1000
+NR_110891	0	0
+NM_001290269	0	0
+NM_153171	0	0
+NM_001145807	0	0
+NM_153539	0	0
+NM_008869	0	0
+NM_001110146	0	0
+NM_024458	0	0
+NM_001159730	0	0
+NR_040603	0	0
+NM_021400	1	7
+NM_001024720	0	0
+NM_197990	349	349
+NM_025819	6926	5296
+NM_022018	25	32
+NM_001039644	217	78
+NM_177756	203	319
+NM_175460	333	138
+NM_016846	27	71
+NM_010877	0	24
+NM_007842	43	63
+NR_108050	0	0
+NR_046171	0	0
+NM_009782	47	63
+NM_011273	954	2136
+NM_021433	2003	2262
+NM_177838	0	0
+NM_001145781	188	300
+NM_001277730	0	0
+NM_028250	8	79
+NM_001039184	229	433
+NM_001134741	3	6
+NM_026683	0	0
+NM_023141	0	0
+NR_027640	1	0
+NM_177644	242	374
+NM_001159968	0	0
+NM_001159965	0	0
+NM_001159967	0	0
+NM_001159966	0	0
+NM_023884	0	0
+NM_001085376	89	30
+NM_011931	10995	17175
+NM_009786	55	26
+NM_001038621	529	1262
+NR_038031	0	0
+NR_038030	0	0
+NM_001162896	0	0
+NM_001039482	112	251
+NM_080844	2371	1696
+NM_001038619	0	0
+NM_172646	0	0
+NM_008030	37	100
+NM_011347	2140	2125
+NM_010629	27	98
+NM_001293313	0	0
+NM_138314	0	0
+NM_001081126	2	0
+NM_028759	133	184
+NM_145512	241	271
+NM_175296	21	1183
+NM_001164306	0	12
+NM_001160049	0	0
+NM_001033344	0	0
+NM_198932	0	0
+NM_011137	0	0
+NM_030724	61	592
+NM_033652	6	76
+NM_175461	0	0
+NM_001160261	0	0
+NM_001160262	0	0
+NM_183355	0	0
+NM_008783	0	0
+NM_001291509	0	0
+NM_001291508	0	0
+NM_023284	8286	12399
+NM_133806	1017	1226
+NM_001081304	3522	2308
+NM_001198955	1615	1879
+NM_001109985	145	48
+NM_027528	2	0
+NM_011070	1429	1293
+NM_013489	0	0
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/fisher.tab	Mon Aug 28 13:25:56 2017 -0400
@@ -0,0 +1,500 @@
+Geneid	2: Oocyte 1: accepted_hits	4: Pronuclei 1: accepted_hits	p.pvalues	p.qvalues
+NM_001195025	0	0	1	1
+NM_028778	0	0	1	1
+NM_001290393	0	0	1	1
+NM_008922	3133	2288	3.43464044020283e-41	4.89681594188918e-40
+NM_175370	4	7	0.550305086197412	1
+NM_001290390	0	0	1	1
+NM_175642	977	1013	0.73622185992851	1
+NM_001290392	0	0	1	1
+NM_027671	0	0	1	1
+NM_178884	0	0	1	1
+NM_199021	332	1436	9.94388878387307e-152	4.96200050315266e-150
+NM_198680	0	0	1	1
+NM_011465	0	0	1	1
+NM_130890	0	0	1	1
+NM_001113393	0	0	1	1
+NM_001113392	0	0	1	1
+NM_001113391	19	88	3.85805937131219e-11	2.37675509417874e-10
+NM_031162	0	0	1	1
+NR_103716	0	0	1	1
+NM_001160371	0	0	1	1
+NM_177834	1	2	1	1
+NM_145381	606	367	2.03973563063845e-17	1.6685706224403e-16
+NM_177173	0	0	1	1
+NM_001160369	0	0	1	1
+NM_001160370	2	9	0.0663808053789194	0.207025136775505
+NM_001289497	0	0	1	1
+NM_027974	0	6	0.0315574986825609	0.104981278950653
+NM_133235	74	0	7.70743067983725e-24	8.18299555157189e-23
+NM_028096	77	2	7.51953165095876e-22	7.65764549760903e-21
+NM_201365	1	0	0.487180835678035	1
+NM_001126046	0	0	1	1
+NM_001081198	31	144	1.44609315021036e-17	1.20347533866442e-16
+NM_001160386	890	1226	7.67992992651892e-10	4.40492532567005e-09
+NM_009418	456	774	1.88962170331458e-16	1.47331442180309e-15
+NM_001114663	369	353	0.205350188340196	0.602763199892693
+NM_001252070	346	463	0.000721229722747277	0.00272646690644615
+NR_045151	0	0	1	1
+NM_001081050	93	35	4.23810487973389e-08	2.15797381121144e-07
+NM_007642	0	0	1	1
+NM_025728	0	0	1	1
+NM_029160	0	0	1	1
+NM_175210	290	346	0.121625356917143	0.363419479650624
+NM_010154	31	23	0.221704114073477	0.639481808801533
+NM_175291	0	0	1	1
+NM_010472	265	1666	2.03104264467209e-233	1.12610031076819e-231
+NM_023314	0	20	2.14887474703225e-06	1.00213878389635e-05
+NM_027456	0	0	1	1
+NM_030184	5	11	0.212514364701721	0.620144257228999
+NM_009183	0	0	1	1
+NM_011066	40	12	4.02537810500803e-05	0.00017022573511856
+NM_145519	32	140	2.5862764748091e-16	1.98546455527653e-15
+NM_173760	8045	6907	1.8325879068154e-36	2.47151720405645e-35
+NM_172851	0	0	1	1
+NM_019933	125	54	1.30109281920789e-08	6.76297204984103e-08
+NM_176957	0	0	1	1
+NM_001081756	0	0	1	1
+NM_181750	1448	1651	0.0266628384791318	0.0892936671213876
+NM_019777	0	0	1	1
+NM_172484	0	0	1	1
+NM_028320	46	70	0.0515714493177683	0.164962520574144
+NM_021605	790	949	0.00615483226156037	0.0217819950249548
+NM_173437	1290	1131	7.28695703860234e-06	3.2758482542906e-05
+NM_001287179	2333	2559	0.150894310882959	0.445540006689919
+NM_175443	0	0	1	1
+NM_007495	0	0	1	1
+NM_022312	0	0	1	1
+NM_001081290	3128	1875	2.2249971323935e-86	5.55136784532179e-85
+NM_013862	79	236	1.07437257606142e-17	9.24330888714909e-17
+NM_001205204	0	0	1	1
+NM_001159525	0	0	1	1
+NM_023041	1	2	1	1
+NM_022563	464	427	0.0477845290978515	0.154834285843038
+NM_011880	0	0	1	1
+NM_026626	31	0	2.07980876187599e-10	1.23550544306681e-09
+NM_001199003	0	0	1	1
+NM_027188	248	359	0.000110641231631144	0.000452540775278205
+NM_001243792	0	1	1	1
+NM_133684	2924	2392	3.38966093808169e-20	3.25277078481301e-19
+NM_001011874	88	28	3.51544698520491e-09	1.97102027597444e-08
+NM_001159750	0	0	1	1
+NM_001159751	3	0	0.115629376515014	0.347584692054169
+NM_001195662	0	1	1	1
+NM_001177795	0	0	1	1
+NM_133826	2265	2743	7.16004962638453e-07	3.40272834625322e-06
+NM_011541	0	0	1	1
+NM_021374	0	0	1	1
+NM_001290372	0	0	1	1
+NM_001244693	0	0	1	1
+NM_001244692	0	0	1	1
+NR_045189	0	0	1	1
+NR_045188	0	0	1	1
+NM_173868	0	0	1	1
+NR_015470	62	217	9.75454644242473e-20	9.18399749956592e-19
+NM_133220	0	0	1	1
+NM_145692	2	6	0.290352459876478	0.804921541546459
+NM_026493	37	10	2.99696454699474e-05	0.000127819257175246
+NR_110971	0	0	1	1
+NM_177547	0	0	1	1
+NM_029525	0	0	1	1
+NM_001033636	11	11	1	1
+NM_001077695	0	0	1	1
+NM_001302702	0	0	1	1
+NM_001081209	0	2	0.500327785998395	1
+NM_172841	76	564	5.4537790032973e-88	1.60084454273256e-86
+NM_172294	0	0	1	1
+NM_001198565	0	0	1	1
+NM_001198566	0	0	1	1
+NM_008678	2	1	0.615637735484249	1
+NM_010164	12	10	0.671901714559592	1
+NM_001252192	0	0	1	1
+NM_001098528	62	55	0.357047830849683	0.984347334773435
+NM_025303	4	6	0.754711865198588	1
+NM_133832	7337	4876	1.01308036780091e-142	4.21272586277212e-141
+NM_001033288	15	8	0.144068159654634	0.427916736117037
+NM_001111272	0	0	1	1
+NM_026456	1919	2756	4.22075212541596e-26	4.89803560600597e-25
+NM_020604	35	27	0.253269070492179	0.710007113346051
+NR_130314	7	1	0.034678455580826	0.113845719308106
+NM_153154	8	18	0.0779921226704073	0.240234995139094
+NM_153179	30	130	4.39220943548333e-15	3.13101786900883e-14
+NM_001077353	0	0	1	1
+NM_001288617	0	0	1	1
+NM_010356	0	0	1	1
+NM_023872	0	0	1	1
+NM_001160139	0	0	1	1
+NM_183018	0	0	1	1
+NM_001012623	0	0	1	1
+NM_053270	0	0	1	1
+NM_001012624	0	0	1	1
+NM_001012625	0	0	1	1
+NR_110972	0	0	1	1
+NM_028534	22	21	0.762721917781111	1
+NM_001290683	0	0	1	1
+NM_007733	111	280	3.34455176601523e-16	2.52868383521454e-15
+NM_007740	0	0	1	1
+NM_026719	1136	514	2.15516220788999e-61	3.70836531633485e-60
+NM_001290691	0	0	1	1
+NM_001033764	0	0	1	1
+NM_001013750	0	0	1	1
+NM_177235	43	42	0.745911161145958	1
+NM_133833	0	0	1	1
+NM_134448	0	0	1	1
+NM_001276764	7	28	0.000581278208632895	0.00221418187868561
+NM_183019	9	2	0.0343617757726414	0.113553153049987
+NM_001277076	12	68	3.25179717952554e-10	1.90899622656852e-09
+NM_001039551	0	0	1	1
+NM_053186	0	0	1	1
+NM_198006	817	1015	0.000406543449914975	0.00158488423052791
+NM_001289766	0	0	1	1
+NM_009539	0	0	1	1
+NM_018872	3195	2254	4.34985298912476e-49	6.78305200491643e-48
+NM_001289765	0	0	1	1
+NM_001290721	0	0	1	1
+NM_001290814	15	332	9.04205916160961e-76	1.87999480068466e-74
+NM_019570	822	545	2.67466821500669e-17	2.15267651498119e-16
+NM_027111	0	0	1	1
+NM_207228	0	0	1	1
+NM_001290720	2	0	0.237344728677233	0.669124404575928
+NR_110976	31	36	0.715293388464529	1
+NM_008719	579	349	6.06669236706449e-17	4.80520554153203e-16
+NM_018775	9019	10143	3.38150776900152e-06	1.54804805204749e-05
+NM_145142	90	9	5.79736921738959e-19	5.25979498086801e-18
+NM_001123382	0	0	1	1
+NM_001252202	0	0	1	1
+NM_001252200	0	0	1	1
+NM_001294171	0	0	1	1
+NM_001252201	1	0	0.487180835678035	1
+NM_008696	0	0	1	1
+NM_008362	0	0	1	1
+NM_010743	0	0	1	1
+NM_001025602	0	0	1	1
+NM_001033289	0	2	0.500327785998395	1
+NR_027826	0	13	0.000256698995061036	0.0010195961497084
+NM_010879	0	0	1	1
+NM_001013025	1593	2602	7.59724771500808e-45	1.11500782640854e-43
+NM_023645	1984	999	7.26734757261621e-86	1.64836656306159e-84
+NM_028450	191	89	6.13762713018285e-11	3.7349706560503e-10
+NM_007737	0	0	1	1
+NM_028599	2428	2519	0.617111589466462	1
+NM_019790	35	105	1.2909047866103e-08	6.76297204984103e-08
+NM_010863	0	0	1	1
+NM_001290982	0	0	1	1
+NM_001161817	0	0	1	1
+NM_001142647	0	0	1	1
+NM_153556	495	503	0.590093047218317	1
+NM_146108	0	0	1	1
+NR_024590	0	0	1	1
+NM_030025	25	53	0.00315162890558066	0.0113961074194547
+NM_172655	0	0	1	1
+NM_001001883	78	147	2.39032697047553e-05	0.000103719405066721
+NM_001081433	13714	20555	6.13594302649696e-244	4.37405081460284e-242
+NM_031179	28668	20953	0	0
+NM_139146	319	229	1.00577044138633e-05	4.44141106417501e-05
+NM_199007	130	41	4.98226470083442e-13	3.31486678095517e-12
+NM_001177867	0	0	1	1
+NM_001164566	0	0	1	1
+NM_025824	1651	1530	0.000326016892914442	0.00128096401231737
+NM_023617	0	0	1	1
+NM_144882	0	0	1	1
+NM_172513	4659	7514	2.25389654339631e-121	8.65149519349813e-120
+NM_001013771	83	56	0.0107049806419357	0.0373551422400415
+NM_001033373	558	310	3.13710421486099e-20	3.06944118277575e-19
+NM_028717	0	0	1	1
+NM_172406	9	346	2.59614755135677e-86	6.16894108631919e-85
+NM_001159948	0	0	1	1
+NM_175200	38	14	0.000440677960081124	0.00170463800062388
+NM_001045513	458	714	3.37138213113651e-11	2.1028996042964e-10
+NM_173444	444	408	0.0505359801146179	0.162693252110931
+NM_139150	0	0	1	1
+NM_007561	3461	3400	0.00413540044650379	0.0148457900921251
+NM_030013	81	461	1.46572272289095e-61	2.61212728115208e-60
+NM_001037725	7010	6779	4.73318607350146e-07	2.27101908718964e-06
+NM_001285463	0	0	1	1
+NM_001077404	0	0	1	1
+NM_001077405	0	0	1	1
+NM_001077406	0	0	1	1
+NM_001077407	0	0	1	1
+NM_010939	0	0	1	1
+NM_001077403	0	0	1	1
+NR_045915	18	26	0.365843978424177	1
+NM_001285937	0	0	1	1
+NM_001285936	0	0	1	1
+NM_001081658	0	0	1	1
+NR_045914	0	2	0.500327785998395	1
+NM_011780	0	2	0.500327785998395	1
+NM_001177600	0	0	1	1
+NM_009952	0	0	1	1
+NM_001097644	8358	3205	0	0
+NM_001039493	1595	1557	0.0349887340201788	0.11411358350372
+NM_001037726	0	0	1	1
+NM_133828	0	0	1	1
+NM_033563	4320	2839	1.0958144252589e-87	3.03784110113439e-86
+NM_011086	1063	819	1.64711002146542e-11	1.04038974773575e-10
+NM_139270	0	0	1	1
+NM_175510	105	116	0.736843026188035	1
+NM_021285	0	0	1	1
+NM_008632	0	0	1	1
+NM_001039934	0	0	1	1
+NM_177645	3	6	0.509101902312159	1
+NM_001122738	0	0	1	1
+NM_001113387	0	0	1	1
+NM_001080809	242	444	1.50608747258379e-12	9.76022920544559e-12
+NM_011770	185	205	0.648418012664926	1
+NM_177164	137	61	8.04305615044474e-09	4.36248371638253e-08
+NM_001276408	0	0	1	1
+NM_001276412	0	0	1	1
+NM_001276411	0	0	1	1
+NM_001276410	0	0	1	1
+NM_001276409	0	0	1	1
+NM_010233	0	0	1	1
+NM_001276413	0	0	1	1
+NM_001005423	109	135	0.223592919821118	0.641223373509988
+NM_008342	0	2	0.500327785998395	1
+NM_001045533	627	848	1.77181465307074e-06	8.3409010554934e-06
+NM_176980	0	7	0.0158403421085777	0.0545126256012432
+NM_001170641	0	0	1	1
+NM_027884	0	0	1	1
+NR_045311	0	0	1	1
+NM_001034060	0	0	1	1
+NM_001289895	0	0	1	1
+NM_029342	121	618	5.71406052461676e-76	1.23970269642772e-74
+NM_009518	0	0	1	1
+NM_001039495	0	0	1	1
+NM_001159884	0	0	1	1
+NM_001159885	0	0	1	1
+NM_001159883	0	0	1	1
+NM_178055	0	0	1	1
+NM_020266	0	0	1	1
+NM_001014974	258	657	1.44901527080536e-36	2.00849616703298e-35
+NM_176972	4705	5534	1.62661141193038e-08	8.36782571704392e-08
+NR_037976	0	0	1	1
+NM_007936	5937	4482	4.11676776757565e-65	7.60839672600093e-64
+NM_008781	0	0	1	1
+NM_026713	0	0	1	1
+NM_001159520	0	0	1	1
+NM_001004173	149	46	3.44855690948509e-15	2.4939563736711e-14
+NM_001285927	0	0	1	1
+NM_009255	23375	13582	0	0
+NM_172849	41	206	6.15604290010093e-26	6.9815122889781e-25
+NM_029777	0	0	1	1
+NM_007735	0	0	1	1
+NM_007734	145	233	5.66606436000516e-05	0.000235613842970215
+NM_172430	0	17	1.6634721170077e-05	7.28133847707756e-05
+NM_133975	22051	37255	0	0
+NM_152915	98	4	8.72976339677448e-26	9.68033763331214e-25
+NM_001003948	2771	2921	0.957494918497635	1
+NM_001081746	17	11	0.256924502620669	0.716230876020747
+NM_013673	30	70	0.00018516086697903	0.000745123166310774
+NM_001144993	0	0	1	1
+NM_001144992	0	0	1	1
+NM_001172157	0	0	1	1
+NM_001081082	0	0	1	1
+NM_008801	252	1115	1.52833227702074e-120	5.44741290166678e-119
+NM_153530	0	0	1	1
+NM_201644	0	0	1	1
+NM_001110192	0	0	1	1
+NM_001205391	0	0	1	1
+NM_019867	0	2	0.500327785998395	1
+NM_001205392	0	0	1	1
+NM_201410	0	0	1	1
+NM_029846	0	0	1	1
+NM_001110212	0	0	1	1
+NM_001111314	0	23	2.78866828318439e-07	1.36426026795001e-06
+NM_146112	0	0	1	1
+NM_010566	0	0	1	1
+NR_122078	0	0	1	1
+NM_201642	0	0	1	1
+NM_201641	0	0	1	1
+NM_001110193	0	0	1	1
+NM_177646	1676	2127	8.96848311468942e-09	4.81212158519357e-08
+NM_145079	0	0	1	1
+NM_029269	0	0	1	1
+NM_139152	0	0	1	1
+NM_029122	7144	9137	4.22803717617408e-36	5.55208039713386e-35
+NM_178119	2	11	0.023011125297134	0.0775848075896613
+NM_001037136	0	0	1	1
+NM_001111312	0	0	1	1
+NM_001033408	7	337	1.27511717522188e-86	3.3488603707143e-85
+NM_001243009	0	0	1	1
+NM_001243008	0	0	1	1
+NM_023046	0	0	1	1
+NM_001039126	0	0	1	1
+NM_207225	338	255	5.49199370856382e-05	0.0002302945260986
+NM_007855	0	0	1	1
+NM_001271884	0	0	1	1
+NM_207031	0	1	1	1
+NM_172463	33	16	0.0098346492224228	0.0345597884647111
+NM_008798	0	0	1	1
+NM_174874	198	623	2.07810359944566e-47	3.14234453370722e-46
+NR_040467	0	0	1	1
+NM_028942	0	0	1	1
+NR_120500	0	64	2.83303172762956e-19	2.61793117053176e-18
+NM_027584	0	0	1	1
+NM_001164233	0	0	1	1
+NM_013626	5	14	0.0651976139111578	0.204613895230615
+NM_011800	468	279	2.35350027532509e-14	1.65407977096792e-13
+NM_173187	0	0	1	1
+NM_009399	0	0	1	1
+NM_013784	1048	1091	0.811604090414923	1
+NM_178779	0	0	1	1
+NM_029349	0	0	1	1
+NM_001160368	0	0	1	1
+NM_009741	2	37	4.01680304456512e-09	2.22709413248666e-08
+NM_133821	96	34	6.28426511846587e-09	3.44598713638953e-08
+NM_198028	0	0	1	1
+NM_001160307	0	0	1	1
+NR_102307	0	0	1	1
+NR_102306	0	0	1	1
+NM_172853	0	31	1.22787550447007e-09	6.96261223557458e-09
+NM_001081386	196	135	0.000141381740564455	0.000573573077574498
+NM_001077425	0	25	7.16506569230139e-08	3.5753677804584e-07
+NM_029709	0	0	1	1
+NM_001293300	0	1	1	1
+NM_001293301	0	0	1	1
+NM_145506	51	74	0.088751460937878	0.271699257717798
+NM_022327	0	0	1	1
+NM_001081125	561	386	1.0139928482173e-10	6.09617387060763e-10
+NM_001113416	1958	2652	1.44706453546825e-17	1.20347533866442e-16
+NM_001172074	3	0	0.115629376515014	0.347584692054169
+NM_133186	0	0	1	1
+NM_001085409	0	0	1	1
+NM_028439	0	0	1	1
+NM_178082	0	0	1	1
+NM_029955	0	0	1	1
+NM_133748	0	0	1	1
+NM_001271532	0	0	1	1
+NM_001271531	0	0	1	1
+NM_001025156	0	0	1	1
+NM_027677	0	0	1	1
+NM_145100	0	0	1	1
+NM_028135	0	47	2.54018064494081e-14	1.76048630809092e-13
+NM_001033041	0	0	1	1
+NM_001285945	0	0	1	1
+NM_145128	0	41	1.4387019292302e-12	9.44621398270881e-12
+NM_027678	0	0	1	1
+NM_177445	331	743	1.28296140404297e-32	1.56145790394498e-31
+NM_026390	7243	6281	4.86292196745426e-30	5.77761443276113e-29
+NM_172485	329	620	2.47916924579126e-18	2.20911688151757e-17
+NM_026976	0	0	1	1
+NM_009581	8	104	2.19954042980037e-21	2.19514134894077e-20
+NM_007799	0	0	1	1
+NM_175382	3678	1498	5.03627751760019e-235	3.14137810160312e-233
+NM_172510	17	40	0.00500305918836176	0.0178323323928037
+NM_001114662	96	89	0.418582308229198	1
+NM_175294	0	0	1	1
+NM_001145804	0	0	1	1
+NM_001081011	107	390	7.40876932817108e-36	9.47942537117274e-35
+NM_001160318	0	0	1	1
+NM_001099276	501	655	0.000257453135998513	0.0010195961497084
+NM_001160317	0	0	1	1
+NM_001160316	0	0	1	1
+NM_182716	0	0	1	1
+NM_172516	4014	2948	5.33175115321359e-51	8.58239943694704e-50
+NM_144530	0	0	1	1
+NM_001276767	0	0	1	1
+NM_001145300	0	0	1	1
+NM_001081307	54	88	0.011624939201388	0.0402836434825877
+NM_053106	0	0	1	1
+NM_001033409	60	74	0.387891180978853	1
+NM_009307	32	96	6.38057383998422e-08	3.21606701631528e-07
+NM_001291015	0	0	1	1
+NM_177642	0	0	1	1
+NM_030676	0	0	1	1
+NM_001159769	0	0	1	1
+NM_011210	0	0	1	1
+NM_001111316	0	0	1	1
+NM_001268286	0	0	1	1
+NM_181347	0	0	1	1
+NM_001166501	0	0	1	1
+NM_133239	208	490	4.21258416929332e-24	4.56973804451601e-23
+NR_040455	0	0	1	1
+NM_001081027	0	0	1	1
+NR_040604	0	0	1	1
+NM_020025	4	0	0.0563321044266021	0.178713719110149
+NM_013835	246	188	0.000908449339122787	0.00338295686733038
+NM_145991	821	1000	0.00193441724417056	0.00709760444736109
+NR_110891	0	0	1	1
+NM_001290269	0	0	1	1
+NM_153171	0	0	1	1
+NM_001145807	0	0	1	1
+NM_153539	0	0	1	1
+NM_008869	0	0	1	1
+NM_001110146	0	0	1	1
+NM_024458	0	0	1	1
+NM_001159730	0	0	1	1
+NR_040603	0	0	1	1
+NM_021400	1	7	0.0710297642097876	0.220148151184373
+NM_001024720	0	0	1	1
+NM_197990	349	349	0.519564231352806	1
+NM_025819	6926	5296	9.79564547122755e-71	1.88001041928559e-69
+NM_022018	25	32	0.508865369520654	1
+NM_001039644	217	78	6.04684457264497e-18	5.29364112587692e-17
+NM_177756	203	319	7.41602080877898e-06	3.30410212819707e-05
+NM_175460	333	138	4.69853525747343e-22	4.88451894474842e-21
+NM_016846	27	71	2.65936860136119e-05	0.000114398701041313
+NM_010877	0	24	1.41321732244732e-07	6.98213310793279e-07
+NM_007842	43	63	0.0989272048843532	0.301004117300563
+NR_108050	0	0	1	1
+NR_046171	0	0	1	1
+NM_009782	47	63	0.216227177458618	0.627310241580525
+NM_011273	954	2136	1.55298456500033e-90	4.84337061209479e-89
+NM_021433	2003	2262	0.0219621561584464	0.0745518090004405
+NM_177838	0	0	1	1
+NM_001145781	188	300	6.82306067358054e-06	3.09518843283336e-05
+NM_001277730	0	0	1	1
+NM_028250	8	79	3.13007447007452e-15	2.29692229495175e-14
+NM_001039184	229	433	2.44670300744058e-13	1.64987135231466e-12
+NM_001134741	3	6	0.509101902312159	1
+NM_026683	0	0	1	1
+NM_023141	0	0	1	1
+NR_027640	1	0	0.487180835678035	1
+NM_177644	242	374	2.73653897453358e-06	1.26438235952987e-05
+NM_001159968	0	0	1	1
+NM_001159965	0	0	1	1
+NM_001159967	0	0	1	1
+NM_001159966	0	0	1	1
+NM_023884	0	0	1	1
+NM_001085376	89	30	1.03530046506338e-08	5.4959035326237e-08
+NM_011931	10995	17175	3.01921101978145e-246	2.5109771647849e-244
+NM_009786	55	26	0.000748314374841077	0.00280758551162178
+NM_001038621	529	1262	2.57158430529302e-61	4.27740189447073e-60
+NR_038031	0	0	1	1
+NR_038030	0	0	1	1
+NM_001162896	0	0	1	1
+NM_001039482	112	251	5.41762264412941e-12	3.46588935823151e-11
+NM_080844	2371	1696	1.1156099981696e-34	1.39172347271658e-33
+NM_001038619	0	0	1	1
+NM_172646	0	0	1	1
+NM_008030	37	100	3.03079923348826e-07	1.46831924030159e-06
+NM_011347	2140	2125	0.0565867086561193	0.178713719110149
+NM_010629	27	98	4.74621994819649e-10	2.75391134203494e-09
+NM_001293313	0	0	1	1
+NM_138314	0	0	1	1
+NM_001081126	2	0	0.237344728677233	0.669124404575928
+NM_028759	133	184	0.0181179792778175	0.0619237784906228
+NM_145512	241	271	0.479203402432331	1
+NM_175296	21	1183	5.96841568668976e-306	5.95647885531638e-304
+NM_001164306	0	12	0.000509510969368339	0.00195573825934462
+NM_001160049	0	0	1	1
+NM_001033344	0	0	1	1
+NM_198932	0	0	1	1
+NM_011137	0	0	1	1
+NM_030724	61	592	1.33758473960323e-104	4.44969856708008e-103
+NM_033652	6	76	6.05582372902713e-16	4.51023289669334e-15
+NM_175461	0	0	1	1
+NM_001160261	0	0	1	1
+NM_001160262	0	0	1	1
+NM_183355	0	0	1	1
+NM_008783	0	0	1	1
+NM_001291509	0	0	1	1
+NM_001291508	0	0	1	1
+NM_023284	8286	12399	3.67036038919774e-143	1.66500894019061e-141
+NM_133806	1017	1226	0.00139054826650452	0.00513987840730191
+NM_001081304	3522	2308	2.29029980995533e-72	4.57143842067083e-71
+NM_001198955	1615	1879	0.00313575635267887	0.0113961074194547
+NM_001109985	145	48	9.35513784370869e-14	6.3948134027543e-13
+NM_027528	2	0	0.237344728677233	0.669124404575928
+NM_011070	1429	1293	8.09329248335086e-05	0.0003337647065448
+NM_013489	0	0	1	1