Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451197,6320TC25.4% 82.1 / 18.7 56T1046T (ACT→ACC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (23/18);  new base C (6/8);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAC‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAC  >  AE015451/197335‑197931
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                                
cAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTc                                                                                                                                                                                                                                                                                                               >  1:827712/1‑301 (MQ=37)
 aGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCa                                                                                                                                                                                                                                                                                                              >  1:870390/1‑301 (MQ=37)
      cACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACACTGTCCAACGGCAAAACCATGAACGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAATTCCAGAACCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCCGCGTc                                                                                                                                                                                                                                                                                                         >  1:615355/1‑301 (MQ=255)
            cGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGCCAGCGTCGATTTCCACACCCCGGCTAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑Tc                                                                                                                                                                                                                                                                                                  >  1:18852/1‑301 (MQ=255)
             gggaACGCCACCGGCGGGAACTTCGCACAGCTGACCCCGAACCCGGCCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATCACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCTAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCg                                                                                                                                                                                                                                                                                                 <  1:520624/299‑1 (MQ=255)
                  cGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTTCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATG‑ACGTCTACAACAACGGCTCGACTGTCAGCGTCACCA‑TCGAGAAc                                                                                                                                                                                                                                                                                            >  1:840596/1‑301 (MQ=255)
                  cGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGAGGCCACCACCGCGACCCTGACGGCGAGCCCGGCGGGCAGCGAAGGCGGCGTCATCACCG‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAAc                                                                                                                                                                                                                                                                                            <  1:121555/301‑1 (MQ=255)
                     cacaggcagtaatccccgtcaactgccgccGAAGCCGACCCCGGCTCAGCCCACGATCACCGACTCGGACGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAACT‑ACACCATGAGCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCTGTGAAGCCGGCAAGCCCCAGCGCAGCCTCGCTTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGcc                                                                                                                                                                                                                                                                                         <  1:489630/273‑1 (MQ=11)
                     cACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGCTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGcc                                                                                                                                                                                                                                                                                         <  1:733480/301‑1 (MQ=37)
                      aCCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGGTTTCCAAAACCCGACGACTG‑ACGTCTACAACACCGATTCGACCGCCTGCGTCACCA‑TCGACACCGCca                                                                                                                                                                                                                                                                                        >  1:266267/1‑301 (MQ=11)
                       ccgacggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCcac                                                                                                                                                                                                                                                                                       <  1:850274/297‑1 (MQ=25)
                           cggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCCAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGACAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATG‑ACGTCTACAACAACGGCTCGCCCGTCCGCGTC‑CCACTCCAGAACGccccggcg                                                                                                                                                                                                                                                                                   >  2:617148/1‑295 (MQ=11)
                            ggcagcttcgactaccggactccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCg                                                                                                                                                                                                                                                                                  <  1:244758/280‑1 (MQ=255)
                            agacgtgtgctcttccgatctccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCg                                                                                                                                                                                                                                                                                  <  1:441320/280‑1 (MQ=255)
                                        caacgcagcccgaacctgtccccggacaaGACCACGCTCAAAGACTCGGTCGATGCCACCCCCGCGACCGTGACGGCTAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAGTCATCACCGTTGAAGCCGGCAAGTCCCAGGTAAGGGTCGGGCGACAGACCACGGCTAAAC‑ACGTCCACAACAACAGGTCGACCCGCAGTGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGc                                                                                                                                                                                                                                                                     <  2:339138/273‑1 (MQ=2)
                                          gCTGAACCCGAAACCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGCTGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCAGCGAATG‑ACGCCTACAACAGCGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCaa                                                                                                                                                                                                                                                                   <  2:827422/301‑1 (MQ=255)
                                           tccgatctgAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAAc                                                                                                                                                                                                                                                                  <  1:705920/293‑1 (MQ=17)
                                           cTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAACG‑GGAGGGCACAAGTCACAACTCCAGCCACCCCAAGTCTGCCCGTACGCCGGCTTCCGTGAGCCAATCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATC‑ACCTGTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAAAGCCACACGCGGTAATTTC‑GAGCAAc                                                                                                                                                                                                                                                                  >  1:129810/1‑301 (MQ=2)
                                            cgACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACt                                                                                                                                                                                                                                                                 <  1:35888/300‑1 (MQ=17)
                                                 ccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACTGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCATCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGATCCAGGGCCGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCCACAAACACGGTTCGACCGCCAGCGCCACCC‑TGCAGAACGCCACAGGccataacggaagagaggcgtgtg                                                                                                                                                                                                                                                             >  2:441175/1‑279 (MQ=11)
                                                   gAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCAACAGGCGGTAATTTC‑Gacaaacggatcgga                                                                                                                                                                                                                                                          >  2:705674/1‑288 (MQ=12)
                                                     aCCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACACAGGGCAGCGTCGCATTCCAGACCCCGGCGAATG‑ACGTCTGCAACCACCGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCACGaa                                                                                                                                                                                                                                                        >  1:517539/1‑301 (MQ=11)
                                                         gaccccGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTT‑ACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCACCTGACCCCGAAcccc                                                                                                                                                                                                                                                    >  1:261030/5‑300 (MQ=11)
                                                          aCGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGCCCCAGGGCAGCGTCGATTCCCATCCCCCGGCGAATC‑ACGTCCATCACACCGGTTCGACTGTCCGCTCCACCA‑TCGTGCACGCCACAGCCTGCAATTTC‑GAGCCACTGACCCCGAACCCgc                                                                                                                                                                                                                                                   >  1:375528/1‑300 (MQ=11)
                                                           ccccGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTT‑ACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGAc                                                                                                                                                                                                                                                  >  1:274813/3‑301 (MQ=255)
                                                            gCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCCACGGCCAAACCATTACCGTTGAAGACGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAATCCGAcc                                                                                                                                                                                                                                                 >  1:342068/1‑300 (MQ=2)
                                                              cgccTCAGACCACGATCAACGCCTCGGTCGATGCCACCACCGCGACCCTCACGGCGATCCCGTCCGTCACCGAAGGCACCGTGATCCCCT‑ACACCGTGATCCTGCGCAATCCTGCCCAGACGACGTTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAATCCGAcccc                                                                                                                                                                                                                                               <  1:782841/298‑4 (MQ=11)
                                                              cGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAACGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATG‑ACGTCTACAACAACGGCTCGACGGTCCGCCTACCCT‑TAGAGAACCCCACGGACGGTAACTTC‑GAACATCTGACTCCGACCCCGAcccc                                                                                                                                                                                                                                               >  2:329755/1‑298 (MQ=11)
                                                               ggCTCAGACCACGATCAACGACTCGATCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTT‑ACACCGTGCACCTGAGCAAACCCGCCCAGACGCCGGTGACCTTGACTCCGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACCCCg                                                                                                                                                                                                                                              <  1:73386/301‑1 (MQ=11)
                                                                  ccAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGCCGGCGTGATCACTT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCt                                                                                                                                                                                                                                           <  1:272954/300‑1 (MQ=255)
                                                                    aGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCTTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCa                                                                                                                                                                                                                                         <  1:449213/301‑1 (MQ=38)
                                                                                            gtGCCACCCCCGCGACCCTGAAGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGa                                                                                                                                                                                                                <  1:450776/300‑1 (MQ=255)
                                                                                            aTGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACCGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACCCCGGCCCAGACTACCATCACCGACTCG‑GTCGa                                                                                                                                                                                                                >  1:263811/1‑301 (MQ=255)
                                                                                                          aCCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGCGATCACTT‑ACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGAACCTGTCCAACGGCCAAACAGTTAGGAATGAAACCGGCAAAAACCAGCGGAACGCCGAGTTCCATCGCCTGCCGAATG‑ACGTTTACAACAACGGTTCGACTGTATGCGTCACCA‑TCGAGAACGCCACAGGCGCGACTTTC‑GAGCATCTGACTCCGAACCCGACACCGGCTGAGACCACAATCAACGACTGG‑GTcactgcgacgtccgctg                                                                                                                                                                                                  >  1:654352/1‑285 (MQ=2)
                                                                                                                           gTCGGTCACCGAATGCGGCGTGATCACCT‑ACACCGTGATACTGAGCAATCCTGCCCAGACGCAGGTAACAGTTATGCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCa                                                                                                                                                                                <  1:227282/301‑1 (MQ=255)
                                                                                                                                                 aTCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGCTCGACTGACGGCGTCACCA‑TCGAGAACGAGACACGCGGTAATATC‑GAGCAACTGACCCTCAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑TCCTCGCCAAGCGTGACCGAAGGCCGCGTGa                                                                                                                                                          >  1:135830/1‑301 (MQ=17)
                                                                                                                                                                                 cccAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGCAATTTC‑GAACAACTGACCCCGAACCCGACCCCGGCCCCGACCACGATCACCGCCCCG‑CTCGACTCCACTACTGCGACGCAAACC‑GCCTCGCCAAGCGTGACCGAAGGCGACGTGATCACCTACACCGTGACCCTGAGCAATCCTGc                                                                                                                           >  1:65471/1‑301 (MQ=18)
                                                                                                                                                                                            gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACCCCGGCCCAGACCACGATCACTGACACG‑GTCGACACCACTACTGCAACCCTGACCAGCACGGCCAA‑CGTGACCGAAGACGACGAGCTCTCCTACACCGTGACCCTGCGCAAACCTGcactgacgacgg                                                                                                                >  1:563195/1‑290 (MQ=14)
                                                                                                                                                                                                 cGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACTGAC‑CGCGTCGACACCACTACTGCGACCCTGCCC‑GCCCAGCCAAGCGTGACCGAAGGCCGGGTGACCACCTACACCGTGACCCTAAGCaagccgcccgagcccccgttaacc                                                                                                           >  1:400394/1‑279 (MQ=17)
                                                                                                                                                                                                              aaCGGGCAAAGCATTACGGTTGAAGCCGGCAAAACCCATGGCGGCGTCGATTTCCAGGCCCCTGCGAATG‑ACGTCTACAACAACGGTTCGACCGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAACAGCTGACCCCGAGCCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCCCTACTGCGACCCTGACC‑GCCTCGCCAAGCGTTACCAAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCa                                                                                              <  1:110066/301‑1 (MQ=17)
                                                                                                                                                                                                                cGGCCTAACCATTTCCGTTGAAGCCAGCAAGAACGAGGGCAGCGCCGATTTCCAGACCCCGGC‑CATGTGCGTCTACGACATCGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAAc                                                                                            <  1:554963/301‑1 (MQ=18)
                                                                                                                                                                                                                                gtgaaaGCCGGCAAGACCCAGGGCAGCGTACATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCAACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCCGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCg                                                                            <  1:140603/297‑1 (MQ=255)
                                                                                                                                                                                                                                 ttggaGCCCGCAAGACTCATGGCCGCGTCGATTTCCAGACCCCGGCCAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGt                                                                           <  1:365694/297‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAGCCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAg                                                                      <  1:830603/301‑1 (MQ=255)
                                                                                                                                                                                                                                        cGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCACACCCTTACCGTTGAAGCg                                                                    >  1:615078/1‑301 (MQ=255)
                                                                                                                                                                                                                                          gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCggg                                                                  <  2:614865/301‑1 (MQ=255)
                                                                                                                                                                                                                                            aaacacccgggCAGCGTCTCTTTCCAGCACCCGACCAATG‑ACTGCTACAACAACGGTTCGACCGTCAGAGTCAACA‑TCGAGAACCCCGCGTGCGGCAATCTCTGA‑CAGCTGCCCCCGAACCCGACGCCGGCTCAGACCTCGATCAACGGCTCG‑GTCGACACCACTACCGCGACCCCGACC‑GCGTCGCCGAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTa                                                                <  2:492912/293‑1 (MQ=2)
                                                                                                                                                                                                                                                              gATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGCAACTTC‑GAGCATCTGACCCCGAACCCGACCCCAGCCCCGACCACGATCACTGACTGC‑ATCGATACCGCTACCGCGACCCTGACG‑TGCACGCCAAGAGTGAACAACGGCGTGGTGGCCGCCTCCTCCATAACCCTGAGCAATCCTGCCCAGACGCCGGTCACCGTGACCTTCACCACCGGCCAAACCATTACCGCTGAACCGGTTAAAAC‑CCAgcgccacgtgc                                             >  2:6287/1‑291 (MQ=9)
                                                                                                                                                                                                                                                                  tCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGGTGAAGCGGGTAAAAC‑CCAGGGCAGCGTCGAttt                                         >  1:242142/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                        cccccgtctaattaCGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑CCGGGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAC‑CCAGGGCAGCGTCGATTTCCAGa                                    <  1:299091/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAC‑CCACGGCAGCGTCGATTTCCAGACCCCGGc                             >  1:515688/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        tCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACTCCGAACCCGACGCCGGCTCAGGACAGGATCACCGACTCG‑GACGACACGAGTACAGAGACCCTGAGC‑GCCCCTCCAATCGAGAACGAAAGCGGCGAGATCCCCTACACCGTGACCCTGAGCAATCCTACCCAGACGCCTGTAACCGTCACCTTGTCCACCGGCTCAACCCGTAACGTTGAAG‑GGGTTAAACGCCAGGACAGctccgctttcacaccccctgacaactgcgc                    >  2:488879/1‑272 (MQ=12)
                                                                                                                                                                                                                                                                                                 acGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAC‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTAcaacaa           >  1:530450/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                    ggtCGACTGTCAGGGTCACCA‑TCGAGAACGCGACAGGCGGTAAGTTC‑GAGCAACTGAACCAGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGGCCACCGGCAACACCGTGACCGTTGCAGCGGGTAAACA‑CCACGCCAGTGTCAATTTCCAGACCCCGGCCAATGACGTCTACAACAACgg        <  2:242055/299‑1 (MQ=18)
                                                                                                                                                                                                                                                                                                         aCTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCTAGACGCCGGTAACCGTGACCTTGTCCAACGGCAAAACCATTACCGTTGAAGCGGGTAAAAC‑CCAGGGCAGCGTCGATTTCCAAACCCCGGCCAATGACGCCCACAACCACgcgccac   >  2:619629/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                                                          cTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGCGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTTACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAC‑CCAGGCCAACGTCCACTTCCAGAACCCGGcccatcaccgctagaaccaacgctcgcc  >  2:179141/1‑275 (MQ=18)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                                
CAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCT‑ACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATG‑ACGTCTACAACAACGGTTCGACTGTCAGCGTCACCA‑TCGAGAACGCCACAGGCGGTAATTTC‑GAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCG‑GTCGACACCACTACTGCGACCCTGACC‑GCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAAC‑CCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAC  >  AE015451/197335‑197931

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