Repository 'bio_hansel'
hg clone https://toolshed.g2.bx.psu.edu/repos/nml/bio_hansel

Changeset 0:9407a9eaad22 (2017-09-27)
Next changeset 1:ba271365095e (2017-10-11)
Commit message:
planemo upload for repository https://github.com/phac-nml/bio_hansel commit c35a4502f95471e22dea2d77f36b92c47424cecc
added:
bio_hansel.xml
test-data/SRR1002850_SMALL.fasta
b
diff -r 000000000000 -r 9407a9eaad22 bio_hansel.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/bio_hansel.xml Wed Sep 27 11:50:41 2017 -0400
[
b'@@ -0,0 +1,252 @@\n+<tool id="bio_hansel" name="Hansel - Salmonella Subtyping" version="0.1.0">\n+    <requirements>\n+        <requirement type="package" version="0.1.0">bio_hansel</requirement>\n+        <requirement type="package" version="17.2.0">attrs</requirement>\n+    </requirements>\n+    <command detect_errors="exit_code"><![CDATA[\n+        ## Variables to be set.\n+        #set outputResultsName = \'results.tab\'\n+        #set outputMatchResultsName = \'match_results.tab\'\n+        #set verboseDebuggingLevel = \'-vvv\'\n+\n+        ## Preparing file input.\n+        #if $single_or_paired.type == "paired":\n+\n+            ln -s \'$single_or_paired.fastq_input1\' fast_q_paired_1.fastq &&\n+            ln -s \'$single_or_paired.fastq_input2\' fast_q_paired_2.fastq &&\n+\n+        #elif $single_or_paired.type == "collection":\n+\n+            ln -s \'$single_or_paired.fastq_input1.forward\' fast_q_paired_1.fastq &&\n+            ln -s \'$single_or_paired.fastq_input1.reverse\' fast_q_paired_2.fastq &&\n+\n+        #elif $single_or_paired.type == "single":\n+\n+            #if $single_or_paired.fastq_input1.is_of_type(\'fastqsanger\') or $single_or_paired.fastq_input1.is_of_type(\'fastq\'):\n+                ln -s \'$single_or_paired.fastq_input1\' fast_q_single.fastq &&\n+            #end if\n+\n+            #if $single_or_paired.fastq_input1.is_of_type(\'fasta\'):\n+                ln -s \'$single_or_paired.fastq_input1\' fast_a_single.fasta &&\n+            #end if\n+\n+        #end if\n+\n+\n+        ## Checking for custom scheme.\n+        #if $type_of_scheme.scheme_type == "custom":\n+            #if $type_of_scheme.scheme_input.is_of_type(\'fasta\'):\n+                ln -s \'$type_of_scheme.scheme_input\' custom_scheme.fasta &&\n+            #end if\n+        #end if\n+\n+\n+        ## Start the actual command here\n+        hansel \n+\n+\n+        ## Select the scheme\n+        -s\n+\n+        #if $type_of_scheme.scheme_type == "heidelberg":\n+            heidelberg\n+        #elif $type_of_scheme.scheme_type == "enteritidis":\n+            enteritidis\n+        #elif $type_of_scheme.scheme_type == "custom":\n+            custom_scheme.fasta\n+        #end if\n+\n+\n+        #if $kmer_min\n+            --min-kmer-freq $kmer_min\n+        #end if\n+\n+        #if $kmer_max\n+            --max-kmer-freq $kmer_max\n+        #end if\n+\n+        ## Adding more parameters to the command.\n+        $verboseDebuggingLevel -t "\\${GALAXY_SLOTS:-1}" -o $outputResultsName -O $outputMatchResultsName \n+\n+\n+        ## Entering the file inputs.\n+        #if $single_or_paired.type == "single":\n+\n+            #if $single_or_paired.fastq_input1.is_of_type(\'fastqsanger\') or $single_or_paired.fastq_input1.is_of_type(\'fastq\'):\n+                fast_q_single.fastq\n+            #end if\n+\n+            #if $single_or_paired.fastq_input1.is_of_type(\'fasta\'):\n+                fast_a_single.fasta\n+            #end if\n+\n+        #else\n+\n+            ##use -p to declare using two files.\n+            -p fast_q_paired_1.fastq  fast_q_paired_2.fastq\n+\n+        #end if\n+    ]]></command>\n+    <inputs>\n+        <conditional name="single_or_paired">\n+            <param name="type" type="select" label="Specify the read type.">\n+                <option value="single">Single-end Data</option>\n+                <option value="paired">Paired-end Data</option>\n+                <option value="collection">Collection Paired-end Data</option>\n+            </param>\n+            <when value="single">\n+                <param name="fastq_input1" type="data" format="fastqsanger, fastq, fasta" label="Single end read file(s)"/>\n+            </when>\n+            <when value="paired">\n+                <param name="fastq_input1" type="data" format="fastqsanger, fastq" label="Forward paired-end read file"/>\n+                <param name="fastq_input2" type="data" format="fastqsanger, fastq" label="Reverse paired-end read file"/>\n+            </when>\n+            <when value="collection">\n+                <param name="fastq_input1" type="data_collection" label='..b'tch    gapopen qstart  qend    sstart  send    evalue  bitscore    qlen    slen    seq coverage    is_trunc    refposition subtype is_pos_tile sample  file_path   scheme\n+\n+    775920-2.2.2.2  NODE_2_length_512016_cov_46.4737_ID_3   100.0   33  0   0   1   33  474875  474907  2.0000000000000002e-11  62.1    33  512016  GTTCAGGTGCTACCGAGGATCGTTTTTGGTGCG   1.0 False   775920  2.2.2.2 \n+\n+    True    SRR1002850  SRR1002850.fasta   heidelberg\n+\n+    negative3305400-2.1.1.1 NODE_3_length_427905_cov_48.1477_ID_5   100.0   33  0   0   1   33  276235  276267  2.0000000000000002e-11  62.1    33  427905  CATCGTGAAGCAGAACAGACGCGCATTCTTGCT   1.0 False   \n+\n+    negative3305400 2.1.1.1 False   SRR1002850  SRR1002850.fasta   heidelberg\n+\n+    negative3200083-2.1 NODE_3_length_427905_cov_48.1477_ID_5   100.0   33  0   0   1   33  170918  170950  2.0000000000000002e-11  62.1    33  427905  ACCCGGTCTACCGCAAAATGGAAAGCGATATGC   1.0 False   \n+\n+    negative3200083 2.1 False   SRR1002850  SRR1002850.fasta   heidelberg\n+\n+    negative3204925-2.2.3.1.5   NODE_3_length_427905_cov_48.1477_ID_5   100.0   33  0   0   1   33  175760  175792  2.0000000000000002e-11  62.1    33  427905  CTCGCTGGCAAGCAGTGCGGGTACTATCGGCGG   1.0 False   \n+\n+    negative3204925 2.2.3.1.5   False   SRR1002850  SRR1002850.fasta   heidelberg\n+\n+    negative3230678-2.2.2.1.1.1 NODE_3_length_427905_cov_48.1477_ID_5   100.0   33  0   0   1   33  201513  201545  2.0000000000000002e-11  62.1    33  427905  AGCGGTGCGCCAAACCACCCGGAATGATGAGTG   1.0 False   \n+\n+    negative3230678 2.2.2.1.1.1 False   SRR1002850  SRR1002850.fasta   heidelberg\n+\n+    negative3233869-2.1.1.1.1   NODE_3_length_427905_cov_48.1477_ID_5   100.0   33  0   0   1   33  204704  204736  2.0000000000000002e-11  62.1    33  427905  CAGCGCTGGTATGTGGCTGCACCATCGTCATTA   1.0 False   \n+\n+    [Next 196 lines omitted.]\n+\n+\n+Analysis of a single FASTQ readset\n+----------------------------------\n+\n+Contents of ``results.tab``:\n+\n+    sample  scheme  subtype all_subtypes    tiles_matching_subtype  are_subtypes_consistent inconsistent_subtypes   n_tiles_matching_all    n_tiles_matching_all_total  n_tiles_matching_positive   n_tiles_matching_positive_total n_tiles_matching_subtype    n_tiles_matching_subtype_total  file_path\n+\n+    SRR5646583  heidelberg  2.2.1.1.1.1 2; 2.2; 2.2.1; 2.2.1.1; 2.2.1.1.1; 2.2.1.1.1.1  1983064-2.2.1.1.1.1; 4211912-2.2.1.1.1.1    True        202 202 20  20  2   2   SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger\n+\n+\n+Contents of ``match_results.tab``:\n+\n+    seq freq    sample  file_path   tilename    is_pos_tile subtype refposition is_kmer_freq_okay   scheme\n+\n+    ACGGTAAAAGAGGACTTGACTGGCGCGATTTGC   68  SRR5646583 SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger    21097-2.2.1.1.1 True    2.2.1.1.1   21097   True    heidelberg\n+\n+    AACCGGCGGTATTGGCTGCGGTAAAAGTACCGT   77  SRR5646583 SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger    157792-2.2.1.1.1    True    2.2.1.1.1   157792  True    heidelberg\n+\n+    CCGCTGCTTTCTGAAATCGCGCGTCGTTTCAAC   67  SRR5646583 SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger    293728-2.2.1.1  True    2.2.1.1 293728  True    heidelberg\n+\n+    GAATAACAGCAAAGTGATCATGATGCCGCTGGA   91  SRR5646583 SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger    607438-2.2.1    True    2.2.1   607438  True    heidelberg\n+\n+    CAGTTTTACATCCTGCGAAATGCGCAGCGTCAA   87  SRR5646583 SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger    691203-2.2.1.1  True    2.2.1.1 691203  True    heidelberg\n+\n+    CAGGAGAAAGGATGCCAGGGTCAACACGTAAAC   33  SRR5646583 SRR5646583_forward.fastqsanger; SRR5646583_reverse.fastqsanger    944885-2.2.1.1.1    True    2.2.1.1.1   944885  True    heidelberg\n+\n+    [Next 200 lines omitted.]\n+    \n+\n+Wrapper written by Matthew Gopez at the Public Health Agency of Canada, National Microbiology Laboratory.\n+    ]]></help>\n+    <citations>\n+    </citations>\n+</tool>\n\\ No newline at end of file\n'
b
diff -r 000000000000 -r 9407a9eaad22 test-data/SRR1002850_SMALL.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/SRR1002850_SMALL.fasta Wed Sep 27 11:50:41 2017 -0400
b
b'@@ -0,0 +1,766 @@\n+>NODE_384_length_56_cov_34125_ID_767\n+CTAGCCTGATTCAGTTTTATCTTTCGTCGGTGAGCGACGGTACTTTCTGACCTCGG\n+>NODE_383_length_56_cov_335_ID_765\n+CAATTGTGATAGCGTCGCGTGTGACGGTATTACAATTACACATCCTGCCCCGAAAA\n+>NODE_382_length_56_cov_201_ID_763\n+CCGTCGCTCACCGACGAAAGATAAAACGGAATCAGGCTAGCGGGGTGTCGGGGTGC\n+>NODE_377_length_57_cov_34097.5_ID_753\n+AAACCGAGGTCAGAAAGTACCGTCGCTCACCGACGAAAGATAAAACTGAATCAGGCT\n+>NODE_376_length_57_cov_200_ID_751\n+GAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAGTGCCGCCAGGCAAATTCTTT\n+>NODE_375_length_58_cov_19568.3_ID_749\n+TAATACCGTCACACGCGACGCTATCACAATTGCCATCTGGTCAGGGGTGCACCCCGAC\n+>NODE_374_length_58_cov_29262.3_ID_747\n+GCACCCCGACACCCCGCTAGCCTGATTCAGTTTTATCTTTCGTCGGTGAGCGACGGTA\n+>NODE_372_length_59_cov_291.25_ID_743\n+TAGCCTGATTCAGTTTTATCTTTCGTCGGTGAGCGACGGGACTTTCTGACCTCGGTTTT\n+>NODE_369_length_60_cov_34761.2_ID_737\n+GCGGTGGCCGAAACCCTGAAATATTCCACGAAGCCGGCTGACATGGTGGCCGATCCTGAA\n+>NODE_367_length_60_cov_167.2_ID_733\n+CCGAGGGCGACGATGACGGTTCACGAATCGCGTTCGAGTGGAAACCCGAGGTCAGAAAGT\n+>NODE_365_length_61_cov_228.5_ID_729\n+TCGCGTTCGAGTGGAAAACCGAGGTCAGAAAGTCCCGTCGCTCACCGACGAAAGATAAAAC\n+>NODE_360_length_62_cov_66.7143_ID_719\n+CAGGGGTGCACCCCGACACCCCGCTAGCCTGATTCAGTTTTATCTTTCGTCGGTGGGCGACG\n+>NODE_357_length_63_cov_202.375_ID_713\n+GGGGTGTGGGGGTGCACCCCTGACCAGATGGCAATTGTGATAGCGTCGCGTGTGACGGTATTA\n+>NODE_356_length_63_cov_142.75_ID_711\n+TGTAATTGTAATAACGTCACACGCGACGCTATCACAATTGCCATCTGGTCAGGGGGGCACCCC\n+>NODE_355_length_64_cov_127.667_ID_709\n+CGTCGCTCCCCGACGAAAGATAAAACTGAATCAGGCTAGCGGGGTGTCGGGGTGCCCCCCTGAC\n+>NODE_352_length_64_cov_30.6667_ID_703\n+CGTCGCTCCCCGACGAAAGATAAAACTGAATCAGGCTAGGGGGGTGTCGGGGTGCACCCCTGAC\n+>NODE_351_length_64_cov_104.444_ID_701\n+ATGATATTTCCGAGGGCGACGATGACGGTTCACGAATCGCGTTCGAGTGGAAACCCGAGGTCAG\n+>NODE_349_length_64_cov_29.3333_ID_697\n+GTCAGGGGGGCACCCCGACACCCCCCTAGCCTGATTCAGTTTTATCTTTCGTCGGGGAGCGACG\n+>NODE_347_length_64_cov_32_ID_693\n+GTCAGGGGGGCACCCCGACACCCCGCTAGCCTGATTCAGTTTTTTCTTTCGTCGGGGAGCGACG\n+>NODE_346_length_65_cov_38729.2_ID_691\n+CTGACCTCGGTTTTCCACTCGAACGCGATTCGTGAACCGTCATCGTCGCCCTCGGAAATATCATC\n+>NODE_342_length_66_cov_28.1818_ID_683\n+CTGACCTCGGGTTTCCACTCGAACGCGATTCGTGAACCGTCATCGTCGCCCTCGGGAATATCATCA\n+>NODE_340_length_67_cov_71.3333_ID_679\n+CCGACGAAAGATAAAACTGAATCAGGCTAGGGGGGTGTCGGGGTGCCCCCCTGACCAGATGGCAATT\n+>NODE_336_length_67_cov_56_ID_671\n+CTTGAATATGGCGGTGAGGGGGGGATTGACTCGCTTTGCTCGCCCCTTCGGGGCAGCCCGTTCGCTG\n+>NODE_335_length_68_cov_32300.2_ID_669\n+CCCCGACACCCCGCTAGCCTGATTCAGTTTTATCTTTCGTCGGTGAGCGACGGTACTTTCTGACCTCG\n+>NODE_333_length_69_cov_126.429_ID_665\n+GGTGTCGGGGTGCCCCCCTGACCAGAGGGCAATTGTGATAGCGTCGCGTGTGACGGTATTACAATTACA\n+>NODE_332_length_69_cov_101.286_ID_663\n+GGTGTCGGGGTGCCCCCCTGACCAGATGGCATTTGTGATAGCGTCGCGTGTGACGGTATTACAATTACA\n+>NODE_327_length_71_cov_31.375_ID_653\n+GCACCCCGACACCCCCCTAGCCTGATTCAGTTTTATCTTTCGTCGGGGAGCGACGGTACTTTCTGACCTCG\n+>NODE_326_length_71_cov_108.125_ID_651\n+CCGAGGTCAGAAAGTCCCGTCGCTCCCCGACGAAAGATAAAACTGAATCAGGCTAGCGGGGTGTCGGGGTG\n+>NODE_324_length_71_cov_35.75_ID_647\n+TGTAATTGTAATACCCTCACACGCGACGCTATCACAATTGCCATCTGGTCAGGGGGGCACCCCGACACCCC\n+>NODE_318_length_74_cov_355.895_ID_635\n+GTAGGCGTTATTGGCGCAGCCAGTTTGGACACGGACAGCGCGCAAAAACCGGAGCGTACACGTAGTACGTGAGG\n+>NODE_316_length_76_cov_59.5714_ID_631\n+CGTCGCTCACCGACGAAAGAAAAAACTGAATCAGGCTAGGGGGGTGTCGGGGTGCACCCCTGACCAGATGGCAATT\n+>NODE_315_length_76_cov_76.8571_ID_629\n+CGTCGCTCACCGACGAAAGATAAAACTGAATCAGGCTAGCGGGGTGTGGGGGTGCACCCCTGACCAGATGGCAATT\n+>NODE_314_length_76_cov_172.333_ID_627\n+AATTGCCATCTGGTCAGGGGTGCACCCCGACACCCCGCTAGCCTGATTCCGTTTTATCTTTCGTCGGTGAGCGACG\n+>NODE_313_length_76_cov_35231.5_ID_625\n+AAACGGCAAAATGATTTACTCGCCGTTTGTGTGGGTTGTTTTGTGGGGTGTCAATGGGGTTGTGGTCTTTTTTGTG\n+>NODE_312_length_76_cov_28.9048_ID_623\n+CGTCGCTCACCGACGAAAGAAAAAACTGAATCAGGCTAGGGGGGTGTCGGGGTGCCCCCCTGACCAGATGGCAATT\n+>NODE_310_length_79_cov_105.5_ID_619\n+CTCAAAAAACCGAGAAACGGCGGGGATTTTCGGGACAGGATGTGTAATTGTAATAACGTCACACGCGACGCTATCACAA\n+>NODE_306_length_82_cov_64.2963_ID_611\n+TATCACAATTGCCATCTGGTCAGGGGGGCACCCCGACACCCCGCTAGCCTGATTCCGTTTTATCTTTCGTCGGTGAGCGACG\n+>NODE_294_length'..b'GGATGGCCGCGTCTTTGAAGCGCTACTGGTTGGCTCCGATTCGCTTACCGATCTGGCGGTGCTGAAGATCAACGCCACTGGCGGGCTGCCTACCATCCCGATTAATACAAAGCGTACACCGCATATTGGCGACGTCGTACTGGCTATCGGCAACCCATATAATCTGGGACAGACCATTACCCAGGGGATCATCAGCGCAACGGGTCGTATCGGCCTGAACCCGACGGGGCGACAGAATTTTCTCCAGACCGACGCCTCGATTAACCACGGTAATTCCGGCGGCGCGCTGGTCAACTCGTTAGGCGAACTGATGGGGATCAACACCCTCTCTTTTGATAAGAGTAACGATGGTGAAACGCCGGAAGGCCTTGGTTTTGCGATTCCCTTCCAGCTAGCCACGAAAATTATGGATAAGCTTATCCGCGACGGTCGTGTGATTCGCGGCTATATCGGTATTGGCGGACGAGAAATCGCGCCGCTGCACGCGCAGCAGGGTAGCGGCATGGACCCGATTCAGGGCATTGTCGTTAATGAAGTGACGCCAAACGGCCCCGCCGCGCTTGCCGGTATTCAGGTTAATGATTTGATTATTTCGGTCAATAATAAACCCGCCGTGTCCGCGCTGGAGACAATGGATCAGGTGGCGGAAATCCGCCCGGGCTCCGTCATTCCGGTCGTGGTAATGCGGGATGATAAGCAACTCACGTTCCAGGTGACGGTGCAGGAATACCCGGCGTCGAACTAAAACGACGCGGCTTTTGCCTGCGTCGTTCGACGGCTAGTCAATAAAGACTGGCCCTGCCGGAAAATGCGCCAGGAAACTCGTGAAGATCAGCAGAGCGGGGAGACTTTTCTCCCCGACGGCTTCTTCACGCAGCGCCTTATCCAGTGCAAGCAGGCTGTGCTCCCCCCCTTCCCCCAAAATGCCGCCCATCCAGAGCAGTCGTTGCGCCTGTTGACGCGCATACTGCGATGCGCTTTTCCAGGCGACCAGATGCAGCAACGTTTGTAACAGCAGGACCCGAGACGCGTGACCATCGCGTTGCAATCGCCTGAATACGGGTAATGCAATACTACGTACCGTCGAGAAACCGCTCGCCGCCTCCCCCTCAATCGGCGATAAACCAACAGCCGCCTGTTTGACGAAATCGCATAACTGATTGTGTGTTAATGGCCGCTCGCCTGCGCTCATATACCCTGCGGCAAAACAGAGTAACCCGATACCCTGAGCCCAATCGGTCGATGGCAGAATCCCCTGCGGAACATTCGGTGCCGACGTCTGTGCGGCACAAAGCCCCTCGTCGTAACACTGTGCAAAACGTTGCTCATATATTGCTGAGTTGATGCCAGACCACGGCAGTAAATGCAGCGCTGGCTCAAGCTCAGTTATCTGATTCACCGTCTCAACAACATAACGCAAGATGGTTTCTCTGGATGGCATGGATGTACCTTGCCGAAACGCAACCATTTCACTCATTACTGTGTCCTCATTTTGCGGGAGCAGAATCACTGCTCCCGATTAGATTGGATTACATTGCCAGCACGTACTTCAGCATCACGCCCGCCGCAATGGCCGAACCGATAACCCCCGCCACGTTCGGGCCCATCGCGTGCATCAGCAGGAAGTTCTGCGGGTCCGACTCCAGGCCCACCTTGTTCGATACGCGGGCCGCCATCGGCACCGCCGACACCCCCGCCGAACCGATAAGCGGGTTGATTTTGTTCTTACTGCACAGGTTCAGCAGCTTCGCCATCAGCACCCCGGCGGCAGTCCCGATACCAAAGGCAATCACACCCAGCAGCAGAATGCCCAGCGTCTGCGGCTGCAGGAACTTGTCCGCCACCAGCTTCGCCCCTACCGACAGCCCGAGGAAAATGGTGACGATGTTGATCAGCCCGTTCTGCACTGTGTCGCTCAGGCGCTCCACCACGCCGCTTTCACGCATCAGATTGCCGAAGCAGAACATCCCCAGCAGCGGCGCGGCGTCCGGCAGCAGCAGCGCCACCAGTAGTAACAGCACTACCGGGAAGAGGATTTTTTCCCGCTTACTCACCGTGCGCAGCTGCACCATGCGGATTTTCCGCTCCGTCTCCGTGGTCAGCGCGCGCATAATCGGCGGCTGGATTAACGGCACCAGCGCCATATACGAGTACGCCGCCACCGCGATGGCCCCCAGCAGCTCCGGCGCCAGTTTGCCCGACAGATAAATGGCCGTCGGACCGTCCGCGCCGCCGATGATGCCTATCGCCGCCGCCTGCGGCAGGGTGAAGGAAATCAGGCCGAAGTAATTCAGCGTCAGCGCCCCCAGCACCGTGGCGAAGATGCCGAACTGCGCCGCCGCGCCGAGCAGCAGGGTACGCGGGTTGGCCAGCAGCGGGCCGAAGTCGGTCATCGCGCCGACGCCCATAAAGATGACCAGCGGCGCGACGCCGGAGCCAATCGCCACTTTATAGAACAGCGCCAGCACCCCCGGCGTGTAGCCCATGTCCACCGCCAGGTTCTCCATCTGGCTCTGAACGGAGGGCAGCGCCAGCGCTAACGCCTCCTTGATGGCGTGCACGTCCGGCGCGCAGTTAAGCTTCGCGGCAATCACCGCCAGCTGCCCGGCGTCGTGGTGCGCCAGCAGGCTTTCCAGCGCGGTCAGCGCCATGCCCGCTTCCGGGATGTTGGAGAGCAGGCCGCCGAAGCCAATCGGCAACAGCAACAACGGCTCGAACTTTTTCGCAATCGCCAGCCACAGCAGCAGCAGGCTGACCAGCAGCATGATGGCCTGGCCTGCGCCAAGGTGCATCAGCCCCATGCCCTGAAGCAGGGCGTTCAGACTTTCCATTCGTCCTCTCCGTTACGCAAGCTGCATCAGGGTGTCGCCCACCGCCACCGCATCCCCGGCTTTCACCGCGATGCCGCGTACCGTCCCGGCCTGCGCGGCGCGGATTTCGGTTTCCATTTTCATGGCTTCCAGAATCAGCAGCACATCGCCTTCGGCCACCGTCTGGCCTTCCGTCGCCACCACTTTCCAGATATTGCCCGCCAGCGGCGCGCTGACCGGGGTGCCCGCGCCCGCCGGGGCGGCGGCAGGCGCCGCGGTCTGAACCGGCGCAGAACTGGCCGTGGCCAGCTGGCTGATATCGCCGCCGTCGCTGACTTTCACCACAAAGGCTTTGCCTTCCACTTCGACGGTGTAGATACCGGAAGCGGCGGGTTTTGCGGCTTTTTCCTCTTTCACGGCAGGTTTCGCGGCTTCCGCCTGCGGCAGTGGCTCAAACGCCGCCGGGTTATGGCGGTTTTCGAGGAATTTCAGCCCGACCTGCGGGAACAGTGCCACGGTGAGGACGTCATCAATGGCGTTCTCTGCAAGCCGGATGCCTTTTTCTGCCGCCTGGCGTTTAATGTCCGCTTCCAGTTGCGCCAGTTCCGGTTTCAGTAAATCCGCCGGGCGGCAGGTGACGGCCTCTTCGCCGTCCAGCACTCTGGCCTGTAACGCCGCGTTGACCGGCGCTGGCGTGTGGCCGTATTCGCCCTTCAGTATCCCGGCGGTTTCTTTGGCGATGGTTT\n+>NODE_216_length_189_cov_23403.1_ID_431\n+CCTGATTCAGTTTTATCTTTCGTCGGTGAGCGACGGTACTTTCTGACCTCGGTTTTTCACTCGAACGCGATTCGTGAACCGTCATCGTCGCCCTCGGAAATATCATCACCGATCACCATGTCGGCGTCAGTTTCCTGATCGAGTTTCAGCACGTCTTTGAGTGCGCCGCCGGTCGCCACAAAGCGGCGC\n+>NODE_222_length_167_cov_18131.6_ID_443\n+CACTCAAAGACGTGCTGAAACTCGATCAGGAAACTGACGCCGACATGGTGATCGGTGATGATATTTCCGAGGGCGACGATGACGGTTCACGAATCGCGTTCGAGTGAAAAACCGAGGTCAGAAAGTACCGTCGCTCCCCGACGAAAGATAAAACTGAATCAGGCTAG\n'