Repository 'isocor'
hg clone https://toolshed.g2.bx.psu.edu/repos/gmat/isocor

Changeset 1:fa183805db31 (2020-05-04)
Previous changeset 0:96f7617b4848 (2019-05-27)
Commit message:
"planemo upload"
modified:
isocor.xml
test-data/Metabolites.dat
test-data/results.log
test-data/results.tsv
b
diff -r 96f7617b4848 -r fa183805db31 isocor.xml
--- a/isocor.xml Mon May 27 03:33:36 2019 -0400
+++ b/isocor.xml Mon May 04 11:38:09 2020 -0400
[
b'@@ -1,22 +1,22 @@\n-<tool id="isocor" name="Isotope Correction for mass spectrometry labeling experiments" version="0.1.2">\n+<tool id="isocor" name="Isotope Correction for mass spectrometry labeling experiments" version="0.1.3">\n     <requirements>\n-        <requirement type="package" version="2.1.3">isocor</requirement>\n+        <requirement type="package" version="2.2.0">isocor</requirement>\n     </requirements>\n     <command detect_errors="exit_code"><![CDATA[\n-        isocorcli -t "$tracer" -M "$input1" -D "$input2" -I "$input3"\n-#if $tracer_purity:\n-    -p $tracer_purity\n+        isocorcli -t "$tracer" -M "$db.input1" -D "$db.input2" -I "$db.input3"\n+#if $tcorrection.tracer_purity:\n+    -p $tcorrection.tracer_purity\n #end if\n-#if $resolution:\n-    -r $resolution\n+#if $HR.resolution:\n+    -r $HR.resolution\n #end if\n-#if $mz:\n-    -m $mz\n+#if $HR.mz:\n+    -m $HR.mz\n #end if \n-#if $resolution_formula_code:\n-    -f $resolution_formula_code\n+#if $HR.resolution_formula_code:\n+    -f $HR.resolution_formula_code\n #end if\n-#if $correct_na_tracer:\n+#if $tcorrection.correct_na_tracer:\n     -n\n #end if\n #if $verbose_log:\n@@ -26,45 +26,81 @@\n     ]]></command>\n     <inputs>\n         <param type="text" name="tracer" argument="-t" label="the isotopic tracer (e.g. \'13C\')" />\n-        <param type="data" name="input1" argument="-M" label="metabolites db" format="tabular" />\n-        <param type="data" name="input2" argument="-D" label="derivatives db" format="tabular" />\n-        <param type="data" name="input3" argument="-I" label="isotopes db" format="csv" />\n-        <param type="integer" optional="true" argument="-r" name="resolution" label="HR only: resolution of the mass spectrometer (e.g. 1e4)" />\n-        <param type="integer" optional="true" argument="-m" name="mz" label="HR only: mz at which resolution is given (e.g. \'400\')" />\n-        <param type="boolean" optional="true" argument="-n" name="correct_na_tracer" label="flag to correct tracer natural abundance" />\n-        <param type="boolean" optional="true" argument="-v" name="verbose_log" label="flag to enable verbose logs" />\n-        <param type="text" optional="true" argument="-p" name="tracer_purity" label="purity vector of the tracer" />\n+        <param type="data" name="input4" label="measurements Data file" format="tabular" help="This file contains the raw MS data for each metabolite of each sample." />\n+        <section name="db" title="Database Files" expanded="true" help="The exact mass and natural abundance of each isotope and the elemental formulas used for correction have to be defined carefully, otherwise the correction will be wrong. IsoCor rely on several flat-files to store this information. Pre-configured files are shipped with IsoCor. For a good start download them on github">\n+            <param type="data" name="input1" argument="-M" label="metabolites db" format="tabular" />\n+            <param type="data" name="input2" argument="-D" label="derivatives db" format="tabular" />\n+            <param type="data" name="input3" argument="-I" label="isotopes db" format="csv" />\n+        </section>\n+        <section name="HR" title="High Resolution Parameters" expanded="true" help="For measurements collected at high or ultrahigh resolution (e.g. on Orbitrap or FT-ICR instruments). Measurements collected at unitary resolution (e.g. on quadrupole instruments) are Low Resolution, don\'t fill this section.">\n+            <param name="resolution_formula_code" optional="true" argument="-f" type="select" label="HR only: spectrometer formula code">\n+                <option value="orbitrap">Orbitrap</option>\n+                <option value="ft-icr">ft-icr</option>\n+                <option value="constant">Constant</option>\n+                <option value="datafile">Datafile</option>\n+            </param>\n+            <param type="integer" optional="true" argument="-r" name="resolution" label="HR only: resolution of the mass spectrometer (e.g. \'60000\' or \'1e4\')" />\n+          '..b'ue">\n+            <param type="boolean" optional="true" argument="-n" name="correct_na_tracer" label="flag to correct tracer natural abundance" help="Correct for the contribution of naturally occurring isotopes of the tracer element at unlabeled positions. This only concerns the tracer element: natural abundance of other elements is always corrected." />\n+            <param type="text" optional="true" argument="-p" name="tracer_purity" label="purity vector of the tracer (e.g. \'0,0,0.99,0.01\')" help="Correct for the contribution of isotopic impurities of the tracer at labeled positions. The isotopic purity is typically obtained from the manufacturer. e.g. \xc2\xb9\xc2\xb3C-substates with purity of 99%, use 0.01 for \xc2\xb9\xc2\xb2C and 0.99 for \xc2\xb9\xc2\xb3C." />\n+        </section>\n+        <param type="boolean" optional="true" argument="-v" name="verbose_log" label="flag to enable verbose logs" help="Useful in case of trouble. Join it to the issue on github." />\n \n-\t<param name="resolution_formula_code" optional="true" argument="-f" type="select" label="HR only: spectrometer formula code">\n-\t\t<option value="orbitrap">Orbitrap</option>\n-  \t\t<option value="ft-icr">ft-icr</option>\n-  \t\t<option value="constant">Constant</option>\n-  \t\t<option value="datafile">Datafile</option>\n-\t</param>\n \n-        <param type="data" name="input4" label="measurements Data file" format="tabular" />\n     </inputs>\n     <outputs>\n-        <data name="output1" label="corrected measurements data file" format="csv" />\n+        <data name="output1" label="corrected measurements data file" format="tabular" />\n         <data name="output2" label="processing log file" format="txt" />\n     </outputs>\n     <tests>\n         <test>\n             <param name="tracer" value="13C"/>\n-            <param name="input1" value="Metabolites.dat"/>\n-            <param name="input2" value="Derivatives.dat"/>\n-            <param name="input3" value="Isotopes.dat"/>\n+            <section name="db">\n+                <param name="input1" value="Metabolites.dat"/>\n+                <param name="input2" value="Derivatives.dat"/>\n+                <param name="input3" value="Isotopes.dat"/>\n+            </section>\n             <param name="input4" value="Data_example.tsv"/>\n             <output name="output1">\n-\t\t<assert_contents>\n-\t\t\t<has_n_columns n="9" />\n-\t\t</assert_contents>\n-\t    </output>\n+                <assert_contents>\n+                    <has_n_columns n="10" />\n+                </assert_contents>\n+\t        </output>\n             <output name="output2">\n-\t\t<assert_contents>\n-\t\t\t<has_text text="- root - INFO -    errors: 0" />\n-\t\t</assert_contents>\n-\t    </output>\n+\t\t        <assert_contents>\n+\t        \t\t<has_text text="- root - INFO -    errors: 0" />\n+\t        \t</assert_contents>\n+\t        </output>\n+        </test>\n+        <test>\n+            <param name="tracer" value="13C"/>\n+            <section name="db">\n+                <param name="input1" value="Metabolites.dat"/>\n+                <param name="input2" value="Derivatives.dat"/>\n+                <param name="input3" value="Isotopes.dat"/>\n+            </section>\n+            <section name="HR">\n+                <param name="mz" value="400"/>\n+                <param name="resolution" value="60000"/>\n+                <param name="resolution_formula_code" value="orbitrap"/>\n+            </section>\n+            <section name="tcorrection">\n+                <param name="correct_na_tracer" value="-n"/>\n+                <param name="tracer_purity" value="0.01,0.99"/>\n+            </section>\n+            <param name="input4" value="Data_example.tsv"/>\n+            <output name="output1">\n+                <assert_contents>\n+                    <has_n_columns n="10" />\n+                </assert_contents>\n+\t        </output>\n+            <output name="output2">\n+\t\t        <assert_contents>\n+\t        \t\t<has_text text="- root - INFO -    errors: 0" />\n+\t        \t</assert_contents>\n+\t        </output>\n         </test>\n     </tests>\n     <help><![CDATA[\n'
b
diff -r 96f7617b4848 -r fa183805db31 test-data/Metabolites.dat
--- a/test-data/Metabolites.dat Mon May 27 03:33:36 2019 -0400
+++ b/test-data/Metabolites.dat Mon May 04 11:38:09 2020 -0400
b
@@ -1,27 +1,27 @@
-name formula charge
-Suc C4H4O3 -1
-Fum C4H3O4 -1
-Mal C4H5O5 -1
-Cit C6H7O7 -1
-aKG C5H5O5 -1
-OA C4H3O5 -1
-Asp C4H6NO4 -1
-Aco C6HO3 -1
-2/3PG C3H6O3 -1
-G3P C3H7O3 -1
-PEP C3H4O3 -1
-E4P C4H8O4 -1
-R1P C5H8O4 -1
-R5P C5H8O4 -1
-G1P C6H10O5 -1
-G6P C6H10O5 -1
-F1P C6H10O5 -1
-F6P C6H10O5 -1
-M6P C6H10O5 -1
-6PG C6H10O6 -1
-S7P C7H11O6 -1
-FBP C6H11O8P -1
-AMP C10H14N5O4 -1
-ADP C10H14N5O4 -1
-ATP C10H16N5O10 -1
-CDP C9H15N3O7P -1
+name formula charge inchi
+Suc C4H4O3 -1 InChI=1S/C4H6O4/c5-3(6)1-2-4(7)8/h1-2H2,(H,5,6)(H,7,8)/p-2
+Fum C4H3O4 -1 InChI=1S/C4H4O4/c5-3(6)1-2-4(7)8/h1-2H,(H,5,6)(H,7,8)/p-2/b2-1+
+Mal C4H5O5 -1
+Cit C6H7O7 -1
+aKG C5H5O5 -1
+OA C4H3O5 -1
+Asp C4H6NO4 -1
+Aco C6HO3 -1
+2/3PG C3H6O3 -1
+G3P C3H7O3 -1
+PEP C3H4O3 -1
+E4P C4H8O4 -1
+R1P C5H8O4 -1
+R5P C5H8O4 -1
+G1P C6H10O5 -1
+G6P C6H10O5 -1
+F1P C6H10O5 -1
+F6P C6H10O5 -1
+M6P C6H10O5 -1
+6PG C6H10O6 -1
+S7P C7H11O6 -1
+FBP C6H11O8P -1
+AMP C10H14N5O4 -1
+ADP C10H14N5O4 -1
+ATP C10H16N5O10 -1
+CDP C9H15N3O7P -1
b
diff -r 96f7617b4848 -r fa183805db31 test-data/results.log
--- a/test-data/results.log Mon May 27 03:33:36 2019 -0400
+++ b/test-data/results.log Mon May 04 11:38:09 2020 -0400
[
b"@@ -1,58 +1,59 @@\n-2019-02-20 15:46:03 - root - INFO - ------------------------------------------------\n-2019-02-20 15:46:03 - root - INFO - Correction process\n-2019-02-20 15:46:03 - root - INFO - ------------------------------------------------\n-2019-02-20 15:46:03 - root - INFO -    data files\n-2019-02-20 15:46:03 - root - INFO -       data file: Data_example.tsv\n-2019-02-20 15:46:03 - root - INFO -       derivatives database: Derivatives.dat\n-2019-02-20 15:46:03 - root - INFO -       metabolites database: Metabolites.dat\n-2019-02-20 15:46:03 - root - INFO -       isotopes database: Isotopes.dat\n-2019-02-20 15:46:03 - root - INFO -    correction parameters\n-2019-02-20 15:46:03 - root - INFO -       isotopic tracer: 13C\n-2019-02-20 15:46:03 - root - INFO -       correct natural abundance of the tracer element: False\n-2019-02-20 15:46:03 - root - INFO -       isotopic purity of the tracer: None\n-2019-02-20 15:46:03 - root - INFO -       mode: low-resolution\n-2019-02-20 15:46:03 - root - INFO -    natural abundance of isotopes\n-2019-02-20 15:46:03 - root - INFO -    {'C': {'mass': [Decimal('12.0'), Decimal('13.003354835')], 'abundance': [0.9893, 0.0107]}, 'H': {'mass': [Decimal('1.0078250322'), Decimal('2.0141017781')], 'abundance': [0.999885, 0.000115]}, 'N': {'mass': [Decimal('14.003074004'), Decimal('15.000108899')], 'abundance': [0.99636, 0.00364]}, 'O': {'mass': [Decimal('15.99491462'), Decimal('16.999131757'), Decimal('17.999159613')], 'abundance': [0.99757, 0.00038, 0.00205]}, 'P': {'mass': [Decimal('30.973761998')], 'abundance': [1.0]}, 'S': {'mass': [Decimal('31.972071174'), Decimal('32.971458910'), Decimal('33.9678670'), Decimal('35'), Decimal('35.967081')], 'abundance': [0.9499, 0.0075, 0.0425, 0.0, 0.0001]}, 'Si': {'mass': [Decimal('27.976926535'), Decimal('28.976494665'), Decimal('29.9737701')], 'abundance': [0.92223, 0.04685, 0.03092]}}\n-2019-02-20 15:46:03 - root - INFO -    IsoCor version: 2.1.0\n-2019-02-20 15:46:03 - root - INFO - ------------------------------------------------\n-2019-02-20 15:46:03 - root - INFO - Constructing correctors for all (metabolite, derivative)...\n-2019-02-20 15:46:03 - root - INFO - ------------------------------------------------\n-2019-02-20 15:46:03 - root - INFO - ('Fum', '') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('OA', '') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('aKG', '') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('G3P', '') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('2/3PG', '') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('E4P', '') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('Fum', 'TMS') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('OA', 'TMS') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('aKG', 'TMS') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('G3P', 'TMS') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('2/3PG', 'TMS') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ('E4P', 'TMS') successfully constructed.\n-2019-02-20 15:46:03 - root - INFO - ------------------------------------------------\n-2019-02-20 15:46:03 - root - INFO - Correcting raw MS data...\n-2019-02-20 15:46:03 - root - INFO - ------------------------------------------------\n-2019-02-20 15:46:03 - root - INFO - Sample_1 - ('Fum', ''): processed\n-2019-02-20 15:46:03 - root - INFO - Sample_2 - ('Fum', ''): processed\n-2019-02-20 15:46:03 - root - INFO - Sample_1 - ('OA', ''): processed\n-2019-02-20 15:46:03 - root - INFO - Sample_2 - ('OA', ''): processed\n-2019-02-20 15:46:03 - root - INFO - Sample_1 - ('aKG', ''): processed\n-2019-02-20 15:46:03 - root - INFO - Sample_2 - ('aKG', ''): processed\n-2019-02-20 15:46:03 - root - INFO - Sample_1 - ('G3P', ''): processed\n-2019-02-20 15:46:04 - root - INFO - Sample_2 - ('G3P', ''): processed\n-2019-02-20 15:46:04 - root - INFO - Sample_1 - ('2/3PG', ''): pr"..b"': [Decimal('12.0'), Decimal('13.003354835')], 'abundance': [0.9893, 0.0107]}, 'H': {'mass': [Decimal('1.0078250322'), Decimal('2.0141017781')], 'abundance': [0.999885, 0.000115]}, 'N': {'mass': [Decimal('14.003074004'), Decimal('15.000108899')], 'abundance': [0.99636, 0.00364]}, 'O': {'mass': [Decimal('15.99491462'), Decimal('16.999131757'), Decimal('17.999159613')], 'abundance': [0.99757, 0.00038, 0.00205]}, 'P': {'mass': [Decimal('30.973761998')], 'abundance': [1.0]}, 'S': {'mass': [Decimal('31.972071174'), Decimal('32.971458910'), Decimal('33.9678670'), Decimal('35'), Decimal('35.967081')], 'abundance': [0.9499, 0.0075, 0.0425, 0.0, 0.0001]}, 'Si': {'mass': [Decimal('27.976926535'), Decimal('28.976494665'), Decimal('29.9737701')], 'abundance': [0.92223, 0.04685, 0.03092]}}\n+2020-05-04 16:44:12 - root - INFO -    IsoCor version: 2.2.0\n+2020-05-04 16:44:12 - root - INFO - ------------------------------------------------\n+2020-05-04 16:44:12 - root - INFO - Constructing correctors for all (metabolite, derivative)...\n+2020-05-04 16:44:12 - root - INFO - ------------------------------------------------\n+2020-05-04 16:44:12 - root - INFO - ('Fum', '') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('OA', '') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('aKG', '') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('G3P', '') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('2/3PG', '') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('E4P', '') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('Fum', 'TMS') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('OA', 'TMS') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('aKG', 'TMS') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('G3P', 'TMS') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('2/3PG', 'TMS') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ('E4P', 'TMS') successfully constructed.\n+2020-05-04 16:44:12 - root - INFO - ------------------------------------------------\n+2020-05-04 16:44:12 - root - INFO - Correcting raw MS data...\n+2020-05-04 16:44:12 - root - INFO - ------------------------------------------------\n+2020-05-04 16:44:12 - root - INFO - Sample_1 - ('Fum', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_2 - ('Fum', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_1 - ('OA', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_2 - ('OA', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_1 - ('aKG', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_2 - ('aKG', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_1 - ('G3P', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_2 - ('G3P', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_1 - ('2/3PG', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_2 - ('2/3PG', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_1 - ('E4P', ''): processed\n+2020-05-04 16:44:12 - root - INFO - Sample_2 - ('E4P', ''): processed\n+2020-05-04 16:44:13 - root - INFO - Sample_3 - ('Fum', 'TMS'): processed\n+2020-05-04 16:44:13 - root - INFO - Sample_3 - ('OA', 'TMS'): processed\n+2020-05-04 16:44:13 - root - INFO - Sample_3 - ('aKG', 'TMS'): processed\n+2020-05-04 16:44:13 - root - INFO - Sample_3 - ('G3P', 'TMS'): processed\n+2020-05-04 16:44:13 - root - INFO - Sample_3 - ('2/3PG', 'TMS'): processed\n+2020-05-04 16:44:13 - root - INFO - Sample_3 - ('E4P', 'TMS'): processed\n+2020-05-04 16:44:13 - root - INFO - ------------------------------------------------\n+2020-05-04 16:44:13 - root - INFO - Correction process summary\n+2020-05-04 16:44:13 - root - INFO - ------------------------------------------------\n+2020-05-04 16:44:13 - root - INFO -    number of samples: 3\n+2020-05-04 16:44:13 - root - INFO -    number of (metabolite, derivative): 12\n+2020-05-04 16:44:13 - root - INFO -    errors: 0\n"
b
diff -r 96f7617b4848 -r fa183805db31 test-data/results.tsv
--- a/test-data/results.tsv Mon May 27 03:33:36 2019 -0400
+++ b/test-data/results.tsv Mon May 04 11:38:09 2020 -0400
b
b'@@ -1,89 +1,89 @@\n-sample\tmetabolite\tderivative\tisotopologue\tarea\tcorrected_area\tisotopologue_fraction\tresiduum\tmean_enrichment\n-Sample_1\tFum\t\t0\t376000.0\t379808.0494972019\t0.41168665351703776\t0.0\t0.2897659902799015\n-Sample_1\tFum\t\t1\t235000.0\t236670.2675765585\t0.2565348221938963\t0.0\t0.2897659902799015\n-Sample_1\tFum\t\t2\t127000.0\t124721.39583086275\t0.13518969421410323\t1.5800125112233282e-17\t0.2897659902799015\n-Sample_1\tFum\t\t3\t143000.0\t142264.79335165164\t0.15420556980234792\t-3.1600250224466564e-17\t0.2897659902799015\n-Sample_1\tFum\t\t4\t40000.0\t39101.348751418525\t0.04238326027261486\t0.0\t0.2897659902799015\n-Sample_2\tFum\t\t0\t376000.0\t379808.0494972019\t0.41168665351703776\t0.0\t0.2897659902799015\n-Sample_2\tFum\t\t1\t235000.0\t236670.2675765585\t0.2565348221938963\t0.0\t0.2897659902799015\n-Sample_2\tFum\t\t2\t127000.0\t124721.39583086275\t0.13518969421410323\t1.5800125112233282e-17\t0.2897659902799015\n-Sample_2\tFum\t\t3\t143000.0\t142264.79335165164\t0.15420556980234792\t-3.1600250224466564e-17\t0.2897659902799015\n-Sample_2\tFum\t\t4\t40000.0\t39101.348751418525\t0.04238326027261486\t0.0\t0.2897659902799015\n-Sample_1\tOA\t\t0\t0.0\t0.0\t\t\t\n-Sample_1\tOA\t\t1\t0.0\t0.0\t\t\t\n-Sample_1\tOA\t\t2\t0.0\t0.0\t\t\t\n-Sample_1\tOA\t\t3\t0.0\t0.0\t\t\t\n-Sample_1\tOA\t\t4\t0.0\t0.0\t\t\t\n-Sample_2\tOA\t\t0\t0.0\t0.0\t\t\t\n-Sample_2\tOA\t\t1\t0.0\t0.0\t\t\t\n-Sample_2\tOA\t\t2\t0.0\t0.0\t\t\t\n-Sample_2\tOA\t\t3\t0.0\t0.0\t\t\t\n-Sample_2\tOA\t\t4\t0.0\t0.0\t\t\t\n-Sample_1\taKG\t\t0\t761000.0\t770757.2822531819\t0.5971059382536738\t-1.9986617443107176e-07\t0.15239309254478398\n-Sample_1\taKG\t\t1\t95100.0\t94408.21105160222\t0.07313807438059411\t-1.1232997022303625e-07\t0.15239309254478398\n-Sample_1\taKG\t\t2\t391000.0\t387815.2020768239\t0.3004405737538707\t3.24340113523371e-05\t0.15239309254478398\n-Sample_1\taKG\t\t3\t39300.0\t37840.9710506128\t0.029315413611861445\t1.0996123584802932e-05\t0.15239309254478398\n-Sample_1\taKG\t\t4\t0.0\t0.0\t0.0\t-0.003157788575981497\t0.15239309254478398\n-Sample_1\taKG\t\t5\t0.0\t0.0\t0.0\t-0.00030802561931763875\t0.15239309254478398\n-Sample_2\taKG\t\t0\t761000.0\t770757.2822531819\t0.5971059382536738\t-1.9986617443107176e-07\t0.15239309254478398\n-Sample_2\taKG\t\t1\t95100.0\t94408.21105160222\t0.07313807438059411\t-1.1232997022303625e-07\t0.15239309254478398\n-Sample_2\taKG\t\t2\t391000.0\t387815.2020768239\t0.3004405737538707\t3.24340113523371e-05\t0.15239309254478398\n-Sample_2\taKG\t\t3\t39300.0\t37840.9710506128\t0.029315413611861445\t1.0996123584802932e-05\t0.15239309254478398\n-Sample_2\taKG\t\t4\t0.0\t0.0\t0.0\t-0.003157788575981497\t0.15239309254478398\n-Sample_2\taKG\t\t5\t0.0\t0.0\t0.0\t-0.00030802561931763875\t0.15239309254478398\n-Sample_1\tG3P\t\t0\t715000.0\t720817.843203928\t0.5037179457196655\t-8.142070347386684e-17\t0.2307972956362909\n-Sample_1\tG3P\t\t1\t557000.0\t560128.1627924398\t0.3914256703849595\t8.142070347386684e-17\t0.2307972956362909\n-Sample_1\tG3P\t\t2\t24800.0\t19465.731702013378\t0.013602935162211748\t-5.088793967116678e-17\t0.2307972956362909\n-Sample_1\tG3P\t\t3\t133000.0\t130583.22551280692\t0.09125344873316321\t1.2213105521080026e-16\t0.2307972956362909\n-Sample_2\tG3P\t\t0\t715000.0\t720817.843203928\t0.5037179457196655\t-8.142070347386684e-17\t0.2307972956362909\n-Sample_2\tG3P\t\t1\t557000.0\t560128.1627924398\t0.3914256703849595\t8.142070347386684e-17\t0.2307972956362909\n-Sample_2\tG3P\t\t2\t24800.0\t19465.731702013378\t0.013602935162211748\t-5.088793967116678e-17\t0.2307972956362909\n-Sample_2\tG3P\t\t3\t133000.0\t130583.22551280692\t0.09125344873316321\t1.2213105521080026e-16\t0.2307972956362909\n-Sample_1\t2/3PG\t\t0\t1430000.0\t1441469.8983039188\t0.26951520374142846\t0.0\t0.350613297990272\n-Sample_1\t2/3PG\t\t1\t2920000.0\t2940779.0438875123\t0.5498447550687304\t0.0\t0.350613297990272\n-Sample_1\t2/3PG\t\t2\t226000.0\t213534.01511096465\t0.03992498466743767\t1.0898270157923124e-17\t0.350613297990272\n-Sample_1\t2/3PG\t\t3\t765000.0\t752597.6842841837\t0.1407150565224034\t2.1796540315846248e-17\t0.350613297990272\n-Sample_2\t2/3PG\t\t0\t1430000.0\t1441469.8983039188\t0.26951520374142846\t0.0\t0.350613297990272\n-Sample_2\t2/3PG\t\t1\t2920000.0\t2940779.0438875123\t0.5498447550687304\t0.0\t0.350613297990272\n-Sample_2\t2/3PG\t\t2\t226000.0\t213534.01511096465\t0.03992498466743767\t1.0898270157923124e-17\t0.3506132'..b'+Sample_1\tE4P\t\t1\t/a(C1+1),(C3+0)\t10200.0\t10307.12055255345\t0.4936185983859252\t-1.0461978931437058e-14\t0.3817239095188967\n+Sample_1\tE4P\t\t2\t/a(C2+1),(C2+0)\t7760.0\t7810.814463659301\t0.37406793373042024\t1.508268629282176e-14\t0.3817239095188967\n+Sample_1\tE4P\t\t3\t/a(C3+1),(C1+0)\t1730.0\t1644.6719684418108\t0.0787650311963941\t-2.3299698911887947e-14\t0.3817239095188967\n+Sample_1\tE4P\t\t4\t/a(C4+1)\t320.0\t254.98554081254488\t0.012211519659909676\t3.423846073733581e-14\t0.3817239095188967\n+Sample_2\tE4P\t\t0\t/a(C4+0)\t854.0\t863.1453281237473\t0.041336917027350754\t2.4329549962951284e-14\t0.3817239095188967\n+Sample_2\tE4P\t\t1\t/a(C1+1),(C3+0)\t10200.0\t10307.12055255345\t0.4936185983859252\t-1.0461978931437058e-14\t0.3817239095188967\n+Sample_2\tE4P\t\t2\t/a(C2+1),(C2+0)\t7760.0\t7810.814463659301\t0.37406793373042024\t1.508268629282176e-14\t0.3817239095188967\n+Sample_2\tE4P\t\t3\t/a(C3+1),(C1+0)\t1730.0\t1644.6719684418108\t0.0787650311963941\t-2.3299698911887947e-14\t0.3817239095188967\n+Sample_2\tE4P\t\t4\t/a(C4+1)\t320.0\t254.98554081254488\t0.012211519659909676\t3.423846073733581e-14\t0.3817239095188967\n+Sample_3\tFum\tTMS\t0\t/a(C4+0)\t376000.0\t425785.0538505561\t0.45766573056460497\t0.0\t0.2626278340043467\n+Sample_3\tFum\tTMS\t1\t/a(C1+1),(C3+0)\t235000.0\t229433.5002431335\t0.24661234478564592\t-3.1600250224466564e-17\t0.2626278340043467\n+Sample_3\tFum\tTMS\t2\t/a(C2+1),(C2+0)\t127000.0\t105318.33031467753\t0.11320404544361702\t-3.1600250224466564e-17\t0.2626278340043467\n+Sample_3\tFum\tTMS\t3\t/a(C3+1),(C1+0)\t143000.0\t141951.9567793512\t0.15258061648002147\t3.1600250224466564e-17\t0.2626278340043467\n+Sample_3\tFum\tTMS\t4\t/a(C4+1)\t40000.0\t27851.853810968023\t0.029937262726110622\t0.0\t0.2626278340043467\n+Sample_3\tOA\tTMS\t0\t/a(C4+0)\t0.0\t0.0\t\t\t\n+Sample_3\tOA\tTMS\t1\t/a(C1+1),(C3+0)\t0.0\t0.0\t\t\t\n+Sample_3\tOA\tTMS\t2\t/a(C2+1),(C2+0)\t0.0\t0.0\t\t\t\n+Sample_3\tOA\tTMS\t3\t/a(C3+1),(C1+0)\t0.0\t0.0\t\t\t\n+Sample_3\tOA\tTMS\t4\t/a(C4+1)\t0.0\t0.0\t\t\t\n+Sample_3\taKG\tTMS\t0\t/a(C5+0)\t761000.0\t864080.4383954712\t0.6634102607174354\t-1.4215099326884681e-05\t0.1302868347875265\n+Sample_3\taKG\tTMS\t1\t/a(C1+1),(C4+0)\t95100.0\t33111.70372452174\t0.02542198969516514\t-7.027619573886684e-05\t0.1302868347875265\n+Sample_3\taKG\tTMS\t2\t/a(C2+1),(C3+0)\t391000.0\t400501.8154913149\t0.3074910645197312\t0.0005087109208496222\t0.1302868347875265\n+Sample_3\taKG\tTMS\t3\t/a(C3+1),(C2+0)\t39300.0\t4788.818975572151\t0.003676685067668318\t0.0011821047313132886\t0.1302868347875265\n+Sample_3\taKG\tTMS\t4\t/a(C4+1),(C1+0)\t0.0\t0.0\t0.0\t-0.013225164964720594\t0.1302868347875265\n+Sample_3\taKG\tTMS\t5\t/a(C5+1)\t0.0\t0.0\t0.0\t-0.000751685958158175\t0.1302868347875265\n+Sample_3\tG3P\tTMS\t0\t/a(C3+0)\t715000.0\t806222.8908729891\t0.5421097940219233\t0.0011462823123024992\t0.20893939381005114\n+Sample_3\tG3P\tTMS\t1\t/a(C1+1),(C2+0)\t557000.0\t555357.5466204376\t0.37342621825203837\t0.002864378095056404\t0.20893939381005114\n+Sample_3\tG3P\tTMS\t2\t/a(C2+1),(C1+0)\t24800.0\t0.0\t0.0\t-0.03321751751634703\t0.20893939381005114\n+Sample_3\tG3P\tTMS\t3\t/a(C3+1)\t133000.0\t125614.40709996362\t0.08446398772603834\t-2.6150701541736555e-12\t0.20893939381005114\n+Sample_3\t2/3PG\tTMS\t0\t/a(C3+0)\t1430000.0\t1613070.9343452677\t0.2934931923774996\t0.0004794245823734611\t0.3239815179271314\n+Sample_3\t2/3PG\tTMS\t1\t/a(C1+1),(C2+0)\t2920000.0\t3153601.731582657\t0.5737879345430535\t0.0011960645935207325\t0.3239815179271314\n+Sample_3\t2/3PG\tTMS\t2\t/a(C2+1),(C1+0)\t226000.0\t0.0\t0.0\t-0.01388900414166337\t0.3239815179271314\n+Sample_3\t2/3PG\tTMS\t3\t/a(C3+1)\t765000.0\t729437.5548178032\t0.1327188730794469\t1.6153633993476813e-12\t0.3239815179271314\n+Sample_3\tE4P\tTMS\t0\t/a(C4+0)\t854.0\t967.700772009107\t0.046150526549923535\t-2.9138586415226663e-06\t0.3493645351728113\n+Sample_3\tE4P\tTMS\t1\t/a(C1+1),(C3+0)\t10200.0\t11473.616488166983\t0.5471871653687257\t-1.4649378122356719e-05\t0.3493645351728113\n+Sample_3\tE4P\tTMS\t2\t/a(C2+1),(C2+0)\t7760.0\t7752.336450775334\t0.369715948921533\t0.00010602215694356356\t0.3493645351728113\n+Sample_3\tE4P\tTMS\t3\t/a(C3+1),(C1+0)\t1730.0\t774.7044932023791\t0.036946359159817824\t0.000251776357076686\t0.3493645351728113\n+Sample_3\tE4P\tTMS\t4\t/a(C4+1)\t320.0\t0.0\t0.0\t-0.002903094538807673\t0.3493645351728113\n \n'