From: martinahansen Date: Fri, 1 Jun 2012 09:43:20 +0000 (+0000) Subject: added illumina 1.8 name scheme support to order pairs X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;h=849409088c8d611089191593ec8ee3bb2a95db45;p=biopieces.git 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 c96f845..e8f7990 100755 --- 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 3e30de4..f9a2a66 100644 --- 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 b0906df..e25e55d 100644 --- 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 +---