Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,571,6610AC50.0% 112.6 / 14.9 76G142G (GGT→GGGgntRgluconate operon transcriptional repressor GntR
Reads supporting (aligned to +/- strand):  ref base A (37/1);  new base C (0/38);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.13e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAT  >  NZ_CP009273/3571579‑3571748
                                                                                    |                                                                                       
ggCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCa                                                                                  >  2:793611/1‑90 (MQ=255)
  cGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACt                                                                                >  2:145337/1‑90 (MQ=255)
  cGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACt                                                                                >  2:629515/1‑90 (MQ=255)
  cGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACt                                                                                >  2:2527139/1‑90 (MQ=255)
   gTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACtt                                                                               >  1:124238/1‑90 (MQ=255)
   gTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACtt                                                                               >  2:617543/1‑90 (MQ=255)
   gTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACtt                                                                               >  2:2037210/1‑90 (MQ=255)
    tgctgcTTAAAATTTATAAAACCCGCCGGCAATACCAAGGCATGGCGTTTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTc                                                                              <  1:537682/90‑1 (MQ=255)
      ctgctAAAAATTTAAAAAACCCGACGGCGATATAAAGGCTTGGCGCCTTGCTGCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCaa                                                                            <  2:211776/90‑1 (MQ=255)
        gctTCAAATTTAAAAAACCGGGCGGAAAAATCAAGGAAGGGCCCTTTGTTGTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAatc                                                                          <  1:2505382/90‑1 (MQ=255)
        gctTAAAATTTATAAAACCCGGGGGCGATATAAAGGCAGGGCGATTTGCTGCCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAatc                                                                          <  2:348275/90‑1 (MQ=255)
        gctTAAAATTTATAAAACCCGACGGCGATACAAAGGCAGGGCGACTTGCTGCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAatc                                                                          <  2:1715200/90‑1 (MQ=255)
         gctaaaaaTTAACAAAACCCGCGGGCAATACAAAGGCAGGGCGTTTTGCTTTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAatc                                                                          <  2:2573614/86‑1 (MQ=255)
          ttCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat                                                                        >  2:33854/1‑90 (MQ=255)
          ttCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat                                                                        >  1:787471/1‑90 (MQ=255)
          ttCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat                                                                        >  1:1510955/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  2:2842343/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  2:1136756/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  1:2183480/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  1:2505461/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  1:422108/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  2:501031/1‑90 (MQ=255)
                  tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg                                                                >  1:660257/1‑90 (MQ=255)
                  aaaaaaaaCCGACGGCAAAAAAAAGGCATGGCGCTTTGCTGCCCATCATTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAgg                                                                <  2:2008524/86‑1 (MQ=255)
                        aacccgccgGCGATAACAGGGCAGGGCGTTTTGCTCCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTgcgc                                                           <  2:2253861/87‑1 (MQ=255)
                        aacccgccgGCAATACAAAGCCAGGGCGCCTTCCTGCCCATCAGCCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTgcgc                                                           <  1:1922804/87‑1 (MQ=255)
                        aacccGGCGGCGATATCAGGGCGGGGCGTTTTGCTCTCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTgcgc                                                           <  1:2819798/87‑1 (MQ=255)
                        aacccGACGGCGATATAAAGGCAGGGCGACTTGCTGCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTgcgc                                                           <  2:887578/87‑1 (MQ=255)
                        aacccGACGGCAATACAAAGGCAGGGCGACTTGCTGTCCATCATTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTgcgc                                                           <  2:1020160/87‑1 (MQ=255)
                         acccGACGGCGATATCAAGGCAGGGCGCCTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg                                                          <  2:1108571/88‑1 (MQ=255)
                          cccgccgGCAATAAAAAGGCATGGCGTCTTGCTGCCCACCAGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                         <  2:294431/89‑1 (MQ=255)
                          cccgccgGCAATAAAAAGGCAGGGCGTCTTGCTCCCCATCATTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                         <  1:586574/89‑1 (MQ=255)
                          cccGGCGGCAATACCAAGGCATGGCGTTTTGCTGCCCACCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                         <  1:458381/89‑1 (MQ=255)
                          cccGGCGGCAAAATCAGGGCAGGGCGCTTTGCTGTCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                         <  1:408052/89‑1 (MQ=255)
                           ccGACGGCGATATCAAGGCGTGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCgcggc                                                        >  1:2432431/1‑90 (MQ=255)
                            cggcGGCGATATAAAGGCAGGGCGTCTTGCTGCCCATCAGTCCCACCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCgcggcg                                                       <  2:2432431/87‑1 (MQ=255)
                              aCGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg                                                     >  2:1026409/1‑90 (MQ=255)
                               cGGCAATATCAGGGCATGGCGTTTTGTTCTCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtgt                                                    <  1:721481/90‑1 (MQ=255)
                                   aatatCAGGGCGTGGCGCCTTGCTTTCCATCAGTCCCCCCGGGGAAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGt                                                <  1:2899258/89‑1 (MQ=255)
                                    atatCAGGGCAGGGCGCTTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTa                                               <  1:1157766/90‑1 (MQ=255)
                                    atatAAAGGCATGGCGTCTTGCTGCCCATCAGTCCCCCCGCGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTa                                               <  2:1049171/90‑1 (MQ=255)
                                       acAAGGCATGGCGCCTTGCTGCCCATCATTCCCCCCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGt                                            <  2:124238/89‑1 (MQ=255)
                                         agggCAGGGCGTCTTGCTGCCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                          <  2:2578839/88‑1 (MQ=255)
                                         agggCAGGGCGCCTTGCTGTCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                          <  2:1912915/88‑1 (MQ=255)
                                         aaGGCGGGGCGCTTTGCTGCCCACCGGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                          <  2:1966715/90‑1 (MQ=255)
                                         aaGGCGGGGCGATTTGTTGTCCCTCATTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                          <  2:1717702/90‑1 (MQ=255)
                                         aaGCCAGGGCGCCTTTCTGTCCATCACTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTc                                          <  1:1811474/90‑1 (MQ=255)
                                             cATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                      >  2:55120/1‑90 (MQ=255)
                                             cATGGCGACTTGCTGCCCATCATTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                      <  1:119326/90‑1 (MQ=255)
                                               tGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                    >  2:2034618/1‑90 (MQ=255)
                                               tGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                    >  2:2027899/1‑90 (MQ=255)
                                               tGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                    >  1:1947241/1‑90 (MQ=255)
                                               tGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg                                    >  2:2732643/1‑90 (MQ=255)
                                                     cTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAg                              >  2:2263485/1‑90 (MQ=255)
                                                      ttGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAgg                             >  2:2096007/1‑90 (MQ=255)
                                                        cctttCCACCAGCCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGcc                           <  1:2203362/86‑1 (MQ=255)
                                                            tCCATCATTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                       <  1:1136756/90‑1 (MQ=255)
                                                            cccATCATTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCg                       <  2:112993/89‑1 (MQ=255)
                                                              catcatTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                     <  2:881007/90‑1 (MQ=255)
                                                              catcatTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                     <  2:489434/90‑1 (MQ=255)
                                                              catcatTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                     <  1:33854/90‑1 (MQ=255)
                                                              cATCAGCCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat                     <  1:250038/90‑1 (MQ=255)
                                                                tCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTGCGTTCGGTGAGGATCAGGCCGTCGatat                   >  1:1903960/1‑90 (MQ=255)
                                                                   gTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTcc                >  1:1645790/1‑90 (MQ=255)
                                                                     tCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg              >  1:2429395/1‑90 (MQ=255)
                                                                     tCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg              >  2:1612819/1‑90 (MQ=255)
                                                                     tCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAg              >  2:2359251/1‑90 (MQ=255)
                                                                      cccccagGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAgg             <  2:943146/87‑1 (MQ=255)
                                                                           acGGGAATACGCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGTTACGTTCGTTGAGGATCAGGCGATCGATATTCCAGGagaga        >  2:2221688/1‑89 (MQ=255)
                                                                             gggAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtgt                                                    >  2:739634/1‑44 (MQ=255)
                                                                             gggAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCAt      >  2:2415301/1‑90 (MQ=255)
                                                                             gggAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCAt      >  1:40720/1‑90 (MQ=255)
                                                                               gAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTGCGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATgg    <  2:1903960/90‑1 (MQ=255)
                                                                                  aTACCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                         <  2:778393/35‑1 (MQ=255)
                                                                                  aTACCCGCCACTTCAATCATCTTTAAGGTGCgcgg                                                         >  1:778393/1‑35 (MQ=255)
                                                                                  aTACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAt  >  1:1630069/1‑90 (MQ=255)
                                                                                    aCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                      >  1:2060424/1‑52 (MQ=255)
                                                                                    aCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTg                                      <  2:2060424/52‑1 (MQ=255)
                                                                                    |                                                                                       
GGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGAT  >  NZ_CP009273/3571579‑3571748

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