added illumina 1.8 name scheme support to order pairs

git-svn-id: http://biopieces.googlecode.com/svn/trunk@1836 74ccb610-7750-0410-82ae-013aeee3265d

diff --git a/bp_bin/order_pairs b/bp_bin/order_pairs
index c96f845f63b14b7a62e5c699fff895562949e626..e8f799052470feb9dd6bfb3eac9a8cd706ea55d5 100755 (executable)
--- a/bp_bin/order_pairs
+++ b/bp_bin/order_pairs
@@ -1,6 +1,6 @@
 #!/usr/bin/env ruby
 
-# Copyright (C) 2007-2011 Martin A. Hansen.
+# Copyright (C) 2007-2012 Martin A. Hansen.
 
 # This program is free software; you can redistribute it and/or
 # modify it under the terms of the GNU General Public License
@@ -24,11 +24,10 @@
 
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> DESCRIPTION <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
-# Order records so SEQ_NAME/1 and SEQ_NAME/2 are in order.
+# Order records with pair end sequence data.
 
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
-
 require 'maasha/biopieces'
 require 'pp'
 
@@ -40,11 +39,19 @@ hash2 = {}
 Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
   input.each_record do |record|
     if record[:SEQ_NAME]
-      name = record[:SEQ_NAME][0 ... -2].to_sym
-      pair = record[:SEQ_NAME][-1].to_i
+      case record[:SEQ_NAME]
+      when /^(.+)\/(\d)$/   # Illumina 1.5
+        name = $1.to_sym
+        pair = $2.to_i
+      when /^(.+) (\d):/    # Illumina 1.8
+        name = $1.to_sym
+        pair = $2.to_i
+      else
+        $stderr.puts "WARNING: Unmatched sequence name: #{record[:SEQ_NAME]}" if options[:verbose]
+      end
 
       if pair == 1
-        if hash2.has_key? name
+        if hash2[name]
           output.puts record
           output.puts hash2[name]
 
@@ -53,7 +60,7 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
           hash1[name] = record
         end
       elsif pair == 2
-        if hash1.has_key? name
+        if hash1[name]
           output.puts hash1[name]
           output.puts record
 
@@ -61,19 +68,16 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
         else
           hash2[name] = record
         end
-      else
-        $stderr.puts "WARNING: Badly formatted name: #{record[:SEQ_NAME]}" if options[:verbose]
       end
     end
   end
 end
 
 if options[:verbose]
-  $stderr.puts "Unmatched names with /1: #{hash1.size}"
-  $stderr.puts "Unmatched names with /2: #{hash2.size}"
+  $stderr.puts "Unmatched names with 1st read pair: #{hash1.size}"
+  $stderr.puts "Unmatched names with 2nd read pair: #{hash2.size}"
 end
 
-
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
 

diff --git a/bp_test/in/order_pairs.in b/bp_test/in/order_pairs.in
index 3e30de4db9fc2302797a745f25f1e9d10755b8f8..f9a2a66b6fb8dcd793074ed44debd2ff927d1e6b 100644 (file)
--- a/bp_test/in/order_pairs.in
+++ b/bp_test/in/order_pairs.in
@@ -38,3 +38,103 @@ SEQ: AANAANANCTNCNNANCGNCGGTTTCGTNANNANNGGTTCGCANCGNNNC
 SEQ_LEN: 50
 SEQ_NAME: ILLUMINA-52179E_0004:2:1:1045:16499#TTAGGC/2
 ---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1958:2175 2:N:0:TAGCTG
