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

Changeset 0:ff827f983bbb (2024-05-03)
Next changeset 1:b820ef95ca20 (2024-05-14)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/fastk commit 57fa4dcee2aaa115df4a5b7173a3a8fb889bc5a3
added:
fastk.xml
macros.xml
test-data/input01.fasta.gz
test-data/test02.tabex.txt
test-data/test03.tabex.txt
b
diff -r 000000000000 -r ff827f983bbb fastk.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/fastk.xml Fri May 03 20:23:04 2024 +0000
[
@@ -0,0 +1,114 @@
+<tool id="fastk_fastk" name="FastK" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="23.2">
+    <description>A k-mer counter for high-quality assembly datasets</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements" />
+    <command detect_errors="exit_code"><![CDATA[
+        mkdir -p outfiles/tmpfiles && 
+        #if $infile.is_of_type("fastq"):
+            #set INPUTFILE="input.fastq"
+        #elif $infile.is_of_type("fastq.gz"):
+            #set INPUTFILE="input.fastq.gz"
+        #else
+            #set INPUTFILE="input."+$infile.ext
+        #end if
+        ln -s '$infile' $INPUTFILE && 
+        FastK $INPUTFILE 
+        -k$kmer_size 
+        #if $sorted_table.sorted_table_option == 'yes_with_default':
+            -t 
+        #elif $sorted_table.sorted_table_option == 'yes_with_custom':
+            -t$sorted_table_cutoff 
+        #end if
+        -T\${GALAXY_SLOTS:-8} -Noutfiles/output -Poutfiles/tmpfiles
+        #if $sorted_table.sorted_table_option == 'yes_with_default':
+            && Tabex outfiles/output.ktab -t${sorted_table.tabex_threshold_for_default} LIST > '$tabex_hist'
+        #elif $sorted_table.sorted_table_option == 'yes_with_custom':
+            && Tabex outfiles/output.ktab -t${sorted_table.tabex_threshold_for_custom} LIST > '$tabex_hist'
+        #end if
+        && tar -c -f fastk.tar ./outfiles/
+    ]]></command>
+    <inputs>
+        <param name="infile" type="data" format="fasta,fasta.gz,fastq,fastq.gz,cram,unsorted.bam,sam" label="Input file"/>   
+        <param name="kmer_size" argument="-k" type="integer" min="5" max="50" value="40" label="Enter desired k-mer size" help="Default: 40" />
+        <conditional name="sorted_table">
+            <param name="sorted_table_option" type="select" label="Sort table" help="Do you want a sorted table of all canonical k-mers and their counts? The sorted table is sorted lexicographically on the k-mer where a &lt; c &lt; g &lt; t.">
+                <option value="no">No</option>
+                <option value="yes_with_default">Yes, Default sorted </option>
+                <option value="yes_with_custom">Yes, Custom sorted </option>
+            </param>
+            <when value="no"/>
+            <when value="yes_with_default">
+                <param name="tabex_threshold_for_default" argument="-t" type="integer" value="5" min="1" label="Tabex count threshold" help="Trim all k-mers with counts less than threshold"/>
+            </when>
+            <when value="yes_with_custom">
+                <param name="sorted_table_cutoff" type="integer" min="2" value="10" label="Enter sorted table cutoff value"/>
+                <param name="tabex_threshold_for_custom" argument="-t" type="integer" value="5" min="1" label="Tabex count threshold" help="Trim all k-mers with counts less than threshold"/>
+            </when>
+        </conditional>
+    </inputs>
+    <outputs>
+        <data name="fastk_out" format="tar" from_work_dir="fastk.tar" label="${tool.name} on ${on_string}: FastK files"/>
+        <data name="fastk_hist_out" format="fastk_hist" from_work_dir="outfiles/output.hist" label="${tool.name} on ${on_string}: FastK hist" />
+        <data name="tabex_hist" format="txt"  label="${tool.name} on ${on_string}: Tabex output">
+            <filter> sorted_table['sorted_table_option'] != 'no' </filter>
+        </data>
+    </outputs>
+    <tests>
+    <!-- TEST 1 -->
+        <test expect_num_outputs="2">
+            <param name="infile" value="input01.fasta.gz"/>
+            <output name="fastk_out" ftype="tar">
+                <assert_contents>
+                    <has_archive_member path="./outfiles/output.hist"/>
+                </assert_contents>
+            </output>
+        </test>
+        <!-- TEST 2 -->
+        <test expect_num_outputs="3">
+            <param name="infile" value="input01.fasta.gz"/>
+            <conditional name="sorted_table">
+                <param name="sorted_table_option" value="yes_with_default"/>
+            </conditional>
+            <output name="fastk_out" ftype="tar">
+                <assert_contents>
+                    <has_archive_member path="./outfiles/output.hist"/>
+                    <has_archive_member path="./outfiles/output.ktab"/>
+                </assert_contents>
+            </output>
+            <output name="tabex_hist" value="test02.tabex.txt"/>
+        </test>
+        <!-- TEST 3 -->
+        <test expect_num_outputs="3">
+            <param name="infile" value="input01.fasta.gz"/>
+            <conditional name="sorted_table">
+                <param name="sorted_table_option" value="yes_with_custom"/>
+                <param name="sorted_table_cutoff" value="5"/>
+            </conditional>
+            <output name="fastk_out" ftype="tar">
+                <assert_contents>
+                    <has_archive_member path="./outfiles/output.hist"/>
+                    <has_archive_member path="./outfiles/output.ktab"/>
+                </assert_contents>
+            </output>
+            <output name="tabex_hist" value="test03.tabex.txt"/>
+        </test>       
+    </tests>
+    <help><![CDATA[
+        FastK is a k‑mer counter that is optimized for processing high quality DNA assembly data sets such as those produced with an Illumina instrument or a PacBio run in HiFi mode.
+
+        The input data can be in CRAM, BAM, SAM, fasta, or fastq files.
+    
+        FastK produces the following outputs:
+
+        1. FastK hist: file in binary format containing histogram information detailing the frequency of occurrence for each k‑mer within the dataset.
+
+        2. A Tabex txt file comprising a table of k‑mer/count pairs, sorted lexicographically on the k‑mer sequence, followingthe order a < c < g < t
+
+        3. A tar file comprising of hidden .ktab files that can be used for downstream FASTK tools.
+
+
+    ]]></help>
+    <expand macro="citations"/>
+</tool>
\ No newline at end of file
b
diff -r 000000000000 -r ff827f983bbb macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml Fri May 03 20:23:04 2024 +0000
b
@@ -0,0 +1,23 @@
+<macros>
+    <token name="@TOOL_VERSION@">1.0.0</token>
+    <token name="@VERSION_SUFFIX@">0</token>
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@TOOL_VERSION@">fastk</requirement>
+        </requirements>
+    </xml>
+    <xml name="citations">
+        <citations>
+            <citation type="bibtex">                                                                   
+                @misc{github,                                                                 
+                author = {Gene Meyers},
+                year = {2020},                                                                 
+                title = {FastK},                                                        
+                publisher = {GitHub},
+                journal = {GitHub repository},                                               
+                url = {https://github.com/thegenemyers/FASTK},                        
+                }                                                                    
+            </citation>
+        </citations>
+    </xml>
+</macros>
\ No newline at end of file
b
diff -r 000000000000 -r ff827f983bbb test-data/input01.fasta.gz
b
Binary file test-data/input01.fasta.gz has changed
b
diff -r 000000000000 -r ff827f983bbb test-data/test02.tabex.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test02.tabex.txt Fri May 03 20:23:04 2024 +0000
b
b'@@ -0,0 +1,311 @@\n+Opening 40-mer table with 502,597 entries\n+      1348: aaaaacgaaaccgttgttgaaatccctggaaagctggttg =     5\n+      3497: aaaacatcttctgcaattccacgaaggcgctctcccattc =     6\n+      4142: aaaacgaaaccgttgttgaaatccctggaaagctggttgc =     5\n+      8271: aaaattagccaccatgccgcgcatctcacacccgcgcatg =     5\n+      9568: aaacaggagaaagactatccccagtaacatcatgcccgta =     5\n+      9898: aaacatcttctgcaattccacgaaggcgctctcccattcg =     6\n+     11316: aaaccgttgttgaaatccctggaaagctggttgcgtgaaa =     5\n+     11609: aaacgaaaccgttgttgaaatccctggaaagctggttgcg =     5\n+     14677: aaagactatccccagtaacatcatgcccgtattccctgcc =     5\n+     23585: aaattagccaccatgccgcgcatctcacacccgcgcatgc =     5\n+     25526: aacaacggtttcgtttttcgctgacgttcagcaagccgct =     5\n+     27763: aacaggagaaagactatccccagtaacatcatgcccgtat =     5\n+     28631: aacatcatgcccgtattccctgccagctttaacacgtccc =     6\n+     30223: aaccagctttccagggatttcaacaacggtttcgtttttc =     5\n+     30777: aacccaccacggatatccgtaaggacagatgcagcaagtg =     5\n+     32167: aaccggaaagacacgacaaagcgccactggcagcagccat =     5\n+     32813: aaccgttgttgaaatccctggaaagctggttgcgtgaaaa =     5\n+     33721: aacgaaaccgttgttgaaatccctggaaagctggttgcgt =     5\n+     34290: aacgagcgggtcatgcagttgagtcagcaggtgcagcact =     5\n+     35527: aacgcgatgctgaacgagcgggtcatgcagttgagtcagc =     7\n+     36570: aacggcatctttagtttgcacgataaccatatcttcaatt =     5\n+     36923: aacgggacaacgcgatgctgaacgagcgggtcatgcagtt =     7\n+     37484: aacggtttcgtttttcgctgacgttcagcaagccgctgct =     5\n+     37751: aacgtcagcgaaaaacgaaaccgttgttgaaatccctgga =     5\n+     39094: aactcaggttttactgagttacttaatcgcgtagaagcct =     5\n+     42359: aagaacggcatctttagtttgcacgataaccatatcttca =     5\n+     44118: aagaggcgatgatgcagcacaggagggacgtgttaaagct =     5\n+     45106: aagatgttttcgcttacgcaacgggacaacgcgatgctga =     6\n+     46997: aagccgctgctctgccggcattccccttatatccgcctct =     7\n+     47380: aagcgaaaacatcttctgcaattccacgaaggcgctctcc =     6\n+     50897: aaggcttctacgcgattaagtaactcagtaaaacctgagt =     5\n+     51296: aaggggaatgccggcagagcagcggcttgctgaacgtcag =     6\n+     52678: aagtaccgttgcaccttcccggttacgctccctccccaga =     5\n+     53910: aagtgctgcacctgctgactcaactgcatgacccgctcgt =     5\n+     57240: aatacgggcatgatgttactggggatagtctttctcctgt =     5\n+     65459: aatgccggcagagcagcggcttgctgaacgtcagcgaaaa =     6\n+     66807: aatggctgctgccagtggcgctttgtcgtgtctttccggt =     5\n+     66876: aatgggagagcgccttcgtggaattgcagaagatgttttc =     6\n+     68705: aattagccaccatgccgcgcatctcacacccgcgcatgca =     5\n+     70117: aattgcagaagatgttttcgcttacgcaacgggacaacgc =     6\n+     73840: acaacgcgatgctgaacgagcgggtcatgcagttgagtca =     7\n+     73929: acaacggtttcgtttttcgctgacgttcagcaagccgctg =     5\n+     79713: acaggagggacgtgttaaagctggcagggaatacgggcat =     7\n+     81836: acatcatgcccgtattccctgccagctttaacacgtccct =     6\n+     84684: accaatggctgctgccagtggcgctttgtcgtgtctttcc =     5\n+     86754: accagctttccagggatttcaacaacggtttcgtttttcg =     5\n+     87548: accatcaccgcccccggtatcccgctcatgctggcgaacc =     6\n+     88541: acccaccacggatatccgtaaggacagatgcagcaagtgc =     5\n+     89933: acccgctcgttcagcatcgcgttgtcccgttgcgtaagcg =     7\n+     92340: accgcccccggtatcccgctcatgctggcgaacccagaca =     6\n+     95583: accgttgcaccttcccggttacgctccctccccagaaggt =     5\n+     95641: accgttgttgaaatccctggaaagctggttgcgtgaaaag =     5\n+     97042: acctgctgactcaactgcatgacccgctcgttcagcatcg =     7\n+     98459: acgaaaccgttgttgaaatccctggaaagctggttgcgtg =     5\n+    101251: acgcaacgggacaacgcgatgctgaacgagcgggtcatgc =     7\n+    102467: acgccatagaggcggatataaggggaatgccggcagagca =     7\n+    103655: acgcgatgctgaacgagcgggtcatgcagttgagtcagca =     7\n+    106812: acggcatctttagtttgcacgataaccatatcttcaattc =     5\n+    107733: acgggacaacgcgatgctgaacgagcgggtcatgcagttg =     7\n+    107912: acgggcatgatgttactggggatagtctttctcctgttta =     5\n+    109184: acggttcaccatgaaccataatatgatttctaagaatacc =     5\n+    109323: acggtttcgtttttcgctgacgttcagcaagccgctgctc =     5\n+    109991: acgtcagcgaaaaacgaaaccgttgttgaaatccctggaa =     5\n+    111273: acgttcagcaagccgctgc'..b'agtcagcaggtgcagcacttga =     5\n+    397948: gatataaggggaatgccggcagagcagcggcttgctgaac =     7\n+    402619: gatgttttcgcttacgcaacgggacaacgcgatgctgaac =     6\n+    402865: gattaccttctggggagggagcgtaaccgggaaggtgcaa =     5\n+    404021: gatttcaacaacggtttcgtttttcgctgacgttcagcaa =     5\n+    407354: gcacaggagggacgtgttaaagctggcagggaatacgggc =     7\n+    409109: gcagagcagcggcttgctgaacgtcagcgaaaaacgaaac =     5\n+    412106: gcatctttagtttgcacgataaccatatcttcaattccaa =     5\n+    412209: gcatgatgttactggggatagtctttctcctgtttaccgc =     5\n+    413124: gccaaaattagccaccatgccgcgcatctcacacccgcgc =     5\n+    413965: gccactcgcggtaaacaggagaaagactatccccagtaac =     5\n+    415288: gccatagaggcggatataaggggaatgccggcagagcagc =     7\n+    416467: gcccccggtatcccgctcatgctggcgaacccagacacgc =     6\n+    421338: gccttcgtggaattgcagaagatgttttcgcttacgcaac =     6\n+    421606: gcgaaaacatcttctgcaattccacgaaggcgctctccca =     6\n+    425114: gcgccttcgtggaattgcagaagatgttttcgcttacgca =     6\n+    428102: gcgggtcatgcagttgagtcagcaggtgcagcacttgagc =     5\n+    432431: gctcgttcagcatcgcgttgtcccgttgcgtaagcgaaaa =     7\n+    432544: gctctgccggcattccccttatatccgcctctatggcgta =     7\n+    432796: gctgaacgtcagcgaaaaacgaaaccgttgttgaaatccc =     5\n+    435138: gctgttacttcccgctctaaagtagaagtaacatgaatca =     5\n+    435931: gcttctacgcgattaagtaactcagtaaaacctgagttaa =     5\n+    436168: gcttgctgaacgtcagcgaaaaacgaaaccgttgttgaaa =     5\n+    436476: gctttccagggatttcaacaacggtttcgtttttcgctga =     5\n+    438332: ggaatgccggcagagcagcggcttgctgaacgtcagcgaa =     6\n+    438435: ggaattgcagaagatgttttcgcttacgcaacgggacaac =     6\n+    439017: ggacgtgttaaagctggcagggaatacgggcatgatgtta =     5\n+    439212: ggagaaagactatccccagtaacatcatgcccgtattccc =     5\n+    440307: ggatataaggggaatgccggcagagcagcggcttgctgaa =     7\n+    441337: ggatttcaacaacggtttcgtttttcgctgacgttcagca =     5\n+    442536: ggcagagcagcggcttgctgaacgtcagcgaaaaacgaaa =     6\n+    443224: ggcatctttagtttgcacgataaccatatcttcaattcca =     5\n+    444964: ggcgatgatgcagcacaggagggacgtgttaaagctggca =     5\n+    447754: ggcttctacgcgattaagtaactcagtaaaacctgagtta =     5\n+    447808: ggcttgctgaacgtcagcgaaaaacgaaaccgttgttgaa =     5\n+    448268: gggaatgccggcagagcagcggcttgctgaacgtcagcga =     6\n+    449217: gggcatgatgttactggggatagtctttctcctgtttacc =     5\n+    451194: gggtcatgcagttgagtcagcaggtgcagcacttgagcga =     5\n+    455477: ggttcaccatgaaccataatatgatttctaagaataccaa =     5\n+    459621: gtattccctgccagctttaacacgtccctcctgtgctgca =     6\n+    459664: gtattcttagaaatcatattatggttcatggtgaaccgta =     5\n+    464117: gtgagagcgctgttacttcccgctctaaagtagaagtaac =     5\n+    465705: gtggaattgcagaagatgttttcgcttacgcaacgggaca =     6\n+    470168: gttgtacgccatagaggcggatataaggggaatgccggca =     5\n+    471622: gtttttcgctgacgttcagcaagccgctgctctgccggca =     6\n+    474188: taaggggaatgccggcagagcagcggcttgctgaacgtca =     6\n+    477666: tactttagagcgggaagtaacagcgctctcacttacaaaa =     5\n+    481177: tattcttagaaatcatattatggttcatggtgaaccgtaa =     5\n+    482917: tcaccgcccccggtatcccgctcatgctggcgaacccaga =     6\n+    484763: tcatcttttcacgcaaccagctttccagggatttcaacaa =     5\n+    486368: tcccgttgcgtaagcgaaaacatcttctgcaattccacga =     6\n+    488979: tcgcgttgtcccgttgcgtaagcgaaaacatcttctgcaa =     6\n+    490132: tcgttcagcatcgcgttgtcccgttgcgtaagcgaaaaca =     7\n+    490368: tctactttagagcgggaagtaacagcgctctcacttacaa =     5\n+    491019: tctgccggcattccccttatatccgcctctatggcgtaca =     7\n+    491604: tcttcatcttttcacgcaaccagctttccagggatttcaa =     5\n+    492432: tgaacgagcgggtcatgcagttgagtcagcaggtgcagca =     5\n+    492674: tgaagtaccgttgcaccttcccggttacgctccctcccca =     5\n+    493518: tgagagcgctgttacttcccgctctaaagtagaagtaaca =     5\n+    494230: tgatgcagcacaggagggacgtgttaaagctggcagggaa =     5\n+    494813: tgcacctgctgactcaactgcatgacccgctcgttcagca =     7\n+    495175: tgcatttacggttcaccatgaaccataatatgatttctaa =     5\n+    497093: tggaattgcagaagatgttttcgcttacgcaacgggacaa =     6\n+    498688: tgtaagtgagagcgctgttacttcccgctctaaagtagaa =     5\n+    499514: tgttgaaatccctggaaagctggttgcgtgaaaagatgaa =     5\n'
b
diff -r 000000000000 -r ff827f983bbb test-data/test03.tabex.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test03.tabex.txt Fri May 03 20:23:04 2024 +0000
b
b'@@ -0,0 +1,311 @@\n+Opening 40-mer table with 310 entries occuring 5-or-more times\n+         0: aaaaacgaaaccgttgttgaaatccctggaaagctggttg =     5\n+         1: aaaacatcttctgcaattccacgaaggcgctctcccattc =     6\n+         2: aaaacgaaaccgttgttgaaatccctggaaagctggttgc =     5\n+         3: aaaattagccaccatgccgcgcatctcacacccgcgcatg =     5\n+         4: aaacaggagaaagactatccccagtaacatcatgcccgta =     5\n+         5: aaacatcttctgcaattccacgaaggcgctctcccattcg =     6\n+         6: aaaccgttgttgaaatccctggaaagctggttgcgtgaaa =     5\n+         7: aaacgaaaccgttgttgaaatccctggaaagctggttgcg =     5\n+         8: aaagactatccccagtaacatcatgcccgtattccctgcc =     5\n+         9: aaattagccaccatgccgcgcatctcacacccgcgcatgc =     5\n+        10: aacaacggtttcgtttttcgctgacgttcagcaagccgct =     5\n+        11: aacaggagaaagactatccccagtaacatcatgcccgtat =     5\n+        12: aacatcatgcccgtattccctgccagctttaacacgtccc =     6\n+        13: aaccagctttccagggatttcaacaacggtttcgtttttc =     5\n+        14: aacccaccacggatatccgtaaggacagatgcagcaagtg =     5\n+        15: aaccggaaagacacgacaaagcgccactggcagcagccat =     5\n+        16: aaccgttgttgaaatccctggaaagctggttgcgtgaaaa =     5\n+        17: aacgaaaccgttgttgaaatccctggaaagctggttgcgt =     5\n+        18: aacgagcgggtcatgcagttgagtcagcaggtgcagcact =     5\n+        19: aacgcgatgctgaacgagcgggtcatgcagttgagtcagc =     7\n+        20: aacggcatctttagtttgcacgataaccatatcttcaatt =     5\n+        21: aacgggacaacgcgatgctgaacgagcgggtcatgcagtt =     7\n+        22: aacggtttcgtttttcgctgacgttcagcaagccgctgct =     5\n+        23: aacgtcagcgaaaaacgaaaccgttgttgaaatccctgga =     5\n+        24: aactcaggttttactgagttacttaatcgcgtagaagcct =     5\n+        25: aagaacggcatctttagtttgcacgataaccatatcttca =     5\n+        26: aagaggcgatgatgcagcacaggagggacgtgttaaagct =     5\n+        27: aagatgttttcgcttacgcaacgggacaacgcgatgctga =     6\n+        28: aagccgctgctctgccggcattccccttatatccgcctct =     7\n+        29: aagcgaaaacatcttctgcaattccacgaaggcgctctcc =     6\n+        30: aaggcttctacgcgattaagtaactcagtaaaacctgagt =     5\n+        31: aaggggaatgccggcagagcagcggcttgctgaacgtcag =     6\n+        32: aagtaccgttgcaccttcccggttacgctccctccccaga =     5\n+        33: aagtgctgcacctgctgactcaactgcatgacccgctcgt =     5\n+        34: aatacgggcatgatgttactggggatagtctttctcctgt =     5\n+        35: aatgccggcagagcagcggcttgctgaacgtcagcgaaaa =     6\n+        36: aatggctgctgccagtggcgctttgtcgtgtctttccggt =     5\n+        37: aatgggagagcgccttcgtggaattgcagaagatgttttc =     6\n+        38: aattagccaccatgccgcgcatctcacacccgcgcatgca =     5\n+        39: aattgcagaagatgttttcgcttacgcaacgggacaacgc =     6\n+        40: acaacgcgatgctgaacgagcgggtcatgcagttgagtca =     7\n+        41: acaacggtttcgtttttcgctgacgttcagcaagccgctg =     5\n+        42: acaggagggacgtgttaaagctggcagggaatacgggcat =     7\n+        43: acatcatgcccgtattccctgccagctttaacacgtccct =     6\n+        44: accaatggctgctgccagtggcgctttgtcgtgtctttcc =     5\n+        45: accagctttccagggatttcaacaacggtttcgtttttcg =     5\n+        46: accatcaccgcccccggtatcccgctcatgctggcgaacc =     6\n+        47: acccaccacggatatccgtaaggacagatgcagcaagtgc =     5\n+        48: acccgctcgttcagcatcgcgttgtcccgttgcgtaagcg =     7\n+        49: accgcccccggtatcccgctcatgctggcgaacccagaca =     6\n+        50: accgttgcaccttcccggttacgctccctccccagaaggt =     5\n+        51: accgttgttgaaatccctggaaagctggttgcgtgaaaag =     5\n+        52: acctgctgactcaactgcatgacccgctcgttcagcatcg =     7\n+        53: acgaaaccgttgttgaaatccctggaaagctggttgcgtg =     5\n+        54: acgcaacgggacaacgcgatgctgaacgagcgggtcatgc =     7\n+        55: acgccatagaggcggatataaggggaatgccggcagagca =     7\n+        56: acgcgatgctgaacgagcgggtcatgcagttgagtcagca =     7\n+        57: acggcatctttagtttgcacgataaccatatcttcaattc =     5\n+        58: acgggacaacgcgatgctgaacgagcgggtcatgcagttg =     7\n+        59: acgggcatgatgttactggggatagtctttctcctgttta =     5\n+        60: acggttcaccatgaaccataatatgatttctaagaatacc =     5\n+        61: acggtttcgtttttcgctgacgttcagcaagccgctgctc =     5\n+        62: acgtcagcgaaaaacgaaaccgttgttgaaatccctggaa =     5\n+        63'..b'agtcagcaggtgcagcacttga =     5\n+       246: gatataaggggaatgccggcagagcagcggcttgctgaac =     7\n+       247: gatgttttcgcttacgcaacgggacaacgcgatgctgaac =     6\n+       248: gattaccttctggggagggagcgtaaccgggaaggtgcaa =     5\n+       249: gatttcaacaacggtttcgtttttcgctgacgttcagcaa =     5\n+       250: gcacaggagggacgtgttaaagctggcagggaatacgggc =     7\n+       251: gcagagcagcggcttgctgaacgtcagcgaaaaacgaaac =     5\n+       252: gcatctttagtttgcacgataaccatatcttcaattccaa =     5\n+       253: gcatgatgttactggggatagtctttctcctgtttaccgc =     5\n+       254: gccaaaattagccaccatgccgcgcatctcacacccgcgc =     5\n+       255: gccactcgcggtaaacaggagaaagactatccccagtaac =     5\n+       256: gccatagaggcggatataaggggaatgccggcagagcagc =     7\n+       257: gcccccggtatcccgctcatgctggcgaacccagacacgc =     6\n+       258: gccttcgtggaattgcagaagatgttttcgcttacgcaac =     6\n+       259: gcgaaaacatcttctgcaattccacgaaggcgctctccca =     6\n+       260: gcgccttcgtggaattgcagaagatgttttcgcttacgca =     6\n+       261: gcgggtcatgcagttgagtcagcaggtgcagcacttgagc =     5\n+       262: gctcgttcagcatcgcgttgtcccgttgcgtaagcgaaaa =     7\n+       263: gctctgccggcattccccttatatccgcctctatggcgta =     7\n+       264: gctgaacgtcagcgaaaaacgaaaccgttgttgaaatccc =     5\n+       265: gctgttacttcccgctctaaagtagaagtaacatgaatca =     5\n+       266: gcttctacgcgattaagtaactcagtaaaacctgagttaa =     5\n+       267: gcttgctgaacgtcagcgaaaaacgaaaccgttgttgaaa =     5\n+       268: gctttccagggatttcaacaacggtttcgtttttcgctga =     5\n+       269: ggaatgccggcagagcagcggcttgctgaacgtcagcgaa =     6\n+       270: ggaattgcagaagatgttttcgcttacgcaacgggacaac =     6\n+       271: ggacgtgttaaagctggcagggaatacgggcatgatgtta =     5\n+       272: ggagaaagactatccccagtaacatcatgcccgtattccc =     5\n+       273: ggatataaggggaatgccggcagagcagcggcttgctgaa =     7\n+       274: ggatttcaacaacggtttcgtttttcgctgacgttcagca =     5\n+       275: ggcagagcagcggcttgctgaacgtcagcgaaaaacgaaa =     6\n+       276: ggcatctttagtttgcacgataaccatatcttcaattcca =     5\n+       277: ggcgatgatgcagcacaggagggacgtgttaaagctggca =     5\n+       278: ggcttctacgcgattaagtaactcagtaaaacctgagtta =     5\n+       279: ggcttgctgaacgtcagcgaaaaacgaaaccgttgttgaa =     5\n+       280: gggaatgccggcagagcagcggcttgctgaacgtcagcga =     6\n+       281: gggcatgatgttactggggatagtctttctcctgtttacc =     5\n+       282: gggtcatgcagttgagtcagcaggtgcagcacttgagcga =     5\n+       283: ggttcaccatgaaccataatatgatttctaagaataccaa =     5\n+       284: gtattccctgccagctttaacacgtccctcctgtgctgca =     6\n+       285: gtattcttagaaatcatattatggttcatggtgaaccgta =     5\n+       286: gtgagagcgctgttacttcccgctctaaagtagaagtaac =     5\n+       287: gtggaattgcagaagatgttttcgcttacgcaacgggaca =     6\n+       288: gttgtacgccatagaggcggatataaggggaatgccggca =     5\n+       289: gtttttcgctgacgttcagcaagccgctgctctgccggca =     6\n+       290: taaggggaatgccggcagagcagcggcttgctgaacgtca =     6\n+       291: tactttagagcgggaagtaacagcgctctcacttacaaaa =     5\n+       292: tattcttagaaatcatattatggttcatggtgaaccgtaa =     5\n+       293: tcaccgcccccggtatcccgctcatgctggcgaacccaga =     6\n+       294: tcatcttttcacgcaaccagctttccagggatttcaacaa =     5\n+       295: tcccgttgcgtaagcgaaaacatcttctgcaattccacga =     6\n+       296: tcgcgttgtcccgttgcgtaagcgaaaacatcttctgcaa =     6\n+       297: tcgttcagcatcgcgttgtcccgttgcgtaagcgaaaaca =     7\n+       298: tctactttagagcgggaagtaacagcgctctcacttacaa =     5\n+       299: tctgccggcattccccttatatccgcctctatggcgtaca =     7\n+       300: tcttcatcttttcacgcaaccagctttccagggatttcaa =     5\n+       301: tgaacgagcgggtcatgcagttgagtcagcaggtgcagca =     5\n+       302: tgaagtaccgttgcaccttcccggttacgctccctcccca =     5\n+       303: tgagagcgctgttacttcccgctctaaagtagaagtaaca =     5\n+       304: tgatgcagcacaggagggacgtgttaaagctggcagggaa =     5\n+       305: tgcacctgctgactcaactgcatgacccgctcgttcagca =     7\n+       306: tgcatttacggttcaccatgaaccataatatgatttctaa =     5\n+       307: tggaattgcagaagatgttttcgcttacgcaacgggacaa =     6\n+       308: tgtaagtgagagcgctgttacttcccgctctaaagtagaa =     5\n+       309: tgttgaaatccctggaaagctggttgcgtgaaaagatgaa =     5\n'