Repository 'nanopolishcomp_freqmethcalculate'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/nanopolishcomp_freqmethcalculate

Changeset 0:80b6e23cb70d (2020-05-12)
Next changeset 1:f08dfdcefdbd (2024-03-13)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/ commit cc7e6263fdb9e59f992aee41fc9d62425f39bf38"
added:
freqmethcalculate.xml
macros.xml
test-data/sample_event.tsv
test-data/sample_meth.tsv
b
diff -r 000000000000 -r 80b6e23cb70d freqmethcalculate.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/freqmethcalculate.xml Tue May 12 17:08:21 2020 -0400
[
@@ -0,0 +1,129 @@
+<?xml version="1.0"?>
+<tool id="nanopolishcomp_freqmethcalculate" name="NanopolishComp: FreqMethCalculate" version="@TOOL_VERSION@+@WRAPPER_VERSION@">
+    <description>calculates methylation frequency at genomic CpG sites</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements"/>
+    <expand macro="version_command"/>
+    <command detect_errors="exit_code"><![CDATA[
+        NanopolishComp Freq_meth_calculate
+            ## required
+            -i '$i'
+            ## optional
+            -d $d
+            #if $f
+                -f '$f'
+            #end if
+            #if $s
+                -s '$s'
+            #end if
+            $strand_specific
+            --min_llr $min_llr
+            -b 'result.bed'
+            -t 'result.tsv'
+            -v
+        ]]></command>
+    <inputs>
+        <param argument="-i" type="data" format="tabular" label="Select nanopolish call methylation file"/>
+        <param argument="-d" type="integer" value="10" min="0" label="Set minimal number of reads covering a site to be reported"/>
+        <param argument="-f" type="data" format="fai" optional="true" label="Select fasta index file" help="Obtained with samtools faidx. Required for coordinate sorting"/>
+        <param argument="-s" type="text" value="" optional="true" label="Set sample ID" help="Used for bed track header."/>
+        <param argument="--strand_specific" type="boolean" truevalue="--strand_specific" falsevalue="" label="Are sites strand specific?"/>
+        <param argument="--min_llr" type="integer" value="2" min="0" label="Set minimal log likelyhood ratio" help="To consider a site significantly methylated or unmethylated."/>
+        <param name="out" type="select" multiple="true" optional="false" label="Select output file(s)">
+            <option value="bed" selected="true">Methylation frequency (BED)</option>
+            <option value="tsv" selected="true">Methylation frequency (Tabular)</option>
+        </param>
+    </inputs>
+    <outputs>
+        <data name="out_bed" format="bed" from_work_dir="result.bed" label="${tool.name} on ${on_string}: Methylation frequency (BED)">
+            <filter>'bed' in out</filter>
+        </data>
+        <data name="out_tsv" format="tabular" from_work_dir="result.tsv" label="${tool.name} on ${on_string}: Methylation frequency (Tabular)">
+            <filter>'tsv' in out</filter>
+        </data>
+    </outputs>
+    <tests>
+        <!-- #1 default -->
+        <test expect_num_outputs="2">
+            <param name="i" value="sample_meth.tsv"/>
+            <output name="out_bed">
+                <assert_contents>
+                    <has_n_lines n="56"/>
+                    <has_line line="track name=Methylation_ itemRgb=On"/>
+                    <has_text_matching expression="tig00000001&#009;201961.+"/>
+                </assert_contents>
+            </output>
+            <output name="out_tsv">
+                <assert_contents>
+                    <has_n_lines n="56"/>
+                    <has_text_matching expression="chromosome&#009;start.+"/>
+                    <has_text_matching expression="tig00000001&#009;201961.+"/>
+                </assert_contents>
+            </output>
+        </test>
+        <!-- #2 -->
+        <test expect_num_outputs="2">
+            <param name="i" value="sample_meth.tsv"/>
+            <param name="d" value="11"/>
+            <param name="s" value="test"/>
+            <param name="strand_specific" value="true"/>
+            <param name="min_llr" value="3"/>
+            <param name="out" value="bed,tsv"/>
+            <output name="out_bed">
+                <assert_contents>
+                    <has_n_lines n="111"/>
+                    <has_line line="track name=Methylation_test itemRgb=On"/>
+                    <has_text_matching expression="tig00000001&#009;200009.+"/>
+                </assert_contents>
+            </output>
+            <output name="out_tsv">
+                <assert_contents>
+                    <has_n_lines n="111"/>
+                    <has_text_matching expression="chromosome&#009;start.+"/>
+                    <has_text_matching expression="tig00000001&#009;201961.+"/>
+                </assert_contents>
+            </output>
+        </test>
+    </tests>
+    <help><![CDATA[
+.. class:: infomark
+
+**What it does**
+
+@WID@
+
+Freq_meth_calculate calculates methylation frequency at genomic CpG sites from the output of nanopolish call-methylation.
+
+**Input**
+
+A nanopolish call methylation tabular output file.
+
+**Output**
+
+*BED file*
+
+Standard genomic BED6. The score correspond to the methylation frequency multiplied by 1000. The file is sorted by coordinates and can be rendered with a genome browser such as IGV
+
+*Tabulated TSV file*
+
+Contrary to the bed file, in the tabulated report, positions are ordered by decreasing methylation frequency.
+
+The file contains the following fields:
+
+- chrom / start / end / strand: Genomic coordinates of the motif or group of motifs in case split_group was not selected.
+- site_id: Unique integer identifier of the genomic position.
+- methylated_reads / unmethylated_reads / ambiguous_reads: Number of reads at a given genomic location with a higher likelyhood of being methylated or unmethylated or with an ambiguous methylation call.
+- sequence: -5 to +5 sequence of the motif or group of motifs in case split_group was not selected.
+- num_motifs: Number of motif in the group.
+- meth_freq: Methylation frequency (out of non anbiguous calls).
+
+.. class:: infomark
+
+**References**
+
+@REFERENCES@
+    ]]></help>
+    <expand macro="citations"/>
+</tool>
\ No newline at end of file
b
diff -r 000000000000 -r 80b6e23cb70d macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml Tue May 12 17:08:21 2020 -0400
[
@@ -0,0 +1,30 @@
+<?xml version="1.0"?>
+<macros>
+    <token name="@TOOL_VERSION@">0.6.11</token>
+    <token name="@WRAPPER_VERSION@">galaxy1</token>
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@TOOL_VERSION@">nanopolishcomp</requirement>
+        </requirements>
+    </xml>
+    <xml name="version_command">
+        <version_command>NanopolishComp --version</version_command>
+    </xml>
+   <xml name="citations">
+        <citations>
+            <citation type="bibtex">@online{NanopolishComp,
+              author = {Adrien Leger },
+              title = {NanopolishComp},
+              year = 2019,
+              url = {https://a-slide.github.io/NanopolishComp/},
+              urldate = {2020-04-29}
+            }</citation>
+        </citations>
+    </xml>
+    <token name="@WID@"><![CDATA[
+NanopolishComp is a companion package for downstream analyses of Nanopolish output files.
+    ]]></token>
+    <token name="@REFERENCES@"><![CDATA[
+More information are available on the `project website <https://a-slide.github.io/NanopolishComp/>`_ and `github <https://github.com/a-slide/NanopolishComp>`_.
+    ]]></token>
+</macros>
\ No newline at end of file
b
diff -r 000000000000 -r 80b6e23cb70d test-data/sample_event.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/sample_event.tsv Tue May 12 17:08:21 2020 -0400
b
b'@@ -0,0 +1,271 @@\n+contig\tposition\treference_kmer\tread_name\tstrand\tevent_index\tevent_level_mean\tevent_stdv\tevent_length\tmodel_kmer\tmodel_mean\tmodel_stdv\tstandardized_level\tsamples\n+chr\t22101\tAGGAA\tf30a0489-ec1e-4f6e-829e-fa6d16dbee2e\tt\t8291\t124.19\t3.490\t0.00664\tAGGAA\t117.46\t3.17\t1.88\t122.692,122.834,120.843,120.558,119.705,122.834,129.235,120.132,126.248,122.977,123.83,127.813,124.541,122.55,122.265,129.947,119.42,133.076,127.528,124.684\n+chr\t22102\tGGAAA\tf30a0489-ec1e-4f6e-829e-fa6d16dbee2e\tt\t8292\t130.37\t3.206\t0.00365\tGGAAA\t121.47\t5.56\t1.42\t131.369,132.08,134.641,127.955,125.822,133.076,132.08,124.257,126.96,132.08,133.787\n+chr\t22102\tGGAAA\tf30a0489-ec1e-4f6e-829e-fa6d16dbee2e\tt\t8293\t121.86\t9.392\t0.00199\tGGAAA\t121.47\t5.56\t0.06\t124.684,132.223,130.8,107.33,126.391,109.748\n+chr\t22101\tAGGAA\te924b30f-22dc-4156-8e5c-ba664abbe7a7\tt\t2723\t118.83\t5.213\t0.00730\tAGGAA\t117.46\t3.17\t0.37\t108.391,113.503,114.639,119.893,121.881,106.687,119.325,116.627,128.555,120.035,120.177,114.781,127.561,122.733,115.349,118.899,118.189,119.609,126.141,116.059,119.609,125.715\n+chr\t22102\tGGAAA\te924b30f-22dc-4156-8e5c-ba664abbe7a7\tt\t2724\t122.61\t5.634\t0.00863\tGGAAA\t121.47\t5.56\t0.17\t122.165,117.337,122.449,111.657,118.331,112.935,119.609,122.307,121.313,142.329,123.869,121.313,120.603,121.455,120.461,134.377,118.473,123.869,122.591,121.739,120.035,124.153,124.863,126.141,127.277,126.141\n+chr\t22102\tGGAAA\te924b30f-22dc-4156-8e5c-ba664abbe7a7\tt\t2725\t124.78\t6.096\t0.01627\tGGAAA\t121.47\t5.56\t0.50\t121.171,121.029,114.781,130.259,129.265,116.485,126.283,110.095,130.827,132.815,128.839,115.065,127.703,131.963,133.667,105.835,132.247,121.881,120.461,126.283,131.395,121.029,110.095,124.721,133.809,131.537,127.987,125.715,120.035,130.685,124.579,123.869,127.419,130.827,131.111,124.721,130.543,120.887,124.863,128.555,126.425,120.461,128.129,127.703,121.455,120.603,125.005,122.165,120.887\n+chr\t22101\tAGGAA\td4988330-6f34-498a-b1e3-6574d8b3dae2\tt\t9924\t115.24\t5.197\t0.01228\tAGGAA\t117.46\t3.17\t-0.59\t104.751,123.941,105.044,118.961,113.248,109.293,114.42,120.865,117.935,108.707,115.005,119.84,118.521,107.974,112.515,109.732,108.853,122.183,106.949,117.349,114.566,108.56,117.789,115.445,119.4,123.209,119.547,114.859,110.318,113.541,123.209,122.037,122.916,108.853,113.834,121.597,118.228\n+chr\t22101\tAGGAA\td4988330-6f34-498a-b1e3-6574d8b3dae2\tt\t9925\t122.43\t2.734\t0.00797\tAGGAA\t117.46\t3.17\t1.32\t119.547,125.406,119.693,121.451,126.871,121.304,124.527,123.502,125.553,117.642,122.916,124.967,124.527,118.228,120.718,120.718,120.865,127.017,127.017,126.138,118.375,120.426,121.011,119.986\n+chr\t22102\tGGAAA\td4988330-6f34-498a-b1e3-6574d8b3dae2\tt\t9926\t125.96\t8.542\t0.00531\tGGAAA\t121.47\t5.56\t0.68\t118.082,116.031,119.693,130.387,119.547,135.66,136.393,111.636,107.242,123.648,133.463,131.705,134.488,134.928,134.342,128.189\n+chr\t22101\tAGGAA\t43000cf1-c697-4e0b-8414-5aed3db85976\tt\t3971\t115.45\t5.142\t0.00863\tAGGAA\t117.46\t3.17\t-0.58\t112.784,113.078,117.2,113.52,118.82,112.342,123.826,104.097,120.439,113.078,122.353,104.981,115.139,115.286,116.611,112.195,106.011,114.109,113.814,117.348,125.592,115.875,112.784,115.728,117.053,127.653\n+chr\t22102\tGGAAA\t43000cf1-c697-4e0b-8414-5aed3db85976\tt\t3972\t129.34\t5.262\t0.00232\tGGAAA\t121.47\t5.56\t1.30\t124.856,124.12,125.151,130.892,141.787,128.095,130.451\n+chr\t22102\tGGAAA\t43000cf1-c697-4e0b-8414-5aed3db85976\tt\t3973\t122.81\t3.290\t0.01527\tGGAAA\t121.47\t5.56\t0.22\t117.642,131.923,125.739,119.998,117.495,116.906,118.084,118.673,122.353,124.709,123.826,119.114,123.678,124.856,119.114,123.973,120.439,121.323,116.464,128.39,117.936,118.673,122.942,123.531,127.065,125.298,118.673,124.414,125.592,125.887,126.917,124.414,121.617,119.703,130.009,124.709,125.592,125.003,121.47,120.881,126.034,124.709,120.881,123.531,124.856,124.414\n+chr\t22102\tGGAAA\t43000cf1-c697-4e0b-8414-5aed3db85976\tt\t3974\t133.48\t2.700\t0.00232\tGGAAA\t121.47\t5.56\t1.98\t136.045,133.542,130.156,136.487,135.604,134.279,128.242\n+chr\t22102\tGGAAA\t43000cf1-c697-4e0b-8414-5aed3db85976\tt\t3975\t117.60\t8.898\t0.00730\tGGAAA\t121.47\t5.56'..b'8.144,121.167,111.26,105.679,103.586,112.097,117.26,125.632,116.004,116.423,124.655,113.213,123.26,116.283,116.144,116.283,114.609,116.562,118.376,127.865,115.027,122.283,100.655,121.027,116.702,117.26,111.539,100.516,107.911,110.702,116.423,109.446,102.888,105.539,107.353,110.981,113.911,105.26,110.702,116.283,110.702,110.144,105.26,113.493,113.632,111.818,115.307,113.213,108.748,110.283,123.4,110.981,112.237,121.307,121.725,113.911,113.911,113.632,105.539,115.586,113.632,104.981,113.353,122.562,120.469,119.493,108.33,113.074,118.097,110.004,110.702,106.655,111.818,122.702,106.097,110.841,109.586,116.283,114.748,106.237,108.19,104.144,106.237\n+chr\t22102\tGGAAA\t6057e4db-745d-407d-98bb-83a8e9a5f0ce\tt\t62\t122.77\t4.401\t0.01726\tGGAAA\t121.47\t5.56\t0.22\t116.562,119.772,117.4,114.609,118.237,121.865,123.26,118.655,124.934,116.702,119.353,127.865,124.795,122.562,121.307,127.027,118.655,124.655,141.121,127.307,120.469,120.609,124.655,121.167,133.446,122.702,121.865,124.934,119.353,126.19,125.493,128.144,125.772,122.841,123.539,118.934,121.586,122.841,123.4,127.446,123.958,123.539,123.26,122.841,117.818,115.027,124.237,124.934,122.004,122.562,116.283,125.632\n+chr\t22102\tGGAAA\t6057e4db-745d-407d-98bb-83a8e9a5f0ce\tt\t63\t123.47\t5.883\t0.02424\tGGAAA\t121.47\t5.56\t0.34\t133.446,128.841,129.121,129.4,131.214,133.307,135.539,133.446,132.33,112.376,129.818,128.841,105.958,133.028,129.958,125.911,118.516,131.214,117.818,126.19,129.121,129.26,115.167,122.981,124.376,118.934,119.353,123.539,118.516,120.748,129.679,124.237,122.841,114.469,115.307,121.586,121.027,131.632,126.19,125.353,114.609,120.888,125.353,134.702,123.12,117.539,119.214,118.655,120.19,124.655,121.446,123.12,128.144,118.516,123.12,116.144,120.888,119.911,120.609,123.12,124.516,114.051,120.19,121.167,115.307,120.33,124.516,125.353,121.586,128.004,115.586,122.144,126.19\n+chr\t22102\tGGAAA\t6057e4db-745d-407d-98bb-83a8e9a5f0ce\tt\t64\t119.90\t8.320\t0.00332\tGGAAA\t121.47\t5.56\t-0.27\t122.981,119.632,105.958,116.144,127.027,128.562,117.818,131.632,105.12,124.097\n+chr\t22102\tGGAAA\t6057e4db-745d-407d-98bb-83a8e9a5f0ce\tt\t65\t108.37\t1.812\t0.01029\tGGAAA\t121.47\t5.56\t-2.25\t104.702,109.167,110.841,109.027,110.004,110.562,109.586,105.818,107.911,109.027,106.655,106.516,109.725,108.748,112.237,106.934,110.562,108.33,111.26,109.167,106.655,106.795,110.144,108.748,105.818,109.865,108.051,107.074,107.213,106.097,106.237\n+chr\t22102\tGGAAA\t6057e4db-745d-407d-98bb-83a8e9a5f0ce\tt\t66\t120.71\t6.942\t0.00232\tGGAAA\t121.47\t5.56\t-0.13\t106.655,121.725,128.144,115.586,119.632,127.307,125.911\n+chr\t22101\tAGGAA\t6b779e7d-26fb-4863-9082-04068823a29d\tt\t11\t125.01\t4.330\t0.00863\tAGGAA\t117.46\t3.17\t1.96\t118.947,115.47,134.526,132.44,123.398,123.815,125.067,121.868,124.65,124.372,125.624,126.737,121.59,126.876,131.466,119.921,121.173,129.797,127.71,122.285,127.015,129.936,124.928,118.252,127.154,125.346\n+chr\t22102\tGGAAA\t6b779e7d-26fb-4863-9082-04068823a29d\tt\t12\t133.43\t2.914\t0.00232\tGGAAA\t121.47\t5.56\t1.77\t132.718,132.579,126.876,135.778,135.639,134.665,135.778\n+chr\t22101\tAGGAA\t484f63e8-a501-4edd-8c0e-fba261ed96e1\tt\t32\t118.35\t4.651\t0.00764\tAGGAA\t117.46\t3.17\t0.25\t115.392,107.698,115.064,123.74,123.249,106.88,115.555,113.264,122.103,120.957,121.776,108.844,122.267,118.502,116.538,120.302,115.392,114.41,126.686,127.505,121.939,122.267,121.776\n+chr\t22102\tGGAAA\t484f63e8-a501-4edd-8c0e-fba261ed96e1\tt\t33\t126.47\t5.067\t0.00697\tGGAAA\t121.47\t5.56\t0.80\t127.505,122.267,128.65,142.891,136.671,128.323,126.359,121.776,123.576,123.085,125.049,128.814,129.305,113.919,126.195,134.543,126.85,125.213,118.502,124.886,121.448\n+chr\t22102\tGGAAA\t484f63e8-a501-4edd-8c0e-fba261ed96e1\tt\t34\t127.41\t5.048\t0.00398\tGGAAA\t121.47\t5.56\t0.95\t135.689,124.886,135.853,133.888,135.525,117.52,123.249,121.121,123.249,127.341,126.686,123.903\n+chr\t22102\tGGAAA\t55fb9fbe-c25a-486b-9fed-b50328c70e5a\tt\t31\t127.10\t7.175\t0.00631\tGGAAA\t121.47\t5.56\t0.89\t127.165,133.781,119.847,105.491,129.277,125.477,131.247,117.736,131.81,133.781,137.299,118.721,131.388,136.173,127.447,126.884,128.01,129.699,123.647\n'
b
diff -r 000000000000 -r 80b6e23cb70d test-data/sample_meth.tsv
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/sample_meth.tsv Tue May 12 17:08:21 2020 -0400
b
b'@@ -0,0 +1,3744 @@\n+chromosome\tstrand\tstart\tend\tread_name\tlog_lik_ratio\tlog_lik_methylated\tlog_lik_unmethylated\tnum_calling_strands\tnum_motifs\tsequence\n+tig00000001\t+\t200009\t200009\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t3.49\t-126.31\t-129.80\t1\t1\tTTCAACGGTAT\n+tig00000001\t+\t200046\t200046\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-1.33\t-98.68\t-97.34\t1\t1\tTGCAGCGCACC\n+tig00000001\t+\t200081\t200081\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-2.51\t-109.65\t-107.14\t1\t1\tAACAACGCCAC\n+tig00000001\t+\t200105\t200122\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-10.74\t-181.29\t-170.55\t1\t3\tCATAACGTGATGCGTAGCAGATCGACCC\n+tig00000001\t+\t200137\t200141\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-11.21\t-164.82\t-153.61\t1\t2\tATTCCCGAGCGTGCT\n+tig00000001\t+\t200157\t200169\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-12.69\t-194.84\t-182.15\t1\t3\tTCTGGCGGGTAAACGCACGGATG\n+tig00000001\t+\t200188\t200202\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-14.53\t-185.43\t-170.90\t1\t3\tATCACCGCAATCCCGCTGGCGGCAG\n+tig00000001\t+\t200224\t200233\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-7.43\t-152.60\t-145.17\t1\t2\tTGCACCGTTAAGCCCGCCAT\n+tig00000001\t+\t200257\t200257\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t0.69\t-99.34\t-100.03\t1\t1\tCTCAACGTGGT\n+tig00000001\t+\t200268\t200268\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-1.03\t-90.55\t-89.52\t1\t1\tTTCACCGCTAT\n+tig00000001\t+\t200307\t200320\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-5.24\t-197.00\t-191.76\t1\t3\tGCTGACGGCGGTATTCAGCGCCAG\n+tig00000001\t+\t200332\t200332\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-0.27\t-94.01\t-93.73\t1\t1\tTCCCACGACAG\n+tig00000001\t+\t200362\t200362\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t1.76\t-143.76\t-145.52\t1\t1\tTTAAGCGGCAG\n+tig00000001\t+\t200379\t200379\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t0.40\t-95.40\t-95.81\t1\t1\tCTGACCGAGCA\n+tig00000001\t+\t200400\t200400\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-5.94\t-109.57\t-103.63\t1\t1\tAGCACCGCCAG\n+tig00000001\t+\t200423\t200423\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-2.93\t-94.91\t-91.98\t1\t1\tCAGGCCGAGAA\n+tig00000001\t+\t200461\t200512\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-26.93\t-416.30\t-389.37\t1\t9\tTCTCACGGTCAGAGACGCCAAGTGCGCGAAAAGTACGCGCTCAACGCCCGTTGTACCGGGAA\n+tig00000001\t+\t200538\t200548\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-7.98\t-147.26\t-139.28\t1\t3\tTCAGCCGCGCCAGCCCGCTGC\n+tig00000001\t+\t200562\t200564\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-1.73\t-90.91\t-89.18\t1\t2\tAAGGACGCGCCAG\n+tig00000001\t+\t200591\t200596\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-6.85\t-149.23\t-142.38\t1\t2\tAAAAACGTGGCGATCA\n+tig00000001\t+\t200615\t200615\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-3.90\t-113.58\t-109.67\t1\t1\tCCTTGCGCAGC\n+tig00000001\t+\t200629\t200635\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-3.60\t-114.95\t-111.35\t1\t2\tCAGAACGCCTCCGCATT\n+tig00000001\t+\t200687\t200699\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-9.56\t-173.23\t-163.67\t1\t3\tCAATACGCCGGTTATCACGCCTG\n+tig00000001\t+\t200742\t200742\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-4.67\t-107.20\t-102.53\t1\t1\tGATACCGGAAG\n+tig00000001\t+\t200773\t200773\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-4.93\t-124.32\t-119.38\t1\t1\tGATTTCGCAAA\n+tig00000001\t+\t200784\t200788\td57afb7d-903e-46cf-a43d-0e17fb0949d8\t-9.11\t-148.71\t-139.60\t1\t2\tATCTGCGGGCGGGGT\n+tig00000001\t+\t200009\t200009\t3221ccdb-0395-4983-8a44-e5558664096f\t-6.47\t-122.75\t-116.27\t1\t1\tTTCAACGGTAT\n+tig00000001\t+\t200046\t200046\t3221ccdb-0395-4983-8a44-e5558664096f\t-1.88\t-87.76\t-85.88\t1\t1\tTGCAGCGCACC\n+tig00000001\t+\t200081\t200081\t3221ccdb-0395-4983-8a44-e5558664096f\t-4.45\t-106.13\t-101.68\t1\t1\tAACAACGCCAC\n+tig00000001\t+\t200105\t200122\t3221ccdb-0395-4983-8a44-e5558664096f\t-10.99\t-256.60\t-245.61\t1\t3\tCATAACGTGATGCGTAGCAGATCGACCC\n+tig00000001\t+\t200137\t200141\t3221ccdb-0395-4983-8a44-e5558664096f\t1.82\t-117.34\t-119.15\t1\t2\tATTCCCGAGCGTGCT\n+tig00000001\t+\t200157\t200169\t3221ccdb-0395-4983-8a44-e5558664096f\t-23.60\t-229.04\t-205.43\t1\t3\tTCTGGCGGGTAAACGCACGGATG\n+tig00000001\t+\t200188\t200202\t3221ccdb-0395-4983-8a44-e5558664096f\t-10.53\t-201.67\t-191.14\t1\t3\tATCACCGCAATCCCGCTGGCGGCAG\n+tig00000001\t+\t200224\t200233\t3221ccdb-0395-4983-8a44-e5558664096f\t3.20\t-128.67\t-131.87\t1\t2\tTGCACCGTTAAGCCCGCCAT\n+tig00000001\t+\t200257\t200257\t3221ccdb-0395-4983-8a44-e5558664096f\t-2.35\t-96.33\t-93.98\t1\t1\tCTCAACGTGGT\n+tig000'..b'b04a38b5225\t0.62\t-143.73\t-144.35\t1\t2\tAGCCTCGGCCCCGGCTG\n+tig00000001\t-\t201761\t201769\tf432cb50-8c47-4c54-825f-ab04a38b5225\t-4.21\t-198.50\t-194.29\t1\t2\tGCCTGCGGATATCCGCTTA\n+tig00000001\t-\t201781\t201783\tf432cb50-8c47-4c54-825f-ab04a38b5225\t0.33\t-137.63\t-137.96\t1\t2\tAGCTGCGCGATAA\n+tig00000001\t-\t201810\t201813\tf432cb50-8c47-4c54-825f-ab04a38b5225\t-5.40\t-135.04\t-129.65\t1\t2\tGATAACGGCGGCAG\n+tig00000001\t-\t201846\t201870\tf432cb50-8c47-4c54-825f-ab04a38b5225\t0.07\t-247.07\t-247.14\t1\t4\tTTTCCCGGTGAGGACGGTTACAGCCGCAGCGAGTC\n+tig00000001\t-\t201889\t201902\tf432cb50-8c47-4c54-825f-ab04a38b5225\t-5.33\t-246.68\t-241.35\t1\t5\tTGGTGCGCGGCGGCGTGGCGAAAC\n+tig00000001\t-\t201922\t201930\tf432cb50-8c47-4c54-825f-ab04a38b5225\t-0.64\t-151.28\t-150.64\t1\t2\tGTCACCGGCTGGCCGCACT\n+tig00000001\t-\t201944\t201950\tf432cb50-8c47-4c54-825f-ab04a38b5225\t0.57\t-132.19\t-132.76\t1\t2\tGCAGGCGCTGCCGGAAG\n+tig00000001\t-\t201961\t201986\tf432cb50-8c47-4c54-825f-ab04a38b5225\t-1.78\t-263.22\t-261.43\t1\t4\tAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAA\n+tig00000001\t-\t201640\t201675\t9a1776fd-9328-408e-9d54-e16535cbe064\t-7.77\t-292.55\t-284.78\t1\t6\tCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCA\n+tig00000001\t-\t201689\t201689\t9a1776fd-9328-408e-9d54-e16535cbe064\t-4.44\t-84.56\t-80.12\t1\t1\tTCTCCCGCACC\n+tig00000001\t-\t201707\t201725\t9a1776fd-9328-408e-9d54-e16535cbe064\t-6.52\t-184.97\t-178.45\t1\t5\tGTTACCGGACAAAACGCCCGCGCCGGTGG\n+tig00000001\t-\t201738\t201744\t9a1776fd-9328-408e-9d54-e16535cbe064\t-3.85\t-130.54\t-126.70\t1\t2\tAGCCTCGGCCCCGGCTG\n+tig00000001\t-\t201761\t201769\t9a1776fd-9328-408e-9d54-e16535cbe064\t-0.66\t-118.76\t-118.10\t1\t2\tGCCTGCGGATATCCGCTTA\n+tig00000001\t-\t201781\t201783\t9a1776fd-9328-408e-9d54-e16535cbe064\t-0.38\t-103.16\t-102.78\t1\t2\tAGCTGCGCGATAA\n+tig00000001\t-\t201810\t201813\t9a1776fd-9328-408e-9d54-e16535cbe064\t2.72\t-123.11\t-125.83\t1\t2\tGATAACGGCGGCAG\n+tig00000001\t-\t201846\t201870\t9a1776fd-9328-408e-9d54-e16535cbe064\t-1.76\t-186.93\t-185.17\t1\t4\tTTTCCCGGTGAGGACGGTTACAGCCGCAGCGAGTC\n+tig00000001\t-\t201889\t201902\t9a1776fd-9328-408e-9d54-e16535cbe064\t-6.80\t-137.98\t-131.18\t1\t5\tTGGTGCGCGGCGGCGTGGCGAAAC\n+tig00000001\t-\t201922\t201930\t9a1776fd-9328-408e-9d54-e16535cbe064\t-6.76\t-136.02\t-129.26\t1\t2\tGTCACCGGCTGGCCGCACT\n+tig00000001\t-\t201944\t201950\t9a1776fd-9328-408e-9d54-e16535cbe064\t-0.88\t-128.19\t-127.31\t1\t2\tGCAGGCGCTGCCGGAAG\n+tig00000001\t-\t201961\t201986\t9a1776fd-9328-408e-9d54-e16535cbe064\t-5.64\t-250.31\t-244.68\t1\t4\tAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAA\n+tig00000001\t-\t201738\t201744\t5c46f248-1989-4a1e-bdf7-265af033510b\t-4.07\t-119.11\t-115.04\t1\t2\tAGCCTCGGCCCCGGCTG\n+tig00000001\t-\t201761\t201769\t5c46f248-1989-4a1e-bdf7-265af033510b\t-11.76\t-179.35\t-167.60\t1\t2\tGCCTGCGGATATCCGCTTA\n+tig00000001\t-\t201781\t201783\t5c46f248-1989-4a1e-bdf7-265af033510b\t-8.43\t-230.67\t-222.24\t1\t2\tAGCTGCGCGATAA\n+tig00000001\t-\t201810\t201813\t5c46f248-1989-4a1e-bdf7-265af033510b\t-13.23\t-143.17\t-129.94\t1\t2\tGATAACGGCGGCAG\n+tig00000001\t-\t201846\t201870\t5c46f248-1989-4a1e-bdf7-265af033510b\t-19.59\t-262.23\t-242.64\t1\t4\tTTTCCCGGTGAGGACGGTTACAGCCGCAGCGAGTC\n+tig00000001\t-\t201889\t201902\t5c46f248-1989-4a1e-bdf7-265af033510b\t-20.86\t-186.25\t-165.38\t1\t5\tTGGTGCGCGGCGGCGTGGCGAAAC\n+tig00000001\t-\t201922\t201930\t5c46f248-1989-4a1e-bdf7-265af033510b\t-7.37\t-156.73\t-149.37\t1\t2\tGTCACCGGCTGGCCGCACT\n+tig00000001\t-\t201944\t201950\t5c46f248-1989-4a1e-bdf7-265af033510b\t-9.38\t-146.62\t-137.23\t1\t2\tGCAGGCGCTGCCGGAAG\n+tig00000001\t-\t201961\t201986\t5c46f248-1989-4a1e-bdf7-265af033510b\t-8.71\t-221.74\t-213.03\t1\t4\tAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAA\n+tig00000001\t-\t201846\t201870\t3ecfda3a-811b-4250-b282-811bffe581d5\t-8.36\t-179.24\t-170.89\t1\t4\tTTTCCCGGTGAGGACGGTTACAGCCGCAGCGAGTC\n+tig00000001\t-\t201889\t201902\t3ecfda3a-811b-4250-b282-811bffe581d5\t-23.12\t-210.39\t-187.27\t1\t5\tTGGTGCGCGGCGGCGTGGCGAAAC\n+tig00000001\t-\t201922\t201930\t3ecfda3a-811b-4250-b282-811bffe581d5\t-10.44\t-199.55\t-189.11\t1\t2\tGTCACCGGCTGGCCGCACT\n+tig00000001\t-\t201944\t201950\t3ecfda3a-811b-4250-b282-811bffe581d5\t-6.80\t-169.92\t-163.12\t1\t2\tGCAGGCGCTGCCGGAAG\n+tig00000001\t-\t201961\t201986\t3ecfda3a-811b-4250-b282-811bffe581d5\t-1.50\t-239.41\t-237.90\t1\t4\tAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAA\n'