Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451195,5590CA32.8% 71.6 / 40.3 58T355T (ACC→ACA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (19/20);  new base A (10/9);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCG  >  AE015451/195259‑195848
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
cGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCAcc                                                                                                                                                                                                                                                                                                           <  1:196367/301‑1 (MQ=255)
 ggtggtGGCGGGCTGAAGGTGAAATGGCCGCAAACGCCAGAGGTTACCCAAGTGACCGATACAATCGACACTACAACCCCGACCCTGGCCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACcg                                                                                                                                                                                                                                                                                                          <  2:515406/301‑1 (MQ=255)
                 cggggACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCa                                                                                                                                                                                                                                                                                          <  1:638620/297‑1 (MQ=255)
                   gTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGCGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                                                                                                                                                                                                                                        >  1:211608/1‑301 (MQ=255)
                         ttGCCGCAAACGCCAGCGGTTACCCAAGTGACCGTTACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAATGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAccc                                                                                                                                                                                                                                                                                  <  1:437989/301‑1 (MQ=255)
                                                   aGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCAGTAACCGTGTCACTTTCCAACACCCAACTTATTACCGTTTACCCCGTTCAACTCCAGGACAGCGTCGATTTCCAGACCCCGGCCACTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGCTAATTTCGACCAACTGACCCCGAATCCGACCC‑‑CGGCTCAcaccacc                                                                                                                                                                                                                                                      >  1:27052/1‑300 (MQ=17)
                                                      gACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGACCAACTGACCCCGAACCCGACCC‑‑CGGCTCAGACCACGA‑Tc                                                                                                                                                                                                                                                  >  2:737865/1‑301 (MQ=255)
                                                      gACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCAAACGGACAGGCCATTACGGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATAACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAACAGCTGACCCCGAACCCGACCC‑‑CGCCTCAGACCACGacc                                                                                                                                                                                                                                                   >  1:192586/1‑299 (MQ=17)
                                                             aCAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGCGCAACCCTGCCCAGACGCCGGCCACCGTGCACCTGTCCAACGGCCAGGTCAACACTGTCGACGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGATCCCGGCCAATGACGCCCACAACAACGGCTCGACCGTCAGCGTCACCATCGACAACGCCCCCGGCGCCTATTTCAAGCACCTCACCCCGAACCCGACCC‑‑CGGCTcggctacggatccgcgacc                                                                                                                                                                                                                                            >  1:631117/1‑283 (MQ=17)
                                                                                   cccTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGACAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGAc                                                                                                                                                                                                                   <  1:233422/301‑1 (MQ=255)
                                                                                      tGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAAACCAGTGCAGCGTAGATATCCAGACCCCGGCCAATGACGTCTACAACAACGGTTAGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGAcccc                                                                                                                                                                                                                >  1:657792/1‑300 (MQ=255)
                                                                                         gggcatgggcAAGCGTGAACGAAGGTGGCGGGATCAACTACACCGTGAAACTGAGCAACCCTGCGCAGACGCCGGTGACCGTAAGCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGGCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGAc                                                                                                                                                                                                             <  2:383666/292‑1 (MQ=17)
                                                                                          cgacTCGCAAAGCGTGAGCGAAAGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACg                                                                                                                                                                                                            <  1:545299/298‑1 (MQ=255)
                                                                                                  cAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCAACACTGTTAGAGCCGGCAAAAGCCAGGGCAGAGATGACGGCCAGAGCCCGGCCAATGAAAACTACAACAACGGTTCGACCGTAAGCGTCACCATCGAGAACGCCAACGGCGGTAATTTCGAGCAACTGAACCCGAAGCCGACAC‑‑CGGCTTCGACCACGA‑TCACAGACCCGGTCCacacaacgaccccaaccatgaagcagaacac                                                                                                                                                                                                      >  2:94006/1‑274 (MQ=12)
                                                                                                      cGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGCCGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCg                                                                                                                                                                                               >  1:168685/1‑301 (MQ=255)
                                                                                                        ggACCGAAGGTGGCGTGATCACCTACCCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGt                                                                                                                                                                                             <  2:168611/300‑1 (MQ=255)
                                                                                                               aGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCGGCGTCCCAATCGACAGCACAACAGGCCGTAATATCGAGCAAGTAACCCCTAATCCGACTC‑‑CGCCTCACATCACGA‑TCAACGACCCAGTCGACACCAC‑‑CACCGCGACCCTCACCGCGAC‑CCCGTCgcgccccgaa                                                                                                                                                                                      >  1:460157/1‑292 (MQ=14)
                                                                                                                 gTGGCATGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAgg                                                                                                                                                                                    <  2:257411/301‑1 (MQ=255)
                                                                                                                              cTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACt                                                                                                                                                                       <  1:452192/301‑1 (MQ=255)
                                                                                                                                        aCCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCACTCGGTCACCGAAGgcggcgggcataataccacactgg                                                                                                                                                             >  1:304270/1‑283 (MQ=14)
                                                                                                                                            cgAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCt                                                                                                                                                         <  1:527533/300‑1 (MQ=255)
                                                                                                                                                         ccggacagcGGTGGCCGTTACCCTGTCCAACGGCCAGGTCATAACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGGCCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGcc                                                                                                                                            <  1:280594/293‑1 (MQ=255)
                                                                                                                                                                          gtccacTGTCCAACGGCCAGGTAACCACTGTTGAAGCCGGCAAAAGCCAGGGCAGCGTCGTTTACCCTACCCCGGCCAATGACGTCTACAACCACGGTTCGACCGTCAGCGTCACCATCGATAACGCCACCGGCGGTAATTTCGAGCACCTGACCCCGAATCCGACAC‑‑CGGCTCAGACCACGA‑TCAACGACTCAGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGt                                                                                                                           <  1:441678/296‑1 (MQ=14)
                                                                                                                                                                           gaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCTGTCACCGAAGGCGGCGTGATCACCTACAGCGTGACCCTGAGCAACCCTGCCCAGACGCCTGTAACCGTg                                                                                                                          >  1:92765/1‑301 (MQ=255)
                                                                                                                                                                             ccctcaccaacgcccgtgTCATCACTGCTGAAGCAGGCAAAACCCAGAGCAGCGTCGATTTCCAGACCCCGACCAATGACGTCTACAACAACGGTTCGCCCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGCATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGaa                                                                                                                        <  1:518685/284‑1 (MQ=14)
                                                                                                                                                                              cgtGTCCAACGGCCAGGTCATCACTGTAGAAGCCGGGAAAACACAGTGCAGCGTCGATGTCCAGACCTCGGCCCATGACGCCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATGTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAAc                                                                                                                       <  1:15539/299‑1 (MQ=14)
                                                                                                                                                                                       gggCCAGGGCAGCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTGCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAGCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑Ac                                                                                                             <  1:439096/300‑1 (MQ=255)
                                                                                                                                                                                          ccAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAACCAGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGAGAG‑CCCGTCGGTCACAGAAGGCGGCGTGAGCACTTACACCGTGACCCGGAGCAACACTGCGCAGACGCCGGTACACGTGAACGTGTACA‑ACggg                                                                                                          >  2:18545/1‑300 (MQ=14)
                                                                                                                                                                                               tcatcaCTGTTTAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAc                                                                                                     <  1:633157/301‑1 (MQ=255)
                                                                                                                                                                                               tcataaCTGTTGAAGCCGGCAAAACCCAGGGCGTCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAc                                                                                                     <  1:94044/301‑1 (MQ=255)
                                                                                                                                                                                                              ccgcaacaaCCCAGAGCAGCGTAGATTTCCTAACCCCAGCCAATGACGTCTACAACAACGGTGCGAACGTCAGCGCCACCATGGAGAGCAGCACCGGTGGTAATGTCCAGAAACTGACACCGAATCCGACGC‑‑CGGCTCAGACCACAA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGc                                                                                     <  1:603627/294‑1 (MQ=14)
                                                                                                                                                                                                               cGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTAGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGAACCTGAGCAACCCTGCCCAGACGACGGTAAACGTGAACCTGTCCA‑ACGGCCACAC‑GCTTACCGTTGAAGCg                                                                                    >  1:11010/1‑301 (MQ=255)
                                                                                                                                                                                                                ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAGCCTGCCCAGACGCCGGGAACGGTGACCCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCgg                                                                                   >  2:567515/1‑301 (MQ=255)
                                                                                                                                                                                                                ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCgg                                                                                   >  1:238345/1‑301 (MQ=255)
                                                                                                                                                                                                                 gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACAAC‑‑TCGCGAGACCCTGACTGCGAA‑CCCGACGGTCAGCAAAGGAGGCGGGATCACTTACACCGGTACCCTGAGCAAACCTGCCCAGACGACGGGAACAGTGAACATGTCCA‑AAGGCCAAAA‑GATTACCGGTGAAGcagg                                                                                  >  1:403206/1‑298 (MQ=11)
                                                                                                                                                                                                                 gCAAGACCCAGGGCAGCATCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAATTCGAGCAACTGACACCGAAAACGACCC‑‑CGGCTCACACCACGA‑TCAACGACTCGCTCCACACCAC‑‑TACCGCGACCCTGACGGCAAG‑CCCGTCGACCACCGCAGGCCGCGTGACCACTTACACTCCGACTCCGAGCAACCCTGCCCAGACGGCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GAGTACCGTTgcaccggg                                                                                  >  1:385606/1‑294 (MQ=11)
                                                                                                                                                                                                                 gCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCggg                                                                                  <  1:249674/301‑1 (MQ=255)
                                                                                                                                                                                                                 gCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAACAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACATACACCGTGACCCTGAGCAATCCTGCCCAGCCGCCGGTAACTGTGAACCAGTccgaaacgcgcagtgtagggaaggagtgagg                                                                                    >  2:317720/1‑272 (MQ=11)
                                                                                                                                                                                                                   aacatccatgtcAGCGTCGACGTCCACACCCCGGCCAACAACGTATACAACAACGGCTCTCTCCTCAGCATCACCATCGAGAACGGCACCGGCGGTAATGTCGAGCAGCTGACACCGACGCCGACGC‑‑CAGCGCTGACCACGA‑TCAACGACACGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTGACCGAAGGCGGCGTGATCACCTACGCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTa                                                                                <  2:745530/290‑1 (MQ=11)
                                                                                                                                                                                                                      gcccactgcaccgtggacttCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAAACCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTaaaa                                                                             <  1:355215/283‑1 (MQ=14)
                                                                                                                                                                                                                      aCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTaaca                                                                             >  1:295697/1‑299 (MQ=255)
                                                                                                                                                                                                                          gtgccgacgtctgcttccaccacccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGATGCCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCCGGCAAAACCCa                                                                         <  1:28078/279‑1 (MQ=11)
                                                                                                                                                                                                                          agggacacGTCGATTTCCAGTACCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTACTTTCGCGCAACTGACCCCTACACCCACCC‑‑CACCTCATACCACGA‑TCAACGCCTTGGTCGACGCTCC‑‑GACCTCTACCCTTCCGACCTG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCa                                                                         <  1:440001/294‑1 (MQ=14)
                                                                                                                                                                                                                              gaGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGATGCCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCCGGCAAAACCCAGGGc                                                                     <  1:269151/300‑1 (MQ=11)
                                                                                                                                                                                                                                  gTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGCGCACCACCGCCGGTATTTTCTACCATCTCCCACCGAATCCGCCTC‑‑TGGCTCATACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGCGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGCTAAAACCCAGGGCAGCg                                                                 >  1:612866/1‑301 (MQ=14)
                                                                                                                                                                                                                                      aTTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACACATCGG‑‑CAGACCACAAGTC‑TGAACTCCAGCCACCGCACTTAATCTCG‑‑TATGCCGGCTTC‑CGCTTCAGACAACAAAGGCGGTGTGATCACTCAAACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCC‑ACCGCCCAAC‑GATTACCGTTGGAGCCGGTACAACCcagcgcagcgctca                                                             >  1:629997/1‑297 (MQ=2)
                                                                                                                                                                                                                                         cccAGACCCCGGCCAAAGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGAttt                                                          <  2:238255/300‑1 (MQ=255)
                                                                                                                                                                                                                                          cctgaccgcggcgaataacgtaaacaacaACGGTTCGACAGTCAGCGTCAACATCGAGAACGCAAGAAGCGGTAATTACGAGCAGCTGACCCCGAATCCGACGC‑‑CGGCTCAGACCACGA‑TCAACAACTCTGTCGACGCAAC‑‑GGACGCGACCCTTGCTGC‑AGTCCCGTCACTGACCTTCGGAGGCGTGCTCATCTACCCCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTc                                                         <  1:593655/278‑1 (MQ=2)
                                                                                                                                                                                                                                           cAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCAGCGGCGGGAATTGGGAGAGAGTGTCCATGAATGCGACCC‑‑CGGCACAGACAACGA‑ACAAAGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑ACCGACGGTCCGCGAAGGCGGCGTTACTACTTACACCCTTACCCTGAGCAACCCTGCCTAGACTCCGGTAACCGTGAACGTGTCCC‑ACAGCCACAC‑GATTCACGTTGCACCGGGAAAAACCCATGGCacactcgatttct                                                        >  2:84146/1‑289 (MQ=9)
                                                                                                                                                                                                                                                 ccGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAGCTCACACCGAATCCGACCC‑‑CGGCTCAGGCCACGA‑GCAATGACCCGGTAGATGCTAC‑‑CACCGCGACCCTGACGGAGAG‑AACGTCGGTCACCGAAGGCGGCGTGATCCCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGATCGTGACCCTGTCCA‑ACGGCCAAAC‑GATTACAGCTGAAGCCGCCAAAACCCAGAGCCGCGTCGACTTCAAGAccc                                                  >  1:40368/1‑301 (MQ=11)
                                                                                                                                                                                                                                                   ggctaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCTAGAACGCCACAGGCGGTAACTTCGAGCAACTGAACCCTAATCCGACCC‑‑CGTCTCAGACTACGA‑TCAACGACTCTGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCg                                                <  1:302258/297‑1 (MQ=14)
                                                                                                                                                                                                                                                     cgaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                              >  1:187347/3‑301 (MQ=255)
                                                                                                                                                                                                                                                              tcgaaaacaacgtggtgaccgtacgCGTCACCATCGAGAACGCCACAAGCGGTAATTTCGCGCAACTGACCCCGAAGCCGACCC‑‑CGGCTCAAACCACGA‑TCAACGACTCGGTTGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGTCTAGACGCTGCTAACAGCAAAGAGGTCCACAC‑GCTATCC‑GATCACCGTGCAAGCCCGTAATACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt                                     <  2:295583/277‑1 (MQ=2)
                                                                                                                                                                                                                                                                  caacaaCGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACGGTTGCAGCGTGTAAAACCCAGGGCGGCTTCGGTTTCAAAAACCCGGCCAATGACGTCTac                                 >  1:77318/1‑301 (MQ=14)
                                                                                                                                                                                                                                                                   aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAATTTCGAGCACCTCACCCCGAATCCGACCG‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACTGTGACCCTGTCCA‑ACGGCCCAAC‑GATTACCGTTGAATCAGGTAAAACCCAGGGCAGCGTCCATTTTCAGACCCCGGCCAATGACGTCTacc                                >  1:187925/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                         ggTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑TACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCAGAGCAACCCTGCGCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGACCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATTGCGTCTACAACAACg                          >  1:726083/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                    cgggacgaCCATCGGGACCGCCACAGGGGGTAAATTCGAGCCACTGACCCCGAATCCCACCC‑‑CGGCTCAGACCACGA‑TCAACCGCCCGGTCGGCAGCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTc               <  2:294439/294‑1 (MQ=14)
                                                                                                                                                                                                                                                                                       cGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGCCCCCGGCCAATGACGTcaaccagaacgacgcacgggtaggg            >  2:546770/1‑277 (MQ=14)
                                                                                                                                                                                                                                                                                       cGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGACCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCAAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGc            >  1:366435/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                 aagaACGCCACCGGCGGTAATGTCGAGCAACTGACCACGAGGCCGACCA‑‑CGGCTCATACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCGTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGGGTGATCACTTACACCGTGACCCTGAGCAACCCTGCTCAAACGCCGGTAACCGTGAACCTGGCCA‑ACAGCCAAACTAATAACCTTTG‑AGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCg  <  1:84183/300‑1 (MQ=17)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
CGGTGGTGGCGGGCTGACGGTGACATTGCCGCAAACGCCAGCGGTTACCCAAGTGACCGATACAATCGACACTACAACCGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGTGGCGTGATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCC‑‑CGGCTCAGACCACGA‑TCAACGACTCGGTCGACACCAC‑‑CACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCA‑ACGGCCAAAC‑GATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCG  >  AE015451/195259‑195848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