Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451196,4320TC27.1% 68.8 / 39.4 62T646T (ACT→ACCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base C (6/10);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCG  >  AE015451/196132‑196729
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                            
ggTCAGCGTCACCATCGAGGACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACTGCGAGCCCGTCGGTCACCGAAGGCGTCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                                                                                                                                                                                                                              <  1:179078/301‑1 (MQ=255)
    aGCGGCACCATCGAGACCGCCACGGGCTGCAATTTCGAAAAACTGACCCCCAACCCGACGCCGCCACAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCCGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCa                                                                                                                                                                                                                                                                                                          <  1:422873/301‑1 (MQ=11)
     gcgccACCATAGAGAGCGCGACGGGCGGCAAGTTCGGACAACGGGACCCGAACCCGACGCCGGGTCAGACCAGGAGCAAAGAATCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATAACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAg                                                                                                                                                                                                                                                                                                         <  2:84259/297‑1 (MQ=11)
      ggttaccctggggaACGCCACGGGCGGGAAAGTAGAGCAAATGACCCCGATCCCGTCGCCGGGTCAGAGACCGAACAACGGCTCGGGCGCCACCACTGGCGTTAACCCGACGGCTATCCCGTCGTTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgt                                                                                                                                                                                                                                                                                                        <  1:493892/289‑1 (MQ=14)
                            ggcagTAAGTTCGTGCAACTGACCCAGAACCCGACGGCGGCTGAGAACAAGATCAACGACTCGGTCGACACCACTGCCGCGACCCGGACGGCGAGCCCGTCGGTCACGGAAGGCGACGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAACCTGTCCAGCGGCCAAACCATTACCGTTGAAGCGGGCAAAAGCAATGACAGCGTCGATTTCCCGACCCCGGCCAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcg                                                                                                                                                                                                                                                                                  <  1:104150/297‑1 (MQ=14)
                            ggCGGTAATTTCGAGCAACTGACCCCGAACCCGCCGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCAcag                                                                                                                                                                                                                                                                                  <  2:352975/301‑3 (MQ=12)
                                  aaTTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAggagaccggacgaggatcctacaggcaagg                                                                                                                                                                                                                                                                            >  2:418806/1‑272 (MQ=11)
                                              cTGACCCCGAACCCGACGGCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCCGTGCCCTAGAAGCCGGCACAACCCAGCGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                                                                                                                                                                                                                <  2:450134/301‑1 (MQ=14)
                                              cTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCCAACCATTACCGTTTAAGACGTCAAAACCCCGGGCATAGTCGATTACGCGACCCCGTCGAATGACGTCTACAACAACGGTTCGACTGTCAGTATCACCATCGAGAACGCCACCGGCGTTAATTTCGAGCAAc                                                                                                                                                                                                                                                                >  1:568276/1‑301 (MQ=14)
                                                   cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTACAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGATCATCTGAcc                                                                                                                                                                                                                                                          >  1:50662/1‑301 (MQ=255)
                                                         cccGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACTGCGAGGCCGTCGGTCACCGACGGCGTCGAGATCACTTACACCGTGACCCTCAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAAt                                                                                                                                                                                                                                                    <  1:632180/301‑1 (MQ=255)
                                                                       aGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGcccc                                                                                                                                                                                                                                      >  1:331123/1‑298 (MQ=255)
                                                                                  aacggtgcgggcGACACCACTACCGCGGCCCGGACGGCGAGCCCGTCGGTGACCGAAGGCGGCGGGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTTACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCTGGCTAATCACACGGCAGCCACCGTGACCAAACCCCTGCCCATAACGGCCACTCTACCGATCTGACGGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCa                                                                                                                                                                                                                           <  2:667530/290‑1 (MQ=2)
                                                                                             cGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCACGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGACTCCGACGCCGGCTCAGACAACGATCAACGACTCgttc                                                                                                                                                                                                                >  1:455417/1‑298 (MQ=255)
                                                                                                   cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGcc                                                                                                                                                                                                          <  1:681047/301‑1 (MQ=40)
                                                                                                    acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACAGCCCGACAATTCGCCTTGTAGCGGGCAAGTCTAGGTCCCGGTGGGCTTCGCAGACTCCGACGAATGACGACTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAACTTCAAGCAACTGACCCCGAATCCGACCACGGCTCCAACCACGATCAACGACTCGGTCGATGccc                                                                                                                                                                                                         >  2:725523/1‑300 (MQ=2)
                                                                                                                                     aCCGAAGGCGTCGTGATCAATTACACCGTGACCCTAAGCAATCCTGCCCAGACGCCGGTGACAGTTACACTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTTAGTGTCACCATCGAGACCGCCACCGCCGATAACTACAATCTACCGACACCCGACCCGCTGCCGATCCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                                                                                                                                                        <  2:330996/301‑1 (MQ=11)
                                                                                                                                                       aCTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTTCAGCGGCCAAACCATGTCCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCa                                                                                                                                                      <  1:156692/301‑1 (MQ=255)
                                                                                                                                                                     cctaagcaagcccttccgatcGCCGGTAACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGAcc                                                                                                                                        <  1:58732/281‑1 (MQ=255)
                                                                                                                                                                           cAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAGCTGACCCCGGATCCGACGCCCGCACATACGACGATCAACGACTCGGCCTATGCTACCATCCCGACCCTGACGGCGAACACGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTCACCCTgcgc                                                                                                                                  >  1:405609/1‑298 (MQ=14)
                                                                                                                                                                              gccTGCCCCGACGCCGGAGTCAAGGACCCTGACCAAAGGTCAAACCAGTACCGTTGAAGACGGCAAAACCCAGGCCAGCTTCGATGGCCAGAGCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAACACGCCACCGGTTTTAAATGCGACGAACTCACCCCAGATCCGCAGCCGGCTCCGACCCGGATCTCCTCCGCGCTCGATGCCACAACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAAt                                                                                                                               <  2:455227/300‑1 (MQ=9)
                                                                                                                                                                                    ccggacgcacggaacaatgccccTGTCAAAAGGCCAAACCATTAGCGTTCAGACCGGCAAAACCCCGGGCTGCGTCGATTTCCAGACCCCGGCAAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                                                                                                                         <  1:426967/281‑1 (MQ=14)
                                                                                                                                                                                       gacgtgtgctcttccgatcttGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAg                                                                                                                      <  1:338955/281‑1 (MQ=255)
                                                                                                                                                                                            cccgacgcgtttccgtgctcAACGGCCAAACGATTACCGTTGAAGCCGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACGGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                                                 <  2:299174/282‑1 (MQ=11)
                                                                                                                                                                                                           tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCGGCCCAGagatcggaagagcgtcgtgt                                                                                                  >  2:338828/1‑282 (MQ=255)
                                                                                                                                                                                                                  cGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                          <  1:725832/301‑1 (MQ=255)
                                                                                                                                                                                                                      cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGCCGAGCCCGTCGGTCACCGAAGGAGGCGGGAACACCTACACCGTGACCCTAGGCAATCCTGCCCCGAcacccgtagacacggccgggtctgacgtag                                                                                       >  2:498626/1‑272 (MQ=11)
                                                                                                                                                                                                                           gaTTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGCACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTCCACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGCGACCCTGTCCACCGGCCAAAcc                                                                                 >  1:328883/2‑301 (MQ=11)
                                                                                                                                                                                                                           ccatACCGTTGAAGCCTGCAATACCCAGGTCACCGTCGATTTCCACACCCCGGCAAATGACGTCGACATCCTCGGTTCAACCGTCAGCGTCACCCTCCAGAACGCCACAGGCGGTAATTTCGAGCGACTGACCCCGTATCCGACTCCTGCTCAGACCACGATCCACGACTCGGTCCATGCTACCGCCGCTACCCTTCCGGCGAGCCCGTCTGTCACCGAATGCGGCGTAATCACCTACACCTTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAcc                                                                                 <  1:80403/298‑1 (MQ=2)
                                                                                                                                                                                                                              cgctcttcggatctgCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCAtt                                                                              <  2:528221/287‑1 (MQ=255)
                                                                                                                                                                                                                                ccgttaaaggaggCAAAACCCAGGGCAGACTCGATTTCCAGGCCCGTGGGAATGACCTCTACCACAACGGTTCGACTGTCAGAGTCAACATCGAGAACGCCACAGGCGGTAATTTCGAACAACTGTCCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCGCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTAAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTAc                                                                            <  2:360696/290‑1 (MQ=11)
                                                                                                                                                                                                                                      aaGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGTATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGa                                                                      >  1:363022/1‑301 (MQ=11)
                                                                                                                                                                                                                                            gCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTagatcggaagagcc                                                                >  1:528442/1‑288 (MQ=255)
                                                                                                                                                                                                                                              aaaaCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCCACACCGTCACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCCTTGAAGCCGGCa                                                              >  1:722378/1‑301 (MQ=11)
                                                                                                                                                                                                                                                      gggCAGCGTCGAGTTCCATACCCCGGCGAACGACGTCTACAACAACGGTTCGACCGGCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAACAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACGCGGTCGATGCCACCACCGCGACCCGGACGGCGAGCACGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                     <  1:382712/301‑1 (MQ=11)
                                                                                                                                                                                                                                                        gCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAggg                                                   <  1:636277/301‑1 (MQ=255)
                                                                                                                                                                                                                                                          aGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCAACCCTGACCTCGAGCCCGTCGGCCACCGATGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTCACAGTGACCCTGTCCAACGGCCAAACCCTTACCGTTGAAGCCGGCAAGACCCAgtgcc                                                 >  1:28997/1‑297 (MQ=11)
                                                                                                                                                                                                                                                           gCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGATGACAGTGACCCGGTCCAACGGCCAAACCATTACCGTTGAAGCCGTCAAGACCCAGGGCac                                                >  1:101211/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                  cttCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAt                                         <  1:57745/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ttCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGACACCACCGCGACCCTGACAGCGAACCTGTCCGTGACCGACGGCGGCGTAATCACCTACACAGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGaat                                        >  1:89898/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                    tCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTAACGTTGAAGCCGGCAAGACCCAGGGCAGCCTCGaata                                       >  2:57721/1‑298 (MQ=11)
                                                                                                                                                                                                                                                                     gtctgACCCCGGCGAAGGACGTCTACAACAACGGTTCGACCGTCCGCGGCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAttt                                       <  2:702519/297‑1 (MQ=11)
                                                                                                                                                                                                                                                                     ccAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGAACCAGGGCAGCGTCGATTTc                                      >  1:493289/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         gACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtttagat                                   >  1:702819/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                          aCCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGa                                  >  1:225578/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                               gggtgaatgacagccacatcctcgtccgaaCTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAcccc                              <  1:506346/272‑1 (MQ=11)
                                                                                                                                                                                                                                                                                 cGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCCACCCTGACGGCGAGCCCGCCGGTCACAGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCggg                           >  1:42947/1‑300 (MQ=11)
                                                                                                                                                                                                                                                                                    aTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAACGCGGCCAGACCCAGGGCAGCgctcattgccaacccccggcgaa                        >  2:444540/1‑279 (MQ=11)
                                                                                                                                                                                                                                                                                        cgtatACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCCAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGGGCAATCCTACCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                    <  1:689009/297‑1 (MQ=11)
                                                                                                                                                                                                                                                                                           cTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACTCTGCCCAACGGCCAAACGATTACAGTTGAAGAGGGCAAACACCAGGCCAGCGTCGacgtaaataacccggcaaaggacgtc                 >  1:219116/1‑276 (MQ=11)
                                                                                                                                                                                                                                                                                              caacaaCGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCGGACGCACGCCACGATCGACGACAGGGTCGATCCCACCACCCCGCGCCTCACGGCGAGCACGTCGGTCACCGAAGGCGGCGTGATCACCTACACCATGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATGACCGTTGAAACCGGCAAGACCCAGGGCAGCGTCCATTTCCAGACCCCAGCGCATGaacttcac              >  2:379740/1‑294 (MQ=14)
                                                                                                                                                                                                                                                                                               aacaacGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCCTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTCCCGTTGAAGCCGGCAAGACCCAGGGCACCGTCGATCTCCAgtccgcacaccatgacgcttacc             >  1:34123/1‑279 (MQ=14)
                                                                                                                                                                                                                                                                                                acaacGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa            >  1:172200/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                acaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGCTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGCGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACCTCTacaa            >  1:351962/1‑301 (MQ=11)
                                                                                                                                                                                                                                                                                                 caacGGTTCGACGGGCAGCGTCACCATCGACAACGCCACAGGAGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTAGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGGGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaac           <  1:104760/301‑1 (MQ=11)
                                                                                                                                                                                                                                                                                                   acGGTTCGACCGTCAGCCTCACAAGCGAGAACGCGACAGGCTGTAATTTCGGGCAACTGACCCCGAAACCTACCCCGGCTCAGACCATGATCAACGGCTCGCTTGATGCCACCACCGCGACCCTGACCGCGAGCCCGTCGGTAACCGAAGGCGGTCTAATCACCTACACCGTGACCCGGAGCAATCCTGGCCATACGCCGGTGACAGTGACCCTGTCAAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGAAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa         <  1:654909/301‑1 (MQ=11)
                                                                                                                                                                                                                                                                                                    cGGTTCGACTGTACGTGTCACCATCGAGAACGCCACCGGCGGTAACTTCGAACAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac        <  2:172126/301‑1 (MQ=12)
                                                                                                                                                                                                                                                                                                     ggTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg       <  1:444730/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     ggTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTCACTCTGAACAATCCTCCCCAGACCCCGGTTCACGTAACCCTGTCCAACGGCCAAACCATTACCGTTGAAGACGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGAGAATGACGTCTCCAACAACg       >  1:4968/1‑301 (MQ=11)
                                                                                                                                                                                                                                                                                                         cGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAGTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc   <  1:708727/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                          ggACTGTCAGTGTAAACATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCGCAGACCACGATCAACGACTCGGTAGATGCCACCACCGCGAACCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGTAATCCTGCCCAGACGCCGTTGACAGTGGTGCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc   <  1:50884/300‑1 (MQ=14)
                                                                                                                                                                                                                                                                                                           gACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAACCCGGCAAGACTCACGGCACCGTCGATCTCCAGCCCCCGGCGACTGACTCCTACAACAACAGTTca  >  1:181743/1‑300 (MQ=14)
                                                                                                                                                                                                                                                                                                           gACTGTCAGTGTCACCAGCGAGAACGCCACCGGCAGGAATGCCGAGCAACTGTCCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCg  <  1:740590/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                            
GGTCAGCGTCACCATCGAGAACGCCACGGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCG  >  AE015451/196132‑196729

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