Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,574,8330TG55.0% 23.9 / 27.0 20G613G (GGT→GGGtkttransketolase
Reads supporting (aligned to +/- strand):  ref base T (0/9);  new base G (11/0);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 5.95e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTA  >  NZ_CP009273/2574746‑2574895
                                                                                       |                                                              
gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGAt                                                              >  2:9954/1‑90 (MQ=255)
gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                                              <  2:93371/90‑1 (MQ=255)
    gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGcc                                                          >  2:295506/1‑90 (MQ=255)
    gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGAAGCAGGGATTGcc                                                          >  2:172873/1‑90 (MQ=255)
    gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGAAGCAGGGATTGcc                                                          >  2:332588/1‑90 (MQ=255)
           tCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACt                                                   <  1:332588/90‑1 (MQ=255)
           tCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAATGCCGATTACt                                                   <  2:382328/90‑1 (MQ=255)
                gCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACa                                             <  2:379336/90‑1 (MQ=255)
                 cccAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGGGGGGATTGCCGATTATTGGTACa                                             >  1:323478/1‑90 (MQ=255)
                         gAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt                                     <  1:349814/90‑1 (MQ=255)
                          aGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGGACAAGTATGTTg                                    >  1:110204/1‑90 (MQ=255)
                                  cGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa                            <  1:372741/90‑1 (MQ=255)
                                         cGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGATTACTGGGACAAAGATGTTGGGCTGGAAGGGGGaa                     >  1:355668/1‑90 (MQ=255)
                                              ttGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTTCCGGTTACTGGGACAAGTATGTTGGGCTGGAAGGGGGAAATGTc                >  1:448218/1‑90 (MQ=255)
                                               tGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCg               <  2:366910/90‑1 (MQ=255)
                                                     cTAACGTTGCGGCTCGCGTGGCGGTGGAAGCCGGGATTGCCGATTTCTGGTACAAGGATGTTGGGTTGTAAGGGGCAATTGTCGGGGTgg         >  1:443488/1‑89 (MQ=255)
                                                     cTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa         >  2:278225/1‑90 (MQ=255)
                                                      tAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGAc        <  1:361183/90‑1 (MQ=255)
                                                           ttGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGtt   <  2:311055/90‑1 (MQ=255)
                                                            tGCGGCTCGCGTGGCGGTGGAAGGAAGGATTGCCGGTTACTTGGAAAAGTATGTTGGGCTGAAAGGGGGAATTGTCGGGGTGACGGGGTa  >  1:349977/1‑90 (MQ=255)
                                                                                       |                                                              
GACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTA  >  NZ_CP009273/2574746‑2574895

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