Mercurial > repos > brenninc > sync_paired_end_reads

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/syncpairs.py	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,111 @@
+"""
+Source: https://raw.githubusercontent.com/mmendez12/sync_paired_end_reads/master/sync_paired_end_reads/syncpairs.py
+"""
+__author__ = 'mickael'
+__author__ = 'mickael'
+
+from Bio import SeqIO
+from itertools import izip
+import argparse
+
+
+def adjust_name(reads1, reads2):
+    for r1, r2 in izip(reads1, reads2):
+        r2.name = r1.description
+        r2.description = r1.description
+        r2.id = r1.description
+        yield r2
+
+
+def remove_space_from_sequence_header(read):
+    """ replaces spaces in a read's name by three underscores.
+    Args:
+        read: A SeqRecord object (see Biopython)
+
+    >>> from Bio.Seq import Seq
+    >>> from Bio.SeqRecord import SeqRecord
+    >>> from Bio.Alphabet import SingleLetterAlphabet
+
+    >>> read = SeqRecord(Seq("AAAAA",SingleLetterAlphabet),\
+                       id="read A", name="read A", description="read A")
+
+    >>> print remove_space_from_sequence_header(read).name
+    read___A
+    """
+
+    read.description = read.description.replace(' ', '___')
+    read.name = read.description
+    read.id = read.description
+    return read
+
+
+def next_matching_read(reads1, reads2):
+    """ return next read2 that matches read2
+    Args:
+        reads1: A generator that contains a SeqRecord (see Biopython)
+
+        reads2: A generator that contains a SeqRecord (see Biopython)
+
+
+    >>> from Bio.Seq import Seq
+    >>> from Bio.SeqRecord import SeqRecord
+    >>> from Bio.Alphabet import SingleLetterAlphabet
+
+    >>> reads1 = []
+    >>> reads2 = []
+
+    >>> reads1.append(SeqRecord(Seq("AAAAA",SingleLetterAlphabet),\
+                       id="read A", name="read A", description="read A"))
+    >>> reads2.append(SeqRecord(Seq("TTTTT",SingleLetterAlphabet),\
+                       id="read A", name="read A", description="read A"))
+
+    >>> reads1.append(SeqRecord(Seq("AAAAA",SingleLetterAlphabet),\
+                       id="read B", name="read B", description="read B"))
+
+    >>> reads1.append(SeqRecord(Seq("AAAAA",SingleLetterAlphabet),\
+                       id="read C", name="read C", description="read C"))
+    >>> reads2.append(SeqRecord(Seq("TTTTT",SingleLetterAlphabet),\
+                       id="read C", name="read C", description="read C"))
+
+    >>> match = [read2 for read2 in next_matching_read(reads1, reads2)]
+    >>> print match[0].name
+    read A
+    >>> print match[1].name
+    read C
+    """
+
+    for read1 in reads1:
+        for read2 in reads2:
+            if read1.name == read2.name:
+                yield read2
+                break
+
+
+def main():
+
+    parser = argparse.ArgumentParser()
+
+    parser.add_argument("reads1",help='modified reads')
+    parser.add_argument("reads2", help='reads to adjust')
+
+    parser.add_argument('reads1_output',  help='output folder and filename. Note that the folder should already exist')
+    parser.add_argument('reads2_output',  help='output folder and filename. Note that the folder should already exist')
+
+    args = parser.parse_args()
+
+    #we'll need to go through the reads1 multiple time and it can be a large file
+    #so it's better to use inline func that return a generator
+    _reads1 = lambda: (rec for rec in SeqIO.parse(args.reads1, 'fastq'))
+    _reads2 = (rec for rec in SeqIO.parse(args.reads2, 'fastq'))
+
+    matching_reads2 = (read2 for read2 in next_matching_read(_reads1(), _reads2))
+    synced_reads2_names = (read2 for read2 in adjust_name(_reads1(), matching_reads2))
+
+    final_reads1 = (remove_space_from_sequence_header(r1) for r1 in _reads1())
+    final_reads2 = (remove_space_from_sequence_header(r2) for r2 in synced_reads2_names)
+
+    SeqIO.write(final_reads1, args.reads1_output, "fastq")
+    SeqIO.write(final_reads2, args.reads2_output, "fastq")
+
+if __name__ == '__main__':
+    main()
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/syncpairs.xml	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,65 @@
+<tool id="syncpairs" name="Sync Fastq Pairs" version="0.0.1">
+    <description>Runes sncy pairs on two fastq files.</description>
+    <requirements>
+        <requirement type="package" version="1.66">biopython</requirement>
+    </requirements>
+    <stdio>
+        <!-- Anything other than zero is an error -->
+        <exit_code range="1:" />
+        <exit_code range=":-1" />
+    </stdio>
+    <command interpreter="python">
+        syncpairs.py $reads1 $reads2 $reads1_output $reads2_output
+    </command>
+    <inputs>
+        <param name="reads1" type="data" format="fastq,fastqsanger" label="Read1 Fastq File" />
+        <param name="reads2" type="data" format="fastq,fastqsanger" label="Read2 Fastq File" />
+    </inputs>
+    <outputs>
+        <data name="reads1_output" format_source="reads1" label="Synced Read 1 from ${on_string}"/>
+        <data name="reads2_output" format_source="reads2" label="Synced Read 2 from ${on_string}"/>
+     </outputs>
+    <tests>
+        <!-- Large dataset in github but not toolshed -->
+        <!--test>
+            <param name="reads1" value="C75_S29_R1_001.fq" ftype="fastq" />
+            <param name="reads2" value="C75_S29_R2_001.fastq" ftype="fastq" />
+            <output name="reads1_output" file="C75_S29_R1_001_extracted.fastq" ftype="fastq" />
+            <output name="reads2_output" file="C75_S29_R2_001_extracted.fastq" ftype="fastq" />
+        </test-->
+        <test>
+            <param name="reads1" value="read1.fq" ftype="fastq" />
+            <param name="reads2" value="read2.fastq" ftype="fastq" />
+            <output name="reads1_output" file="read1_extracted.fastq" ftype="fastq" />
+            <output name="reads2_output" file="read2_extracted.fastq" ftype="fastq" />
+        </test>
+    </tests>
+    <help>
+<![CDATA[
+Description taken from source: https://github.com/mmendez12/sync_paired_end_reads
+
+sync_paired_end_reads
+
+sync_paired_end_reads is a python tool to synchronise paired-end reads when reads1 or reads2 were modified.
+
+When working with paired-end sequencing data, it is common to filter out reads that do not pass basic quality controls. This leads to pairs that are not synced anymore. This tool streams reads1 and search for the associated read2 in reads2. Additionally it synchronises the sequence identifiers of the reads so if a software modified the sequence identifier of the reads1 then same identifiers will be used for reads2. Finally it replaces all space characters by an arbitrary '___' pattern.
+
+This tool was mainly developed to process the output of tagdust2 when ran in single-end mode which appends the UMIs found in the raw sequences to the sequence identifier.
+Install
+
+]]>
+    </help>
+    <citations>
+        <citation type="bibtex">
+            @misc{
+                sync_paired_end_reads,
+                author = {Mickaël Mendez},
+                title = {sync_paired_end_reads on Github},
+                publisher = {GitHub},
+                journal = {GitHub repository},
+                url = {https://github.com/mmendez12/sync_paired_end_reads},
+                commit = {afd9bf092ed55a65f90f28f40d6230c5fe849223}
+            }
+        </citation>
+    </citations>
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/read1.fq	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,328 @@
+@NB500968:8:H5CFGAFXX:1:11101:15939:1071 1:N:0:TCCTGAGC+TACTCCTT;FP:23586;RQ:34.63
+GGCACGTCAGATGAGTATAAGAGACAGCAAATAGGTATAGGGGGCACGTCAGATGAGTAA
++
+EAEEEEEEA///A/EAEEE/EAEEEEEEEEEEAEEE//E/AA/E6AE/EE/<E////A//
+@NB500968:8:H5CFGAFXX:1:11101:1121:1139 1:N:0:TCGTGAGC+TACTCCTT;FP:27828;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGGTATTCCCTATAGTCTACATC
++
+EEEEEEEEEEEEEEEAAEEEEAEEEAEEE/EEE/EE<AAEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:7529:1175 1:N:0:TCCTGAGC+TACTCCTT;FP:17358;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGCAACCCTAAAAATATCCA
++
+EEEE/EEEEEEA/EA<EEEE/AEEEE//E//6A/AE/A/AAE/
+@NB500968:8:H5CFGAFXX:1:11101:5052:1198 1:N:0:TCCTGAGC+TACTCCTT;FP:38799;RQ:34.63
+GGTTTTTAGTGGCATTCATTGTGTGGATTCCTCTGTTGCTGCCCCTGGGCTAACAAGGGT
++
+EEEEEEE6EEEEAEEEEEEEEEEEAEEEEEEEEEEEEEAAE/AEEEEEAAAEEEE/EAEE
+@NB500968:8:H5CFGAFXX:1:11101:15805:1205 1:N:0:TCCTGAGC+TACTCCTT;FP:44542;RQ:34.63
+GAGGGCCAGGCAGCTCAAAGCCTGAAGGAGGTACTTCAAGGAGGGTCCAGAGAAGCTGGG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEE<EAAEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:11834:1217 1:N:0:TCCTGAGC+TACTCCTT;FP:56795;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGGATTCCCCTATATAGCGACCC
++
+EEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEEEEA<EAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:20105:1230 1:N:0:TCCTGAGC+TACTCCTT;FP:2265;RQ:34.63
+GTTGTTTCATGTTGCAATTTTTAGGTTACTTTTGAGTTGCTTAGGAATTGTTTTATTTAA
++
+EEEEEEEEEE/E/AEEEEE<E</EEAEEEEEEAAEAEEEAAE//EA//</E6EAAEA/6/
+@NB500968:8:H5CFGAFXX:1:11101:14524:1264 1:N:0:TCCTGAGC+TACTCCTT;FP:7937;RQ:2.65
+CCCAACTTAATATTATAAAATTATCATAAAAACAACTACGAGCCCACGATACAACTGATA
++
+//EEEE<///A/AAA/A////E/<///E/6/E/E/E6///E//6E<//6/AA//E//A//
+@NB500968:8:H5CFGAFXX:1:11101:14981:1265 1:N:0:TCCTGAGC+TACTCCTT;FP:26333;RQ:34.63
+GTCTGGTTAATTCCGATAACGAACGAGACTCTGGCATGCTAACTAGTTACGCGACCCCCG
++
+EEEEEEEEEEEEEAEEEEEEEEEAAEEEEEEEA<EEEEEEEEE6/E66A<A/<AA<A//6
+@NB500968:8:H5CFGAFXX:1:11101:5538:1286 1:N:0:TCCTGAGC+TACTCCTT;FP:5079;RQ:34.62
+AGTGAACTAACTAGAAAACTTCTAACTAAAAGAATTACAGCTAGAAACCCCGAAACCAAA
++
+AAEEEEEEEEAE/AEEEEEEEAEEEEE/AE6/EEA<E/EA6EAAEEE6//6<<EEA/<<A
+@NB500968:8:H5CFGAFXX:1:11101:11473:1333 1:N:0:TCCTGAGC+TACTCCTT;FP:48246;RQ:35.03
+GACAGCGTCAGATGTGTATAAGAGACAGGGTATC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:2141:1369 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:19474:1422 1:N:0:TCCTGAGC+TACTCCTT;FP:5233;RQ:35.03
+GGAGCGTCAGATGTGTATAAGAGACAGGTACCC
++
+EEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12411:1441 1:N:0:TCCTGAGC+TACTCCTT;FP:46039;RQ:34.83
+GGGAGCGTCAGATGTGTATAAGAGACAGTCTCCCTATATTTTGGC
++
+//AEEE/EAEAEE/AEAE<EEEEAE/EAE<E/E</EEEE/6E/EE
+@NB500968:8:H5CFGAFXX:1:11101:5889:1458 1:N:0:TCCTGAGC+TACTCCTT;FP:52002;RQ:34.83
+GAGCGTCAGATGTGTATAAGAGACAGCGTATCCCTATACCCCTTCGCTG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:8662:1469 1:N:0:TCCTGAGC+TACTCCTT;FP:62135;RQ:35.03
+GGAGCGTCAGATGTGTATAAGAGACAGGTTATC
++
+EEEEEEEEEEEEEAEAE<AEAEEEEEEEEEE<E
+@NB500968:8:H5CFGAFXX:1:11101:9851:1471 1:N:0:TCCTGAGC+TACTCCTT;FP:43356;RQ:3.17
+CCCCTGTCTCTTATACCNNN
++
+EEEEEA/AA/EEEEE//###
+@NB500968:8:H5CFGAFXX:1:11101:5425:1488 1:N:0:TCCTGAGC+TACTCCTT;FP:13390;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGACTCCCCTATATGGACACA
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:21998:1542 1:N:0:TCCTGAGC+TACTCCTT;FP:15180;RQ:35.09
+GGCGTCAGATGTGTATAAGAGACAGGG
++
+EEEEEEEAEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:2814:1550 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:8024:1581 1:N:0:TCCTGAGC+TACTCCTT;FP:38311;RQ:34.63
+GTGACCCGCCGCTTCAGGGACCCGCTTCTCTCATCTTCCCCAAAGTTCTGGACAGATGAT
++
+EEEEEEEAEEEEEAEAEEEEEEEEEEAEEEEAEEEAEAEEEEEEEEAEEEEEEEEEE/A<
+@NB500968:8:H5CFGAFXX:1:11101:8843:1628 1:N:0:TCCTGAGC+TACTCCTT;FP:44455;RQ:34.63
+GGCGTTCAGCCACCCGAGATTGAGCAATAACAGGTCTGTGATGCCCTTAGATGTCCGGGG
++
+EEEEEEEEEAEAEEEEEEEEEEE6/AEEEEAAEE6EA<EEEEE<EEEEEEEEE/AAAAEA
+@NB500968:8:H5CFGAFXX:1:11101:23474:1673 1:N:0:TCCTGAGC+TACTCCTT;FP:2968;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGAGGCGCCCTATAGAATATAT
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEE/EEEEE
+@NB500968:8:H5CFGAFXX:1:11101:1482:1752 1:N:0:TCCTGAGC+TACTCCTT;FP:20529;RQ:35.03
+GGGAGCGTCAGATGTGTATAAGAGACAGTTTTACCC
++
+E/EEEEEEEEEE</EEEEEEEEEEEEEEEEE<EAEA
+@NB500968:8:H5CFGAFXX:1:11101:21720:1786 1:N:0:TCCTGAGC+TACTCCTT;FP:16257;RQ:35.14
+GGGAGCGTCAGATGTGTATAAGAGAC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:9012:1809 1:N:0:TCCTGAGC+TACTCCTT;FP:4167;RQ:34.83
+GGGAGCGTCAGATGTGTATAAGAGACAGGTACCCCCTATATCATAGGC
++
+EEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12567:1818 1:N:0:TCCTGAGC+TACTCCTT;FP:61471;RQ:34.63
+GAGCCGCCTGGATACCGCAGCTAGGAATAATGGAATAGGACCGCGGTTCTATTTTGTTGG
++
+EEEEE/EA/EEEAEEEEEA/EEAEEEEE/EEEA/EEAEAEE/A<E/EEEAA<<A6AA<//
+@NB500968:8:H5CFGAFXX:1:11101:20110:1819 1:N:0:TCCTGAGC+TAATCCTT;FP:64559;RQ:34.63
+GGATTCTTGAAGTATAACAGGAGAATCAAATAATTAAAGAGGGGGCTGGAGAGATAGCT
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:24842:1887 1:N:0:TCCTGAGC+TACTCCTT;FP:13772;RQ:35.14
+GGGAGCGTCAGATGTGTATAAGAG
++
+EEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:16709:1904 1:N:0:TCCTGAGC+TACTCCTT;FP:5671;RQ:34.73
+GTTCCATGCGTCAATGTTTCTTCTGGATACAGCAGGCCTGTTCCCGGGAGTCAATGC
++
+EEEEEAEEEEEEEEEEAEEEEEEEEEEEAEE6<EEEEEAAA<EAEEEEEEEEEEE<E
+@NB500968:8:H5CFGAFXX:1:11101:12189:1919 1:N:0:TCCTGAGC+TACTCCTT;FP:53204;RQ:34.62
+ACACCCTTTATAGTGATGGATTGGATCTGTTTGAAGGCAAAGTCCTTCCCAGACTTATCA
++
+EEEEEEE/EEE/EEAEEEEEEEEEEE/EE<EAAAE/E<</EA///A</E/A/AAA<A/</
+@NB500968:8:H5CFGAFXX:1:11101:26750:1960 1:N:0:TCCTGAGC+TACTCCTT;FP:578;RQ:35.03
+GGGAGCGTCAGATGTGTATAAGAGACAGCTGACCC
++
+EEEAEEEAEEAEAEAEE/EEEEAEAAEEEEE/EEE
+@NB500968:8:H5CFGAFXX:1:11101:21600:1964 1:N:0:TCCTGAGC+TACTCCTT;FP:38123;RQ:18.05
+GAACGAGACTCTGGCATGCTAACTAGTTACGCGACCCCCGAGCGGTCGGCGTCCCCCAA
++
+EEEEEEAEE/AAAEEEEEEEEEEEEEE<EEEAEEEEAAAEAAEEEE<EA<A<EA<<E<E
+@NB500968:8:H5CFGAFXX:1:11101:11625:2006 1:N:0:TCCTGAGC+TACTCCTT;FP:4030;RQ:35.03
+GACAGCGTCAGATGTGTATAAGAGACAGCAACCACC
++
+AEA66E/<AEAEE/E/EEAEEEEAEEEEEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10474:2017 1:N:0:TCCTGAGC+TACTCCTT;FP:16506;RQ:35.03
+GGAGCGTCAGATGTGTATAAGAGACAGTCGTCC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:17381:2018 1:N:0:TCCTGAGC+TACTCCTT;FP:64754;RQ:35.18
+AACGTCAGATGTGTATAAGN
++
+EEE/EEEAEEEEEAE</EA#
+@NB500968:8:H5CFGAFXX:1:11101:9520:2027 1:N:0:TCCTGAGC+TACTCCTT;FP:57159;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGCTACATCCCTATAAGGATGTAC
++
+EEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEE/EEEEE<AEEEA
+@NB500968:8:H5CFGAFXX:1:11101:23302:2031 1:N:0:TCCTGAGC+TACTCCTT;FP:9055;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGATCGTCCCTATAGCTCCAATC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEAEE
+@NB500968:8:H5CFGAFXX:1:11101:19677:2089 1:N:0:TCCTGAGC+TACTCCTT;FP:14776;RQ:34.63
+GGCGCCAGTCCTCCGATTGACTGAGTCGCCCGGGTACCCGTGTATCCAATAAACCCTCTT
++
+EEEEEEEAEEEEEEEAE/EEEEEEEEAEEE//E/EEE/A<AAAAAA<AEAE/EE/AAA<6
+@NB500968:8:H5CFGAFXX:1:11101:26158:2095 1:N:0:TCCTGAGC+TACTCCTT;FP:23868;RQ:34.83
+GGGAGCGTCAGATGTGTATAAGAGACAGCCACTCCCTATAGGATTTCC
++
+EEEEEEEEEAEEEEEEEEEEEEEEEEEEA<EEEEAEEA<6EEEEE<EA
+@NB500968:8:H5CFGAFXX:1:11101:7883:2101 1:N:0:TCCGGAGC+TACTCCTT;FP:59501;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGCATCACCCTATAGGCTGTCT
++
+EEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEE<EEEEEEAE
+@NB500968:8:H5CFGAFXX:1:11101:24496:2131 1:N:0:TCCTGAGC+TACTCCTT;FP:247;RQ:35.19
+GGGAGCGTCAGATGTGTATAAG
++
+EEEEEEEEEEE/EEEEEEE/AE
+@NB500968:8:H5CFGAFXX:1:11101:5814:2162 1:N:0:TCCTGAGC+TACTCCTT;FP:5751;RQ:34.63
+GCCCTTGGCTTCACCACACTTTCTACTGGAAATGTGCCTTTTCCTCCATTGCTGCTAC
++
+EEEEEEEEEEEEEEEAAEEEAEEEEEA/EEEEEE<EEEAEAAEEAEEE<AAEAEEAAE
+@NB500968:8:H5CFGAFXX:1:11101:5796:2184 1:N:0:TCCTGAGC+TATTCCTT;FP:4250;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGGCTGTCCCTATAGACTGACT
++
+EEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEAEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:23937:2223 1:N:0:TCCTGAGC+TACTCCTT;FP:26022;RQ:34.63
+GCTTCTTCTGAAGAGCTCCACCTGTGGTCTCTCCACCTCTCGAGGCTCTTCACTAGCCT
++
+EEEEAEA/AEEEEEEEA/</</</A//EAEEE<AE66A666/A/E6E/AAAAA////E/
+@NB500968:8:H5CFGAFXX:1:11101:25267:2246 1:N:0:TCCTGAGC+TACTCCTT;FP:4173;RQ:34.83
+GAGCGTCAGATGTGTATAAGAGACAGGATGCCCCTATATATAGAGGCTG
++
+E6EEEEAAEEEEEEEAEEEEEEEEEEEEEEE6EEEEEEEEEAEAEEAEA
+@NB500968:8:H5CFGAFXX:1:11101:19822:2253 1:N:0:TCCTGAGC+TACTCCTT;FP:41486;RQ:34.63
+GACGTGCAAATCGGTCGTCCGACCTGGGTATAGGGGCGAAAGACTAATCGAACCATCTAG
++
+EEEEEEEEEEEEEAEEEEEEEEEEEEAEEEEEAEAEE/EEEEEEEEE<<EEEEAEEEA//
+@NB500968:8:H5CFGAFXX:1:11101:6670:2262 1:N:0:TCCTGAGC+TACTCCTT;FP:14539;RQ:3.28
+GGAGGCCAAATTGTGCTGATTTTCCTGTAGCTGCGATTAATAGGCCTATAAGTGGAATTA
++
+EEEEEEEEEEAEAEEAEEEEEEEE/EAEEEAEEAEEAE/EEAE/</E/EA66A<//<EE/
+@NB500968:8:H5CFGAFXX:1:11101:14858:2304 1:N:0:TCCTGAGC+TACTCCTT;FP:4202;RQ:34.63
+GTCTCCGGCTTGTGTAGACTCTGCCTGGAAGTGAAGTTGTGGACTCCCGACCTCCAGTC
++
+EEEEEEEEEEEEEEEEEEEEEA/E/AEEEEEAEEEEAE/<AAEEE<EAEEE<EEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10463:2345 1:N:0:GCCTGAGC+TACTCCTT;FP:24337;RQ:34.81
+TCCAGCGTCAGATGTGTATAAGAGACAGAGTTCCCCTATAGAGAGACG
++
+EEEEEEEEEEE6EEEEA/A//E/EEEEEEEEEEEEEEAA/AAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:24405:2351 1:N:0:TCCTGAGC+TACTCCTT;FP:5580;RQ:35.18
+TCTGTCTCTTATACCNNNNN
++
+EEEEEEAEEEEEEEE#####
+@NB500968:8:H5CFGAFXX:1:11101:6508:2370 1:N:0:TCCTGAGC+TACTCCTT;FP:64422;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGCCCAGCCCTATACCGTCTAGCTG
++
+EEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:21730:2488 1:N:0:TCGTGAGC+TACTCCTT;FP:57284;RQ:35.18
+AGCGTCAGATGTGTATCNNN
++
+EEEEEEEEEEEEEAEAE###
+@NB500968:8:H5CFGAFXX:1:11101:5499:2506 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:20021:2516 1:N:0:TCCTGGGC+TACTCCTT;FP:28351;RQ:34.83
+GGAGCGTCAGATGAGTATAAGAGACAGTCATCCCCTATAGTCGTTCA
++
+EAEEAAEEAAE///E/EAAE/EAEE//EEE/EEEEEE/EE/EA6EEE
+@NB500968:8:H5CFGAFXX:1:11101:4827:2579 1:N:0:TCCTGAGC+TACTCCTT;FP:8626;RQ:34.83
+GACAGCGTCAGATGTGTATAAGAGACAGGGAATCCCTATAGACTTTGC
++
+EEEEEEEEE/E<AAAEEEEEEEE/E/EEEEEEEEEEEAEEAEAAEEE/
+@NB500968:8:H5CFGAFXX:1:11101:25015:2619 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:11758:2662 1:N:0:TCCTGAGC+TACTCCTT;FP:32308;RQ:34.83
+GACAGCGTCAGATGTGTATAAGAGACAGAGATGCCCTATAGACCGAAC
++
+EEEEEE/EAE/EEE/EEEAEE6E6EEEEEAAEA/E/AEEEA/E6EA/E
+@NB500968:8:H5CFGAFXX:1:11101:25905:2738 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:23131:2747 1:N:0:TCCTGAGC+TACTCCTT;FP:61606;RQ:34.61
+CAGACGTTCGAATGGGTCGTCGCCGCCACGGGGGGCGTGCGATCGGCCCGAGGTTATCTA
++
+EEEEEEAEAEEEE/EEEAEEEEA/EEEEAEA/EA<EAAAEEAE<AE/A<<EEEA<<EEA<
+@NB500968:8:H5CFGAFXX:1:11101:20373:2751 1:N:0:TCCTGATC+TACTCCTT;FP:283;RQ:34.93
+GGAGCGTCAGATGTGTATAAGAGACAGCCCTATACAGTGATC
++
+EEEEEEEEE<EEEEAEAEEEEEEEEEEEEEE/EEEEEE<EAA
+@NB500968:8:H5CFGAFXX:1:11101:21204:2795 1:N:0:TCCTGAGC+TACTCCTT;FP:16279;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGCTCCACCCTATAGGCCCTAACTG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10694:2817 1:N:0:TCCTGAGC+TACTCCTT;FP:10036;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGCCCTTCCCTATACGCTTGGCC
++
+EEEEEEEEEEEEE<EEEEAEEEEEEEEEE/A<E<EEEA/EEEAEEEE/
+@NB500968:8:H5CFGAFXX:1:11101:8407:2853 1:N:0:TCCTGAGC+TACTCCTT;FP:21171;RQ:34.62
+AGAGGGACAAGTGGCGTTCAGCCACCCGAGATTGAGCAATAACAGGTCTGTGATGCCCTT
++
+/EEEEEEAEEEEEAEEEEEEEEEAEAEEAAEEEE/EEEAAAEA//EEEAEAAA/EAEAE6
+@NB500968:8:H5CFGAFXX:1:11101:17699:2880 1:N:0:TCCTGAGC+TACTCCTT;FP:52710;RQ:34.62
+CTTTGGAATCTTCACAAGCAGAACCAGTCAAAACCCCGGCAAGCACAAAGAGACTCTCCA
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEE6AEAEAEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:2261:2941 1:N:0:TCCTGAGC+TACTCCTT;FP:2596;RQ:34.63
+GGCAGCATCAGATGTGTATAAGAGACAGAAGGAGCAGAGAGCGGCCAGCAACCGATGAGA
++
+EEEEEEEEEEEEEEAEEEAEEEEE///////EA//E/E/A/////E<EE//</E<<////
+@NB500968:8:H5CFGAFXX:1:11101:11672:2944 1:N:0:TCCTGAGC+AACTCCTT;FP:32501;RQ:34.63
+GGGTGAAATACTGTGAGGAACCCAGAAGTGAGTGGAAAGATAACTGGCTTTTTTTTTTTT
++
+AEEEEEEEEEEEAAAEEEEEAA<E/EEEEEEEAAEAEEEE<EEE<</EE6AA/6/<66//
+@NB500968:8:H5CFGAFXX:1:11101:15968:2945 1:N:0:TCCTGAGC+TACTCCTT;FP:8706;RQ:34.63
+GCAAATCGGTCGTCCGACCTGGGTATAGGGGCGAAAGGCTAATCGAACCATCTAGTAG
++
+EEEEEEEAEEEEEEEEEEEEEEEEAEEEEEEAEEEEE<EEEEAAEAEEAE<AEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:25690:2959 1:N:0:TCCTGCGC+TACTCCTT;FP:27004;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGGCCTCCCCTATAGCTTCGGA
++
+AE6EEEEEEEEEAAEEEEEEEEEEEEEEAEE/E/EAEEEEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:8090:2972 1:N:0:TCCTGAGC+TACTCCTT;FP:12902;RQ:34.81
+CAACGTCAGATGTGTATAAGAGACAGGGACCCCCTATAATTTATGCC
++
+EEEEEEEEEEEEEEEEEAAEEEEEEEEEAEE<EEEEEAE<EEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:22125:3141 1:N:0:TCCTGAGC+TACTCCTT;FP:65197;RQ:35.07
+AGCGTCAGATGTGTATAAGAGACAGCTC
++
+EEEEEEEEEEEEEEEEEEEEEEEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:6707:3162 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:24232:3202 1:N:0:TCCTGAGC+TACTCCTT;FP:19483;RQ:34.63
+GTTGGTTTTCGGAACTGAGGCCATGATTAAGATGGACGGCCGGGGGCATTCGTATTGCG
++
+EEEEEEEEEEAEEEEEEE/E<EEEEEA</EEE/A/EE<E666<EEAEEEE/6/A66/66
+@NB500968:8:H5CFGAFXX:1:11101:16936:3211 1:N:0:TCCTGAGC+TACTCCTT;FP:7401;RQ:34.83
+GAGCGTCAGATGTGTATAAGAGACAGTGCGTCCCTATAAAAGTGAACTG
++
+EEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEA/EEEAEAEEE
+@NB500968:8:H5CFGAFXX:1:11101:6680:3212 1:N:0:TCCTGAGC+TACTCCTT;FP:62506;RQ:34.81
+TCCAGCGTCAGATGTGTATAAGAGACAGCTCCTCCCTATAAG
++
+EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:6132:3257 1:N:0:TCCTGAGC+TACTCCTT;FP:7557;RQ:34.62
+CGTCAGATGTGTATAAGAGACAGCATTCCCCTATACACGCCTCCTGTCTCTTATAC
++
+EEEEAEEEEEEEEEEEEEEAEAEEAEEEEAEEEEEEEEEEEEEEEEAE/EEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:3840:3285 1:N:0:TCCTGAGC+TACTCCTT;FP:16156;RQ:34.93
+GGGAGCGTCAGATGTGTATAAGAGACAGGTATACCCTATAG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:2858:3312 1:N:0:TCCTGAGC+TACTCCTT;FP:14452;RQ:35.12
+AACGTCAGATGTGTATAAGAG
++
+EE/AE/E/EAAEEEEEEEAEE
+@NB500968:8:H5CFGAFXX:1:11101:3641:3398 1:N:0:TCCTGAGC+TACTCCTT;FP:36051;RQ:35.14
+GCAGCGTCAGATGTGTATAAGAG
++
+EEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:20540:3401 1:N:0:TCCTGAGC+TACTCCTT;FP:48603;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGATCAACCCTATATGCCCTAACTG
++
+E/E///E/E</A/EEA<E/E/EEE/E</E///A/E6<E/A/<E/EEE/
+@NB500968:8:H5CFGAFXX:1:11101:19438:3463 1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:4772:3477 1:N:0:TCCTGAGC+AACTCCTT;FP:61788;RQ:3.33
+TAACTGTCTCTTATACACCN
++
+EEEEEEEAEEEAEEE<EEE#
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/read1_extracted.fastq	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,328 @@
+@NB500968:8:H5CFGAFXX:1:11101:15939:1071___1:N:0:TCCTGAGC+TACTCCTT;FP:23586;RQ:34.63
+GGCACGTCAGATGAGTATAAGAGACAGCAAATAGGTATAGGGGGCACGTCAGATGAGTAA
++
+EAEEEEEEA///A/EAEEE/EAEEEEEEEEEEAEEE//E/AA/E6AE/EE/<E////A//
+@NB500968:8:H5CFGAFXX:1:11101:1121:1139___1:N:0:TCGTGAGC+TACTCCTT;FP:27828;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGGTATTCCCTATAGTCTACATC
++
+EEEEEEEEEEEEEEEAAEEEEAEEEAEEE/EEE/EE<AAEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:7529:1175___1:N:0:TCCTGAGC+TACTCCTT;FP:17358;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGCAACCCTAAAAATATCCA
++
+EEEE/EEEEEEA/EA<EEEE/AEEEE//E//6A/AE/A/AAE/
+@NB500968:8:H5CFGAFXX:1:11101:5052:1198___1:N:0:TCCTGAGC+TACTCCTT;FP:38799;RQ:34.63
+GGTTTTTAGTGGCATTCATTGTGTGGATTCCTCTGTTGCTGCCCCTGGGCTAACAAGGGT
++
+EEEEEEE6EEEEAEEEEEEEEEEEAEEEEEEEEEEEEEAAE/AEEEEEAAAEEEE/EAEE
+@NB500968:8:H5CFGAFXX:1:11101:15805:1205___1:N:0:TCCTGAGC+TACTCCTT;FP:44542;RQ:34.63
+GAGGGCCAGGCAGCTCAAAGCCTGAAGGAGGTACTTCAAGGAGGGTCCAGAGAAGCTGGG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEE<EAAEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:11834:1217___1:N:0:TCCTGAGC+TACTCCTT;FP:56795;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGGATTCCCCTATATAGCGACCC
++
+EEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEEEEA<EAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:20105:1230___1:N:0:TCCTGAGC+TACTCCTT;FP:2265;RQ:34.63
+GTTGTTTCATGTTGCAATTTTTAGGTTACTTTTGAGTTGCTTAGGAATTGTTTTATTTAA
++
+EEEEEEEEEE/E/AEEEEE<E</EEAEEEEEEAAEAEEEAAE//EA//</E6EAAEA/6/
+@NB500968:8:H5CFGAFXX:1:11101:14524:1264___1:N:0:TCCTGAGC+TACTCCTT;FP:7937;RQ:2.65
+CCCAACTTAATATTATAAAATTATCATAAAAACAACTACGAGCCCACGATACAACTGATA
++
+//EEEE<///A/AAA/A////E/<///E/6/E/E/E6///E//6E<//6/AA//E//A//
+@NB500968:8:H5CFGAFXX:1:11101:14981:1265___1:N:0:TCCTGAGC+TACTCCTT;FP:26333;RQ:34.63
+GTCTGGTTAATTCCGATAACGAACGAGACTCTGGCATGCTAACTAGTTACGCGACCCCCG
++
+EEEEEEEEEEEEEAEEEEEEEEEAAEEEEEEEA<EEEEEEEEE6/E66A<A/<AA<A//6
+@NB500968:8:H5CFGAFXX:1:11101:5538:1286___1:N:0:TCCTGAGC+TACTCCTT;FP:5079;RQ:34.62
+AGTGAACTAACTAGAAAACTTCTAACTAAAAGAATTACAGCTAGAAACCCCGAAACCAAA
++
+AAEEEEEEEEAE/AEEEEEEEAEEEEE/AE6/EEA<E/EA6EAAEEE6//6<<EEA/<<A
+@NB500968:8:H5CFGAFXX:1:11101:11473:1333___1:N:0:TCCTGAGC+TACTCCTT;FP:48246;RQ:35.03
+GACAGCGTCAGATGTGTATAAGAGACAGGGTATC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:2141:1369___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:19474:1422___1:N:0:TCCTGAGC+TACTCCTT;FP:5233;RQ:35.03
+GGAGCGTCAGATGTGTATAAGAGACAGGTACCC
++
+EEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12411:1441___1:N:0:TCCTGAGC+TACTCCTT;FP:46039;RQ:34.83
+GGGAGCGTCAGATGTGTATAAGAGACAGTCTCCCTATATTTTGGC
++
+//AEEE/EAEAEE/AEAE<EEEEAE/EAE<E/E</EEEE/6E/EE
+@NB500968:8:H5CFGAFXX:1:11101:5889:1458___1:N:0:TCCTGAGC+TACTCCTT;FP:52002;RQ:34.83
+GAGCGTCAGATGTGTATAAGAGACAGCGTATCCCTATACCCCTTCGCTG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:8662:1469___1:N:0:TCCTGAGC+TACTCCTT;FP:62135;RQ:35.03
+GGAGCGTCAGATGTGTATAAGAGACAGGTTATC
++
+EEEEEEEEEEEEEAEAE<AEAEEEEEEEEEE<E
+@NB500968:8:H5CFGAFXX:1:11101:9851:1471___1:N:0:TCCTGAGC+TACTCCTT;FP:43356;RQ:3.17
+CCCCTGTCTCTTATACCNNN
++
+EEEEEA/AA/EEEEE//###
+@NB500968:8:H5CFGAFXX:1:11101:5425:1488___1:N:0:TCCTGAGC+TACTCCTT;FP:13390;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGACTCCCCTATATGGACACA
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:21998:1542___1:N:0:TCCTGAGC+TACTCCTT;FP:15180;RQ:35.09
+GGCGTCAGATGTGTATAAGAGACAGGG
++
+EEEEEEEAEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:2814:1550___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:8024:1581___1:N:0:TCCTGAGC+TACTCCTT;FP:38311;RQ:34.63
+GTGACCCGCCGCTTCAGGGACCCGCTTCTCTCATCTTCCCCAAAGTTCTGGACAGATGAT
++
+EEEEEEEAEEEEEAEAEEEEEEEEEEAEEEEAEEEAEAEEEEEEEEAEEEEEEEEEE/A<
+@NB500968:8:H5CFGAFXX:1:11101:8843:1628___1:N:0:TCCTGAGC+TACTCCTT;FP:44455;RQ:34.63
+GGCGTTCAGCCACCCGAGATTGAGCAATAACAGGTCTGTGATGCCCTTAGATGTCCGGGG
++
+EEEEEEEEEAEAEEEEEEEEEEE6/AEEEEAAEE6EA<EEEEE<EEEEEEEEE/AAAAEA
+@NB500968:8:H5CFGAFXX:1:11101:23474:1673___1:N:0:TCCTGAGC+TACTCCTT;FP:2968;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGAGGCGCCCTATAGAATATAT
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEE/EEEEE
+@NB500968:8:H5CFGAFXX:1:11101:1482:1752___1:N:0:TCCTGAGC+TACTCCTT;FP:20529;RQ:35.03
+GGGAGCGTCAGATGTGTATAAGAGACAGTTTTACCC
++
+E/EEEEEEEEEE</EEEEEEEEEEEEEEEEE<EAEA
+@NB500968:8:H5CFGAFXX:1:11101:21720:1786___1:N:0:TCCTGAGC+TACTCCTT;FP:16257;RQ:35.14
+GGGAGCGTCAGATGTGTATAAGAGAC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:9012:1809___1:N:0:TCCTGAGC+TACTCCTT;FP:4167;RQ:34.83
+GGGAGCGTCAGATGTGTATAAGAGACAGGTACCCCCTATATCATAGGC
++
+EEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12567:1818___1:N:0:TCCTGAGC+TACTCCTT;FP:61471;RQ:34.63
+GAGCCGCCTGGATACCGCAGCTAGGAATAATGGAATAGGACCGCGGTTCTATTTTGTTGG
++
+EEEEE/EA/EEEAEEEEEA/EEAEEEEE/EEEA/EEAEAEE/A<E/EEEAA<<A6AA<//
+@NB500968:8:H5CFGAFXX:1:11101:20110:1819___1:N:0:TCCTGAGC+TAATCCTT;FP:64559;RQ:34.63
+GGATTCTTGAAGTATAACAGGAGAATCAAATAATTAAAGAGGGGGCTGGAGAGATAGCT
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:24842:1887___1:N:0:TCCTGAGC+TACTCCTT;FP:13772;RQ:35.14
+GGGAGCGTCAGATGTGTATAAGAG
++
+EEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:16709:1904___1:N:0:TCCTGAGC+TACTCCTT;FP:5671;RQ:34.73
+GTTCCATGCGTCAATGTTTCTTCTGGATACAGCAGGCCTGTTCCCGGGAGTCAATGC
++
+EEEEEAEEEEEEEEEEAEEEEEEEEEEEAEE6<EEEEEAAA<EAEEEEEEEEEEE<E
+@NB500968:8:H5CFGAFXX:1:11101:12189:1919___1:N:0:TCCTGAGC+TACTCCTT;FP:53204;RQ:34.62
+ACACCCTTTATAGTGATGGATTGGATCTGTTTGAAGGCAAAGTCCTTCCCAGACTTATCA
++
+EEEEEEE/EEE/EEAEEEEEEEEEEE/EE<EAAAE/E<</EA///A</E/A/AAA<A/</
+@NB500968:8:H5CFGAFXX:1:11101:26750:1960___1:N:0:TCCTGAGC+TACTCCTT;FP:578;RQ:35.03
+GGGAGCGTCAGATGTGTATAAGAGACAGCTGACCC
++
+EEEAEEEAEEAEAEAEE/EEEEAEAAEEEEE/EEE
+@NB500968:8:H5CFGAFXX:1:11101:21600:1964___1:N:0:TCCTGAGC+TACTCCTT;FP:38123;RQ:18.05
+GAACGAGACTCTGGCATGCTAACTAGTTACGCGACCCCCGAGCGGTCGGCGTCCCCCAA
++
+EEEEEEAEE/AAAEEEEEEEEEEEEEE<EEEAEEEEAAAEAAEEEE<EA<A<EA<<E<E
+@NB500968:8:H5CFGAFXX:1:11101:11625:2006___1:N:0:TCCTGAGC+TACTCCTT;FP:4030;RQ:35.03
+GACAGCGTCAGATGTGTATAAGAGACAGCAACCACC
++
+AEA66E/<AEAEE/E/EEAEEEEAEEEEEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10474:2017___1:N:0:TCCTGAGC+TACTCCTT;FP:16506;RQ:35.03
+GGAGCGTCAGATGTGTATAAGAGACAGTCGTCC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:17381:2018___1:N:0:TCCTGAGC+TACTCCTT;FP:64754;RQ:35.18
+AACGTCAGATGTGTATAAGN
++
+EEE/EEEAEEEEEAE</EA#
+@NB500968:8:H5CFGAFXX:1:11101:9520:2027___1:N:0:TCCTGAGC+TACTCCTT;FP:57159;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGCTACATCCCTATAAGGATGTAC
++
+EEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEEE/EEEEE<AEEEA
+@NB500968:8:H5CFGAFXX:1:11101:23302:2031___1:N:0:TCCTGAGC+TACTCCTT;FP:9055;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGATCGTCCCTATAGCTCCAATC
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEAEEEEEEAEE
+@NB500968:8:H5CFGAFXX:1:11101:19677:2089___1:N:0:TCCTGAGC+TACTCCTT;FP:14776;RQ:34.63
+GGCGCCAGTCCTCCGATTGACTGAGTCGCCCGGGTACCCGTGTATCCAATAAACCCTCTT
++
+EEEEEEEAEEEEEEEAE/EEEEEEEEAEEE//E/EEE/A<AAAAAA<AEAE/EE/AAA<6
+@NB500968:8:H5CFGAFXX:1:11101:26158:2095___1:N:0:TCCTGAGC+TACTCCTT;FP:23868;RQ:34.83
+GGGAGCGTCAGATGTGTATAAGAGACAGCCACTCCCTATAGGATTTCC
++
+EEEEEEEEEAEEEEEEEEEEEEEEEEEEA<EEEEAEEA<6EEEEE<EA
+@NB500968:8:H5CFGAFXX:1:11101:7883:2101___1:N:0:TCCGGAGC+TACTCCTT;FP:59501;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGCATCACCCTATAGGCTGTCT
++
+EEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEE<EEEEEEAE
+@NB500968:8:H5CFGAFXX:1:11101:24496:2131___1:N:0:TCCTGAGC+TACTCCTT;FP:247;RQ:35.19
+GGGAGCGTCAGATGTGTATAAG
++
+EEEEEEEEEEE/EEEEEEE/AE
+@NB500968:8:H5CFGAFXX:1:11101:5814:2162___1:N:0:TCCTGAGC+TACTCCTT;FP:5751;RQ:34.63
+GCCCTTGGCTTCACCACACTTTCTACTGGAAATGTGCCTTTTCCTCCATTGCTGCTAC
++
+EEEEEEEEEEEEEEEAAEEEAEEEEEA/EEEEEE<EEEAEAAEEAEEE<AAEAEEAAE
+@NB500968:8:H5CFGAFXX:1:11101:5796:2184___1:N:0:TCCTGAGC+TATTCCTT;FP:4250;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGGCTGTCCCTATAGACTGACT
++
+EEEEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEAEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:23937:2223___1:N:0:TCCTGAGC+TACTCCTT;FP:26022;RQ:34.63
+GCTTCTTCTGAAGAGCTCCACCTGTGGTCTCTCCACCTCTCGAGGCTCTTCACTAGCCT
++
+EEEEAEA/AEEEEEEEA/</</</A//EAEEE<AE66A666/A/E6E/AAAAA////E/
+@NB500968:8:H5CFGAFXX:1:11101:25267:2246___1:N:0:TCCTGAGC+TACTCCTT;FP:4173;RQ:34.83
+GAGCGTCAGATGTGTATAAGAGACAGGATGCCCCTATATATAGAGGCTG
++
+E6EEEEAAEEEEEEEAEEEEEEEEEEEEEEE6EEEEEEEEEAEAEEAEA
+@NB500968:8:H5CFGAFXX:1:11101:19822:2253___1:N:0:TCCTGAGC+TACTCCTT;FP:41486;RQ:34.63
+GACGTGCAAATCGGTCGTCCGACCTGGGTATAGGGGCGAAAGACTAATCGAACCATCTAG
++
+EEEEEEEEEEEEEAEEEEEEEEEEEEAEEEEEAEAEE/EEEEEEEEE<<EEEEAEEEA//
+@NB500968:8:H5CFGAFXX:1:11101:6670:2262___1:N:0:TCCTGAGC+TACTCCTT;FP:14539;RQ:3.28
+GGAGGCCAAATTGTGCTGATTTTCCTGTAGCTGCGATTAATAGGCCTATAAGTGGAATTA
++
+EEEEEEEEEEAEAEEAEEEEEEEE/EAEEEAEEAEEAE/EEAE/</E/EA66A<//<EE/
+@NB500968:8:H5CFGAFXX:1:11101:14858:2304___1:N:0:TCCTGAGC+TACTCCTT;FP:4202;RQ:34.63
+GTCTCCGGCTTGTGTAGACTCTGCCTGGAAGTGAAGTTGTGGACTCCCGACCTCCAGTC
++
+EEEEEEEEEEEEEEEEEEEEEA/E/AEEEEEAEEEEAE/<AAEEE<EAEEE<EEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10463:2345___1:N:0:GCCTGAGC+TACTCCTT;FP:24337;RQ:34.81
+TCCAGCGTCAGATGTGTATAAGAGACAGAGTTCCCCTATAGAGAGACG
++
+EEEEEEEEEEE6EEEEA/A//E/EEEEEEEEEEEEEEAA/AAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:24405:2351___1:N:0:TCCTGAGC+TACTCCTT;FP:5580;RQ:35.18
+TCTGTCTCTTATACCNNNNN
++
+EEEEEEAEEEEEEEE#####
+@NB500968:8:H5CFGAFXX:1:11101:6508:2370___1:N:0:TCCTGAGC+TACTCCTT;FP:64422;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGCCCAGCCCTATACCGTCTAGCTG
++
+EEEEEEEEEEEEEAEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:21730:2488___1:N:0:TCGTGAGC+TACTCCTT;FP:57284;RQ:35.18
+AGCGTCAGATGTGTATCNNN
++
+EEEEEEEEEEEEEAEAE###
+@NB500968:8:H5CFGAFXX:1:11101:5499:2506___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:20021:2516___1:N:0:TCCTGGGC+TACTCCTT;FP:28351;RQ:34.83
+GGAGCGTCAGATGAGTATAAGAGACAGTCATCCCCTATAGTCGTTCA
++
+EAEEAAEEAAE///E/EAAE/EAEE//EEE/EEEEEE/EE/EA6EEE
+@NB500968:8:H5CFGAFXX:1:11101:4827:2579___1:N:0:TCCTGAGC+TACTCCTT;FP:8626;RQ:34.83
+GACAGCGTCAGATGTGTATAAGAGACAGGGAATCCCTATAGACTTTGC
++
+EEEEEEEEE/E<AAAEEEEEEEE/E/EEEEEEEEEEEAEEAEAAEEE/
+@NB500968:8:H5CFGAFXX:1:11101:25015:2619___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:11758:2662___1:N:0:TCCTGAGC+TACTCCTT;FP:32308;RQ:34.83
+GACAGCGTCAGATGTGTATAAGAGACAGAGATGCCCTATAGACCGAAC
++
+EEEEEE/EAE/EEE/EEEAEE6E6EEEEEAAEA/E/AEEEA/E6EA/E
+@NB500968:8:H5CFGAFXX:1:11101:25905:2738___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:23131:2747___1:N:0:TCCTGAGC+TACTCCTT;FP:61606;RQ:34.61
+CAGACGTTCGAATGGGTCGTCGCCGCCACGGGGGGCGTGCGATCGGCCCGAGGTTATCTA
++
+EEEEEEAEAEEEE/EEEAEEEEA/EEEEAEA/EA<EAAAEEAE<AE/A<<EEEA<<EEA<
+@NB500968:8:H5CFGAFXX:1:11101:20373:2751___1:N:0:TCCTGATC+TACTCCTT;FP:283;RQ:34.93
+GGAGCGTCAGATGTGTATAAGAGACAGCCCTATACAGTGATC
++
+EEEEEEEEE<EEEEAEAEEEEEEEEEEEEEE/EEEEEE<EAA
+@NB500968:8:H5CFGAFXX:1:11101:21204:2795___1:N:0:TCCTGAGC+TACTCCTT;FP:16279;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGCTCCACCCTATAGGCCCTAACTG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10694:2817___1:N:0:TCCTGAGC+TACTCCTT;FP:10036;RQ:34.83
+GGAGCGTCAGATGTGTATAAGAGACAGCCCTTCCCTATACGCTTGGCC
++
+EEEEEEEEEEEEE<EEEEAEEEEEEEEEE/A<E<EEEA/EEEAEEEE/
+@NB500968:8:H5CFGAFXX:1:11101:8407:2853___1:N:0:TCCTGAGC+TACTCCTT;FP:21171;RQ:34.62
+AGAGGGACAAGTGGCGTTCAGCCACCCGAGATTGAGCAATAACAGGTCTGTGATGCCCTT
++
+/EEEEEEAEEEEEAEEEEEEEEEAEAEEAAEEEE/EEEAAAEA//EEEAEAAA/EAEAE6
+@NB500968:8:H5CFGAFXX:1:11101:17699:2880___1:N:0:TCCTGAGC+TACTCCTT;FP:52710;RQ:34.62
+CTTTGGAATCTTCACAAGCAGAACCAGTCAAAACCCCGGCAAGCACAAAGAGACTCTCCA
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<EEEEEEEEEEE6AEAEAEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:2261:2941___1:N:0:TCCTGAGC+TACTCCTT;FP:2596;RQ:34.63
+GGCAGCATCAGATGTGTATAAGAGACAGAAGGAGCAGAGAGCGGCCAGCAACCGATGAGA
++
+EEEEEEEEEEEEEEAEEEAEEEEE///////EA//E/E/A/////E<EE//</E<<////
+@NB500968:8:H5CFGAFXX:1:11101:11672:2944___1:N:0:TCCTGAGC+AACTCCTT;FP:32501;RQ:34.63
+GGGTGAAATACTGTGAGGAACCCAGAAGTGAGTGGAAAGATAACTGGCTTTTTTTTTTTT
++
+AEEEEEEEEEEEAAAEEEEEAA<E/EEEEEEEAAEAEEEE<EEE<</EE6AA/6/<66//
+@NB500968:8:H5CFGAFXX:1:11101:15968:2945___1:N:0:TCCTGAGC+TACTCCTT;FP:8706;RQ:34.63
+GCAAATCGGTCGTCCGACCTGGGTATAGGGGCGAAAGGCTAATCGAACCATCTAGTAG
++
+EEEEEEEAEEEEEEEEEEEEEEEEAEEEEEEAEEEEE<EEEEAAEAEEAE<AEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:25690:2959___1:N:0:TCCTGCGC+TACTCCTT;FP:27004;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGGCCTCCCCTATAGCTTCGGA
++
+AE6EEEEEEEEEAAEEEEEEEEEEEEEEAEE/E/EAEEEEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:8090:2972___1:N:0:TCCTGAGC+TACTCCTT;FP:12902;RQ:34.81
+CAACGTCAGATGTGTATAAGAGACAGGGACCCCCTATAATTTATGCC
++
+EEEEEEEEEEEEEEEEEAAEEEEEEEEEAEE<EEEEEAE<EEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:22125:3141___1:N:0:TCCTGAGC+TACTCCTT;FP:65197;RQ:35.07
+AGCGTCAGATGTGTATAAGAGACAGCTC
++
+EEEEEEEEEEEEEEEEEEEEEEEAEEEE
+@NB500968:8:H5CFGAFXX:1:11101:6707:3162___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:24232:3202___1:N:0:TCCTGAGC+TACTCCTT;FP:19483;RQ:34.63
+GTTGGTTTTCGGAACTGAGGCCATGATTAAGATGGACGGCCGGGGGCATTCGTATTGCG
++
+EEEEEEEEEEAEEEEEEE/E<EEEEEA</EEE/A/EE<E666<EEAEEEE/6/A66/66
+@NB500968:8:H5CFGAFXX:1:11101:16936:3211___1:N:0:TCCTGAGC+TACTCCTT;FP:7401;RQ:34.83
+GAGCGTCAGATGTGTATAAGAGACAGTGCGTCCCTATAAAAGTGAACTG
++
+EEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEA/EEEAEAEEE
+@NB500968:8:H5CFGAFXX:1:11101:6680:3212___1:N:0:TCCTGAGC+TACTCCTT;FP:62506;RQ:34.81
+TCCAGCGTCAGATGTGTATAAGAGACAGCTCCTCCCTATAAG
++
+EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:6132:3257___1:N:0:TCCTGAGC+TACTCCTT;FP:7557;RQ:34.62
+CGTCAGATGTGTATAAGAGACAGCATTCCCCTATACACGCCTCCTGTCTCTTATAC
++
+EEEEAEEEEEEEEEEEEEEAEAEEAEEEEAEEEEEEEEEEEEEEEEAE/EEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:3840:3285___1:N:0:TCCTGAGC+TACTCCTT;FP:16156;RQ:34.93
+GGGAGCGTCAGATGTGTATAAGAGACAGGTATACCCTATAG
++
+EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:2858:3312___1:N:0:TCCTGAGC+TACTCCTT;FP:14452;RQ:35.12
+AACGTCAGATGTGTATAAGAG
++
+EE/AE/E/EAAEEEEEEEAEE
+@NB500968:8:H5CFGAFXX:1:11101:3641:3398___1:N:0:TCCTGAGC+TACTCCTT;FP:36051;RQ:35.14
+GCAGCGTCAGATGTGTATAAGAG
++
+EEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:20540:3401___1:N:0:TCCTGAGC+TACTCCTT;FP:48603;RQ:34.81
+AGCGTCAGATGTGTATAAGAGACAGATCAACCCTATATGCCCTAACTG
++
+E/E///E/E</A/EEA<E/E/EEE/E</E///A/E6<E/A/<E/EEE/
+@NB500968:8:H5CFGAFXX:1:11101:19438:3463___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+NNNNNNNNNNNNNNNNNNNN
++
+####################
+@NB500968:8:H5CFGAFXX:1:11101:4772:3477___1:N:0:TCCTGAGC+AACTCCTT;FP:61788;RQ:3.33
+TAACTGTCTCTTATACACCN
++
+EEEEEEEAEEEAEEE<EEE#
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/read2.fastq	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,401 @@
+@NB500968:8:H5CFGAFXX:1:11101:15939:1071 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:1121:1139 2:N:0:TCGTGAGC+TACTCCTT
+GATGTAGACTATAGGGAATACCNNNNNNNNNNNNN
++
+A6AAAEEEEEEEEEEEAEE/E<#############
+@NB500968:8:H5CFGAFXX:1:11101:7529:1175 2:N:0:TCCTGAGC+TACTCCTT
+GATGTGTAAAAGAGACAGAGGATAATTATAGGGATCCTGTCTCTAAAACACATCTGACGCTCCCTATACATAATA
++
+/AAAAE<A/E//EEE/EE//EEEEAEEEEEEEEEE/A/<EEE<<///<AE</AE//AAA//E/////EEE<6/A/
+@NB500968:8:H5CFGAFXX:1:11101:5052:1198 2:N:0:TCCTGAGC+TACTCCTT
+GGATTGAGGCTAGAAACAGCTGCTTTTCTCACTTTCCCAGCAGCAAGACATCAATCTGTTGAGACCTCTCATTTT
++
+AA6AAEEAEEEEEEEEAAEEEEE///AEEA/<E/E<A//EEEEAE/EE//AE///E</EAEEA<EEEAEE/E/E/
+@NB500968:8:H5CFGAFXX:1:11101:15805:1205 2:N:0:TCCTGAGC+TACTCCTT
+GCGTCAGATGTGGGTAGTCAGAAAAGCAGATGGATGGGGGCCTTGCCCAACCAGATCCAGGCCGAGTTCACCT
++
+AAAAAEEAEEEEEEEEEEEEEEA//EEEEEEEEEEEEEEEEEEAEE/E//EEEEAEAE/EEEEE<EEEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:11834:1217 2:N:0:TCCTGAGC+TACTCCTT
+GGGTCGCTATATAGGGGAATCCNNNNNNNNNNNNN
++
+AAAAAEEEEEEEEEEEEEAEEE#############
+@NB500968:8:H5CFGAFXX:1:11101:20105:1230 2:N:0:TCCTGAGC+TACTCCTT
+GTACACGTATTACTTCTGTTTGGTGAATTAAATAAAACAATTCCTAAGCAACTCAAAAGTAACCTAAAAATTGCA
++
+AA/A//EEE/EAEEEEEEE/EEEAEE6EEE<EEA<6E//E/EEAEE//EE/EEE<EAEE6</EEEEE/A/EEE/A
+@NB500968:8:H5CFGAFXX:1:11101:14524:1264 2:N:0:TCCTGAGC+TACTCCTT
+GTATAAGAGACAGTTGGGACTTATAGGGTAAGTCN
++
+AAAAAEEEAEEEEAEEE/EEEEEAE<EE/EEEEE#
+@NB500968:8:H5CFGAFXX:1:11101:14981:1265 2:N:0:TCCTGAGC+TACTCCTT
+GGACTTAATCAACGCAAGCTTATGACCCGCACTTACTGGGAATTCCTCGTTCATGGGGAATAATTGCAATCCCCG
++
+AA6AAEE/EE//EEE6EEEEEEEE/AAEEAEEEE/EEEEEAEEEE/EEEEEEEEEEEE/AEE<EEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:5538:1286 2:N:0:TCCTGAGC+TACTCCTT
+ATCACCAAGCTCGTTAGGCTTTTCACCTCTACCTAAAAATCTTCTCACTATTTTGCCACATAGACGAGTTGATT
++
+AAA/AEE/EEEEEEEAEEEEEEEE/EEEEE/EEEA</EEEEEEEEAAAEAEEEEE<EAE6EEEEEEEEEEEEAE
+@NB500968:8:H5CFGAFXX:1:11101:11473:1333 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:2141:1369 2:N:0:TCCTGAGC+TACTCCTT
+GATGTGTATAAGAAACAGACCACATATATAGGGGGGAGCGTCAGATGTGTATAAGAGACAGCCCCCCCCAATAGA
++
+AAAAAEEEEEEEAEEEEE/E/EAEEEEEAEEEA/EEEEEEEAA//AEEEE/</AA/////////AAEE////AA/
+@NB500968:8:H5CFGAFXX:1:11101:24482:1380 2:N:0:TCCTGAGC+TACTCCTT
+CATAGACACTAGACAGTAGGACAGACACAGATTAGTTGCTGGCGAGCTTACCTCGCGGTGGTGGGTCGGTGGTC
++
+AA/AAE//A<EA/E//E/E66//////A/A/E//E<AEEAA/A/EE/EEEA6<E/AAEE//6///6/A/E/E6/
+@NB500968:8:H5CFGAFXX:1:11101:2671:1393 2:N:0:TCCTGAGC+TACTCCTT
+GTATAAGAGACAGAATCCCCTTATAGGGGATTCNN
++
+A///A///AA/E//EE/EEE/EEEEAAEEAEEE##
+@NB500968:8:H5CFGAFXX:1:11101:19474:1422 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:12411:1441 2:N:0:TCCTGAGC+TACTCCTT
+GACAAAATATAGGGGGACTGTCTGTTATAGACGTCGGACGCGCGCCCCTATAAGGAATAGGGGTCTCTGAGAGA
++
+A/////E////////A///E/E///EEE///A//A/A//<///////E/////AEA//////////////////
+@NB500968:8:H5CFGAFXX:1:11101:5889:1458 2:N:0:TCCTGAGC+TACTCCTT
+CAGCGAAGGGGTATAGGGATACGNNNNNNNNNNNN
++
+A/AAA//EEEEA/A/EEE////E############
+@NB500968:8:H5CFGAFXX:1:11101:8662:1469 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:9851:1471 2:N:0:TCCTGAGC+TACTCCTT
+GTATAAGAGACAGGGGCCCGTTAAAGGGCCCNNNN
++
+/A6<//EA//E/E/A/EEE/E/E/E/6EEE/####
+@NB500968:8:H5CFGAFXX:1:11101:7850:1487 2:N:0:TCCTGAGC+TACTCCTT
+CCGCTCCATTGAGCAGTTACCTAGAATCATGAGGGCCACGTCCATTTCCTCCTGGGGAATTTGTCCGAAGAT
++
+AAA<AA///EE/EE/E///EE//E///A//E/EEEEA/AEA/<//6AA</6//EEEEA//<EEE/EE//AA/
+@NB500968:8:H5CFGAFXX:1:11101:5425:1488 2:N:0:TCCTGAGC+TACTCCTT
+GTATAAGAGACAGTGTGTCCATATAGGGGAGTCNN
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEE##
+@NB500968:8:H5CFGAFXX:1:11101:21998:1542 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:2814:1550 2:N:0:TCCTGAGC+TACTCCTT
+CCGTCAGATGAGTATAAGAGACAGGTAAGACATAAAGAGCGTAGTTGTCTCTAATACACATCTGACGCTGCCGA
++
+/AAA6AEEE//E/EAAEE6/EEEEEEEE//EE/E/EE/AE/<A//</AE/E////AE<A/<A/6/E<EEAAEEE
+@NB500968:8:H5CFGAFXX:1:11101:8024:1581 2:N:0:TCCTGAGC+TACTCCTT
+GAACTATGACTGTGAGCAGGCATTGTGGTTCTGCCCAAGCTACATCCAGATCACCCAGGCCTCTGGGTCTTTTGG
++
+AAAA6EAEE/E/EEEEEEEAEEEEEEEAEEEEAEEEEEAE6EEE6EEEAA/EEEEEE/A/E/AEAEE<EEEEA/<
+@NB500968:8:H5CFGAFXX:1:11101:8843:1628 2:N:0:TCCTGAGC+TACTCCTT
+GGGCAGGGACTTAATCAACGCAAGCTTATGACCCGCACTTACTGGGAATTCCTCGTTCATGGGGAATAATTGCAA
++
+AAAAAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEE<EEEEEEEEEEEE<AA/E6EEAA6
+@NB500968:8:H5CFGAFXX:1:11101:23474:1673 2:N:0:TCCTGAGC+TACTCCTT
+GATGTGTATAAGAGACAGATATAGACTATAGGGCGCCT
++
+A/AA////A//EA/6<EE6/EEE6////A//EE/E6//
+@NB500968:8:H5CFGAFXX:1:11101:13010:1719 2:N:0:TCCTGAGC+TACTCCTT
+GAAAAGACGTATAGGGCTTAGCTGTCTCTTATGCACGTCTGACGCTGTCCGCTATATAATCACACTGTCGT
++
+A///////E////EEE////E//E////////////////E//E//E///////////////////E////
+@NB500968:8:H5CFGAFXX:1:11101:15935:1725 2:N:0:GCCTGAGC+TACTCCTT
+GCGTGGTATAGGGGAGGGCTGGCTGTTATGCGCATCTGACGGTGCCGCTATATTAGGCCTGTTTCGGATGGGCCG
++
+A6AAAEA///EEAE//////E/<//////////////E//EE/E////////////<<//E//////////////
+@NB500968:8:H5CFGAFXX:1:11101:1482:1752 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:21720:1786 2:N:0:TCCTGAGC+TACTCCTT
+GTCTCTTATACACATCTGACGCTCCCCCCTATAGTTCAAAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:9012:1809 2:N:0:TCCTGAGC+TACTCCTT
+GTGTATAAGAGACAGGCCTATGATATAGGGGGTAC
++
+AAAAAEEEEEEEEEEEEEAEE/EEEAE/AAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12567:1818 2:N:0:TCCTGAGC+TACTCCTT
+GGTCATGGGAATAACGCCGCCGCATCGCCAGTCGGCATCGTTTATGGTCGGAACTACGACGGTATCTGATCGT
++
+AA/AAEEAEEE//AEAEEEEEEEEEEAEEEEEEEE/EAEE<EE<EA/<E<EEEEE6EEEE/EEEEEE/EEAEE
+@NB500968:8:H5CFGAFXX:1:11101:20110:1819 2:N:0:TCCTGAGC+TAATCCTT
+GAAGAGGGTATCAGATCCCTTGGATCTGGAGTCGCAGATGGCTGCGAGCCACCATGTGGGTTTTGAGAATTGAAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:24842:1887 2:N:0:TCCTGAGC+TACTCCTT
+CTCATAAACACATCAGGCGGACCCCGCAAAATAAAGGATCTATCTCTGATAAAGAACTGACGCTGCCGACGAAA
++
+A///AE/6EE///////E//6E/AE/<////EE/A/A/A////EAA///////////E/EE/E/////A//E/A
+@NB500968:8:H5CFGAFXX:1:11101:12108:1890 2:N:0:TCCTGAGC+TACTCCTT
+GTACAAGCCAGATGCCGAAGTTCAAGGACAGCCAGGGATAGACACAGAGAAACCTTGTCTTGGGATCCCTATAA
++
+AAAAAEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEAEEEEEEAEEEEEEE/EEAEEEEEEEEAEEEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:16709:1904 2:N:0:TCCTGAGC+TACTCCTT
+CGGATCGCTTCTTCTGATTGAGTTGCCACTCTTTGAAATGCGGCACGGCCTCGTCCACGTTGACGATGCCCAG
++
+AAAAAEEEEEEEAEEEEEEEEEEEEEEEEEEE/EEEEAEEEEEEEEEEEEEEEEE//EEEEEEEAEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12189:1919 2:N:0:TCCTGAGC+TACTCCTT
+GTCTATAGGTATCAGTTTCACTGTAATACCTTTCCTTCCATGGCTAAGTCTGGGAAGGACTTTGCCTTCAAACAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<AEEEEEEEEE<E<EAEAEEEEEEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:26750:1960 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:21600:1964 2:N:0:TCCTGAGC+TACTCCTT
+GCAATCCCCGATCCCCATCACGAATGGGGTTCAACGGGTTACCCGCGCCTGCCGGCGTAGGGTAGGCACACG
++
+AAAAAEEEEEEAEEEEE/E/EEE6/EEEEEEEE/EEEEEEAEEEEEEAEEEAEEEEEEEEEAEE<EAEAEEE
+@NB500968:8:H5CFGAFXX:1:11101:6961:1984 2:N:0:TCCTGAGC+TACTCCTT
+CTACTGGCAGGATCAGCCAGGTAAGGAGGCGCGAGCGAAGCCGTCGGGGGTGGGACAGGACCAGGCGGGACAGAG
++
+A/6A/EEA/EE/<A///A/EEE//EE6/EEEEE/E6E//6//E//EEEEE/EEA///EE////EEEEEA/E/E/E
+@NB500968:8:H5CFGAFXX:1:11101:20646:1998 2:N:0:TCCTGAGC+TACTCCTT
+CGTCAGATGTGTATAAGAAACAGAATCTGATTATAGGGGACTCTTAGCGGTGGATCACTCGGCTCGTGCGTCGAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:11625:2006 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:10474:2017 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:17381:2018 2:N:0:TCCTGAGC+TACTCCTT
+CTTATACACATCTGACGTTCCCTATACTAATAAAN
++
+A/AAAAE6E////E/AEAE/AEEAEA/EA/EEEE#
+@NB500968:8:H5CFGAFXX:1:11101:9520:2027 2:N:0:TCCTGAGC+TACTCCTT
+GTACATCCTTATAGGGATGTAGCNNNNNNNNNNNN
++
+AAAAAEEEEAEEEEEEEEAEEEE############
+@NB500968:8:H5CFGAFXX:1:11101:23302:2031 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:19677:2089 2:N:0:TCCTGAGC+TACTCCTT
+CTCCTCTAAGGGGGGTGCACCAAAGAGTCCAAAACGATCGGGATTTTTGGACTCAGGTCGGGCCACAAAAACGG
++
+AAAAAA/EEEEEEEEEEE/EEEAEEEEEEE/AAAEA/<EEEEE/E</<<E/EEEE/E6EEAEAE/EEE///EEA
+@NB500968:8:H5CFGAFXX:1:11101:26158:2095 2:N:0:TCCTGAGC+TACTCCTT
+GATGTGTATAAGAGACAGGGAAATCCTATAGGGAGTGG
++
+AAAAAEEEEEEEEEEEAEEEAEAEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:7883:2101 2:N:0:TCCGGAGC+TACTCCTT
+GTATAAGAGACAGAGACAGCCTATAGGGTGATGCN
++
+AAA/AEEE/EE/EAAEE//EEA/E/EAAEE/E6<#
+@NB500968:8:H5CFGAFXX:1:11101:19922:2125 2:N:0:TCCTGAGC+TACTCCTT
+ATCCTTGGCAAGATGACCAGACTCACCACAGCGGTAGCAGATGTCAGGGAGAGCCGAGGAAACAAACTGGAACC
++
+/////EE/A//E/EE/E/EE//EA/EE/A/EEEEA/E//E/AE/AEEEEE</A/EEEEAA///////EEE//E/
+@NB500968:8:H5CFGAFXX:1:11101:24496:2131 2:N:0:TCCTGAGC+TACTCCTT
+CTTATACACATCTGACGCTCCCCCCTATAAGAATTTT
++
+AAAAAAEAAE//EEEEEEEE/EEEEEEEAAAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:5814:2162 2:N:0:TCCTGAGC+TACTCCTT
+CTCCTACAGGACAGAAGGACTCTGCTGCCTGTTGCCAAGGAAGGCAGAGCTAGCTGAACTGCGTGTAAAAACGGG
++
+AAA/AEA/EEEEEEEEEEE/EAEEAE<EEA/EAEAA<EEEE/EEEEAEAEEEAEEAEAEEAAAEAEEE/E/<A<<
+@NB500968:8:H5CFGAFXX:1:11101:5796:2184 2:N:0:TCCTGAGC+TATTCCTT
+GTGTATAAGAGACAGAGTCAGTCTATAGGGACAGC
++
+AAAAAEEEEEEEEEEEEEEEEEAEEEEAEAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:23937:2223 2:N:0:TCCTGAGC+TACTCCTT
+GAGATGGTCAGTGCAGAGTGATCGCAGCTACCGGAAAAGAGTTAGCTAGCAATGGGAGCTGAATGAGGCCGTGGT
++
+AAAAAEEEEEEEEEE6EEE6EEEAEEEEE/EEEE/E/E/EEEAE<AAEEA/EEEAAEEEAAE/AA//EEEE/AEA
+@NB500968:8:H5CFGAFXX:1:11101:25267:2246 2:N:0:TCCTGAGC+TACTCCTT
+CAGCCTCTATATATAGGGGCATCCNNNNNNNNNNN
++
+A/AAAAEEEEEEEEEAEEEE6EAE###########
+@NB500968:8:H5CFGAFXX:1:11101:19822:2253 2:N:0:TCCTGAGC+TACTCCTT
+GGCTTCTTACCCATTTAAAGTTTGAGAATAGGTTGAGATCGTTTCGGCCCCAAGACCTCTAATCATTCGCTTTA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEE<
+@NB500968:8:H5CFGAFXX:1:11101:6670:2262 2:N:0:TCCTGAGC+TACTCCTT
+CCTCTATAACCGCATCGGAGACATCGGATTCATTTTAGCTATAGTTTGATTTTCCCTAAACATAAGCTCATGAGA
++
+AAAA/////EEEE///EE/E/<//E<EE/////<AE/AA<A6<EE/EE6<<E////<///</////A/E//AAE/
+@NB500968:8:H5CFGAFXX:1:11101:14935:2297 2:N:0:TCCTGAGC+TACTCCTT
+CGTCAGATGTGTATAAGAGACAGGTGTCTGCTATAGGGCGAAT
++
+AAAA/EA//EEE/<A/E/AEE/AE/EA/EEEEEEE/EEEE///
+@NB500968:8:H5CFGAFXX:1:11101:14858:2304 2:N:0:TCCTGAGC+TACTCCTT
+GCACTGCAGAGGTCTCTAGACAGCACTCTCTTCCTCTCCATGGTCCTCTCTTCCTCCAGCGCTGAGTACTGCTG
++
+AAAAAEEAEEEEEEEEEE6/E/EE/EEEEEEEE/EEEEEEEEEAEEEE/EEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10463:2345 2:N:0:GCCTGAGC+TACTCCTT
+GTGTATAAGAGACAGCGTCTCTGTATAGGGGAAGT
++
+AAAA6E66E/E/A/A/E//EE//A///EEEE///E
+@NB500968:8:H5CFGAFXX:1:11101:24405:2351 2:N:0:TCCTGAGC+TACTCCTT
+GTATAAGAGACAGACCCTATATATAGGGTCNNNNN
++
+AAAAAEEEEEEEEEEEEEEEAEAEEEEEEE#####
+@NB500968:8:H5CFGAFXX:1:11101:6508:2370 2:N:0:TCCTGAGC+TACTCCTT
+CAGCTAGACGGTATAGGGCTGGGCNNNNNNNNNNN
++
+A/AAA/E/EEE//A/EEEAAEEEA###########
+@NB500968:8:H5CFGAFXX:1:11101:13250:2472 2:N:0:TCCTGAGC+TACTCCTT
+CTATGAAATCCCTGTTTGGATTTACAAATAAATTGGACTCCTTGGAGCCTGAAGAGGAAGATGCTGTGCTGAAGG
++
+AAAAAEEAAEEEAAEEEEEEAEEEEEE//EEEAE/<EEEEAEEEEEEEEEAEAEEEEE6AEEEEEEAEEEEAEEA
+@NB500968:8:H5CFGAFXX:1:11101:21730:2488 2:N:0:TCGTGAGC+TACTCCTT
+ATACACATCTGACGCTCCCTATATGTAAAGANNNN
++
+AA/A///E</E/EEAE////////AEA//E/####
+@NB500968:8:H5CFGAFXX:1:11101:5499:2506 2:N:0:TCCTGAGC+TACTCCTT
+GATTTGCCTTATAGGGAGCGTCAGATTTGTATAAGAGACAGCTCCACCCTATACTCGCCGC
++
+A//A<<AA<////A/E/E6<///A//////////////6///<//////////////////
+@NB500968:8:H5CFGAFXX:1:11101:20021:2516 2:N:0:TCCTGGGC+TACTCCTT
+CTAAAAGAGACAGTGCACGACTTAAGGGGATGAGTGTCTCGTATACACGTCTGACGCTCCCGCTATAAAACCACG
++
+//////E/E/E/E////EE/EE///E/EA//////E/E/////E/A///A//E/<<////A/E<///////E//E
+@NB500968:8:H5CFGAFXX:1:11101:4827:2579 2:N:0:TCCTGAGC+TACTCCTT
+ATGTGTATAAGAGACAGTCGAAGTCTATAGGGATTCC
++
+AAAAAE/EA/E/E///E/A/A/EE//EE/E/E/E/AE
+@NB500968:8:H5CFGAFXX:1:11101:25015:2619 2:N:0:TCCTGAGC+TACTCCTT
+GCGTCAGATGTGAGCGTCAGATGTGTATAAGAGACAGGATCCCCCTATAGCAAGAAG
++
+AAAA/AAA<//A/E6A6EEE//A/6//</A6E/AE6/EE<EAEEEA</AEE///AAE
+@NB500968:8:H5CFGAFXX:1:11101:11758:2662 2:N:0:TCCTGAGC+TACTCCTT
+GTGTATAAGAGACAGGTTCGGTCTATAGGGCATCT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEAEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:25905:2738 2:N:0:TCCTGAGC+TACTCCTT
+GTCAGATGTGTATAAGAGACAGTGTGTGAGTATAGGGAAACC
++
+/AAAAEE/AEAA////AEAA/AAEE/EEEAE/E//AE/E/EA
+@NB500968:8:H5CFGAFXX:1:11101:23131:2747 2:N:0:TCCTGAGC+TACTCCTT
+GAATAAGAGACAGATGAGCAGTATAGGGAGCATATTATTGTCTCAAAGATTAAGCCATGCATGTCTAAGTACTAA
++
+A/AAAEEEEAEEAE<EEEEE/E6E///EE/EEE/6//6E/EEAAEEE/EAE<E//6AAAEEE/EE6AEE6EE//A
+@NB500968:8:H5CFGAFXX:1:11101:20373:2751 2:N:0:TCCTGATC+TACTCCTT
+GATGTGTATAAGAGACAGGATCACTTTATAGGCCN
++
+AAAAAEE6AAEE/E<EEEEAE/E6</A/EAEA/E#
+@NB500968:8:H5CFGAFXX:1:11101:8485:2792 2:N:0:TCCTGAGC+TACTCCTT
+GCGGGAACAGGTCCGGGGGCTGTCGGCCAGCGGCCGGGAGAGACTCGGGGGCGGCAGCGGGGGGAGGCGCGGTG
++
+AA/AA////EE///E//E///E///E///<//////EE/E////////A///////E///AE/E//</E//E//
+@NB500968:8:H5CFGAFXX:1:11101:21204:2795 2:N:0:TCCTGAGC+TACTCCTT
+CAGTTAGGGCCTATAGGGTGGAGCNNNNNNNNNNN
++
+AAA<AEEEEEEEEEEEEEEEEEEE###########
+@NB500968:8:H5CFGAFXX:1:11101:10694:2817 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:8407:2853 2:N:0:TCCTGAGC+TACTCCTT
+CACTTACTGGGAATTCCTCGTTCATGGGGAATAATTGCAATCCCCGATCCCCATCACGGATGGGGTTCAACGGGT
++
+AAA6AEEEEEEE/EE/EEEEEEEEEEEEEEA6EEEEEEEE<EEEE</EEEAEEEEAEEEEEEEEAEEEE/EAEEA
+@NB500968:8:H5CFGAFXX:1:11101:17699:2880 2:N:0:TCCTGAGC+TACTCCTT
+CATCTGACGCAGCCGACGAAAAAACATGTGTAGATCTCTGACTTCTCAGTATCATTAAATAAAAAAAATTGATGG
++
+AAAAAEEEEE/EEEEEEE6EEA//<<EEEE//E//EEEAE///E/6///E/EA/EA//////66A//<6A///AA
+@NB500968:8:H5CFGAFXX:1:11101:2261:2941 2:N:0:TCCTGAGC+TACTCCTT
+CTCTTATACACATCTGATGCTGCCCCCTATATGCTCCTT
++
+A/AAAEAEE/E/AEEEEEEEEEEEEEE/E/A<EEEEE//
+@NB500968:8:H5CFGAFXX:1:11101:11672:2944 2:N:0:TCCTGAGC+AACTCCTT
+AAAACCCTAACTCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
++
+A/AA<//6/6//////66666666/6///AEE//EE/</AE</E/<EEEAEEEEEA////<EEE/E<E<A///<6
+@NB500968:8:H5CFGAFXX:1:11101:15968:2945 2:N:0:TCCTGAGC+TACTCCTT
+GATGTGGCGTCGGCATCGGGCGCCTTAACCCGGCGTTCGGTTCATCCCGCAGCGCCAGTTCTGCTTACCAAAAGT
++
+AAAA/EEEEEEEEEAEEEEEEE/EEEEAEEEEEEEEEEEEEAEEE<<EEA<AEE//AAA<EAA/EEEEEEEEEE<
+@NB500968:8:H5CFGAFXX:1:11101:25690:2959 2:N:0:TCCTGCGC+TACTCCTT
+GTGTATAAGAGACAGTCCGAAGCTATAGGGGCGGC
++
+A/AA/EEEEEE/6EEE//E//6EEEAAE/EE/EEE
+@NB500968:8:H5CFGAFXX:1:11101:8090:2972 2:N:0:TCCTGAGC+TACTCCTT
+GGCATAAATTATAGGGGGTCCCNNNNNNNNNNNNN
++
+AAAAAEEEEEEEE6EAEEEEAE#############
+@NB500968:8:H5CFGAFXX:1:11101:23310:3069 2:N:0:TCCTGAGC+TACTCCTT
+GAAAGGAACTGTGACAGTTAGTTGCATGGTCATCACTGGAGACAAACCCTCTGGATGCCGGCATTTCTAAGCTGA
++
+AAAAAEA/AEEEE/EAEEA/EEEEEEEEE<EEA/EEEEEAEEE666E<AEA/AEEEEE/EEEAE/EE<</EEAE/
+@NB500968:8:H5CFGAFXX:1:11101:9254:3110 2:N:0:TCCTGAGC+TACTCCTT
+GTTAGGCTCTTTATAGGGGATCCCNNNNNNNNNNN
++
+AA//AE//A6<//E/EEEE/EA/6###########
+@NB500968:8:H5CFGAFXX:1:11101:22125:3141 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:6707:3162 2:N:0:TCCTGAGC+TACTCCTT
+GCGTAAGAAAAATATAAAACACACTTAAACTCAATAAGTCCATAAATACTATTAAACACAAAAAAAGCATCCAAA
++
+///A/A/E</</EA<///A//<//E/EE//6EA/<//A/</6E//6//<///<//A6<<A/AA//////////</
+@NB500968:8:H5CFGAFXX:1:11101:4099:3187 2:N:0:TCCTGAGC+TACTCCTT
+GATTAGTCAAAAGGCCCATGGTTGCGATGTGGCAGAGAAATCACTTGTGGGGAGACCTAGTAACAGGATTCTGGA
++
+/AAAAE/E/EE66/EEE///EE//E/E///EEEAEEEA///A/EE/6E//A/EEEE///EEEEEEEAAE<<AA<E
+@NB500968:8:H5CFGAFXX:1:11101:24232:3202 2:N:0:TCCTGAGC+TACTCCTT
+TGATTAATGAAAACATTCTTGGCAAATGCTTTCGCTCTGGTCCGTCTTGCGCCGGTCCAAGAATTTCACCTCTAG
++
+AA/AAE/AA/E6/E/EEA<AAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEA6AAEE/E/EEEEA<EAEE
+@NB500968:8:H5CFGAFXX:1:11101:16936:3211 2:N:0:TCCTGAGC+TACTCCTT
+CAGTTCACTTTTATAGGGACGCACNNNNNNNNNNN
++
+AAAAAEEEEEEEEEEEEEEEEEEE###########
+@NB500968:8:H5CFGAFXX:1:11101:6680:3212 2:N:0:TCCTGAGC+TACTCCTT
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:15733:3232 2:N:0:TCCTGAGC+TACTCCTT
+GTGTATAAGAGACAGTCCCCAAATATAGGGGACNN
++
+AAAAAEEEEEEEEEEEEEEEEEEEEE<EEEEEE##
+@NB500968:8:H5CFGAFXX:1:11101:3009:3242 2:N:0:TCCTGAGC+TACTCCTT
+GAAAGAAAGCAAGGACCAGGAGGTTTTCTGGTGGGTTGGTTGGATTTTTTTGAAACAAGAGTCTATGTGGCC
++
+AA/A6E6EE/E/E///EAEEA/E/EEEEAEEEE/EE/AA<E//E/<EEEEEEE//<</AEAEEEA/EEE/E/
+@NB500968:8:H5CFGAFXX:1:11101:6132:3257 2:N:0:TCCTGAGC+TACTCCTT
+GTATAAGAGACAGGAGGCGTGTATAGGGGAATGCN
++
+AAAAAEEEEEEEEEEEEEEEEEEEAEEEE6EEEE#
+@NB500968:8:H5CFGAFXX:1:11101:3840:3285 2:N:0:TCCTGAGC+TACTCCTT
+CAATAGGGTATACCTGTCACTAATACACATCAGACGCTCCCCCCTATATAGTAAAT
++
+A/AAAE//EE/EEAE//E/EE/EEEEEEEEE/EEEAEEEE<A<A/</AE<EEAE<6
+@NB500968:8:H5CFGAFXX:1:11101:2858:3312 2:N:0:TCCTGAGC+TACTCCTT
+CTCTTATACACATCTGACGTTCCCTATATGAGTCAT
++
+/AAAAEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:3641:3398 2:N:0:TCCTGAGC+TACTCCTT
+CTCTTATACACATCTGACGCTGCCCCTATAATGATATC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:20540:3401 2:N:0:TCCTGAGC+TACTCCTT
+CAGTTAGGGCAAAAAGGGTTGATNNNNNNNNNNNN
++
+AAAA//EEAE6/E//AAA//EEE############
+@NB500968:8:H5CFGAFXX:1:11101:18503:3442 2:N:0:TCCTGAGC+TACTCCGT
+GTATAAGAGACAGTAATCCCTTATAGGGATTANNN
++
+AAAAAEAEEEEEAEA<EEAEEEEEEEEEEEE<###
+@NB500968:8:H5CFGAFXX:1:11101:19438:3463 2:N:0:TCCTGAGC+TACTCCTT
+GCGTCGAGATGTGTATAAGAGACAGCCGAGCTCTATAGGGCAAGT
++
+AAAAAEEEEEEEEEEE<AEEEEEEEEEEEEEEEEEEEEEAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:4772:3477 2:N:0:TCCTGAGC+AACTCCTT
+GTTTATAATAGACACATAACAATAAAAGAGTAACAATCACAAAAAAAAATAAGACTCTCACGACGAAACACCTTA
++
+A//A/E//E/A/E/A/E//A<E/////E/E//AE////E/AE</////EEA/E/A/E//AE/E///<<///A<A/
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/read2_extracted.fastq	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,328 @@
+@NB500968:8:H5CFGAFXX:1:11101:15939:1071___1:N:0:TCCTGAGC+TACTCCTT;FP:23586;RQ:34.63
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:1121:1139___1:N:0:TCGTGAGC+TACTCCTT;FP:27828;RQ:34.83
+GATGTAGACTATAGGGAATACCNNNNNNNNNNNNN
++
+A6AAAEEEEEEEEEEEAEE/E<#############
+@NB500968:8:H5CFGAFXX:1:11101:7529:1175___1:N:0:TCCTGAGC+TACTCCTT;FP:17358;RQ:34.81
+GATGTGTAAAAGAGACAGAGGATAATTATAGGGATCCTGTCTCTAAAACACATCTGACGCTCCCTATACATAATA
++
+/AAAAE<A/E//EEE/EE//EEEEAEEEEEEEEEE/A/<EEE<<///<AE</AE//AAA//E/////EEE<6/A/
+@NB500968:8:H5CFGAFXX:1:11101:5052:1198___1:N:0:TCCTGAGC+TACTCCTT;FP:38799;RQ:34.63
+GGATTGAGGCTAGAAACAGCTGCTTTTCTCACTTTCCCAGCAGCAAGACATCAATCTGTTGAGACCTCTCATTTT
++
+AA6AAEEAEEEEEEEEAAEEEEE///AEEA/<E/E<A//EEEEAE/EE//AE///E</EAEEA<EEEAEE/E/E/
+@NB500968:8:H5CFGAFXX:1:11101:15805:1205___1:N:0:TCCTGAGC+TACTCCTT;FP:44542;RQ:34.63
+GCGTCAGATGTGGGTAGTCAGAAAAGCAGATGGATGGGGGCCTTGCCCAACCAGATCCAGGCCGAGTTCACCT
++
+AAAAAEEAEEEEEEEEEEEEEEA//EEEEEEEEEEEEEEEEEEAEE/E//EEEEAEAE/EEEEE<EEEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:11834:1217___1:N:0:TCCTGAGC+TACTCCTT;FP:56795;RQ:34.81
+GGGTCGCTATATAGGGGAATCCNNNNNNNNNNNNN
++
+AAAAAEEEEEEEEEEEEEAEEE#############
+@NB500968:8:H5CFGAFXX:1:11101:20105:1230___1:N:0:TCCTGAGC+TACTCCTT;FP:2265;RQ:34.63
+GTACACGTATTACTTCTGTTTGGTGAATTAAATAAAACAATTCCTAAGCAACTCAAAAGTAACCTAAAAATTGCA
++
+AA/A//EEE/EAEEEEEEE/EEEAEE6EEE<EEA<6E//E/EEAEE//EE/EEE<EAEE6</EEEEE/A/EEE/A
+@NB500968:8:H5CFGAFXX:1:11101:14524:1264___1:N:0:TCCTGAGC+TACTCCTT;FP:7937;RQ:2.65
+GTATAAGAGACAGTTGGGACTTATAGGGTAAGTCN
++
+AAAAAEEEAEEEEAEEE/EEEEEAE<EE/EEEEE#
+@NB500968:8:H5CFGAFXX:1:11101:14981:1265___1:N:0:TCCTGAGC+TACTCCTT;FP:26333;RQ:34.63
+GGACTTAATCAACGCAAGCTTATGACCCGCACTTACTGGGAATTCCTCGTTCATGGGGAATAATTGCAATCCCCG
++
+AA6AAEE/EE//EEE6EEEEEEEE/AAEEAEEEE/EEEEEAEEEE/EEEEEEEEEEEE/AEE<EEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:5538:1286___1:N:0:TCCTGAGC+TACTCCTT;FP:5079;RQ:34.62
+ATCACCAAGCTCGTTAGGCTTTTCACCTCTACCTAAAAATCTTCTCACTATTTTGCCACATAGACGAGTTGATT
++
+AAA/AEE/EEEEEEEAEEEEEEEE/EEEEE/EEEA</EEEEEEEEAAAEAEEEEE<EAE6EEEEEEEEEEEEAE
+@NB500968:8:H5CFGAFXX:1:11101:11473:1333___1:N:0:TCCTGAGC+TACTCCTT;FP:48246;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:2141:1369___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+GATGTGTATAAGAAACAGACCACATATATAGGGGGGAGCGTCAGATGTGTATAAGAGACAGCCCCCCCCAATAGA
++
+AAAAAEEEEEEEAEEEEE/E/EAEEEEEAEEEA/EEEEEEEAA//AEEEE/</AA/////////AAEE////AA/
+@NB500968:8:H5CFGAFXX:1:11101:19474:1422___1:N:0:TCCTGAGC+TACTCCTT;FP:5233;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:12411:1441___1:N:0:TCCTGAGC+TACTCCTT;FP:46039;RQ:34.83
+GACAAAATATAGGGGGACTGTCTGTTATAGACGTCGGACGCGCGCCCCTATAAGGAATAGGGGTCTCTGAGAGA
++
+A/////E////////A///E/E///EEE///A//A/A//<///////E/////AEA//////////////////
+@NB500968:8:H5CFGAFXX:1:11101:5889:1458___1:N:0:TCCTGAGC+TACTCCTT;FP:52002;RQ:34.83
+CAGCGAAGGGGTATAGGGATACGNNNNNNNNNNNN
++
+A/AAA//EEEEA/A/EEE////E############
+@NB500968:8:H5CFGAFXX:1:11101:8662:1469___1:N:0:TCCTGAGC+TACTCCTT;FP:62135;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:9851:1471___1:N:0:TCCTGAGC+TACTCCTT;FP:43356;RQ:3.17
+GTATAAGAGACAGGGGCCCGTTAAAGGGCCCNNNN
++
+/A6<//EA//E/E/A/EEE/E/E/E/6EEE/####
+@NB500968:8:H5CFGAFXX:1:11101:5425:1488___1:N:0:TCCTGAGC+TACTCCTT;FP:13390;RQ:34.83
+GTATAAGAGACAGTGTGTCCATATAGGGGAGTCNN
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEE##
+@NB500968:8:H5CFGAFXX:1:11101:21998:1542___1:N:0:TCCTGAGC+TACTCCTT;FP:15180;RQ:35.09
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:2814:1550___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+CCGTCAGATGAGTATAAGAGACAGGTAAGACATAAAGAGCGTAGTTGTCTCTAATACACATCTGACGCTGCCGA
++
+/AAA6AEEE//E/EAAEE6/EEEEEEEE//EE/E/EE/AE/<A//</AE/E////AE<A/<A/6/E<EEAAEEE
+@NB500968:8:H5CFGAFXX:1:11101:8024:1581___1:N:0:TCCTGAGC+TACTCCTT;FP:38311;RQ:34.63
+GAACTATGACTGTGAGCAGGCATTGTGGTTCTGCCCAAGCTACATCCAGATCACCCAGGCCTCTGGGTCTTTTGG
++
+AAAA6EAEE/E/EEEEEEEAEEEEEEEAEEEEAEEEEEAE6EEE6EEEAA/EEEEEE/A/E/AEAEE<EEEEA/<
+@NB500968:8:H5CFGAFXX:1:11101:8843:1628___1:N:0:TCCTGAGC+TACTCCTT;FP:44455;RQ:34.63
+GGGCAGGGACTTAATCAACGCAAGCTTATGACCCGCACTTACTGGGAATTCCTCGTTCATGGGGAATAATTGCAA
++
+AAAAAEEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEE<EEEEEEEEEEEE<AA/E6EEAA6
+@NB500968:8:H5CFGAFXX:1:11101:23474:1673___1:N:0:TCCTGAGC+TACTCCTT;FP:2968;RQ:34.81
+GATGTGTATAAGAGACAGATATAGACTATAGGGCGCCT
++
+A/AA////A//EA/6<EE6/EEE6////A//EE/E6//
+@NB500968:8:H5CFGAFXX:1:11101:1482:1752___1:N:0:TCCTGAGC+TACTCCTT;FP:20529;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:21720:1786___1:N:0:TCCTGAGC+TACTCCTT;FP:16257;RQ:35.14
+GTCTCTTATACACATCTGACGCTCCCCCCTATAGTTCAAAT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:9012:1809___1:N:0:TCCTGAGC+TACTCCTT;FP:4167;RQ:34.83
+GTGTATAAGAGACAGGCCTATGATATAGGGGGTAC
++
+AAAAAEEEEEEEEEEEEEAEE/EEEAE/AAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12567:1818___1:N:0:TCCTGAGC+TACTCCTT;FP:61471;RQ:34.63
+GGTCATGGGAATAACGCCGCCGCATCGCCAGTCGGCATCGTTTATGGTCGGAACTACGACGGTATCTGATCGT
++
+AA/AAEEAEEE//AEAEEEEEEEEEEAEEEEEEEE/EAEE<EE<EA/<E<EEEEE6EEEE/EEEEEE/EEAEE
+@NB500968:8:H5CFGAFXX:1:11101:20110:1819___1:N:0:TCCTGAGC+TAATCCTT;FP:64559;RQ:34.63
+GAAGAGGGTATCAGATCCCTTGGATCTGGAGTCGCAGATGGCTGCGAGCCACCATGTGGGTTTTGAGAATTGAAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEEEEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:24842:1887___1:N:0:TCCTGAGC+TACTCCTT;FP:13772;RQ:35.14
+CTCATAAACACATCAGGCGGACCCCGCAAAATAAAGGATCTATCTCTGATAAAGAACTGACGCTGCCGACGAAA
++
+A///AE/6EE///////E//6E/AE/<////EE/A/A/A////EAA///////////E/EE/E/////A//E/A
+@NB500968:8:H5CFGAFXX:1:11101:16709:1904___1:N:0:TCCTGAGC+TACTCCTT;FP:5671;RQ:34.73
+CGGATCGCTTCTTCTGATTGAGTTGCCACTCTTTGAAATGCGGCACGGCCTCGTCCACGTTGACGATGCCCAG
++
+AAAAAEEEEEEEAEEEEEEEEEEEEEEEEEEE/EEEEAEEEEEEEEEEEEEEEEE//EEEEEEEAEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:12189:1919___1:N:0:TCCTGAGC+TACTCCTT;FP:53204;RQ:34.62
+GTCTATAGGTATCAGTTTCACTGTAATACCTTTCCTTCCATGGCTAAGTCTGGGAAGGACTTTGCCTTCAAACAG
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE<AEEEEEEEEE<E<EAEAEEEEEEEEEEA
+@NB500968:8:H5CFGAFXX:1:11101:26750:1960___1:N:0:TCCTGAGC+TACTCCTT;FP:578;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:21600:1964___1:N:0:TCCTGAGC+TACTCCTT;FP:38123;RQ:18.05
+GCAATCCCCGATCCCCATCACGAATGGGGTTCAACGGGTTACCCGCGCCTGCCGGCGTAGGGTAGGCACACG
++
+AAAAAEEEEEEAEEEEE/E/EEE6/EEEEEEEE/EEEEEEAEEEEEEAEEEAEEEEEEEEEAEE<EAEAEEE
+@NB500968:8:H5CFGAFXX:1:11101:11625:2006___1:N:0:TCCTGAGC+TACTCCTT;FP:4030;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:10474:2017___1:N:0:TCCTGAGC+TACTCCTT;FP:16506;RQ:35.03
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:17381:2018___1:N:0:TCCTGAGC+TACTCCTT;FP:64754;RQ:35.18
+CTTATACACATCTGACGTTCCCTATACTAATAAAN
++
+A/AAAAE6E////E/AEAE/AEEAEA/EA/EEEE#
+@NB500968:8:H5CFGAFXX:1:11101:9520:2027___1:N:0:TCCTGAGC+TACTCCTT;FP:57159;RQ:34.81
+GTACATCCTTATAGGGATGTAGCNNNNNNNNNNNN
++
+AAAAAEEEEAEEEEEEEEAEEEE############
+@NB500968:8:H5CFGAFXX:1:11101:23302:2031___1:N:0:TCCTGAGC+TACTCCTT;FP:9055;RQ:34.81
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:19677:2089___1:N:0:TCCTGAGC+TACTCCTT;FP:14776;RQ:34.63
+CTCCTCTAAGGGGGGTGCACCAAAGAGTCCAAAACGATCGGGATTTTTGGACTCAGGTCGGGCCACAAAAACGG
++
+AAAAAA/EEEEEEEEEEE/EEEAEEEEEEE/AAAEA/<EEEEE/E</<<E/EEEE/E6EEAEAE/EEE///EEA
+@NB500968:8:H5CFGAFXX:1:11101:26158:2095___1:N:0:TCCTGAGC+TACTCCTT;FP:23868;RQ:34.83
+GATGTGTATAAGAGACAGGGAAATCCTATAGGGAGTGG
++
+AAAAAEEEEEEEEEEEAEEEAEAEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:7883:2101___1:N:0:TCCGGAGC+TACTCCTT;FP:59501;RQ:34.83
+GTATAAGAGACAGAGACAGCCTATAGGGTGATGCN
++
+AAA/AEEE/EE/EAAEE//EEA/E/EAAEE/E6<#
+@NB500968:8:H5CFGAFXX:1:11101:24496:2131___1:N:0:TCCTGAGC+TACTCCTT;FP:247;RQ:35.19
+CTTATACACATCTGACGCTCCCCCCTATAAGAATTTT
++
+AAAAAAEAAE//EEEEEEEE/EEEEEEEAAAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:5814:2162___1:N:0:TCCTGAGC+TACTCCTT;FP:5751;RQ:34.63
+CTCCTACAGGACAGAAGGACTCTGCTGCCTGTTGCCAAGGAAGGCAGAGCTAGCTGAACTGCGTGTAAAAACGGG
++
+AAA/AEA/EEEEEEEEEEE/EAEEAE<EEA/EAEAA<EEEE/EEEEAEAEEEAEEAEAEEAAAEAEEE/E/<A<<
+@NB500968:8:H5CFGAFXX:1:11101:5796:2184___1:N:0:TCCTGAGC+TATTCCTT;FP:4250;RQ:34.83
+GTGTATAAGAGACAGAGTCAGTCTATAGGGACAGC
++
+AAAAAEEEEEEEEEEEEEEEEEAEEEEAEAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:23937:2223___1:N:0:TCCTGAGC+TACTCCTT;FP:26022;RQ:34.63
+GAGATGGTCAGTGCAGAGTGATCGCAGCTACCGGAAAAGAGTTAGCTAGCAATGGGAGCTGAATGAGGCCGTGGT
++
+AAAAAEEEEEEEEEE6EEE6EEEAEEEEE/EEEE/E/E/EEEAE<AAEEA/EEEAAEEEAAE/AA//EEEE/AEA
+@NB500968:8:H5CFGAFXX:1:11101:25267:2246___1:N:0:TCCTGAGC+TACTCCTT;FP:4173;RQ:34.83
+CAGCCTCTATATATAGGGGCATCCNNNNNNNNNNN
++
+A/AAAAEEEEEEEEEAEEEE6EAE###########
+@NB500968:8:H5CFGAFXX:1:11101:19822:2253___1:N:0:TCCTGAGC+TACTCCTT;FP:41486;RQ:34.63
+GGCTTCTTACCCATTTAAAGTTTGAGAATAGGTTGAGATCGTTTCGGCCCCAAGACCTCTAATCATTCGCTTTA
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEAEEEEEEEE<
+@NB500968:8:H5CFGAFXX:1:11101:6670:2262___1:N:0:TCCTGAGC+TACTCCTT;FP:14539;RQ:3.28
+CCTCTATAACCGCATCGGAGACATCGGATTCATTTTAGCTATAGTTTGATTTTCCCTAAACATAAGCTCATGAGA
++
+AAAA/////EEEE///EE/E/<//E<EE/////<AE/AA<A6<EE/EE6<<E////<///</////A/E//AAE/
+@NB500968:8:H5CFGAFXX:1:11101:14858:2304___1:N:0:TCCTGAGC+TACTCCTT;FP:4202;RQ:34.63
+GCACTGCAGAGGTCTCTAGACAGCACTCTCTTCCTCTCCATGGTCCTCTCTTCCTCCAGCGCTGAGTACTGCTG
++
+AAAAAEEAEEEEEEEEEE6/E/EE/EEEEEEEE/EEEEEEEEEAEEEE/EEEEEEEEEEEEEEEEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:10463:2345___1:N:0:GCCTGAGC+TACTCCTT;FP:24337;RQ:34.81
+GTGTATAAGAGACAGCGTCTCTGTATAGGGGAAGT
++
+AAAA6E66E/E/A/A/E//EE//A///EEEE///E
+@NB500968:8:H5CFGAFXX:1:11101:24405:2351___1:N:0:TCCTGAGC+TACTCCTT;FP:5580;RQ:35.18
+GTATAAGAGACAGACCCTATATATAGGGTCNNNNN
++
+AAAAAEEEEEEEEEEEEEEEAEAEEEEEEE#####
+@NB500968:8:H5CFGAFXX:1:11101:6508:2370___1:N:0:TCCTGAGC+TACTCCTT;FP:64422;RQ:34.83
+CAGCTAGACGGTATAGGGCTGGGCNNNNNNNNNNN
++
+A/AAA/E/EEE//A/EEEAAEEEA###########
+@NB500968:8:H5CFGAFXX:1:11101:21730:2488___1:N:0:TCGTGAGC+TACTCCTT;FP:57284;RQ:35.18
+ATACACATCTGACGCTCCCTATATGTAAAGANNNN
++
+AA/A///E</E/EEAE////////AEA//E/####
+@NB500968:8:H5CFGAFXX:1:11101:5499:2506___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+GATTTGCCTTATAGGGAGCGTCAGATTTGTATAAGAGACAGCTCCACCCTATACTCGCCGC
++
+A//A<<AA<////A/E/E6<///A//////////////6///<//////////////////
+@NB500968:8:H5CFGAFXX:1:11101:20021:2516___1:N:0:TCCTGGGC+TACTCCTT;FP:28351;RQ:34.83
+CTAAAAGAGACAGTGCACGACTTAAGGGGATGAGTGTCTCGTATACACGTCTGACGCTCCCGCTATAAAACCACG
++
+//////E/E/E/E////EE/EE///E/EA//////E/E/////E/A///A//E/<<////A/E<///////E//E
+@NB500968:8:H5CFGAFXX:1:11101:4827:2579___1:N:0:TCCTGAGC+TACTCCTT;FP:8626;RQ:34.83
+ATGTGTATAAGAGACAGTCGAAGTCTATAGGGATTCC
++
+AAAAAE/EA/E/E///E/A/A/EE//EE/E/E/E/AE
+@NB500968:8:H5CFGAFXX:1:11101:25015:2619___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+GCGTCAGATGTGAGCGTCAGATGTGTATAAGAGACAGGATCCCCCTATAGCAAGAAG
++
+AAAA/AAA<//A/E6A6EEE//A/6//</A6E/AE6/EE<EAEEEA</AEE///AAE
+@NB500968:8:H5CFGAFXX:1:11101:11758:2662___1:N:0:TCCTGAGC+TACTCCTT;FP:32308;RQ:34.83
+GTGTATAAGAGACAGGTTCGGTCTATAGGGCATCT
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEAEAEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:25905:2738___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+GTCAGATGTGTATAAGAGACAGTGTGTGAGTATAGGGAAACC
++
+/AAAAEE/AEAA////AEAA/AAEE/EEEAE/E//AE/E/EA
+@NB500968:8:H5CFGAFXX:1:11101:23131:2747___1:N:0:TCCTGAGC+TACTCCTT;FP:61606;RQ:34.61
+GAATAAGAGACAGATGAGCAGTATAGGGAGCATATTATTGTCTCAAAGATTAAGCCATGCATGTCTAAGTACTAA
++
+A/AAAEEEEAEEAE<EEEEE/E6E///EE/EEE/6//6E/EEAAEEE/EAE<E//6AAAEEE/EE6AEE6EE//A
+@NB500968:8:H5CFGAFXX:1:11101:20373:2751___1:N:0:TCCTGATC+TACTCCTT;FP:283;RQ:34.93
+GATGTGTATAAGAGACAGGATCACTTTATAGGCCN
++
+AAAAAEE6AAEE/E<EEEEAE/E6</A/EAEA/E#
+@NB500968:8:H5CFGAFXX:1:11101:21204:2795___1:N:0:TCCTGAGC+TACTCCTT;FP:16279;RQ:34.81
+CAGTTAGGGCCTATAGGGTGGAGCNNNNNNNNNNN
++
+AAA<AEEEEEEEEEEEEEEEEEEE###########
+@NB500968:8:H5CFGAFXX:1:11101:10694:2817___1:N:0:TCCTGAGC+TACTCCTT;FP:10036;RQ:34.83
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:8407:2853___1:N:0:TCCTGAGC+TACTCCTT;FP:21171;RQ:34.62
+CACTTACTGGGAATTCCTCGTTCATGGGGAATAATTGCAATCCCCGATCCCCATCACGGATGGGGTTCAACGGGT
++
+AAA6AEEEEEEE/EE/EEEEEEEEEEEEEEA6EEEEEEEE<EEEE</EEEAEEEEAEEEEEEEEAEEEE/EAEEA
+@NB500968:8:H5CFGAFXX:1:11101:17699:2880___1:N:0:TCCTGAGC+TACTCCTT;FP:52710;RQ:34.62
+CATCTGACGCAGCCGACGAAAAAACATGTGTAGATCTCTGACTTCTCAGTATCATTAAATAAAAAAAATTGATGG
++
+AAAAAEEEEE/EEEEEEE6EEA//<<EEEE//E//EEEAE///E/6///E/EA/EA//////66A//<6A///AA
+@NB500968:8:H5CFGAFXX:1:11101:2261:2941___1:N:0:TCCTGAGC+TACTCCTT;FP:2596;RQ:34.63
+CTCTTATACACATCTGATGCTGCCCCCTATATGCTCCTT
++
+A/AAAEAEE/E/AEEEEEEEEEEEEEE/E/A<EEEEE//
+@NB500968:8:H5CFGAFXX:1:11101:11672:2944___1:N:0:TCCTGAGC+AACTCCTT;FP:32501;RQ:34.63
+AAAACCCTAACTCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
++
+A/AA<//6/6//////66666666/6///AEE//EE/</AE</E/<EEEAEEEEEA////<EEE/E<E<A///<6
+@NB500968:8:H5CFGAFXX:1:11101:15968:2945___1:N:0:TCCTGAGC+TACTCCTT;FP:8706;RQ:34.63
+GATGTGGCGTCGGCATCGGGCGCCTTAACCCGGCGTTCGGTTCATCCCGCAGCGCCAGTTCTGCTTACCAAAAGT
++
+AAAA/EEEEEEEEEAEEEEEEE/EEEEAEEEEEEEEEEEEEAEEE<<EEA<AEE//AAA<EAA/EEEEEEEEEE<
+@NB500968:8:H5CFGAFXX:1:11101:25690:2959___1:N:0:TCCTGCGC+TACTCCTT;FP:27004;RQ:34.81
+GTGTATAAGAGACAGTCCGAAGCTATAGGGGCGGC
++
+A/AA/EEEEEE/6EEE//E//6EEEAAE/EE/EEE
+@NB500968:8:H5CFGAFXX:1:11101:8090:2972___1:N:0:TCCTGAGC+TACTCCTT;FP:12902;RQ:34.81
+GGCATAAATTATAGGGGGTCCCNNNNNNNNNNNNN
++
+AAAAAEEEEEEEE6EAEEEEAE#############
+@NB500968:8:H5CFGAFXX:1:11101:22125:3141___1:N:0:TCCTGAGC+TACTCCTT;FP:65197;RQ:35.07
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:6707:3162___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+GCGTAAGAAAAATATAAAACACACTTAAACTCAATAAGTCCATAAATACTATTAAACACAAAAAAAGCATCCAAA
++
+///A/A/E</</EA<///A//<//E/EE//6EA/<//A/</6E//6//<///<//A6<<A/AA//////////</
+@NB500968:8:H5CFGAFXX:1:11101:24232:3202___1:N:0:TCCTGAGC+TACTCCTT;FP:19483;RQ:34.63
+TGATTAATGAAAACATTCTTGGCAAATGCTTTCGCTCTGGTCCGTCTTGCGCCGGTCCAAGAATTTCACCTCTAG
++
+AA/AAE/AA/E6/E/EEA<AAEEAEEEEEEEEEEEEEEEEEEEEEEEEEEEEAEEEA6AAEE/E/EEEEA<EAEE
+@NB500968:8:H5CFGAFXX:1:11101:16936:3211___1:N:0:TCCTGAGC+TACTCCTT;FP:7401;RQ:34.83
+CAGTTCACTTTTATAGGGACGCACNNNNNNNNNNN
++
+AAAAAEEEEEEEEEEEEEEEEEEE###########
+@NB500968:8:H5CFGAFXX:1:11101:6680:3212___1:N:0:TCCTGAGC+TACTCCTT;FP:62506;RQ:34.81
+NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
++
+###################################
+@NB500968:8:H5CFGAFXX:1:11101:6132:3257___1:N:0:TCCTGAGC+TACTCCTT;FP:7557;RQ:34.62
+GTATAAGAGACAGGAGGCGTGTATAGGGGAATGCN
++
+AAAAAEEEEEEEEEEEEEEEEEEEAEEEE6EEEE#
+@NB500968:8:H5CFGAFXX:1:11101:3840:3285___1:N:0:TCCTGAGC+TACTCCTT;FP:16156;RQ:34.93
+CAATAGGGTATACCTGTCACTAATACACATCAGACGCTCCCCCCTATATAGTAAAT
++
+A/AAAE//EE/EEAE//E/EE/EEEEEEEEE/EEEAEEEE<A<A/</AE<EEAE<6
+@NB500968:8:H5CFGAFXX:1:11101:2858:3312___1:N:0:TCCTGAGC+TACTCCTT;FP:14452;RQ:35.12
+CTCTTATACACATCTGACGTTCCCTATATGAGTCAT
++
+/AAAAEEAEEEEEEEEEEEEEEEEEAEEEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:3641:3398___1:N:0:TCCTGAGC+TACTCCTT;FP:36051;RQ:35.14
+CTCTTATACACATCTGACGCTGCCCCTATAATGATATC
++
+AAAAAEEEEEEEEEEEEEEEEEEEEEEEEAEEEEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:20540:3401___1:N:0:TCCTGAGC+TACTCCTT;FP:48603;RQ:34.81
+CAGTTAGGGCAAAAAGGGTTGATNNNNNNNNNNNN
++
+AAAA//EEAE6/E//AAA//EEE############
+@NB500968:8:H5CFGAFXX:1:11101:19438:3463___1:N:0:TCCTGAGC+TACTCCTT;FP:0;RQ:1.56
+GCGTCGAGATGTGTATAAGAGACAGCCGAGCTCTATAGGGCAAGT
++
+AAAAAEEEEEEEEEEE<AEEEEEEEEEEEEEEEEEEEEEAEEEEE
+@NB500968:8:H5CFGAFXX:1:11101:4772:3477___1:N:0:TCCTGAGC+AACTCCTT;FP:61788;RQ:3.33
+GTTTATAATAGACACATAACAATAAAAGAGTAACAATCACAAAAAAAAATAAGACTCTCACGACGAAACACCTTA
++
+A//A/E//E/A/E/A/E//A<E/////E/E//AE////E/AE</////EEA/E/A/E//AE/E///<<///A<A/
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml	Mon May 09 05:15:01 2016 -0400
@@ -0,0 +1,6 @@
+<?xml version="1.0"?>
+<tool_dependency>
+    <package name="biopython" version="1.66">
+        <repository changeset_revision="8433ee4531ff" name="package_biopython_1_66" owner="biopython" prior_installation_required="True" toolshed="https://toolshed.g2.bx.psu.edu" />
+    </package>
+</tool_dependency>