Repository 'nanofilt'
hg clone https://toolshed.g2.bx.psu.edu/repos/leomrtns/nanofilt

Changeset 0:ff1cc05e789a (2019-05-14)
Next changeset 1:a2c318b04f49 (2019-05-16)
Commit message:
planemo upload
added:
nanofilt.xml
test-data/input.fastq
test-data/output.fastq
b
diff -r 000000000000 -r ff1cc05e789a nanofilt.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/nanofilt.xml Tue May 14 05:57:15 2019 -0400
[
@@ -0,0 +1,107 @@
+<tool id="nanofilt" name="NanoFilt"  version="0.1.0">
+  <description>Filtering and trimming of long read sequencing data</description>
+  <requirements>
+    <requirement type="package" version="2.3.0">nanofilt</requirement>
+  </requirements>
+  <command detect_errors="exit_code"><![CDATA[
+    NanoFilt 
+    #if str($filter.quality) != ""
+    -q ${filter.quality} 
+    #end if
+    #if str($filter.minlength) != ""
+    --length ${filter.minlength} 
+    #end if
+    #if str($filter.maxlength) != ""
+    --maxlength ${filter.maxlength}
+    #end if
+    #if str($crop.head) != ""
+    --headcrop ${crop.head}
+    #end if
+    #if str($crop.tail) != ""
+    --tailcrop ${crop.tail}
+    #end if
+    #if $summary.selector
+    --readtype "${summary.readtype}"
+    --summary "${summary.summary_file}"
+    #end if
+    --logfile $logfile  
+    "$input1" > "$output1"
+    ]]></command>
+  <inputs>
+    <param type="data" name="input1" format="fastqsanger" label="uncompressed fastq file" />
+    <section name="filter" title="Options for filtering reads based on length or quality">  
+      <param name="quality"  type="integer" label="Filter on a minimum average read quality score" value="30" min="0" max="255" help="(-q)" optional="true" />
+      <param name="maxlength" type="integer" label="Filter on a maximum read length" help="(--maxlength)" optional="true" />
+      <param name="minlength" type="integer" label="Filter on a minimum read length" help="(--length)" optional="true" />
+    </section>
+    <section name="crop" title="Options for trimming reads">  
+      <param name="head" type="integer" label="Trim n nucleotides from start of read" help="(--headcrop)" optional="true" />
+      <param name="tail" type="integer" label="Trim n nucleotides from end of read" help="(--tailcrop)" optional="true" />
+    </section>
+    <conditional name="summary">
+      <param name="selector" type="boolean" checked="false" label="Use albacore or guppy summary file for quality scores"/>
+      <when value="true">
+        <param type="data" name="summary_file" format="tabular" label="Summary file generated by albacore or guppy."/>
+        <param name="readtype" argument="--readtype" type="select" label="Which read type to extract information about from summary">
+          <option value="1D" selected="true">1D</option>
+          <option value="2D">2D</option>
+          <option value="1D2">1D2</option>
+        </param>
+      </when>
+      <when value="false">
+        <param name="minGC" type="float" label="Sequences must have GC content higher than this" min="0.0" max="1.0" help="(--minGC)" optional="true" />
+        <param name="maxGC" type="float" label="Sequences must have GC content lower than this" min="0.0" max="1.0" help="(--maxGC)" optional="true" />
+      </when>
+    </conditional>
+  </inputs>
+  <outputs>
+    <data name="output1" format="fastq" />
+    <data name="logfile" format="tabular" label="Log file" from_work_dir="logfile.txt" />
+  </outputs>
+  <tests>
+    <test>
+      <param name="input1" value="input.fastq"/>
+      <output name="output1" file="output.fastq"/>
+    </test>
+  </tests>
+  <help><![CDATA[
+    usage: NanoFilt [-h] [-v] [--logfile LOGFILE] [-l LENGTH] [--maxlength MAXLENGTH] 
+    [-q QUALITY] [--minGC MINGC] [--maxGC MAXGC] [--headcrop HEADCROP] [--tailcrop TAILCROP] 
+    [-s SUMMARY] [--readtype {1D,2D,1D2}]
+
+    Perform quality and/or length and/or GC filtering of (long read) fastq data.           Reads on stdin.
+
+    EXAMPLES:
+
+    gunzip -c reads.fastq.gz | NanoFilt -q 10 -l 500 --headcrop 50 | minimap2 genome.fa - | samtools sort -O BAM -@24 -o alignment.bam -
+
+    gunzip -c reads.fastq.gz | NanoFilt -q 12 --headcrop 75 | gzip > trimmed-reads.fastq.gz
+
+
+    ]]></help>
+  <citations>
+    <citation type="bibtex">
+      @misc{githubnanofilt,
+      title = {nanofilt},
+      publisher = {GitHub},
+      journal = {GitHub repository},
+      url = {https://github.com/wdecoster/nanofilt}
+      }
+      @article{10.1093/bioinformatics/bty149,
+      author = {De Coster, Wouter and D’Hert, Svenn and Schultz, Darrin T and Cruts, Marc and Van Broeckhoven, Christine},
+      title = "{NanoPack: visualizing and processing long-read sequencing data}",
+      journal = {Bioinformatics},
+      volume = {34},
+      number = {15},
+      pages = {2666-2669},
+      year = {2018},
+      month = {03},
+      abstract = "{Here we describe NanoPack, a set of tools developed for visualization and processing of long-read sequencing data from Oxford Nanopore Technologies and Pacific Biosciences.The NanoPack tools are written in Python3 and released under the GNU GPL3.0 License. The source code can be found at https://github.com/wdecoster/nanopack, together with links to separate scripts and their documentation. The scripts are compatible with Linux, Mac OS and the MS Windows 10 subsystem for Linux and are available as a graphical user interface, a web service at http://nanoplot.bioinf.be and command line tools.Supplementary data are available at Bioinformatics online.}",
+      issn = {1367-4803},
+      doi = {10.1093/bioinformatics/bty149},
+      url = {https://doi.org/10.1093/bioinformatics/bty149},
+      eprint = {http://oup.prod.sis.lan/bioinformatics/article-pdf/34/15/2666/25230836/bty149.pdf},
+      }
+    </citation>
+  </citations>
+</tool>
b
diff -r 000000000000 -r ff1cc05e789a test-data/input.fastq
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/input.fastq Tue May 14 05:57:15 2019 -0400
b
b'@@ -0,0 +1,73000 @@\n+@NB501061:3:HF573AFXX:1:11101:3103:1113 1:N:0:CGTACTAG+NTAGAGAG\n+TTCCTATAACTGGTTGACTTGAACCAACAGTTGCTATTAAATTTGCTCGACCATCTGTAACTGTATCAATA\n++\n+AAAAAEEEEEAEEEEE6EEEEAAEEEEEEAEAEEEEEEEEAEEEEEEEEEEEEEEEAEEEEEEEEEEEAEE\n+@NB501061:3:HF573AFXX:1:11101:2518:1121 1:N:0:CGTACTAG+NTAGAGAG\n+GATTTACTTATGGACGATATGGATGGTGTAGAAGCAACTACTGAAATAAAAAAAGATTTACCTCAAATTAAAGTAGTCATGTTAACAAGCTTTATAGAGGATAAAGAAGTTTATCGTGCACTTGATTCTGGAGTAGATAGTTATATTTTAA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAA6EEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEAEEEEEAEEEEEEE/EEAEEEEEEEA<EEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAE\n+@NB501061:3:HF573AFXX:1:11101:22477:1123 1:N:0:CGTACTAG+NTAGAGAG\n+GGTCTATTACGTAAAGGTGAAGGTGACATGGTCATGGAGCAATTCGGTGAAGGTTTTGATATGAATATTATTCGTGTTAATGCTCAAGAACGCTTCATGAGTAAACTAAAAGGTGTGTCTGATCCAGAACGTAAACGTAAAATCATTGGTA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEAEEEEEEAEEEEEEE/EE/EEEE/EEEEEEE/AEEEEEEEEEEEEEEAEEEEEEEEEEEEEEAEEEEE<AEEEAEEEEEEEEEEAE\n+@NB501061:3:HF573AFXX:1:11101:24857:1142 1:N:0:CGTACTAG+NTAGAGAG\n+GTGCTGTATATCTTGTAAATCTATTCTCTAATATTAAAACTCCAGAGAATCTTAAACATATCAAAGAGCCATACGATGAGCACACAGATATCCACTTGATGAAAGCGATAAGCGAAAGTGACACAGTCATTCTAGCTTATGGCGCCTATG\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEAEEEEEAEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEE/E<AAEEEEEEA/<<AAAAEEAEEEEEAEAAEEAE/EEEAA/EAEEEE/EEAEA/E<A<A<EA<<EA<E<EE/EEE<EEAE//E\n+@NB501061:3:HF573AFXX:1:11101:12352:1144 1:N:0:CGTACTAG+ATAGAGAG\n+TAGCTGGCTTAAACCAAAACGATCAGCCTCTTCTATGATTAAAGTATTAGCATTTGGTACATCTACACCTGTTTCAATAATTGTAGTCGTTACTAAAATATCGTACTCGTGATTAATAAAGCTTAACATTGTTTCCTCTAAATCACGTTCA\n++\n+AAAAAEEEEEEA6EEEE6EEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEEE6EEEEEAEEEEEEEEEEEEEAEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEAEEEA\n+@NB501061:3:HF573AFXX:1:11101:15400:1147 1:N:0:CGTACTAG+ATAGAGAG\n+TATGTACAGTATGACAGTATTTAGTTTAATGAACTTTTGTCTATTGTCATTTCCATTAGATTTGATTGAAGAGCGAAATCAAGGTTGGTATAAGCGGTTAATGGTGACCCCATTATCATCATTTCAATATTATTTAGTTAAAATTAGTAAG\n++\n+AAAAAEEEEEEEEEEEEEEEEEE6EEEE<EAEAEEEEEEEEEEEEEEEEEEEE/EEEAEEEEEEEAEEEEAEE66EEEEEEEEEEEE/EEEEEAA<EEEEEEEEEEA/<AEEE/EE<E<<E/EEEEE/EEEE<EAEEEEEEEEEEAEEAE/\n+@NB501061:3:HF573AFXX:1:11101:2858:1149 1:N:0:CGTACTAG+ATAGAGAG\n+ACTAAAGATGAGCGTTTAGATGCAGTGATGAGTGCAGGAGGACTGCAAGAATGTGGTAATTCTCAAAATAGTGTGAATGCATGTCCTAAAGGTATACCATTGACAACGTCAAATGCCGCAATGAATAGAGAAACTTCTTTCCATATGTTAA\n++\n+/AAA66E//EE6EE/AEEEEE<///EA//EEE/E/EEEAAE//EEAA/A6E/6A//E/EEE<EAEEEEE//EE/E/EEEEA</EEE/EAEE/AE/EEEEE/EAEE//<AEE///</<<E<E/EEE/AA/<AAAA/EA/E/A<EEAAEE///\n+@NB501061:3:HF573AFXX:1:11101:4491:1150 1:N:0:CGTACTAG+ATAGAGAG\n+ATGCATTGTAATCGAAACATACCCCTTGTATAGACTTCATCGGATCGATTAACACCTATCACTTGAGTCACAGCACTTGGAATAGCATTCATAGATAAAATAACCCCTAGTGCTACGACAGGATATATTTGTTGATAAGCATATAAACCGT\n++\n+<AAAAEEEEEE6AE6EE6EEE/EEEEEEEEEEE6EEEE/EEEEEEEEEE6EEE6EEEEAEE6EEE/EEE6E6E6EAEEEE6EA<EEAEAEE/AAEEAEAEEEEAEEEE/AEEAAEEEEEEEEEEEEEEEEAEE/A<EEEEEEE</E/AEE/\n+@NB501061:3:HF573AFXX:1:11101:12967:1151 1:N:0:CGTACTAG+ATAGAGAG\n+TAGTAAGACTATACGGTTTGTCTTTTTAGAAAATTCTGATAGAATAAAGTTATATCCGATAGAAAAACTAAGTATTGATAAATTGTAAAACTAAAAGGAGTGTTTTGGATGGCACAAAAACCTGTAGATTATGTTACACAAATTATTGGGA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEE<AEEAEEEEEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEE<EEEEEEEEEEEEEEEEEEAEEEAEEEEEE\n+@NB501061:3:HF573AFXX:1:11101:24302:1152 1:N:0:CGTACTAG+ATAGAGAG\n+ATTGTATAATTATAAAGAAATAAGAACGCCAATTTTTGAAAGTACTGAACTTTTTGCAAGAGGCGTGGGGGATTCTACTGAGGTTGATCAAAAGGAAATGTATACATTTAAAGATAAAGGGGATCGTAGTTTAACATTACGTCCTGAAGGA\n++\n+/AA/A///EA/AEEE6A//6EEE6E6A/////E/EEEEEEEE/EE/EE6/EA/A</E/E/E//E////AEE//<6A///////EEA///E//E//E/EA//A//AEE/E//E6EEAA///E/AAA</E///A/E/EE//</E666<//EE/\n+@NB501061:3:HF573AFXX:1:11101:6728:1157 1:N:0:CGTACTAG+ATAGAGAG\n+TACTAAAAGGTAGCGATATTTTAGGAACGAATAACGCGCTTAGTCACGCTTTTAGTGGTTTTTCAACCTCAGCCGTAGCACTTGTAGCTGCAGCATTATTTTTAGCAGTAGCTATGCAGGAAACCAATTTACATAAACGACTTGCATTA\n++\n+AAAAAEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEE/EEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEE<EEEEEEEEEEE/A6AEEEEEE\n+@NB501061:3:HF573AFXX:1:'..b'NB501061:3:HF573AFXX:1:11106:25404:2014 1:N:0:CGTACTAG+ATAGAGAG\n+TCCGACGATATCTTCTTCTGATAAACCTTCCAAATACTTAATTAATGCGCGAATGGAATTTCCGTGGGCCGAAACAAGTACAGTTTTATCATCTAATAGGTGTTGTGAAATATGATCAGTCCAGAACGGAACAACTCGTTCTAATGTAT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEAEEEAEEEE/A/EEEEEEEEEEEEEAAE\n+@NB501061:3:HF573AFXX:1:11106:20809:2016 1:N:0:CGTACTAG+ATAGAGAG\n+CTTTATTAATATGTAAAAATAAGTGTTTTATTTGCACGCTAACTTGGCTTTTCAACATAAATACTTTTTGAGTAAAAAGGAGGAACTTTTATGGATTATGCACATTTAAATCTCGAGCACTTTTTTGCAAGAAATGATGATTTGGATATTA\n++\n+AAAAAEEEEEEEEEEEEEEEE6EEEAAEEEEEEAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEAEEE6EEEEEEAEEEEEEEE<EEEE/\n+@NB501061:3:HF573AFXX:1:11106:2456:2024 1:N:0:CCTACTAG+ATAGAGAG\n+GAATATCGAGAGAGTGAATATTGCTTCGCACTTATCAAAGTAGTAACAGTCATCATATTAATAGCAAATGGAATCTTAACTATTTTACGCATTATGGGTGGAGAATTTGAAGCAAATGATACCTTTAGAAAAGGAGAAGGCCCAATAGAAG\n++\n+A6/AAA/EEEE6EE6E6E/AEAE/EE/EEE///EEEA//E/6EE/EE66E/A/EAE6/E6EEEA/EE///E/E/EAEE///EEE//E/E/E//A///EAE<E/EE//EE//A/E///A///</E//A//6<EAEEE///E/6///////<A\n+@NB501061:3:HF573AFXX:1:11106:15329:2025 1:N:0:CGTACTAG+ATAGAGAG\n+CCTTATACAACTTCACAAACTCTTCTTCATTCATGGCAAGCTTTGGAAATAACATATTATATTTTTTTTCGATATCTTTCTGTAACTTCACTACTTTAGGCAACTCTATATGCCACTTTAATAATTGCGTTTTTTGATTAAATAAATTT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6E/EEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEE<EEEEEEEEEEEAAEE/EEEEA\n+@NB501061:3:HF573AFXX:1:11106:20334:2026 1:N:0:CGTACTAG+ATAGAGAG\n+CTATAATGCCAAAGATCCACAAAAACATGCACTGTTTCCTAAATATGATTTCTTTAGAGCAGATACTGACTATGCTGACATTGCAAAATTTTTAGGACTCAAAGGTAATACAACTGAAGAATTAGTGGATGCTCTAGCTAATGCGGTGTAT\n++\n+AAAAAEEEEAEEEEEEEEEEEEEE6EEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/AEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEAEEEEEEAEEEEEEEAEEEEE/EEEEEE<EEEAE/EEEE<AA<EEEE\n+@NB501061:3:HF573AFXX:1:11106:2427:2031 1:N:0:CGTACTAG+ATAGAGAG\n+ACATCCGATAGTGCATCTACATCAAAGAGTGTAAGTGAGTCAAACAGCGCCAGTACATCGGTAAGTGGCTCACCAAGTACAAGCATTTCAGATTCATCAAGCACGTCCACATCAATGAGTACATCTGAAACTTTCACTTCTCAATCTCCTA\n++\n+6AAAAEEAEEEAEEEEEEEE6EEEEEEEEEEEEAEEEEAAEEEEA/EEAEEEEE/EEEEEEEEEEEEEEEEE/<E/AEEAAEEAEEAEEEE/EAEEEA/<<EEE<EEEEAEAAA/E<EA</AEEEEA/EE///AE//E<EEE///<E6AA/\n+@NB501061:3:HF573AFXX:1:11106:10496:2032 1:N:0:CGTACTAG+ATAGAGAG\n+CCATTAGGTCGACTAGGAACACCGGATGAAGTCGGGAAACTTGTAGCATTTTTAGCTTCCGATGATAGTTCATTTATAACTGGTGAAACTATTCGTATAGATGGTGGCGTTATGGCTTATACATGGCCTGGAGAAATGTTAAGTGATGAAA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEE6EEEEE/EEEEEEEE/\n+@NB501061:3:HF573AFXX:1:11106:15681:2048 1:N:0:CGTACTAG+ATAGAGAG\n+ATATATGACTGCTGAAGAAACAAAACAATTTAATGAAGGTGTATCAGGTGATTTCGTTGGCATAGGTGCTGAAATGCAAAAGAAAAATGAACAGATAAGTGTTACTAGTCCAATGAAGGATTCACCAGCAGAAAAAGCTGGTATTCAAC\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEEE<EEEEEEEEEEEEAEEEEEAEEEE/<E\n+@NB501061:3:HF573AFXX:1:11106:11089:2054 1:N:0:CGTACTAG+ATAGAGAG\n+CTAAAGGACTAAATGCAGCTCCAAAACCAGCCCAAGCATTACCAACTAAATTTAGTATTGTATCATTTGGATGCCAAGCAATCGTAATAGCAACTATCGCAACAAGTAAAACTGATAAGCGTCCAATCAAAACAAACTCTTTTTGGTGTGA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEA6EAEEEE/EAEEEEEEEEEEEEEE/AEEE/EEEEEE/EE//EEAE//EEE6EEE<AEAEEE<EEE/\n+@NB501061:3:HF573AFXX:1:11106:18404:2057 1:N:0:CGTACTAG+ATAGAGAG\n+TTATAAGACACATACAATTTTTAATATATCGATTAACTAAAGGTGAGTATTTAGAAGCACAGGAGAATTTTATATCGATGATAAAAACTATGTATC\n++\n+AAAAAEEEEEEEEAEEEEEEAEEEEEEEAA<EEEEEEE/EEEAEEEEEEEEEEEEEA<E<EEAEEEEEEEEEEEE//EEAEEEEEEEEAEEEEEE<\n+@NB501061:3:HF573AFXX:1:11106:22800:2060 1:N:0:CGTACTAG+ATAGAGAG\n+CAAGTATATCATGCTAAAGTAAAACGTATTGAAAAGTATGGTGCTTTCGTTGAATTGTTCCCTGGTAAAGACGCGTTATTACACATTTCTCAAATTTCACAAGAAAGAATTAATAAAGTAGAAGATGTTCTTAAAATTTGAGATACAATTG\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEE/E/<EAEEEEEAEEEEEAEEEEEEEEEAEEEEEEAEE/EEEAEEEEEEEE/EE/EE/EEEAEEEEEEEE\n'
b
diff -r 000000000000 -r ff1cc05e789a test-data/output.fastq
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output.fastq Tue May 14 05:57:15 2019 -0400
b
b'@@ -0,0 +1,27872 @@\n+@NB501061:3:HF573AFXX:1:11101:3103:1113 1:N:0:CGTACTAG+NTAGAGAG\n+TTCCTATAACTGGTTGACTTGAACCAACAGTTGCTATTAAATTTGCTCGACCATCTGTAACTGTATCAATA\n++\n+AAAAAEEEEEAEEEEE6EEEEAAEEEEEEAEAEEEEEEEEAEEEEEEEEEEEEEEEAEEEEEEEEEEEAEE\n+@NB501061:3:HF573AFXX:1:11101:2518:1121 1:N:0:CGTACTAG+NTAGAGAG\n+GATTTACTTATGGACGATATGGATGGTGTAGAAGCAACTACTGAAATAAAAAAAGATTTACCTCAAATTAAAGTAGTCATGTTAACAAGCTTTATAGAGGATAAAGAAGTTTATCGTGCACTTGATTCTGGAGTAGATAGTTATATTTTAA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAA6EEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEAEEEEEAEEEEEEE/EEAEEEEEEEA<EEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEAAAE\n+@NB501061:3:HF573AFXX:1:11101:12352:1144 1:N:0:CGTACTAG+ATAGAGAG\n+TAGCTGGCTTAAACCAAAACGATCAGCCTCTTCTATGATTAAAGTATTAGCATTTGGTACATCTACACCTGTTTCAATAATTGTAGTCGTTACTAAAATATCGTACTCGTGATTAATAAAGCTTAACATTGTTTCCTCTAAATCACGTTCA\n++\n+AAAAAEEEEEEA6EEEE6EEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEEE6EEEEEAEEEEEEEEEEEEEAEEEE/EEEEEEEEEEEEEEEEEEEEEEEEEEAEEEA\n+@NB501061:3:HF573AFXX:1:11101:12967:1151 1:N:0:CGTACTAG+ATAGAGAG\n+TAGTAAGACTATACGGTTTGTCTTTTTAGAAAATTCTGATAGAATAAAGTTATATCCGATAGAAAAACTAAGTATTGATAAATTGTAAAACTAAAAGGAGTGTTTTGGATGGCACAAAAACCTGTAGATTATGTTACACAAATTATTGGGA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEE<AEEAEEEEEEEEAEEEEEEEEEEEEEEEEEEEEAEEEEEEE<EEEEEEEEEEEEEEEEEEAEEEAEEEEEE\n+@NB501061:3:HF573AFXX:1:11101:9265:1161 1:N:0:CGTACTAG+ATAGAGAG\n+CAATTAAATGACCGCCGTCATCATCTGGTAATCGGTCTTCTCTTCCTACTGTTCGTTGCGCATATGAATTACGTCCGCCATCTAGAAGTGATAGGTCATCCGCATATACTTCTTTAATACGGCCTTTGTGATC\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEE</EEEEEEEEEEEEEEEE<EEEAEEEEAEEEEEEEAAAEEE<EEEEEEEEE\n+@NB501061:3:HF573AFXX:1:11101:15298:1176 1:N:0:CGTACTAG+ATAGAGAG\n+GCTCATCATTGTACTGGATTTTTTAAAACTACATTGCCTGGCTTTTATAGTGTAGTTAAAGATATTTTTACACCGAACTATGAAAGGAAAGAGAGTTAGAATATGTTGCAACCTTTGATACATGGACTATTATTAGCTTTAGGTTTAATTT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEE<EEEEEEE/EEEEEEEEEEAEEEEEEEEEEEEEE<AEEEEEEEEAEEEE\n+@NB501061:3:HF573AFXX:1:11101:21083:1190 1:N:0:CGTACTAG+ATAGAGAG\n+CTTTTACTTATTCTAATTTACATTGATTAAAACACCGATTTTATTCTATCTTTAATTTAGCCATGATGTATTACCTCTTAGCTATAGCTTTAATTAGATGGCATTTCATCCTCAAACTCATCTCAATTGATTATGTTTGTTGAATTAAAAA\n++\n+AAAAAEEEEEEEEEEEEEEEEAEEEEEEEEEAEEEAEAEEEEAEEEEEEEEEEEEEEEEEEEAEEEEEEAEEEEAEEEAEEAA<EEEEEEAEEEEAE/EEEEEEEEEEEE6AEEEEEEAAEEE<EEEE<<EEAEAAEEEEEEAAEE/AAE6\n+@NB501061:3:HF573AFXX:1:11101:21687:1270 1:N:0:CGTACTAG+ATAGAGAG\n+CAGTAAGCATGATGATAGGTGTATCTTTATGTTCTCTAAGTTTAGAGGCAACTTCAATTCCATCCATTTCAGGAAGCATTAAATCTAAGAGGATGCACGCATAGTTGTTTTCAAGAGCAAGCTCGTAGGCTTCTTTACCATCTCTTGCTT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEEEEEEEEEEEE/EAEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEE\n+@NB501061:3:HF573AFXX:1:11101:6075:1283 1:N:0:CGTACTAG+ATAGAGAG\n+ACTTCCCACTTGCTTCAGGTATACTAGCTGGCAAGTACGATGAAAATCAAACATTTGAGGATCATCGTGCAGAACGTCCTGATTTTCAACCAGACGTATTCAAAGATAATGTGCGTCGTGTAAAAGCATTAACTGATATTGCAGAGGCACA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAAEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEE/<EEEAE<EEE<EEEEEEEEEEEEEEAEEE\n+@NB501061:3:HF573AFXX:1:11101:3339:1291 1:N:0:CGTACTAG+ATAGAGAG\n+ATTTCTGAATGATGATTTTATATTACCATCGAAAATAGATATTTGTGTTGAGCGCAAAAAGTTATTGGAAATTATCAGGGTGATTCCTGAAGAATTTACGATTCACTATTATGATAATCAGTTATATGAAAGATTAAGAGAATCTTTAT\n++\n+AAAAAEEAEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAAEE<EEEEEEEEEEEEAEEEEAE<EEEEEEEAEEEEEEEEEEAEEEAEEEEE<EEEE/EA6EEEEE\n+@NB501061:3:HF573AFXX:1:11101:22451:1312 1:N:0:CGTACTAG+ATAGAGAG\n+AATGTGATGTTAGTAGATTTAAATGAAAATAATCAAATTCAAGCTGAAAATGTGATTTTTATAGGCTTAATATGCAAAAATGCTGTTTGGCATGCTGTTGTATTTGAAATCGGTAGAGGTTACACTCGTGAATTATCTATCCTGGATATTA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEE<EEEEEEEEEEEEEEAAEEEEEEEEAEEA66AE/EEEEE/\n+@NB501061:3:HF573AFXX:1:11101:16981:1322 1:N:0:CGTACTAG+ATAG'..b'EEEEEEEEEEEEEEEEEEEEEEAEEAEEEEEEEEAEEEEEEEEEEEEEEEE\n+@NB501061:3:HF573AFXX:1:11106:20291:1959 1:N:0:CGTACTAG+ATAGAGAG\n+ATTTAAAGAATTATGGTTACCAATCATAACAATTTGTTTCATCTTAGCAATTTCAAAAATCACAACATATGGTGGTTTAAGTAATGCTATGGGACAAGGCATCTCAAAAGCAGGAAGCGTATTCCCAATATTATCACCAATCCTTGGTTGG\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEAEEEEEEAEEEEEEEEEEEEEEEEAEAAEEEEEEEEEEEEE\n+@NB501061:3:HF573AFXX:1:11106:20230:1964 1:N:0:CGTACTAG+ATAGAGAG\n+CAGCAGTAGGTGGTTGGAAAATCATCAAAACTGTTGGCGGAAATATTATGAAAATTAGACCTGCTAACGGTGCAGCCGCTGATATTTCTTCAGCATTAACGATTTTCGTAGCGTCATCATTACATTTCCCATTATCTACTACACACGTTGT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEE<EEEEEEEEEAAEEEAEEEEEEE\n+@NB501061:3:HF573AFXX:1:11106:18598:1965 1:N:0:CGTACTAG+ATAGAGAG\n+GCTTAAGTACGGTTGGTACAATTAAAGATGCTGACAACTGGACTTTAGATAGTGGATTCTACAAAACGCCTTCTAAATACAGCTTAGGTGATTATGTATGGTACGATTCTAATAAAGATGGTTACCAAGACAATGATGAAAAAGGTATTAA\n++\n+AAAAAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAEEEEEEEEEE6AAEEEAEEAEEEEEAEEEEEEEEEEE\n+@NB501061:3:HF573AFXX:1:11106:17916:1985 1:N:0:CGTACTAG+ATAGAGAG\n+TCCCAAGCAGGCGTTAAAATACAACTCATCATTCGTGGTATATGTTGTCTTAAACCAGGCATTCCAGGTATTAGCGAAAATATAGAGGTTGTTAGTATAGTGGGTCGTTTGCTTGAACATTCACGTATTTACTACTTCCATAATAATAGTG\n++\n+AAAAAAEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEAEEAEEE<EEEEEAA<EE/EEEEEEE<EEEEEEEEEEEEEAEEEEEE\n+@NB501061:3:HF573AFXX:1:11106:25172:1992 1:N:0:CGTACTAG+ATAGAGAG\n+TTTTATTAATGAATTCATTGAGTGTTTGACTTTTATTAAAATAAGTTTCTAAATTTCTTATATAAAAC\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE\n+@NB501061:3:HF573AFXX:1:11106:11380:2005 1:N:0:CGTACTAG+ATAGAGAG\n+GACCAGTTAAATATGCAACTAATGTTGCAGTAGAAATAACCGCAGCTAGTGTACCCCAGTCACGGAACAATGTCACCATAACCATGCTGATAATTGCATTAAATGCAATGGCAACACGTGGGATATTATACTTTTTATTAATTTTACCTAA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEA6EEE\n+@NB501061:3:HF573AFXX:1:11106:25404:2014 1:N:0:CGTACTAG+ATAGAGAG\n+TCCGACGATATCTTCTTCTGATAAACCTTCCAAATACTTAATTAATGCGCGAATGGAATTTCCGTGGGCCGAAACAAGTACAGTTTTATCATCTAATAGGTGTTGTGAAATATGATCAGTCCAGAACGGAACAACTCGTTCTAATGTAT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEAEEEAEEEE/A/EEEEEEEEEEEEEAAE\n+@NB501061:3:HF573AFXX:1:11106:20809:2016 1:N:0:CGTACTAG+ATAGAGAG\n+CTTTATTAATATGTAAAAATAAGTGTTTTATTTGCACGCTAACTTGGCTTTTCAACATAAATACTTTTTGAGTAAAAAGGAGGAACTTTTATGGATTATGCACATTTAAATCTCGAGCACTTTTTTGCAAGAAATGATGATTTGGATATTA\n++\n+AAAAAEEEEEEEEEEEEEEEE6EEEAAEEEEEEAEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEE/EEEEEEEEEEEEEEEAEEE6EEEEEEAEEEEEEEE<EEEE/\n+@NB501061:3:HF573AFXX:1:11106:15329:2025 1:N:0:CGTACTAG+ATAGAGAG\n+CCTTATACAACTTCACAAACTCTTCTTCATTCATGGCAAGCTTTGGAAATAACATATTATATTTTTTTTCGATATCTTTCTGTAACTTCACTACTTTAGGCAACTCTATATGCCACTTTAATAATTGCGTTTTTTGATTAAATAAATTT\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE6E/EEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEE<EEEEEEEEEEEAAEE/EEEEA\n+@NB501061:3:HF573AFXX:1:11106:10496:2032 1:N:0:CGTACTAG+ATAGAGAG\n+CCATTAGGTCGACTAGGAACACCGGATGAAGTCGGGAAACTTGTAGCATTTTTAGCTTCCGATGATAGTTCATTTATAACTGGTGAAACTATTCGTATAGATGGTGGCGTTATGGCTTATACATGGCCTGGAGAAATGTTAAGTGATGAAA\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEAEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEE6EEEEE/EEEEEEEE/\n+@NB501061:3:HF573AFXX:1:11106:15681:2048 1:N:0:CGTACTAG+ATAGAGAG\n+ATATATGACTGCTGAAGAAACAAAACAATTTAATGAAGGTGTATCAGGTGATTTCGTTGGCATAGGTGCTGAAATGCAAAAGAAAAATGAACAGATAAGTGTTACTAGTCCAATGAAGGATTCACCAGCAGAAAAAGCTGGTATTCAAC\n++\n+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEE6EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEEE<EEEEEEEEEEEEAEEEEEAEEEE/<E\n'