+SEQ: CTACCTGTGCAAAAAGCACCAAGCAGAGCCATGCTNAACATGGACAGGCGAGCTGGTGTCATAAACATTTCCTGNTTTTANTCTATNGTAGNTANTAA
+SEQ_LEN: 98
+SCORES: bbbeeeecggfeghfggfhihhieh]ffbeghiiiBOYcfghighfghghaedfeeV^abcccce`bdbdc```BKZ^aaBKW^^aBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2412:2210 2:N:0:TAGCTG
+SEQ: CCGGGCTTGATGCCCTGGCATTGGAGCAATACGAGTTGTTGCTGGAAGAACAAGCCTATCCGTTTAGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT
+SEQ_LEN: 98
+SCORES: ^^_eee^ccgeggiidfdd_ggfhhhfhifhhfgfggffgdeghe_ebfgedfXZcffdggdT`^aeeee^`bcacc\^Z^a[]Q[`b^]^aaBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2110:2236 2:N:0:TAGCTG
+SEQ: CCATCACGCAATCTGCTAACACCGACAGATAACGCTCGCGCCAACGCAGATCTGTATTGCGGGCATCCTGCATATTACTGTTAGGAATAGATCGGAAG
+SEQ_LEN: 98
+SCORES: ZZac^cd`gf`eeifhh`dfbaafagha`g`_`e`a^cZ`edS\deedcc`^_]]bccRTZQW[T^bbb__bbbb__bc]bc_R^bc]_]bbba[[^R
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2150 1:N:0:TAGCTG
+SEQ: CATAGGCTTTTTTAACGGCCCAGCCAAAGGCGATCAGGCCAAATAGNANTACCCACCAGTACTGTTGCATNGTCTCAGAAATGNNCAATACAAATTGG
+SEQ_LEN: 98
+SCORES: _[_cccccgcaeeaeJb`b^dd_cg]cZedbfU_eghhhddfg[f]BLBLTZ^abbb^aT]b_`]`b`Y`BKT^^ac_YYY_`BBJGQW_b]bbbbbb
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1832:2218 1:N:0:TAGCTG
+SEQ: GCTTTGACATAGTCGCTCCAGAATTGCCAGCTAGGGTTAGCTTGGCAACTGCAGCGACGTAATGTGCTGTGGCAGATCAATTTATCTGTTTTGAATCA
+SEQ_LEN: 98
+SCORES: ^P^PJ\Y`eea`e[daYdecggadgdXJIYVbdc`efg_cdedI^aXIO^abeb\eL_daQU^_V]``]UGTZ\^BBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2155:2169 2:N:0:TAGCTG
+SEQ: AACAAACTGTTTTGCAGAATGTTATGAATGGTCCGNGCGCNTNAACCCAGATCGGAAGAGCGTCGTGTAGGGAANGANGGNNGNNCNCGGNNNNNGNC
+SEQ_LEN: 98
+SCORES: a__eeecegcgegihhifhhfcgdghfffiibgheBPPYcBOBOXaffffgfhhhffa`bfe\^cR]Z]`baZ`BBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2037:2173 1:N:0:TAGCTG
+SEQ: CTGGGTGCCACTGTCGCCACCCTGTTAAGTCGTGAACCCGTTACCATCAGCACTGATGCCAACGTCTGGCAGGCCGCCCAGCTAATGACCGAGCAAAA
+SEQ_LEN: 98
+SCORES: _[_cccceeegggffhiiihdfdefhhafeff_efdfgfhifhghhfbgbgdggfdhihhegggeeeebdZ^aaa^_aaX``b]b_bcccXXaV_ccb
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2284:2226 1:N:0:TAGCTG
+SEQ: CGCCTTTTCATCAGCAACAGACAGTTACGACGCATATCGCCACCCAAATGTTTGCCGAGGCTAAACGGGTATTGAAACGCGGTGGCCGCTTGCGTATT
+SEQ_LEN: 98
+SCORES: ___eeeeegagffgghdefb_e_eggia`\egfdhhdf`gdfghihiiib_`gghf``cbcabcbc^^_KW`c]Y_bbcc[a[[_^aa]cccc]__ab
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2412:2210 1:N:0:TAGCTG
+SEQ: AAACGGATAGGCTTGTTCTTCCAGCAACAACTCGTATTGCTCCAATGCCAGGGCATCAAGCCCGGAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
+SEQ_LEN: 98
+SCORES: _bbaceceggeccgeZddghaegffiiWdghfff``ffgfacfcffhiihhfi\_e\egfcggde]acdccc^aa_`[`aca_a^b_`bbbbbbbS_b
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2284:2226 2:N:0:TAGCTG
+SEQ: GGCCAAATAACCGCTGCAGTTGCTGCTGATAATTCAGGTGACGGTTAGCCACAATACGCAAGCGGCCACCGCGTTTCAATACCCGTTTAGCCTCGGCA
+SEQ_LEN: 98
+SCORES: bbbeeeeegfgggiehihhef_gfgfhfghdffhihhfXc^egc_cgabgffghhfhfhfg`d^cba_ac^\_W]abcb_bccc]X`bQ]`b`^_a__
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2037:2173 2:N:0:TAGCTG
+SEQ: AGCCTGGCCCGGATATCCTTGTCAGTGATAATGCCNGTCAGATGCTGCTTGGCATCTAAAATAAGCAAGGATGANACCTTNTGNTCNGTCNNTANCTG
+SEQ_LEN: 98
+SCORES: b__eeeeeggggcgghiiiihefggcgeffagdbfBOY^cffhfgbgfhdggXegghagghfffeggggeXV^_BKZ\_cBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1958:2175 1:N:0:TAGCTG
+SEQ: TACCTAAAATAATTAATATCTACAATAGAATAAAATCAGGAAATGTTTATGACACCAGCTCGCCTGTCCATGTTCAGCATGGCTCTGCTTGGTGCTTT
+SEQ_LEN: 98
+SCORES: __acceeegcgggiihihifiiiiiiiiiiefihhhhchhdhiafcghidffghhiihidfggighedhdhibdggegdgeeeedddcddccbcccab
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2155:2169 1:N:0:TAGCTG
+SEQ: GGGTTCAGGCGCCCGGACCATTCATAACATTCTGCAAAACAGTTTGTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTAGCTGATCTCGTATG
+SEQ_LEN: 98
+SCORES: ___J`cccgZccggghhhhffhbghhfeegffdghhfchhhfaffgbb]e`_c\a_\degdddbbb``^aY]bbb_bbbbbbb`bbbb_]`R]^[WR[
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2175 1:N:0:TAGCTG
+SEQ: GTGCCAAAGAGGTACTGATGCTTTGTTGTTGTGGTTTAGCCGCACATTCAGCCTCGCTGACGGCAGCTAACTGCACATCTTTGGGGATATCCACCAGA
+SEQ_LEN: 98
+SCORES: ^PWYP`R^cRcccgdghS`d_ebc[`ebffddb[ffffagf[O^_Uc`SabeffeegffW\\[XXR^bTGT__]b]_Y_]Zbbba][abb_cbba_aa
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2175 2:N:0:TAGCTG
+SEQ: GGTGGATATCCCCAAAGATGTGCAGTTAGCTGCCGNCAGCGAGGCTGAATGTGCGGCTAAACCACAACAACAAANCATCANTACCTNTTTGNCANAGA
+SEQ_LEN: 98
+SCORES: ___cc\c`eggcghad`acabgfdg`g^daagaefBPYOafY_fYaa]e\aMWW__[W^Z_RZ^`aa__`Q[`^BBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2139:2216 2:N:0:TAGCTG
+SEQ: ATGTGTCGATGCTGGATAACGAATGGGGTATGTCGGTGTATCCGCTGGGGGCTGGTCATGAAGGGGTGGGGCATATTGCAGCTGTGGGTTCGGGGGCC
+SEQ_LEN: 98
+SCORES: JJ\caa`^Ye`gW^[J[Q[`YJJ`dfa^PYbfX^^_H^^eh]fb`aefGZ_BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2139:2216 1:N:0:TAGCTG
+SEQ: CCACCTGCCCGGTTTTAAGATGACTGACCGCCGAACCCACAGCTGCAATATGCCCCACCACTTCATGACCAGCCACCAGCGGATACGCCGACATACCC
+SEQ_LEN: 98
+SCORES: _Z_cccccaY_^cghhbagg[e]de][^ZcghfY_efgfbfdQ^^\bdZaaRMZZT\FRT\LLYYZTGT]_[^a^[^aBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2110:2236 1:N:0:TAGCTG
+SEQ: ATTCCTAACAGTAATATGCAGGATGCCCGCAATACAGATCTGCGTTGGCGCGAGCGTTATCTGTCGGTGTTAGCAGATTTCGTGATGGAGATCGGAAG
+SEQ_LEN: 98
+SCORES: ___cc`Sc^J[ecK``f]cY^RPbcRYJQQY^YccacfehhheU^^^eXFS_BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2150 2:N:0:TAGCTG
+SEQ: CGGCCNNNACCCAATTTGNNTTGGGCATTTCTGANNNNNNNCNACNNTACTGNNGGGTANTGCTNNNTGNNNNGNNNNNNNNNNNNNNNNNNNNNNNN
+SEQ_LEN: 98
+SCORES: WZ^ceBBBQQQQbffgggBBSQ`_fRZdegffdhBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1832:2218 2:N:0:TAGCTG
+SEQ: GGTTATCGATCTGGAAAAAGCAACTAAACCTAAAGCTAAACCACGTAGCGCCGGGTAAATGATTCAAAACAGATAAATTGATCTGCCACAGCACATTA
+SEQ_LEN: 98
+SCORES: ^VYPJQ`c^JJ[b[efg^dHJ`aa`adXd_ZXXbIIIY[af_H^aWHWPZ[`gggFFZ^bd_Z]Zb_]ba\^ZGY_`TZ``cc[bbR]]^aaXQ[bbb
+---

diff --git a/bp_test/out/order_pairs.out.1 b/bp_test/out/order_pairs.out.1
index b0906df71803d7d658f8f765fb9ee1f76bc5ceb8..e25e55d02e1d10db14690bde01e2b33669599089 100644 (file)
--- a/bp_test/out/order_pairs.out.1
+++ b/bp_test/out/order_pairs.out.1
@@ -28,3 +28,103 @@ SEQ: AANAANANCTNCNNANCGNCGGTTTCGTNANNANNGGTTCGCANCGNNNC
 SEQ_LEN: 50
 SEQ_NAME: ILLUMINA-52179E_0004:2:1:1045:16499#TTAGGC/2
 ---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2412:2210 1:N:0:TAGCTG
+SEQ: AAACGGATAGGCTTGTTCTTCCAGCAACAACTCGTATTGCTCCAATGCCAGGGCATCAAGCCCGGAGATCGGAAGAGCACACGTCTGAACTCCAGTCA
+SEQ_LEN: 98
+SCORES: _bbaceceggeccgeZddghaegffiiWdghfff``ffgfacfcffhiihhfi\_e\egfcggde]acdccc^aa_`[`aca_a^b_`bbbbbbbS_b
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2412:2210 2:N:0:TAGCTG
+SEQ: CCGGGCTTGATGCCCTGGCATTGGAGCAATACGAGTTGTTGCTGGAAGAACAAGCCTATCCGTTTAGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT
+SEQ_LEN: 98
+SCORES: ^^_eee^ccgeggiidfdd_ggfhhhfhifhhfgfggffgdeghe_ebfgedfXZcffdggdT`^aeeee^`bcacc\^Z^a[]Q[`b^]^aaBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2284:2226 1:N:0:TAGCTG
+SEQ: CGCCTTTTCATCAGCAACAGACAGTTACGACGCATATCGCCACCCAAATGTTTGCCGAGGCTAAACGGGTATTGAAACGCGGTGGCCGCTTGCGTATT
+SEQ_LEN: 98
+SCORES: ___eeeeegagffgghdefb_e_eggia`\egfdhhdf`gdfghihiiib_`gghf``cbcabcbc^^_KW`c]Y_bbcc[a[[_^aa]cccc]__ab
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2284:2226 2:N:0:TAGCTG
+SEQ: GGCCAAATAACCGCTGCAGTTGCTGCTGATAATTCAGGTGACGGTTAGCCACAATACGCAAGCGGCCACCGCGTTTCAATACCCGTTTAGCCTCGGCA
+SEQ_LEN: 98
+SCORES: bbbeeeeegfgggiehihhef_gfgfhfghdffhihhfXc^egc_cgabgffghhfhfhfg`d^cba_ac^\_W]abcb_bccc]X`bQ]`b`^_a__
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2037:2173 1:N:0:TAGCTG
+SEQ: CTGGGTGCCACTGTCGCCACCCTGTTAAGTCGTGAACCCGTTACCATCAGCACTGATGCCAACGTCTGGCAGGCCGCCCAGCTAATGACCGAGCAAAA
+SEQ_LEN: 98
+SCORES: _[_cccceeegggffhiiihdfdefhhafeff_efdfgfhifhghhfbgbgdggfdhihhegggeeeebdZ^aaa^_aaX``b]b_bcccXXaV_ccb
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2037:2173 2:N:0:TAGCTG
+SEQ: AGCCTGGCCCGGATATCCTTGTCAGTGATAATGCCNGTCAGATGCTGCTTGGCATCTAAAATAAGCAAGGATGANACCTTNTGNTCNGTCNNTANCTG
+SEQ_LEN: 98
+SCORES: b__eeeeeggggcgghiiiihefggcgeffagdbfBOY^cffhfgbgfhdggXegghagghfffeggggeXV^_BKZ\_cBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1958:2175 1:N:0:TAGCTG
+SEQ: TACCTAAAATAATTAATATCTACAATAGAATAAAATCAGGAAATGTTTATGACACCAGCTCGCCTGTCCATGTTCAGCATGGCTCTGCTTGGTGCTTT
+SEQ_LEN: 98
+SCORES: __acceeegcgggiihihifiiiiiiiiiiefihhhhchhdhiafcghidffghhiihidfggighedhdhibdggegdgeeeedddcddccbcccab
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1958:2175 2:N:0:TAGCTG
+SEQ: CTACCTGTGCAAAAAGCACCAAGCAGAGCCATGCTNAACATGGACAGGCGAGCTGGTGTCATAAACATTTCCTGNTTTTANTCTATNGTAGNTANTAA
+SEQ_LEN: 98
+SCORES: bbbeeeecggfeghfggfhihhieh]ffbeghiiiBOYcfghighfghghaedfeeV^abcccce`bdbdc```BKZ^aaBKW^^aBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2155:2169 1:N:0:TAGCTG
+SEQ: GGGTTCAGGCGCCCGGACCATTCATAACATTCTGCAAAACAGTTTGTTAGATCGGAAGAGCACACGTCTGAACTCCAGTCACTAGCTGATCTCGTATG
+SEQ_LEN: 98
+SCORES: ___J`cccgZccggghhhhffhbghhfeegffdghhfchhhfaffgbb]e`_c\a_\degdddbbb``^aY]bbb_bbbbbbb`bbbb_]`R]^[WR[
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2155:2169 2:N:0:TAGCTG
+SEQ: AACAAACTGTTTTGCAGAATGTTATGAATGGTCCGNGCGCNTNAACCCAGATCGGAAGAGCGTCGTGTAGGGAANGANGGNNGNNCNCGGNNNNNGNC
+SEQ_LEN: 98
+SCORES: a__eeecegcgegihhifhhfcgdghfffiibgheBPPYcBOBOXaffffgfhhhffa`bfe\^cR]Z]`baZ`BBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2175 1:N:0:TAGCTG
+SEQ: GTGCCAAAGAGGTACTGATGCTTTGTTGTTGTGGTTTAGCCGCACATTCAGCCTCGCTGACGGCAGCTAACTGCACATCTTTGGGGATATCCACCAGA
+SEQ_LEN: 98
+SCORES: ^PWYP`R^cRcccgdghS`d_ebc[`ebffddb[ffffagf[O^_Uc`SabeffeegffW\\[XXR^bTGT__]b]_Y_]Zbbba][abb_cbba_aa
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2175 2:N:0:TAGCTG
+SEQ: GGTGGATATCCCCAAAGATGTGCAGTTAGCTGCCGNCAGCGAGGCTGAATGTGCGGCTAAACCACAACAACAAANCATCANTACCTNTTTGNCANAGA
+SEQ_LEN: 98
+SCORES: ___cc\c`eggcghad`acabgfdg`g^daagaefBPYOafY_fYaa]e\aMWW__[W^Z_RZ^`aa__`Q[`^BBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2139:2216 1:N:0:TAGCTG
+SEQ: CCACCTGCCCGGTTTTAAGATGACTGACCGCCGAACCCACAGCTGCAATATGCCCCACCACTTCATGACCAGCCACCAGCGGATACGCCGACATACCC
+SEQ_LEN: 98
+SCORES: _Z_cccccaY_^cghhbagg[e]de][^ZcghfY_efgfbfdQ^^\bdZaaRMZZT\FRT\LLYYZTGT]_[^a^[^aBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2139:2216 2:N:0:TAGCTG
+SEQ: ATGTGTCGATGCTGGATAACGAATGGGGTATGTCGGTGTATCCGCTGGGGGCTGGTCATGAAGGGGTGGGGCATATTGCAGCTGTGGGTTCGGGGGCC
+SEQ_LEN: 98
+SCORES: JJ\caa`^Ye`gW^[J[Q[`YJJ`dfa^PYbfX^^_H^^eh]fb`aefGZ_BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2110:2236 1:N:0:TAGCTG
+SEQ: ATTCCTAACAGTAATATGCAGGATGCCCGCAATACAGATCTGCGTTGGCGCGAGCGTTATCTGTCGGTGTTAGCAGATTTCGTGATGGAGATCGGAAG
+SEQ_LEN: 98
+SCORES: ___cc`Sc^J[ecK``f]cY^RPbcRYJQQY^YccacfehhheU^^^eXFS_BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2110:2236 2:N:0:TAGCTG
+SEQ: CCATCACGCAATCTGCTAACACCGACAGATAACGCTCGCGCCAACGCAGATCTGTATTGCGGGCATCCTGCATATTACTGTTAGGAATAGATCGGAAG
+SEQ_LEN: 98
+SCORES: ZZac^cd`gf`eeifhh`dfbaafagha`g`_`e`a^cZ`edS\deedcc`^_]]bccRTZQW[T^bbb__bbbb__bc]bc_R^bc]_]bbba[[^R
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2150 1:N:0:TAGCTG
+SEQ: CATAGGCTTTTTTAACGGCCCAGCCAAAGGCGATCAGGCCAAATAGNANTACCCACCAGTACTGTTGCATNGTCTCAGAAATGNNCAATACAAATTGG
+SEQ_LEN: 98
+SCORES: _[_cccccgcaeeaeJb`b^dd_cg]cZedbfU_eghhhddfg[f]BLBLTZ^abbb^aT]b_`]`b`Y`BKT^^ac_YYY_`BBJGQW_b]bbbbbb
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:2326:2150 2:N:0:TAGCTG
+SEQ: CGGCCNNNACCCAATTTGNNTTGGGCATTTCTGANNNNNNNCNACNNTACTGNNGGGTANTGCTNNNTGNNNNGNNNNNNNNNNNNNNNNNNNNNNNN
+SEQ_LEN: 98
+SCORES: WZ^ceBBBQQQQbffgggBBSQ`_fRZdegffdhBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1832:2218 1:N:0:TAGCTG
+SEQ: GCTTTGACATAGTCGCTCCAGAATTGCCAGCTAGGGTTAGCTTGGCAACTGCAGCGACGTAATGTGCTGTGGCAGATCAATTTATCTGTTTTGAATCA
+SEQ_LEN: 98
+SCORES: ^P^PJ\Y`eea`e[daYdecggadgdXJIYVbdc`efg_cdedI^aXIO^abeb\eL_daQU^_V]``]UGTZ\^BBBBBBBBBBBBBBBBBBBBBBB
+---
+SEQ_NAME: HWI-ST575:107:C0HE6ACXX:5:1101:1832:2218 2:N:0:TAGCTG
+SEQ: GGTTATCGATCTGGAAAAAGCAACTAAACCTAAAGCTAAACCACGTAGCGCCGGGTAAATGATTCAAAACAGATAAATTGATCTGCCACAGCACATTA
+SEQ_LEN: 98
+SCORES: ^VYPJQ`c^JJ[b[efg^dHJ`aa`adXd_ZXXbIIIY[af_H^aWHWPZ[`gggFFZ^bd_Z]Zb_]ba\^ZGY_`TZ``cc[bbR]]^aaXQ[bbb
+---