Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,621,9790TA33.9% 97.8 / 19.4 62T107S (ACT→TCT) pdeFcyclic‑guanylate‑specific phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base T (30/11);  new base A (0/21);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.78e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATT  >  NZ_CP009273/2621733‑2622226
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
cccGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTaa                                                                                                                                                                                                                                                        <  1:487514/249‑1 (MQ=255)
    aaTGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc                                                                                                                                                                                                                                                   <  2:726295/250‑1 (MQ=255)
            aGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAAcc                                                                                                                                                                                                                                          >  2:153574/1‑250 (MQ=255)
                gCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                                                                                                                                      <  2:488357/250‑1 (MQ=255)
                gCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGcc                                                                                                                                                                                                                                       <  1:387457/249‑1 (MQ=255)
                  acaGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGcccc                                                                                                                                                                                                                                     >  2:680288/1‑249 (MQ=255)
                   caGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                                                                                                                                    >  2:219574/1‑249 (MQ=255)
                               gATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCACCAACAGAGaatgat                                                                                                                                                                                                                        >  2:555211/1‑249 (MQ=255)
                                     aTTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataatgat                                                                                                                                                                                                                  >  2:622949/1‑249 (MQ=255)
                                      ttACCCACCAGTAAGGCCTGATAATTAATTAAAGTCCCGAGGTAAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAAGGCAGAAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGAAAATTAAACGGGTACCACCATGTGAACCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa                                                                                                                                                                                                                <  1:134232/250‑1 (MQ=255)
                                         cccACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGataata                                                                                                                                                                                                              >  2:739285/1‑249 (MQ=255)
                                                 tAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                                                                                                                                     >  1:693143/1‑250 (MQ=255)
                                                   aGGCATGATAATTAATTAAGGCCCGGAATAAAAAAGCCATGACACCAACCAGATGTTCTCTGCCCGCCAGTAACCCAACAAAGGAAGCAAAAGGAAAAAGAATCAGAAACTGCGTAGAAGGAAACAAAATCTACAGGAATTTACGCGGGGAAATAAAACGGAAACCGCCATATAAAACGTTGTTACGCCGGGGATTAAAACCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                                                                                                                   <  2:518583/250‑1 (MQ=255)
                                                       ctaaTAATTGATTAAAGCCCAGAGGTTAAAAGGCATGACACCGCCCAGATTTTCTCTGCTCGCAAGTAATCCTAAAAATGCAGAAAACTGAAAAAGTATCAGAACCAGCGTTGCAGAAAACACAATCTGCAAAAATATACGCTGGGAAATTAAACGGGTACCACCATGTGAAACGTTGTTCCGCGGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAAt                                                                                                                                                                                               <  1:137477/247‑1 (MQ=255)
                                                          aTAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGa                                                                                                                                                                                             >  1:394163/1‑249 (MQ=255)
                                                          aTAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGa                                                                                                                                                                                             >  1:540105/1‑249 (MQ=255)
                                                           tAATTGATTAAAGTCCCGAGGATAAAAGGCATGACACCGACCAGATTTTCTCTACTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACTCTGTAACCGACCCAATGGAGAATGGTCATGATAGTAAACTTAGCCGTAACTGa                                                                                                                                                                                             >  2:192385/1‑248 (MQ=255)
                                                             tttGATAAAATCCACGAGTTAAAAGGGCATGACACCGACAAGTGTTCCTCTGCACGCAAGTAATCCAAAAAAGGAAAAAAATGGAAAAAGAAACAAAAAAAGCGTTGCAGGAAAAAAAATCGGCAAAAAAAAACGCGGGAAAAAAAAACGGGTATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCTGGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                                                                                                                                         <  2:630941/249‑1 (MQ=255)
                                                              ttGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                                                                                                                                        >  2:170668/1‑250 (MQ=255)
                                                                 aTTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACCGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                                                                                                                     >  2:336759/1‑250 (MQ=255)
                                                                      actccgtagttaaaaaGGCAGCACACCGACAAGATTTTCTCGGTCCGCCAGAAATGCTACAAATGAAGAAAACTGGAAAAGTATCAGAAACATGGTGGCAGGAAACACAATCGGAAATAAAATGCGCTGGGAAATACAACGCTTACCCCCAGGTAAAACTTTTTAACGCGGGGGATAAAAGCCTTATTATCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                                                                                                                <  2:710479/238‑1 (MQ=255)
                                                                                  aaaaGGCATGACCCCGACCAGATTTTCTCTGCCCCCCAGTAATCCTAAAAAGGAAGCAAACTGAAAAAGTATCAAAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATCCGCTGGAAAATAAAACGGGTATCCCCATGTGAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                                                                                                                                    <  1:112534/250‑1 (MQ=255)
                                                                                  aaaaGACATGACCCCCACCATATTTCCTTCGCTCCCAAGTAATCCTAAAAAGCAAGAAAACTGAAAAATATAAAAAAACGGGTTGGCAGAAAACACAATCTCCAAGAATATACGGTGGGAAATAACCCGGTTATCCCCATTTAAAATTTTTTTCCGCGGGGAAAAAAAACCCCTATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                                                                                                                                    <  2:608972/250‑1 (MQ=255)
                                                                                     aGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                                                                                                                                 >  2:334396/1‑250 (MQ=255)
                                                                                     aGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                                                                                                                                 >  2:228033/1‑250 (MQ=255)
                                                                                     aGGCATAACCCCGACCATATTTTCTCGGCTCGCCAGTAATCCTAAAAATGAAGAAACCGGGAAAAGTATCAAAAACAGCTTTGCAGAAAACAAAATCTGCCAGAATATACGCGGGGAAATAAAACGGGTATCCCCATGTAAAACGTTGTTACGCCGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                                                                                                                                 <  1:114625/250‑1 (MQ=255)
                                                                                            acacCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                                                                                                                                          >  1:403867/1‑250 (MQ=255)
                                                                                             cacCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                                                                                                                                                         >  2:586893/1‑250 (MQ=255)
                                                                                                 gACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGAc                                                                                                                                                      >  2:646676/1‑249 (MQ=255)
                                                                                                  aCCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACa                                                                                                                                                     >  2:89943/1‑249 (MQ=255)
                                                                                                           tccttGCTCGACAGAAATCCAACAAACAAAACAAACTGAAAAATTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGGGGGAAAATAAAACGGGTATCCCCATGTAAAACTTTGTTACGCCGGGGAGAAAACACCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                                                                                                                            <  1:105755/245‑1 (MQ=255)
                                                                                                               ggCTCGCCAGTAACCCAAAAAATGAAGAAAACGGAAAAAGTATCAGAACCAGCGTTGCAGAAAAACAATTCGGCAAAAATATACGCTGGAAAATAAAACGGGTATCACAATGTAAAACGTTGTTCCGCCGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTaa                                                                                                                                          <  1:604106/246‑1 (MQ=255)
                                                                                                                                    aTGCAGAAAACTAAAAAAGTAACAGAAACAGCTTGGCAGAAAACACAATCTGCGAGAATATACGGTGGAAAATAAAACGGGTACCACCAGGTAAAATTTTGTTCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCa                                                                                                                  <  2:624742/250‑1 (MQ=255)
                                                                                                                                            aacggaaaaaGTATAAGAAACAGCGTTGCAGAAAACACAACCTGCAAGAATATACGCGGGAAAATAAACCGGTTCCCCCCATGTGAAATTTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaa                                                                                                          <  2:307295/244‑1 (MQ=255)
                                                                                                                                              cgtgaaatagaacAGAAACAACGTTAGAAGGAACCCAATCCGGCAAGAATATACGTGGGAAAATAAACGGGGTACCTCATGGTAAAACTTGTTTCCGCGGGAAATAAAGCCCCTTGAATCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa                                                                                                         <  2:157001/238‑1 (MQ=255)
                                                                                                                                                              aaCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGatca                                                                                        >  1:624046/1‑250 (MQ=255)
                                                                                                                                                                  gCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAg                                                                                    >  2:618912/1‑250 (MQ=255)
                                                                                                                                                                    gTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTg                                                                                  >  2:345705/1‑250 (MQ=255)
                                                                                                                                                                             aaaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAgg                                                                           >  2:473937/1‑249 (MQ=255)
                                                                                                                                                                                cacaATCTGCCAGAATATCCGCTGGGAAATTAAACGGGTATCACCATGTAAAACTTTGTTACGCGGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGa                                                                       <  1:205897/250‑1 (MQ=255)
                                                                                                                                                                                 acaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGa                                                                       >  2:174193/1‑249 (MQ=255)
                                                                                                                                                                                  caAGCTGCCAAAATAACCGAGGGAAAATAAAACGGGTACCACCAGGTGAAACGTTTTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATa                                                                     <  2:117855/250‑1 (MQ=255)
                                                                                                                                                                                              gAATATACGGAGGAAAATAAAGCGGGTACCCCCATGGAAAACGTTGTTACGCGGGGAAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaaa                                                           <  1:406151/250‑1 (MQ=255)
                                                                                                                                                                                              gAATATACGCGGGAAAATTAACCGGGTACCACCATGTGAAACTTTGTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGaaaa                                                           <  1:716938/250‑1 (MQ=255)
                                                                                                                                                                                                    aCGCTGGAAAATTAAACGGTTATCCCCATGAAAAACTTTGTTACGCGGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCa                                                     <  1:85758/250‑1 (MQ=255)
                                                                                                                                                                                                      gCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCAcc                                                   >  1:285891/1‑250 (MQ=255)
                                                                                                                                                                                                        tGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCg                                                  >  1:683333/1‑249 (MQ=255)
                                                                                                                                                                                                         gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCgg                                                 >  1:154942/1‑249 (MQ=255)
                                                                                                                                                                                                          ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCgg                                                 >  2:89033/1‑248 (MQ=255)
                                                                                                                                                                                                          ggAAATAAAACGGGAATCCCCAGTGAAAATGTTGTTACGCCGGGAATAAAAGACCCGGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaa                                               <  1:197910/250‑1 (MQ=255)
                                                                                                                                                                                                           aaaaTAAACCGGGTATCCCCAGGTGAAACGTTGTTACGCCGGGGAAAAAAACCCCTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaag                                              <  2:671945/249‑1 (MQ=255)
                                                                                                                                                                                                            aaaTTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaag                                              >  2:530262/1‑249 (MQ=255)
                                                                                                                                                                                                            aaaTAAAAGGGTTACCACCAGGGAAACGGTGGTGACGCGGGGGAGAAAAGACCTGGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGgaaga                                             <  1:229179/250‑1 (MQ=255)
                                                                                                                                                                                                                aaagcGGGTAACCCCATGAGAAACGTTGTTACGCGGGGAAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAata                                         <  2:464881/246‑1 (MQ=255)
                                                                                                                                                                                                                 aaaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAatat                                        >  2:622324/1‑250 (MQ=255)
                                                                                                                                                                                                                   aCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATgc                                      <  1:345705/250‑1 (MQ=255)
                                                                                                                                                                                                                     gggTACCCCCATGAGAAACGTTGTTACGCCGGGGAATAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATgcgc                                    <  1:438871/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggTATCACCATGGAAAACGTTGTAACGCGGGGGAGAAAAGACCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGcgcg                                   <  2:615642/250‑1 (MQ=255)
                                                                                                                                                                                                                             ccATGAGAAAGTTTTTTACGCGGGGGAAAAAAATCCCCTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTAt                            <  2:342976/250‑1 (MQ=255)
                                                                                                                                                                                                                                     aaaCGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAAtgt                      <  1:176567/249‑1 (MQ=255)
                                                                                                                                                                                                                                          tttgttACGCCGGGGATAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCg                 <  2:698680/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          gTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGAtt  >  1:738066/1‑249 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
CCCGAATGATGAAGTAGCACAGCGGGACACCGATCAGATTACCCACCAGTAAGGCCTGATAATTGATTAAAGTCCCGAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATCATGCTGATCATCAGTGCCAGAGGCAGATAAAAGAGAAAAACCTCACCGGAAGAAATATGCGCGAAAGTATTAATGTGGGCGAAAATGGGCAAGATT  >  NZ_CP009273/2621733‑2622226

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: