Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,832,7180TG52.5% 68.4 / 77.2 61intergenic (‑95/‑165)ascG/ascFDNA‑binding transcriptional regulator AscG/PTS cellobiose/arbutin/salicin transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base T (13/16);  new base G (32/0);  total (45/16)
Fisher's exact test for biased strand distribution p-value = 3.35e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-07
Rejected as consensus: Frequency below/above cutoff threshold.

TTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATT  >  NZ_CP009273/2832636‑2832803
                                                                                  |                                                                                     
tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACaa                                                                                <  1:1622876/90‑1 (MQ=255)
tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTTACaa                                                                                >  2:1902126/1‑90 (MQ=255)
        cAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTTAGGAGACCCGGTTCACAAATATaaaa                                                                        >  1:1172901/1‑90 (MQ=255)
              aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGaa                                                                  <  2:1377208/90‑1 (MQ=255)
              aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTTAGGGGACCGGGTTCACAAATATAAAAAATGaa                                                                  >  1:172932/1‑90 (MQ=255)
               cGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                 <  1:1035019/90‑1 (MQ=255)
                  aGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGGGACCGGGTTCACAAATATAAAAAATGAACAAt                                                              >  2:1532612/1‑90 (MQ=255)
                  aGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAAt                                                              >  2:1838056/1‑90 (MQ=255)
                  aGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTTAGGTGACCGGGGTCACAAATATAAAAAATGAACAAt                                                              >  1:1176459/1‑90 (MQ=255)
                  aGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGACCGGTGTCACAAATATAAAAAAAGAACAAt                                                              >  1:1238332/1‑90 (MQ=255)
                   gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAtt                                                             <  1:1237116/90‑1 (MQ=255)
                    tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGTGACCGGGTTCACAAATATAAAAAATGAAGAATTc                                                            >  2:106036/1‑90 (MQ=255)
                    tgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGTTTCACAAATATAAAAAATGAACAATTc                                                            >  2:1664362/1‑90 (MQ=255)
                      tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                          <  1:1454280/90‑1 (MQ=255)
                      tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                          >  2:249358/1‑90 (MQ=255)
                        cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActc                                                        <  2:1330838/90‑1 (MQ=255)
                        cTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCCGGGGACCGGGTGCACAAAAATAAAAAATGGACAAATCActc                                                        >  1:1250450/1‑90 (MQ=255)
                         tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAATGAACAATTCActct                                                       >  2:1588091/1‑90 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTTAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTTActctc                                                      >  1:1970737/1‑90 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctc                                                      >  2:489317/1‑90 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActctc                                                      <  2:2015327/90‑1 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCActctc                                                      >  2:690951/1‑90 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCActctc                                                      >  2:458821/1‑90 (MQ=255)
                          gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAAAGAACAATTCActctc                                                      >  1:527280/1‑90 (MQ=255)
                          gCACCGTTTGCCAATTGTCCTTGCACAATCGGCGGGGAAAAAATTTAGGGGGCCGGGTTCACAAAAATAAAAAAAGAACAAATTActctc                                                      >  1:1868313/1‑90 (MQ=255)
                             ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTCAGGTGACCGGGTTTACAAAAATAAAAAAAGAAAAATTCACTCTCTTg                                                   >  1:1633992/1‑90 (MQ=255)
                             ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTg                                                   >  2:179590/1‑90 (MQ=255)
                             ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCttt                                                   >  2:1594905/1‑89 (MQ=255)
                             ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTGTGGGGACCGGGGGCCCAAATATAAAAAATGAAGCATTCACTCTCTTg                                                   >  1:909034/1‑90 (MQ=255)
                             ccGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGTGACCGGGGTCACAAATATAAAAAATGAACAATTCACTCTCTTg                                                   >  1:1682314/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTTAGGGGACCGGGTTCACAAAAATAAAAAATGAACAATTCACTCTCTTGCt                                                 >  1:1928227/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTCAGGTGGCCGGTTTCACAAATATAAAAAAAGAACAATTCAATCTCTTggt                                                 >  2:1859870/1‑88 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGGAAAAAATTCAGGTGACCGGGGTCACAAATATAAAAAAAGAACAATTCACTCTCTTGCt                                                 >  1:1769328/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTTAGGGGACCGGTTTCACAAATATAAAAAAAGAACAATTTACTCTCTTGCt                                                 >  2:1188005/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                          >  1:644087/1‑81 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                 >  1:758283/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                 >  2:1293297/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                 >  2:449524/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                 >  2:67134/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATAATTAGGGGGCCGGGGTCCCAAATAAAAAAAAATTACAAATCACTTTCTTGCt                                                 >  1:1273466/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTTAGGTGACCGGGTTCACAAATATAAAAAAAGAACAATTCACTCTCTTTCt                                                 >  1:1594251/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGGCCGGGTTCACAAATATAAAAAAAGAACAAATCACTCTCTTGCt                                                 >  1:1506455/1‑90 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAAATTCAGGGGACCGGGTTCACAAAAATAAAAAAAAAACAATTCACTCTCTTggt                                                 >  2:882488/1‑88 (MQ=255)
                               gTTTGCCGATTGTCCTTGCACAATCCGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCAc                                                          <  2:644087/81‑1 (MQ=255)
                               gTTTGCCGATTGTCCTTGCAAAATCGGCGGGAAAAATATTCAGGTGGCCGGGGCGGCAAATATAAAAAATTAACAATTTTCTCTCTTggt                                                 >  2:1205371/1‑88 (MQ=255)
                                  tGCCGATTGTCCTTGCACAATCGGCGGGGAAAATATTCAGGAGACCGGGTTCACAAAAATAAAAAATGAACAAATCAACCTCCTCCTTAt                                              >  2:1391674/1‑90 (MQ=255)
                                      gATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAg                                          <  1:1589068/90‑1 (MQ=255)
                                             cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                   <  1:306545/90‑1 (MQ=255)
                                             cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                   <  1:1902126/90‑1 (MQ=255)
                                              ttGCACAATCGGCGGGAAAAATATTCAGGGGACCGGGTTCACAAATATAAAAAAAGAACAATTCACTCTCTTGCTTTTTTTGTGACAACt                                  >  1:1122131/1‑90 (MQ=255)
                                               tGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTa                                 <  2:1633992/90‑1 (MQ=255)
                                                    aaTCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTg                                        <  2:1646953/78‑1 (MQ=255)
                                                    aaTCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTg                                        >  1:1646953/1‑78 (MQ=255)
                                                     aTCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATg                           >  2:824696/1‑90 (MQ=255)
                                                     aTCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATg                           >  2:1680544/1‑90 (MQ=255)
                                                     aTCGGCGGGAAAAATATTCAGGTGACCGGGTTCACAAATATAAAAAATGAAAAATTCACTCTCTTGCTTATTTTGTGACAAATATTCATg                           >  1:172373/1‑90 (MQ=255)
                                                     aTCGGCGGGAAAAATATTCAGGGGGCCGGGTTCACAAATATAAAAAAAGAAAAATTTACTTTCTTGCTTATTTAATGACAAATATTTATg                           >  1:1882621/1‑90 (MQ=255)
                                                                     ttCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCCCTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGt           <  2:2003535/90‑1 (MQ=255)
                                                                     ttCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGt           <  1:1293297/90‑1 (MQ=255)
                                                                              ccGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAAtt  >  2:463341/1‑90 (MQ=255)
                                                                              ccGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAAtt  <  2:172932/90‑1 (MQ=255)
                                                                                  |                                                                                     
TTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATT  >  NZ_CP009273/2832636‑2832803

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