Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,098,3370TG52.6% 19.4 / 24.4 19G286G (GGT→GGGpfkA6‑phosphofructokinase
Reads supporting (aligned to +/- strand):  ref base T (0/9);  new base G (10/0);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.08e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: Frequency below/above cutoff threshold.

TTACGACCGTATTCTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCTATCGAAAACATGAAGCGTCCGTTCAAAGGTGACTGGCTGGACTGCGC  >  NZ_CP009273/4098256‑4098426
                                                                                 |                                                                                         
ttACGACCGTATTCTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTTGGGATTCAgaa                                                                                   >  1:202925/1‑90 (MQ=255)
 tACGACCGTATTCTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGGATCCAgaac                                                                                  >  1:214743/1‑90 (MQ=255)
        gTATTCTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGGCGTTGTGTAGGGATCCAGAAAGAACAGc                                                                           >  2:367105/1‑90 (MQ=255)
         tATTCTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCt                                                                          <  1:353053/90‑1 (MQ=255)
             cTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGtt                                                                      <  2:584246/90‑1 (MQ=255)
               ggCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGGCGGTTTGTTGGGATCCACAAAGAAAAGCTGGGTcc                                                                    >  2:122398/1‑89 (MQ=255)
                 cTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTATGTCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTcacc                                                                  <  1:194332/90‑1 (MQ=255)
                   tCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTcaccac                                                                <  2:477933/90‑1 (MQ=255)
                         aTGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGGAGGGATCCAGAAAGAACAGCTGGGTCACCACGAcatc                                                          >  1:395641/1‑90 (MQ=255)
                            ggCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACatcatc                                                       <  1:208933/90‑1 (MQ=255)
                               gcTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGGTGTGGAGGGATCCACAACGAACAGCTGGGTCACCACCACATCATCGAc                                                    >  2:355955/1‑90 (MQ=255)
                                  tACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCt                                                 <  1:514286/90‑1 (MQ=255)
                                  tACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGGGGTTGTTTTGGGATCCACAAAGGAAACCCGGTTTACCACGACATCCTCGACGCt                                                 >  2:448944/1‑90 (MQ=255)
                                     gCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGGATCCAGAACGAACAGCTGGTTCCCCACCACATCATCGGCGCTGTc                                              >  1:145312/1‑90 (MQ=255)
                                       tATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCTATCGa                                            <  1:448944/90‑1 (MQ=255)
                                         tCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCTATCGaaa                                          <  1:34938/90‑1 (MQ=255)
                                                                 gcggTCGTTGTGTAGGGATCCAGAAAGAACAGCTGGGTCACCCCGACATCATCGACGCTATTGAAAAAATGAAGCGGCCGGTCAAAgggg                  >  1:385097/1‑88 (MQ=255)
                                                                    gTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCTATCGAAAACATGAAGCGTCCGTTCAAAGGTGAct               <  2:214743/90‑1 (MQ=255)
                                                                     tCGTTGTGTAGGGATCCAGAACGAACAACTGGGTTACCCCCACATCATCGACGCTATCGAAAAAAAGAAGCGGCCGTTCAAAAGGGActg              >  2:71317/1‑90 (MQ=255)
                                                                                 tATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCTATCGAAAACATGAAGCGTCCGTTCAAAGGTGACTGGCTGGACTgcgc  <  2:395641/90‑1 (MQ=255)
                                                                                 |                                                                                         
TTACGACCGTATTCTGGCTTCCCGTATGGGCGCTTACGCTATCGATCTGCTGCTGGCAGGTTACGGCGGTCGTTGTGTAGGTATCCAGAACGAACAGCTGGTTCACCACGACATCATCGACGCTATCGAAAACATGAAGCGTCCGTTCAAAGGTGACTGGCTGGACTGCGC  >  NZ_CP009273/4098256‑4098426

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