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

Changeset 3:915637bd1d1a (2024-11-12)
Previous changeset 2:f355085dd2aa (2024-01-11)
Commit message:
planemo upload for repository https://github.com/mesocentre-clermont-auvergne/galaxy-tools/tree/master/tools/polypolish commit f3b5a8f978923c5bbf663b0a99025245458e6ee0
modified:
macro.xml
polypolish.xml
test-data/polished.fasta
b
diff -r f355085dd2aa -r 915637bd1d1a macro.xml
--- a/macro.xml Thu Jan 11 16:44:56 2024 +0000
+++ b/macro.xml Tue Nov 12 09:56:40 2024 +0000
[
@@ -1,48 +1,43 @@
-
 <macros>
-  <token name="@TOOL_VERSION@">0.5.0</token>
-  <token name="@VERSION_SUFFIX@">2</token>
-  <token name="@PROFILE@">21.05</token>
-
-  <token name="@THREADS@"><![CDATA[
+    <token name="@TOOL_VERSION@">0.6.0</token>
+    <token name="@VERSION_SUFFIX@">0</token>
+    <token name="@PROFILE@">23.02</token>
+    <token name="@THREADS@"><![CDATA[
       ##compute the number of ADDITIONAL threads to be used by samtools (-@)
       addthreads=\${GALAXY_SLOTS:-1} && (( addthreads-- )) &&
   ]]></token>
-
-  <xml name="version_command">
-     <version_command><![CDATA[polypolish -V]]></version_command>
-  </xml>
-  <xml name="xrefs">
-    <xrefs>
-      <xref type='bio.tools'>Polypolish</xref>
-    </xrefs>
-  </xml>
-  <xml name="requirements">
-    <requirements>
-      <requirement type="package" version="@TOOL_VERSION@">polypolish</requirement>
-      <requirement type="package" version="1.19">samtools</requirement>
-    </requirements>
-  </xml>
-  <xml name="citations">
-      <citations>
-          <citation type="doi">10.1371/journal.pcbi.1009802</citation>
-      </citations>
-  </xml>
-  <!-- Filter option of polypolish -->
-  <xml name="filter_option">
-      <conditional name="insert_filter">
-          <param name="filter_select" type="select" label="Filter by insert size ? [recommanded]"
-                 help="Exclude some alignments based on their insert size. It reduce the number of excessive alignments,
-                 particularly near the edges of repeat sequences">
-              <option value="filter" selected="true">Filter by insert size</option>
-              <option value="non_filter"> No filtering step</option>
-          </param>
-          <when value="filter">
-              <param argument="--low" type="float" min="0" value="0.1" label="Low percentile threshold" help="Select the lower value to remove [default: 0.1]"/>
-              <param argument="--high" type="float" min="0" value="99.9" label="High percentile threshold" help="Select the lower value to remove [default: 99.9]"/>
+    <xml name="version_command">
+        <version_command><![CDATA[polypolish -V]]></version_command>
+    </xml>
+    <xml name="xrefs">
+        <xrefs>
+            <xref type="bio.tools">Polypolish</xref>
+        </xrefs>
+    </xml>
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@TOOL_VERSION@">polypolish</requirement>
+            <requirement type="package" version="1.19.2">samtools</requirement>
+        </requirements>
+    </xml>
+    <xml name="citations">
+        <citations>
+            <citation type="doi">10.1371/journal.pcbi.1009802</citation>
+        </citations>
+    </xml>
+    <!-- Filter option of polypolish -->
+    <xml name="filter_option">
+        <conditional name="insert_filter">
+            <param name="filter_select" type="select" label="Filter by insert size ? [recommanded]" help="Exclude some alignments based on their insert size. It reduce the number of excessive alignments,                  particularly near the edges of repeat sequences">
+                <option value="filter" selected="true">Filter by insert size</option>
+                <option value="non_filter"> No filtering step</option>
+            </param>
+            <when value="filter">
+                <param argument="--low" type="float" min="0" value="0.1" label="Low percentile threshold" help="Select the lower value to remove [default: 0.1]"/>
+                <param argument="--high" type="float" min="0" value="99.9" label="High percentile threshold" help="Select the lower value to remove [default: 99.9]"/>
+            </when>
+            <when value="non_filter">
           </when>
-          <when value="non_filter">
-          </when>
-      </conditional>
-  </xml>
+        </conditional>
+    </xml>
 </macros>
b
diff -r f355085dd2aa -r 915637bd1d1a polypolish.xml
--- a/polypolish.xml Thu Jan 11 16:44:56 2024 +0000
+++ b/polypolish.xml Tue Nov 12 09:56:40 2024 +0000
[
b'@@ -5,9 +5,9 @@\n     <macros>\n         <import>macro.xml</import>\n     </macros>\n-    <expand macro=\'xrefs\'/>\n-    <expand macro="requirements" />\n-    <expand macro="version_command" />\n+    <expand macro="xrefs"/>\n+    <expand macro="requirements"/>\n+    <expand macro="version_command"/>\n     <command detect_errors="aggressive"><![CDATA[\n         ln -s \'$input.fasta_file\' input_data &&\n         #*======================================\n@@ -21,7 +21,7 @@\n             #else\n                 echo "${input.sam_data_type.single_sam} not a sam/bam file"\n             #end if\n-            polypolish input_data input_sam > \'$polished_fasta\'\n+            polypolish polish input_data input_sam > \'$polished_fasta\'\n         #*======================================\n                     For paired SAM\n         ======================================*#\n@@ -41,15 +41,15 @@\n                 echo "${input.sam_data_type.single_sam} not a sam/bam file"\n             #end if\n             #if $input.sam_data_type.insert_filter.filter_select == \'filter\'\n-                polypolish_insert_filter.py --low \'$input.sam_data_type.insert_filter.low\'\n+                polypolish filter --low \'$input.sam_data_type.insert_filter.low\'\n                                             --high \'$input.sam_data_type.insert_filter.high\'\n                                             --in1 sample_R1.sam\n                                             --in2 sample_R2.sam\n                                             --out1 \'filtered_1.sam\'\n                                             --out2 \'filtered_2.sam\' &&\n-                polypolish input_data \'filtered_1.sam\' \'filtered_2.sam\' >  $polished_fasta\n+                polypolish polish input_data \'filtered_1.sam\' \'filtered_2.sam\' >  $polished_fasta\n             #else\n-                polypolish input_data sample_R1.sam sample_R2.sam  >  $polished_fasta\n+                polypolish polish input_data sample_R1.sam sample_R2.sam  >  $polished_fasta\n             #end if\n         #*======================================\n             For multiple single-end SAM\n@@ -65,7 +65,7 @@\n                     echo "${single_sam} is not a sam/bam file"\n                 #end if\n             #end for\n-            polypolish input_data single_collection/*.sam > \'$polished_fasta\'\n+            polypolish polish input_data single_collection/*.sam > \'$polished_fasta\'\n         #*======================================\n             For multiple paired-end SAM\n         ======================================*#\n@@ -83,7 +83,7 @@\n                     ln -s \'$paired_sam.reverse\' \'paired_collection/reverse_input$(value)$(paired_sam.reverse.element_identifier).sam\' &&\n                 #end if\n                 #if $input.sam_data_type.insert_filter.filter_select == \'filter\'\n-                    polypolish_insert_filter.py --low \'$input.sam_data_type.insert_filter.low\'\n+                    polypolish filter --low \'$input.sam_data_type.insert_filter.low\'\n                                                 --high \'$input.sam_data_type.insert_filter.high\'\n                                                 --in1 \'paired_collection/forward_input$(value)$(paired_sam.forward.element_identifier).sam\'\n                                                 --in2 \'paired_collection/reverse_input$(value)$(paired_sam.reverse.element_identifier).sam\'\n@@ -95,9 +95,9 @@\n                     Filtering option\n         ======================================*#\n             #if $input.sam_data_type.insert_filter.filter_select == \'filter\'\n-                polypolish input_data paired_collection/*_filtered.sam > \'$polished_fasta\'\n+                polypolish polish input_data paired_collection/*_filtered.sam > \'$polished_fasta\'\n             #else\n-                polypolish input_data paired_collection/*.sam > \'$polished_fasta\'\n+                polypolish polish input_data paired_collection/*.sam > \'$polished_fasta\'\n             #end if\n         #end if\n         #*===================================='..b'an" truevalue="true" falsevalue="false" label="Keep per base information file"/>\n         </section>\n@@ -369,29 +364,29 @@\n             <section name="input">\n                 <param name="fasta_file" value="contigs.fa"/>\n                 <conditional name="sam_data_type">\n-                  <param name="sam_selector" value="multiple_paired"/>\n-                  <param name="paired_collection">\n-                      <collection type="list:paired">\n-                          <element name="paired_1">\n-                              <collection type="paired">\n-                                  <element name="forward" value="aligned_test_file/alignement_R1.bam" ftype="unsorted.bam"/>\n-                                  <element name="reverse" value="aligned_test_file/alignement_R2.bam" ftype="unsorted.bam"/>\n-                              </collection>\n-                          </element>\n-                          <element name="paired_2">\n-                              <collection type="paired">\n-                                  <element name="forward" value="aligned_test_file/alignement_R1_bis.bam" ftype="unsorted.bam"/>\n-                                  <element name="reverse" value="aligned_test_file/alignement_R2_bis.bam" ftype="unsorted.bam"/>\n-                              </collection>\n-                          </element>\n-                          <element name="paired_3">\n-                              <collection type="paired">\n-                                  <element name="forward" value="aligned_test_file/alignement_R1_ter.bam" ftype="unsorted.bam"/>\n-                                  <element name="reverse" value="aligned_test_file/alignement_R2_ter.bam" ftype="unsorted.bam"/>\n-                              </collection>\n-                          </element>\n-                      </collection>\n-                  </param>\n+                    <param name="sam_selector" value="multiple_paired"/>\n+                    <param name="paired_collection">\n+                        <collection type="list:paired">\n+                            <element name="paired_1">\n+                                <collection type="paired">\n+                                    <element name="forward" value="aligned_test_file/alignement_R1.bam" ftype="unsorted.bam"/>\n+                                    <element name="reverse" value="aligned_test_file/alignement_R2.bam" ftype="unsorted.bam"/>\n+                                </collection>\n+                            </element>\n+                            <element name="paired_2">\n+                                <collection type="paired">\n+                                    <element name="forward" value="aligned_test_file/alignement_R1_bis.bam" ftype="unsorted.bam"/>\n+                                    <element name="reverse" value="aligned_test_file/alignement_R2_bis.bam" ftype="unsorted.bam"/>\n+                                </collection>\n+                            </element>\n+                            <element name="paired_3">\n+                                <collection type="paired">\n+                                    <element name="forward" value="aligned_test_file/alignement_R1_ter.bam" ftype="unsorted.bam"/>\n+                                    <element name="reverse" value="aligned_test_file/alignement_R2_ter.bam" ftype="unsorted.bam"/>\n+                                </collection>\n+                            </element>\n+                        </collection>\n+                    </param>\n                 </conditional>\n             </section>\n             <section name="options">\n@@ -429,7 +424,7 @@\n             <output name="debug_file" value="debug_file_test_2.tsv"/>\n         </test>\n     </tests>\n-  <help><![CDATA[\n+    <help><![CDATA[\n     **What it does**\n     Polypolish is a tool for polishing genome assemblies with short reads.\n     Polypolish uses SAM/BAM files where each read has been aligned to all possible locations (not just a single best location).\n'
b
diff -r f355085dd2aa -r 915637bd1d1a test-data/polished.fasta
--- a/test-data/polished.fasta Thu Jan 11 16:44:56 2024 +0000
+++ b/test-data/polished.fasta Tue Nov 12 09:56:40 2024 +0000
b
b'@@ -1,8 +1,8 @@\n->contig00001_polypolish\n+>contig00001 len=5513 cov=11.7 corr=0 origname=NODE_1_length_5513_cov_11.741181_pilon sw=shovill-spades/1.1.0 date=20220302 polypolish\n CAGCATCGGACTCAGATAGTAATCCACGCTCTTTTAAAATGTCAACAAGAGAATCTCTACCATGAACAAAATGTGACTCATATCTAAACCAGTCCTTGACGAACGTGCCAAGCATATTAAGCCACTTCTCCTCATCCAACGCGTCAGTTTTTGACAGAATCGTTAGTTGATGGCGAAAGGTCGCAAAGTAAGAGCTTCTCGAGCTGCGCAAGGATAGGTCGAATTTTCTCATTTTCCGCCAGCAGTCCACTTCGATTTAATTCGTAAACAAGCAGTAGTAATTCCTGCTTTATCAAGATAATTTTTCGACTCATCAGAAATATCCGAAAGTGTTAACTTCTGCGTCATGGAAGCGATAAAACTCTGCAGGTTGGATACGCCAATCATTTTTATCGAAGCGCGCATAAATTTGAGCAGATTTGTCGTCACAGGTTGCGCCGCCAAAACGTCGGCTACAGTAACTTTTCCCAGCCTCAATCTCATCTCTCTTTTTGCGTTCTGCTTCAATATCTGGTTGAACGGCGTCGCGTCGTAACCCAGCTTGGTAAGTTGGATTAAGCACTCCGTGGACAGATTTGTCATTGTGAGCATTTTCATCCCGAAGTTGCGGCTCATTCTGATTCTGAACAGCTTCTTGGGAAGTAGCGACAGCTTGGTTTTTAGTGAGTTGTTCCATTCTTTAGCTCCTAGACCTTTAGCAGCAAGGTCCATATCTGACTTTTTGTTAACGTATTTAGCCACATAGAAACCAACAGCCATATAACTGGTAGCTTTAAGCGGCTCACCTTTAGCATCAACAGGCCACAACCAACCAGAACGTGAAAAAGCGTCCTGCGTGTAGCGAACTGCGATGGGCATACTGTAACCATAAGGCCACGTATTTTGCAAGCTATTTAACTGGCGGCGATTGCGTACCCGACGACCAAAATTAGGGTCAACGCTACCTGTAGGAAGTGTCCGCATAAAGTGCACCGCATGGAAATGAAGACGGCCATTAGCTGTACCATACTCAGGCACACAAAAATACTGATAGCAGTCGGCGTGTGAATCATTAGCCTTGCGACCCTCGGCAGCAAGAACCATACGACCAATATCACGAAAATAGTCACGCAAAGCATTGGGATTATCATAAAACGCCTCTAATCGGTCGTCAGCCAACGTGAGAGTGTCAAAAACGATAAACCAACCATCAGCATGAGCCTGTCGCATTGCATTCATCAAACGCTGAATAGCAAAGCCTCTACGCGATTTCATAGTGGAGGCCTCCAGCAATCTTGAACACTCATCCTTAATACCTTTCTTTTTGGGGTAATTATACTCATCGCGAATATCCTTAAGAGGGCGTTCAGCAGCCAGCTTGCGGCAAAACTGCGTAACCGTCTTCTCGTTCTCTAAAAACCATTTTTCGTCCCCTTCGGGGCGGTGGTCTATAGTGTTATTAATATCAAGTTGGGGGAGCACATTGTAGCATTGTGCCAATTCATCCATTAACTTCTCAGTAACAGATACAAACTCATCACGAACGTCAGAAGCAGCCTTATGGCCGTCAACATACATATCACCATTATCGAACTCAACGCCCTGCATACGAAAAGACAGAATCTCTTCCAAGAGCTTGATGCGGTTATCCATCTGCTTATGGAAGCCAAGCATTGGGGATTGAGAAAGAGTAGAAATGCCACAAGCCTCAATAGCAGGTTTAAGAGCCTCGATACGCTCAAAGTCAAAATAATCAGCGTGACATTCAGAAGGGTAATAAGAACGAACCATAAAAAAGCCTCCAAGATTTGGAGGCATGAAAACATACAATTGGGAGGGTGTCAATCCTGACGGTTATTTCCTAGACAAATTAGAGCCAATACCATCAGCTTTACCGTCTTTCCAGAAATTGTTCCAAGTATCGGCAACAGCTTTATCAATACCATGAAAAATATCAACCACACCAGAAGCAGCATCAGTGACGACATTAGAAATATCCTTTGCAGTAGCGCCAATATGAGAAGAGCCATACCGCTGATTCTGCGTTTGCTGATGAACTAAGTCAACCTCAGCACTAACCTTGCGAGTCATTTCTTTGATTTGGTCATTGGTAAAATACTGACCAGCCGTTTGAGCTTGAGTAAGCATTTGGCGCATAATCTCGGAAACCTGCTGTTGCTTGGAAAGATTGGTGTTTTCCATAATAGACGCAACGCGAGCAGTAGACTCCTTCTGTTGATAAGCAAGCATCTCATTTTGTGCATATACCTGGTCTTTCGTATTCTGGCGTGAAGTCGCCGACTGAATGCCAGCAATCTCTTTTTGAGTCTCATTTTGCATCTCGGCAATCTCTTTCTGATTGTCCAGTTGCATTTTAGTAAGCTCTTTTTGATTCTCAAATCCGGCGTCAACCATACCAGCAGAGGAAGCATCAGCACCAGCACGCTCCCAAGCATTAAGCTCAGGAAATGCAGCAGCAAGATAATCACGAGTATCCTTTCCTTTATCAGCGGCAGACTTGCCACCAAGTCCAACCAAATCAAGCAACTTATCAGAAACGGCAGAAGTGCCAGCCTGCAACGTACCTTCAAGAAGTCCTTTACCAGCTTTAGCCATAGCACCAGAAACAAAACTAGGGACGGCCTCATCAGGGTTAGGAACATTAGAGCCTTGAATGGCAGATTTAATACCAGCATCACCCATGCCTACAGTATTGTTATCGGTAGCAAGCACATCACCTTGAATGCCACCGGAGGCGGCTTTTTGACCGCCTCCAAACAATTTAGACATGGCGCCACCAGCAAGAGCAGAAGCAATACCGCCAGCAATAGCACCAAACATAAATCACCTCACTTAAGTGGCTGGAGACAAATAATCTCTTTAATAACCTGATTCAGCGAAACCAATCCGCGGCATTTAGTAGCGGTAAAGTTAGACCAAACCATGAAACCAACATAAACATTATTGCCCGGCGTACGGGGAAGGACGTCAATAGTCACACAGTCCTTGACGGTATAATAACCACCATCATGGCGACCATCCAAAGGATAAACATCATAGGCAGTCGGGAGGGTAGTCGGAACCGAAGAAGACTCAAAGCGAACCAAACAGGCAAAAAATTTAGGGTCGGCATCAAAAGCAATATCAGCACCAACAGAAACAACCTGATTAGCGGCGTTGACAGATGTATCCATCTGAATGCAATGAAGAAAACCACCATTACCAGCATTAACCGTCAAACTATCAAAATATAACGTTGACGATGTAGCTTTAGGTGTCTGTAAAACAGGTGCCGAAGAAGCTGGAGTAACAGAAGTGAGAACCAGCTTATCAGAAAAAAAGTTTGAATTATGGCGAGAAATAAAAGTCTGAAACATGATTAAACTCCTAAGCAGAAAACCTACCGCGCTTCGCTTGGTCAACCCCTCAGCGGCAAAAATTAAAATTTTTACCGCTTCGGCGTTATAACCTCACACTCAATCTTTTATCACGAAGTCATGATTGAATCGCGAGTGGTCGGCAGATTGCGATAAACGGTCACATTAAATTTAACCTGACTATTCCACTGCAACAACTGAACGGACTGGAAACACTGGTCATAATCATGGTGGCGAATAAGTACGCGTTCTTGCAAATCACCAGAAGGCGGTTCCTGAATGAATGGGAAGCCTTCAAGAAGGTGATAAGCAGGAGAAACATACGAAGGCGCATAACGATACCACTGACCCTCAGCAATCTTAAACTTCTTAGACGAATCACCAGAACGGAAAACATCCTTCATAGAAATTTCACGCGGCGGCAAGTTGCCATACAAAACAGGGTCGCCAGCAATATCGGTATAA'..b'TGATGCAGCTACTTACTATGGAAATGGTCTTTATTGTAACAAAGAAAAATGTTGGGTAGATTGGAATCAAGCTAAAGGAGAAATTGGAAAAATTATTGTTAATGGTTGGGTTAATCATGGTCCATGGGCACCTAGAAGGTAGTGGTAGGGAGGTTATTATAATGGATTTTACTAAAGAAGAAAAACTTTTAAATGCAATTAGTAAAGTATACAATGAAGCAACTATAGATGACTATCCTGACTTAAAAGAAAAGCTCTTTCTTTATTCTAAAGAAATCAGTGAGGGAAAAAGTGTTGGTGAAGTTAGTATGAAATTAAGTAGTTTTCTTGGAAGATATATTTTAAAACATAAATTTGGATTACCTAAATCTTTAATAGAATTACAAGAAATTGTTAGTAAGGAATCTCAAGTATATAGAGGATGGGCTTCTATTGGTATTTGGAGTTAATTCTTTGGAATAAAATGAAAATGGGACTGAAGAATCAATTCTGATTCATCAGTCCCATAAATTATAGAGCCGAAAAATTATCATATTTATAGATTTATTAATAGTAAATGGCATAAGAATCTTTGAACAACGAAAAAAATGTAGCTCTACCTAATTATTGATTTGCAACAATTTCACCGAAGGTGAAACAAGCGAAACTTTACAAAAAGTGGCTAGCATTTGCTAGTCTTTTTTGTATTGTGGAGAGTGG\n->contig00004_polypolish\n+>contig00004 len=3138 cov=2.3 corr=1 origname=NODE_4_length_3138_cov_2.261375_pilon sw=shovill-spades/1.1.0 date=20220302 polypolish\n TTTTTTTCCGATATAATTTGCTGGTAAATCTATCACAGCCAATGAGAGTCCGGGATTGCTTTCTTCCTGATAATCTTTCCGTTGGTATAAAATCTTTTCGCCTAGATATACAGTAATCCATTGATGGTTTCCTTTGATTAGTAAATTTGCTTCGTCCAACGGTTTTTCAAGTGTCGTTTCCATAATGAATGTTTCGTCTTTATCCACGGTAGGAATCGTACGAGTGTACTTGCTGTCAAACTTAAAATCAGGATCATCTTTTGTATAGAAAACCCACGAATCTTTTATTTCTTTCATTTGATTGCTGTCCGCAGGAGAAAAAATAAAATAAAAGAAATAAGCAATTCCTGTTAAAAATAACAATGATAATATAATGATTATGCTCCACTTTATACGTAAATAAAAAGATTCTTTTTTATTATATATGTTATTAAAAAAAGACATAAAGGTATTCCTTTGAGAGTATAAAATATTTTGTGTAAATGAAAAAATCCATACAAAAAAGGAAGTCCCTTCTGTAGAATAAAGTTAACGACAACCAATTCACAGAAAAGAGGACTTCCCTATGAATGATTTTACTACAGAAATTGTGCAAACTCTAGTCACTAAAGGCGATTTAAATGAATTATTCCGTTCGCACTTAGAAAAAGCGATAAACACACTCCTACGGACTGAATTAACGGCTTTTTTAGATTACGAAAAATATGATCGCACTGGTTTTAATTCAGGTAATTCGAGAAACGGTTCTTACTTTCGATCAATCAAAACCGAATATGGTGAATTAACATTGGAAATACCTAGAGATCGTAATGGTGAGTTTAAACAACAAACTTTACCAGCCTACAAAAGAACAAACGATACATTGGAAACCACTATTATCCATTTATTCGAAAAAGGTGTTACGATGTCTGAAATTGCTGATTTGATCGAAAAAATGTACGGTCATCACTATACTCCACAAACCATGTCCAACATGACTAAAGTTCTGACTGAAGAAGTAAATGCCTTTAAATCCAGAGCCTTAAATGATAAGTATGTCGCTATTTTTATGGACGTTACTTACATTCCACTAAAACGTCAAACTGTATCCAAAGAAGCGATTTATATTGCCATTGGTATACGAGAAGACGGCACTAAAGAAGTACTGAGTTATGCGATTGCTCCAACTGAATCAACATACGTTTGGAATGAGCTGCTACAGGATATTAACTCCAGAGGAGTTCAAGAAGTCTTGCTTTTTATTACGGACGGCTTAAAAGGCATGAAAGATACTATCCATCAAATTTATCCTAAAGCAAAATATCAGCATTGTTGTATCCATGTATCTCGTAACATCGCTCATAAAGTACGTGTCAAAGACCGAAAAGAAATCTGTGATGACTTTAAGGCTGTTTATCAAGCTAACTCAAAAGAAGAAGCGAATACCTTCTTATCCGGCATGATTGAGAAATGGAAGAAAAACTATCCTAAAGTGACGCAGTCACTCATAGAAAACCAAGACTTATTAACTTTTTATGATTTTCCACCTAGCATTCGTAGAACCATTTACTCAACCAATCTAATCGAGTCTTTCAATAAGCAAATTAAAAGATACAGCCGTAGAAAAGAGCAGTTTCAAAATGAAGAATCACTAGAACGCTTTCTAGTCAGCATTTTTGATACATACAATCAAAAATTTCTAAACAGAAGCCATAAAGGTTTTCAACAGGTAACCGATACATTAGTTTCAATGTTTACTGAGTAACTAATTATTTTGCAGGAGGACAATTTATTTACACAAAATTATTGACGCTCCCTAACAAAATCCCACCTACTTTATCTTATAACTTGTTTTTTTGCATAAAAACTAAATAAAGAGTGCGGCATTTCGCACTTTTTTATTATATATTATTTGTCAAATTAAGCGAGACAAAACTTTCCATTTACGTAACTCAAAAAAACTTCTAATGGTGTTTGGTAATTTAGTGATTTTCTAGGGATATGGTTTCTTTTGGACGCAACGGATGAGATGAATCCTTGATTGACTTGGTGGAAGTCCATTTCTTTTGGTAGTCCATCTTTTCGGAGGAGTCCGTTTGAATGTTCATTTAATCCCCGTTGGGAAGGCGTTCCTGGGTCTGCGAAATAAATATCAATATCATTCGTATTACTAATACTTTTCCAATTAGAGAATTCTTTTCCACAATCAAAGGTAATTGATTTGAAAAGATTTTTGGGTACGGATTGAAGCCATTCATTAATCGAATTTTCAATATCAACTGCCTTACGGCCTTCTGGTTTCAAGGCGATAATGGCTTTTGAAAGTCGCTCAACGAGTGTGATGACGGCACTTTTATGGTGGATGCCGACAATCGTATCACCCTCTAAATGTCCAAATTCTTCTTCGAAAACGACATAATCTTTTTTACGTTCAGAAATATTTCTTTTGAATGCCTGTTTTCCACGTTTTTCCTTATAACCATTTGGTTTTCGTTTTCCTTGCATCGGTAAATGTAAAACATTGAATTCTCCCGTCTTAAACCGGCGATAGAGGGTACTTACTGAACAAGAGAAGGTTCGCTCTGCACGGCCAATAATGACGTCTGGGGTCCAACCTTCAGTTGATCTTTTTGCAATGTATTCTTTTTCCTCAGCAGGAAAAATGATTTCGGTTCTCCCACAACGTCGCTTATTTTCTTTGTATTGTTCAAAGTATTCAAGTGCTGTTCCACCACATTTGAGAAAGTTATAGACATTGTAAATTGTTTGGCGTCCACGTTTAAGCTGCTTCGCAACGATAGCAACCGGAGTTTCTTGATGAAAATATGCTTCTATCATTACAAGCTCGTTTGGTGTAAGATGGGTATAAGTCATTTATGTTCACTCTCCTTGTATGCTTTAGCGGGTATTACAATTGGAGTGTAACATAAATGGCTTTTTTATTTGTCTCGCTTAATTATACAAACGGCGTTGTTAAATAAATAGGAACAAAAAAATTAAATTATTTAATACACACATTAGGGATCATTCATTGAAGGATTAGAAGACAACGAGGATGGCGAAAAATATGATTACAATAGAAAAGTCGTTCTAAACCAATCAGTTCCTAGTATCACACAACGCAAGCCAGTCGATTTGGTGAGATTAGAGGAGTGAAAG\n'