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

Changeset 11:c59d48774d03 (2024-03-05)
Previous changeset 10:65b93b623c77 (2022-02-03) Next changeset 12:d60c3f704da0 (2024-08-13)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/fastp commit 61d055a8abacf009e8afe2d72377e6769ad08f42
modified:
fastp.xml
macros.xml
test-data/bwa-mem-fastq2.fq
added:
test-data/bwa-mem-fastq2_too_long.fq
b
diff -r 65b93b623c77 -r c59d48774d03 fastp.xml
--- a/fastp.xml Thu Feb 03 17:20:32 2022 +0000
+++ b/fastp.xml Tue Mar 05 08:41:43 2024 +0000
[
b'@@ -1,10 +1,11 @@\n-<tool id="fastp" name="fastp" version="@WRAPPER_VERSION@+galaxy0">\n-    <description>- fast all-in-one preprocessing for FASTQ files</description>\n+<tool id="fastp" name="fastp" version="@TOOL_VERSION@+galaxy0" profile="23.1">\n+    <description>fast all-in-one preprocessing for FASTQ files</description>\n     <macros>\n         <import>macros.xml</import>\n     </macros>\n+    <expand macro="biotools" />\n     <requirements>\n-        <requirement type="package" version="@WRAPPER_VERSION@">fastp</requirement>\n+        <requirement type="package" version="@TOOL_VERSION@">fastp</requirement>\n     </requirements>\n     <version_command>fastp -v</version_command>\n     <command detect_errors="exit_code"><![CDATA[\n@@ -24,19 +25,19 @@\n     #set $in2_name = re.sub(\'[^\\w\\-\\s]\', \'_\', str("%s_%s" % ($single_paired.paired_input.name, "R2"))) + $ext\n     #set out1 = $output_paired_coll.forward\n     #set out2 = $output_paired_coll.reverse\n-    ln -s \'$in1\' \'$in1_name\' &&\n-    ln -s \'$in2\' \'$in2_name\' &&\n+    ln -sf \'$in1\' \'$in1_name\' &&\n+    ln -sf \'$in2\' \'$in2_name\' &&\n #else\n     #if $in1.is_of_type(\'fastq.gz\')\n         #set ext = \'.fastq.gz\'\n     #end if\n \n     #set $in1_name = re.sub(\'[^\\w\\-\\s]\', \'_\', str($in1.element_identifier)) + $ext\n-    ln -s \'$in1\' \'$in1_name\' &&\n+    ln -sf \'$in1\' \'$in1_name\' &&\n \n     #if str($single_paired.single_paired_selector) == \'paired\':\n         #set $in2_name = re.sub(\'[^\\w\\-\\s]\', \'_\', str("%s_R2" % $in2.element_identifier)) + $ext\n-        ln -s \'$in2\' \'$in2_name\' &&\n+        ln -sf \'$in2\' \'$in2_name\' &&\n     #end if\n #end if\n \n@@ -209,7 +210,6 @@\n #end if\n ]]></command>\n     <inputs>\n-\n         <conditional name="single_paired">\n             <param name="single_paired_selector" type="select" label="Single-end or paired reads">\n                 <option value="single" selected="true">Single-end</option>\n@@ -229,8 +229,8 @@\n                 </expand>\n                 <expand macro="global_trimming_options_paired" />\n             </when>\n-        <when value="paired_collection">\n-                <param name="paired_input" type="data_collection" collection_type="paired" format="fastq,fastq.gz" label="Select paired collection(s)"/>\n+            <when value="paired_collection">\n+                <param name="paired_input" type="data_collection" format="fastq,fastq.gz" label="Select paired collection(s)" collection_type="paired"/>\n                 <expand macro="adapter_trimming_options">\n                     <expand macro="adapter_sequence" read_number="2"/>\n                 </expand>\n@@ -255,7 +255,7 @@\n             <section name="length_filtering_options" title="Length filtering options" expanded="True">\n                 <param name="disable_length_filtering" argument="-L" type="boolean" truevalue="-L" falsevalue="" checked="false" label="Disable length filtering" help="Length filtering is enabled by default. If this option is specified, length filtering is disabled."/>\n                 <param name="length_required" argument="-l" type="integer" optional="true" label="Length required" help="Reads shorter than this value will be discarded. Default is 15."/>\n-                <param name="length_limit" argument="--length_limit" type="integer" optional="true" label="Maximum length" help="Reads longer than this value will be discarded. Default is 0 and means no limitation."/>\n+                <param argument="--length_limit" type="integer" optional="true" label="Maximum length" help="Reads longer than this value will be discarded. Default is 0 and means no limitation."/>\n             </section>\n \n             <section name="low_complexity_filter" title="Low complexity filtering options" expanded="True">\n@@ -287,17 +287,17 @@\n                     <option value="-x">Enable polyX tail trimming</option>\n                 </param>\n                 <when value="-x">\n-                    <param name="poly_x_min_len" argument="--poly_x_min_len" type="integer" optional="true" label="PolyX minimum length"\n-    '..b'   <param name="in1" ftype="fastq" value="R1.fq"/>\n             <param name="single_paired_selector" value="single"/>\n@@ -396,7 +396,7 @@\n             </section>\n             <output name="out1" ftype="fastq" file="out3.fq"/>\n         </test>\n-        <!-- Ensure paired-end fastq works -->\n+        <!-- 6. Ensure paired-end fastq works -->\n         <test expect_num_outputs="3">\n             <param name="in1" ftype="fastq" value="bwa-mem-fastq1.fq"/>\n             <param name="in2" ftype="fastq" value="bwa-mem-fastq2.fq"/>\n@@ -404,7 +404,7 @@\n             <output name="out1" ftype="fastq" file="out_bwa1.fq"/>\n             <output name="out2" ftype="fastq" file="out_bwa2.fq"/>\n         </test>\n-        <!-- Ensure paired-end UMI processing of Read 1 works -->\n+        <!-- 7. Ensure paired-end UMI processing of Read 1 works -->\n         <test expect_num_outputs="3">\n             <param name="in1" ftype="fastq" value="bwa-mem-fastq1.fq"/>\n             <param name="in2" ftype="fastq" value="bwa-mem-fastq2.fq"/>\n@@ -417,7 +417,7 @@\n             <output name="out1" ftype="fastq" file="out_bwa_umi_read1_1.fq"/>\n             <output name="out2" ftype="fastq" file="out_bwa_umi_read1_2.fq"/>\n         </test>\n-        <!-- Ensure paired-end UMI processing of Read 2 works -->\n+        <!-- 8. Ensure paired-end UMI processing of Read 2 works -->\n         <test expect_num_outputs="3">\n             <param name="in1" ftype="fastq" value="bwa-mem-fastq1.fq"/>\n             <param name="in2" ftype="fastq" value="bwa-mem-fastq2.fq"/>\n@@ -430,7 +430,7 @@\n             <output name="out1" ftype="fastq" file="out_bwa_umi_read2_1.fq"/>\n             <output name="out2" ftype="fastq" file="out_bwa_umi_read2_2.fq"/>\n         </test>\n-        <!-- Ensure JSON report output works -->\n+        <!-- 9. Ensure JSON report output works -->\n         <test expect_num_outputs="2">\n             <param name="in1" ftype="fastqsanger" value="R1.fq"/>\n             <param name="single_paired_selector" value="single"/>\n@@ -443,7 +443,7 @@\n                 </assert_contents>\n             </output>\n         </test>\n-        <!-- Ensure polyG trimming works -->\n+        <!-- 10. Ensure polyG trimming works -->\n         <test expect_num_outputs="2">\n             <param name="in1" ftype="fastq.gz" value="R1.fq.gz"/>\n             <param name="single_paired_selector" value="single"/>\n@@ -451,7 +451,7 @@\n             <param name="poly_g_min_len" value="10"/>\n             <output name="out1" ftype="fastq.gz" decompress="True" file="out1.fq.gz"/>\n         </test>\n-        <!-- Ensure polyX trimming works -->\n+        <!-- 11. Ensure polyX trimming works -->\n         <test expect_num_outputs="2">\n             <param name="in1" ftype="fastq.gz" value="R1.fq.gz"/>\n             <param name="single_paired_selector" value="single"/>\n@@ -460,13 +460,24 @@\n             <param name="poly_x_min_len" value="10"/>\n             <output name="out1" ftype="fastq.gz" decompress="True" file="out1.fq.gz"/>\n         </test>\n+        <!-- 12. Test fastq files with different length -->\n+        <test expect_exit_code="255" expect_failure="true">\n+            <param name="single_paired_selector" value="paired_collection"/>\n+            <param name="paired_input">\n+                <collection type="paired">\n+                    <element name="forward" value="bwa-mem-fastq1.fq" ftype="fastqsanger" />\n+                    <element name="reverse" value="bwa-mem-fastq2_too_long.fq" ftype="fastqsanger" />\n+                </collection>\n+            </param>\n+        </test>\n     </tests>\n     <help><![CDATA[\n .. class:: infomark\n \n **What it does**\n \n-fastp_ is a tool designed to provide fast all-in-one preprocessing for FASTQ files. This tool is developed in C++ with multithreading supported to afford high performance.\n+fastp_ is a tool designed to provide fast all-in-one preprocessing for FASTQ files. This tool is developed in C++ with multithreading supported to\n+afford high performance.\n \n *Features*\n \n'
b
diff -r 65b93b623c77 -r c59d48774d03 macros.xml
--- a/macros.xml Thu Feb 03 17:20:32 2022 +0000
+++ b/macros.xml Tue Mar 05 08:41:43 2024 +0000
b
@@ -1,16 +1,18 @@
 <macros>
-    <token name="@WRAPPER_VERSION@">0.23.2</token>
-
+    <token name="@TOOL_VERSION@">0.23.4</token>
+    <xml name="biotools">
+        <xrefs>
+            <xref type="bio.tools">
+                fastp
+            </xref>
+        </xrefs>
+    </xml>
     <xml name="adapter_trimming_options">
         <section name="adapter_trimming_options" title="Adapter Trimming Options" expanded="False">
             <param name="disable_adapter_trimming" argument="-A" type="boolean" truevalue="-A" falsevalue="" checked="false" label="Disable adapter trimming" help="Adapter trimming is enabled by default. If this option is specified, adapter trimming is disabled."/>
             <expand macro="adapter_sequence" read_number="1" />
             <yield />
         </section>
-        <section name="global_trimming_options" title="Global Trimming Options" expanded="False">
-            <param name="trim_front1" argument="-f" type="integer" optional="true" label="Trim front for input 1" help="Trimming how many bases in front for read1, default is 0."/>
-            <param name="trim_tail1" argument="-t" type="integer" optional="true" label="Trim tail for input 1" help="Trimming how many bases in tail for read1, default is 0."/>
-        </section>
     </xml>
 
     <xml name="global_trimming_options">
@@ -29,7 +31,7 @@
     </xml>
 
     <xml name="adapter_sequence" token_read_number="1">
-        <param name="adapter_sequence@READ_NUMBER@" argument="--adapter_sequence@READ_NUMBER@" type="text" optional="true" label="Adapter sequence for input @READ_NUMBER@" help="The adapter for read@READ_NUMBER@. For SE data, if not specified, the adapter will be auto-detected. For PE data, this is used if R1/R2 are found not overlapped.">
+        <param argument="--adapter_sequence@READ_NUMBER@" type="text" optional="true" label="Adapter sequence for input @READ_NUMBER@" help="The adapter for read@READ_NUMBER@. For SE data, if not specified, the adapter will be auto-detected. For PE data, this is used if R1/R2 are found not overlapped.">
             <sanitizer>
                 <valid>
                     <add value="A"/>
@@ -46,7 +48,7 @@
     </xml>
 
     <xml name="poly_g_min_len">
-        <param name="poly_g_min_len" argument="--poly_g_min_len" type="integer" optional="true" label="PolyG minimum length"
+        <param argument="--poly_g_min_len" type="integer" optional="true" label="PolyG minimum length"
             help="The minimum length to detect polyG in the read tail. 10 by default."/>
     </xml>
 </macros>
b
diff -r 65b93b623c77 -r c59d48774d03 test-data/bwa-mem-fastq2.fq
--- a/test-data/bwa-mem-fastq2.fq Thu Feb 03 17:20:32 2022 +0000
+++ b/test-data/bwa-mem-fastq2.fq Tue Mar 05 08:41:43 2024 +0000
b
@@ -393,8 +393,4 @@
 @M01368:8:000000000-A3GHV:1:1114:2404:13066/2
 ATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCTATTATTTATCGCACCTACGTTCAATATTACAGGCGACCATACTTACTAAAGTGTGTTAATTAATTAATGCTTGTAGGACTGTCTCTTATACACATT
 +
-CCCCCFFFFFFCGGGGGGGGGGHHHHHHHHHFFHHHHGGGGGHFFFHHFHHHHHHHHHHHHHHHGFEGGGHGEDFCDFHGHFG@@DGGHHHHHHGGGGCGGGGGEHGGCGBB?CF99EGFGGFGG?D9CFFFF/BBFFFFFEF9BFFAFFFFEFFFFFFFFFFFFFFFFFFFFF.FFBBFFFFFFFFFFFF-9;;;BFFFFFB9BFBFBFABFFEFFFFFFFFFF::BFFBFFFF.9//;FFFFF/BFFB/
-@M01368:8:000000000-A3GHV:1:1114:9184:6959/2
-AAAGTGAACTGTATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCCTGTCTCTTA
-+
-CCCCBFFFFFFFGGGGGGGGGGHHHHHHHHHHHHHHHHHHHHHHHHHHGHHHHHHEHIHHGGGGHHHHHHHHHHHHHGHHHHHHHHGGGGGHHFHHHHHBGHHHHHHHHHHHHHHHHHGHHHHHGGGGGHHHHGHHHHHHHHHHHHHHHHHGHGHGHHGGGGCFFFFFFFFFFFFFFFFFFFFFFFFFF.CFFFFAF=D=EAEFFF0B:0AF-DAFBFFFFFFFFFBFFFFFFFFFFBFFFEFF9B900B0
+CCCCCFFFFFFCGGGGGGGGGGHHHHHHHHHFFHHHHGGGGGHFFFHHFHHHHHHHHHHHHHHHGFEGGGHGEDFCDFHGHFG@@DGGHHHHHHGGGGCGGGGGEHGGCGBB?CF99EGFGGFGG?D9CFFFF/BBFFFFFEF9BFFAFFFFEFFFFFFFFFFFFFFFFFFFFF.FFBBFFFFFFFFFFFF-9;;;BFFFFFB9BFBFBFABFFEFFFFFFFFFF::BFFBFFFF.9//;FFFFF/BFFB/
\ No newline at end of file
b
diff -r 65b93b623c77 -r c59d48774d03 test-data/bwa-mem-fastq2_too_long.fq
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/bwa-mem-fastq2_too_long.fq Tue Mar 05 08:41:43 2024 +0000
b
b'@@ -0,0 +1,400 @@\n+@M01368:8:000000000-A3GHV:1:1101:6911:8255/2\n+TCGCCTTACCGCTACTCACCCACGGCGGCCATCAGCCGATACTAAGTTTGGGGTATGGTGGGGGGGATAATGAATTAGGTTGTGGGGGAGGGTTTGTGGTTGAGAGAAACACAAAAAACAATCTTATATATGGGTAGTCGTTTTGTATTGGTTTTTTGTTTTGTTTGTGTTTTGAGTGTCGGTTTAGTTCGGTGTACTAGGGGGGGTGGATGGGGTCGGCTGGTGAGGGGGTCTTAGTGTATTGAGTGTGG\n++\n+1>11111@11111A111A100000000////011110//>>/12@1@22B/////1@>21/>>/-----9/;////9////--;-;-;-----;--------9/-/-///9-;-------9//////9/////-//-/9-;-;9--/////99-;--9-:-;----/---/-----////---9-/////--;A-//////---------9/-----;-----/-/-----;--;//////////9;///-\n+@M01368:8:000000000-A3GHV:1:1101:14518:9998/2\n+CATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCTATTATTTATCGCACCTACGTTCAATATTACAGGCGAACATACTTACTAAAGTGTGTTAATTAATTAATGCTTGTAGGACATAATAATAA\n++\n+CCCCCFCCCCCFGGGGGGGGGGHHHHHHHHHHHHHHHHGFHHHHGGGGGHGFHHHHHHHHHHHHHHHHHHHGHGGEHGGGGCGGGHHGGCGGGGGHHGHHHGGGGGGGG.BFFFGAGADFGAFDGFGGCFFF;DDFFFFFFFFFFFFFFFFFFEFFFFFFFFFFFBFFFFFFFFFFFFFFFFFFF09FFFE00;BE@;DABBFFFFFBBFB00;F:9;FFBFFF9BFFFFFFFFFFFFF90/::BFFFBF0\n+@M01368:8:000000000-A3GHV:1:1101:18422:19051/2\n+CTACAAGCATTAATTAATTAACACACTTTAGTAAGTATGTTCGCCTGTAATATTGAACGTAGGTGCGATAAATAATAGGATGAGGCAGGAATCAAAGACAGATACTGCGACATAGGGTGCTCCGGCTCCAGCGTCTCGCAATGCTATCGCGTGCACACCCCCCAGACGAAAATACCAAATGCATGGAGAGCTCCCGTGAGTGGTTAATAGGGGGATAGACCTGTGATCCATCGTGATGTCTTATTTAAGGG\n++\n+BCCCCFFCFFFFGGGGGGGGGGHHHGHHHHHHHHHHHHHHHHGGGGHHHHHHHHEHHHHHHHGGHHGGHGGHHHHHHHGHGGHHHGGGGGHGHHHHGGGHFHFHHHHHGGGGGHBFFCGDHHHGGGGGGHGGGGGGHHGCGGGFGHHBGGGGGFFFHEGGGGGCDCCE@EFGHHHHFHEGHGFFHHGB;ECBFGGGEFEFFGF0AFGFGFFG.;;DFFFFFFFFFF090BFFFE?FEFBBFBFFFB990BF\n+@M01368:8:000000000-A3GHV:1:1101:25545:21098/2\n+GCCCATAACACTTGGGGGTAGCTAAAGTGAACTGTATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTTCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCAACCTATGTCGCAGTATCTGTCTTTGAT\n++\n+3AAA34BFFFFFGCGGCGGGCGHHHGHHHHCGHFGHFHGFGGCDGGBD5FGFFFHFHGGGHAEBFGHHHEHGGFHBGGFHHHGGGEEFHHHFHFFHDHGHGHHHHFFHHGGFAEGFFFFFBGHHHFFHHHHHHHFHFHHFHGFHGHGA/<@EFHEHH1GGHFFFHG1>=FGHFHGF.GE//-<BB?9.9BFFFB==AA/FFFFFA-@BA.;;D?F9FAB;---./99BFFFBBAE-.9B/BFB9F/9BFBB\n+@M01368:8:000000000-A3GHV:1:1101:5446:12248/2\n+ATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCTATTATTTATCGCACCTACGTT\n++\n+CCDDDCCCCFFFGGGGGGGGGGHHHHHGHHHHHHHHGHHHHHHGHHHHGGGGHHHHHHHHHHHHHHHHHHHHHHGGGGGHHHHHHHGFGHHHHHHHHHHHHHHHHHGHHHHHGGGGGHHHHHHHHHHHHHHHGHHHHHHGHGHGHHGGGGGGGGGGGFFFFFFFFFFFFFFAFFFFFFFFFFFFFBFC?.EAFFFFFFFFAFFFFFFFFFFFFFFFEFFFF0FFFFFFF0BFFFFFFFFF?.BDCFFEDA/\n+@M01368:8:000000000-A3GHV:1:1101:5861:6452/2\n+TTATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCTATTATTTATCGCACCTACG\n++\n+BCCCCFBBCCCFGGGGGGGGGGHHHHHHHHHHHHHHHHGHHHHHHHHHHHGGGGHHHHHHHHGHHHHHHHHHHHHHGGHGGHHHHHHHGHHHHHHHHHHHHHHHHHHHGGHHHHGGGGGGEHGHHHHHHHHHHHGHHHHHHGHGHGHGGGGCGGGGGFFFFFFFFFFFFFFFFFFFFFFFF?FFEAF:DAEF9DEFFFFFFF-A.BFFFFFEFFFFEFFFFFFFF0FFEB0FFFFBFFFFFFADEFCF0/.\n+@M01368:8:000000000-A3GHV:1:1102:10403:6021/2\n+GGATGAGGCAGGAATCAAAGACAGATACTGCGACATAGGGGGCTCCGGCTCCAGCGTCTCGCAATGCTATCGCGGGCACACCCCCCAGACGAAAATACCAAATGCATGGAGAGCTCCCGTGAGAGGTTAATAGGGGGATAGACCTGTGATCCATCGTGATGTCTTATTTAAGGGGAACGTGTGGGCTATTTAGGTTTTATGACGCTGAAGTAGGAACCAGATGTCGGATACAGTTGACTTTAGGTAACCCA\n++\n+AA?A34BF>AADEBGGGGFFCGGHHA5DFGAE?AAEDFHG2AEE1FBE00A1BGE0?E1FFGEEEGB4F43?EE/>///F??//@BCBAC<<BAGHB?11F/FGFGGBFBGBD.11><.<<<C0<A</<<=0GFHC@-.;FF09BFFGB00B0;FE.C/:CCF0;0C0;BBBFF@.9AF.:9EFF.;.BFF/FF/;/./BFB/..;;ABF/9FF//.;AAF9BBBA;-A.B/9///9/BF////;///:9.\n+@M01368:8:000000000-A3GHV:1:1102:10677:23253/2\n+AATTAACACACTTTAGTAAGTATGTTCGCCTGTAATATTG'..b'HGGGGGGGGGFEACFFFFFFFFFFFFFFFDFFFAFFFFBFAB@EFAAEFFFFF.ACF.BBFFEBFFFEB;FFFFFFFA/BFBFBBBFFFBFFFFFED.>DFFF.\n+@M01368:8:000000000-A3GHV:1:1113:5741:16959/2\n+GTATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCTATTATTTATCGCACCTACG\n++\n+BBBBCFCCCCCFGGGGGGGGGGHHHHHHHGHHHHHHHHGHHHHGHDGHHHGGGGHHHHHHHHHHHHHHHGHHHHHGGGEGGHHHHFHHHHHHHHHHHHHHHHHHHHHHGHHGHHGGGGGHHHHHHHHHHHGHHHFHHHHHHGHGHGHGGGGCGGFGGFFFFFFFFFFFFFBFFFFFFFFFFFFFFAFFFFEAEFFFFFFFFFFF9BFFBFFFFFFFFFFFFFFBFFFFFFFFFFFFFFFFFFADAB-/BF.\n+@M01368:8:000000000-A3GHV:1:1114:10130:11959/2\n+ATCAGAGACAGATACTGCGACATAGGGTGCTCCGGCTCCAGCGTCTCGCAATGCTATCGCGTGCACACCCCCCAGACGAAAATACCAAATGCATGGAGAGCTCCCGTGAGTGGTTAATAGGGTGATAGACCTGTGATCCATCGTGATGTCTTATTTAAGGGGAACGCTGTCTCTTATACACATCTGACGCTGCCGACGAGCGATCTAGTGGAGAGTTCTGTGGGGCCGGGGATCCTTATAAAAAAAATAGA\n++\n+BCCDDFFFFFFFGGGGGGGGGGGHHHHHGGHHHGGGCGGHGHGGGGHHGGGGHHHHHHGGGEGGHHHFGGGGG?E1FE?/EEHHHHHGHHGHHHHGHFHGHGHHGDGGFG2FF2?GHHHHHGCCCFHGHGHHHHGHHFEHHFGHHGHH<1=DGHHHGHHGHGAGAEEDG.CGCGHC0CGBFHGFBBF0ABDDEFF@?--:BB@.;:BF;0.0:0.0:000:BF.-.------.9/;000::0;0:--.000\n+@M01368:8:000000000-A3GHV:1:1114:14540:5315/2\n+CACTTGGGGGTAGCTAAAGTGAACTGTATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCC\n++\n+CCCCCFFCCDCCGGGGGGGGGGHHHHHHHHHGGGGGHHHHHHHHHHHHHHHHHHHGHHHHHHHGHHHHFHFHHHGGGGGHHHHHHHHHHHHHGHHHHHHHHGGGGGGHHHFHHHGHHHHHHHHHHHHHHHHHHGHGHHHGGGFGHHHHHHHHFHHHHHF?1FHHGHGHGHGHHGGFFFFDBFBE;BCC.:BFFFFFFFFFFFFFF;AFFFFF-=-.AEDEFFFFF..9A;9FFFF0FFFFE00FFF0:BA.\n+@M01368:8:000000000-A3GHV:1:1114:15066:16302/2\n+TTATTATTATGTCCTACAAGCATTAATTAATTAACACACTTTAGTAAGTATGTTCGCCTGTAATATTGAACGTAGGTGCGATAAATAATAGGATGAGGCAGGAATCAAAGACAGATACTGCGACATAGGGTGCTCCGGCTCCAGCGTCTCGCAATGCTATCGCGTGCACACCCCCCAGACGAAAATACCAAATGCATGGAGAGCTCCCGTGAGTGGTTAATAGGGGGATAGACCTGTGATCCATCGTGTGG\n++\n+CCCCCFFFFFFFGGGGGGGGGGHHHHHHHHHHHHHHHDGHHHHFHHGHHGHHHHHGGGGEHHHHFHHFF5FHHFEGHHHGDHGGHGHGFGGGEHFGHHGGGGGHGHHHGHHFHHB3FGHHFGGGG?GFFHCCEBGFFECCDFEGFCFGCHHGFDDHHHGHHCFGGGGGFBFDGFG?-:..AFG.-C0C009;:00;00:9/:CEFFF?AE::9;9?0:FEF0;0..--./00::B/090000000;A....\n+@M01368:8:000000000-A3GHV:1:1114:16639:15258/2\n+TTATTATTATGTCCTACAAGCATTAATTAATTAACACACTTTAGTAAGTATGTTCGCCTGTAATATTGAACGTAGGTGCGATAAATAATAGGATGAGGCAGGAATCAAAGACAGATACTGCGACATAGGGTGCTCCGGCTCCAGCGTCTCGCAATGCTATCGCGTGCGCACCCCCCAGACGAAAATACCAAATGCATGGAGAGCTCCCGTGAGTGGTTAATAGGGGGATAGACCTGTGATCCATCGTGATG\n++\n+CDCCCFFFFFFFGGGGGGGGGGHHHHHHHHHHHHHHHHGHHHHHHHHGHHHHHHGGGGGHGIHHHIH5DEGHHHF?FGHGGHGGHEGGHFHHGHGEHHGGGGGFFFGHFBG2GHEBGHHGHGGEG/GFGABEDFGHEED?GGHHFFGGGCFEGD/GFHFFGEFGCGG?CC??D-EF@EEEFGCDDBBFGGGEBBFFF09090A.BFGA.9CCA0;EBAB00BBFF.@-./;BB;BFFF0:00099AAFFF0\n+@M01368:8:000000000-A3GHV:1:1114:2404:13066/2\n+ATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCTATTATTTATCGCACCTACGTTCAATATTACAGGCGACCATACTTACTAAAGTGTGTTAATTAATTAATGCTTGTAGGACTGTCTCTTATACACATT\n++\n+CCCCCFFFFFFCGGGGGGGGGGHHHHHHHHHFFHHHHGGGGGHFFFHHFHHHHHHHHHHHHHHHGFEGGGHGEDFCDFHGHFG@@DGGHHHHHHGGGGCGGGGGEHGGCGBB?CF99EGFGGFGG?D9CFFFF/BBFFFFFEF9BFFAFFFFEFFFFFFFFFFFFFFFFFFFFF.FFBBFFFFFFFFFFFF-9;;;BFFFFFB9BFBFBFABFFEFFFFFFFFFF::BFFBFFFF.9//;FFFFF/BFFB/\n+@M01368:8:000000000-A3GHV:1:1114:9184:6959/2\n+AAAGTGAACTGTATCCGACATCTGGTTCCTACTTCAGGGTCATAAAACCTAAATAGCCCACACGTTCCCCTTAAATAAGACATCACGATGGATCACAGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGGGTGTGCACGCGATAGCATTGCGAGACGCTGGAGCCGGAGCACCCTATGTCGCAGTATCTGTCTTTGATTCCTGCCTCATCCCTGTCTCTTA\n++\n+CCCCBFFFFFFFGGGGGGGGGGHHHHHHHHHHHHHHHHHHHHHHHHHHGHHHHHHEHIHHGGGGHHHHHHHHHHHHHGHHHHHHHHGGGGGHHFHHHHHBGHHHHHHHHHHHHHHHHHGHHHHHGGGGGHHHHGHHHHHHHHHHHHHHHHHGHGHGHHGGGGCFFFFFFFFFFFFFFFFFFFFFFFFFF.CFFFFAF=D=EAEFFF0B:0AF-DAFBFFFFFFFFFBFFFFFFFFFFBFFFEFF9B900B0\n'