Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,154,0680TG52.0% 12.1 / 21.3 26G248G (GGT→GGGptsGPTS glucose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base G (13/0);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 2.02e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-04
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAG  >  NZ_CP009273/1153982‑1154147
                                                                                      |                                                                               
aTGCAGATTGGTGATTACAGCAACGCAGGAGGGCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                                                              <  1:992764/90‑1 (MQ=255)
aTGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCCACTGCGGGTgaa                                                                              >  2:827274/1‑90 (MQ=38)
  gCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACt                                                                            <  1:629778/90‑1 (MQ=255)
  gCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGaaaat                                                                            >  2:92868/1‑88 (MQ=255)
              aTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGGTGGCtt                                                                >  1:1459036/1‑90 (MQ=255)
              aTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCCACTGCGGGAAAACTGTCTGGTGGCtt                                                                >  2:1379366/1‑90 (MQ=255)
              aTACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAACTGGCTGGGGGCtt                                                                >  1:776143/1‑90 (MQ=255)
               tACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGGGGGCTTc                                                               >  1:152948/1‑90 (MQ=255)
               tACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCCACCGCGGGGAAACTGTCTGGGGGCTTc                                                               >  2:591436/1‑90 (MQ=255)
                acacCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTcc                                                              <  2:1751589/90‑1 (MQ=255)
                 cacCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCt                                                             <  1:827274/90‑1 (MQ=255)
                   ccAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGGCCCCACTGCGGGGAAACTATCTGGTGGCTTccctt                                                           >  2:588519/1‑87 (MQ=255)
                      aCGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAAATGTCTGGTGGCTTCCTTTTCa                                                        >  2:343322/1‑90 (MQ=255)
                       cGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaa                                                       <  1:1707927/90‑1 (MQ=255)
                           gcagGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCCACTGCGGGGAAAATGTCTGGGGGGTTCCTGTTTAAAATg                                                   >  2:1707927/1‑90 (MQ=255)
                             aggtcaggtTTTCCACGGCGACATTACGCGTTATATGGCGGGGGACCCCACTGCGGGGAAACTGTCTTGGGGCTTCCCGTTCAAAATGTa                                                 >  1:1011870/1‑90 (MQ=255)
                              ggtcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGGGACCCGACTGCGGGGAAACTGTCTGGGGGGTTCCTGTTTAAAATGTaa                                                >  2:1509343/1‑89 (MQ=255)
                                tcaggtTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGGAAACTGTCTGGGGGCTTCCTGTTCAAAATGTACgg                                              >  1:1311828/1‑90 (MQ=255)
                                              gCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCg                                <  1:567158/90‑1 (MQ=255)
                                                                   cGGGTGACCCGACTGCGGGGAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTg           >  2:463203/1‑90 (MQ=255)
                                                                       tGACCCGACTGCGGGTAAACTGTCTGTTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa       <  1:656806/90‑1 (MQ=255)
                                                                       tGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa       <  2:67809/90‑1 (MQ=255)
                                                                       tGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa       <  1:21834/90‑1 (MQ=255)
                                                                       tGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaa       <  1:1639385/90‑1 (MQ=255)
                                                                        gACCCCACTGCGGGGAAAATGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTaaa      >  2:1045931/1‑90 (MQ=255)
                                                                            cGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAg  >  2:144156/1‑90 (MQ=255)
                                                                                      |                                                                               
ATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAG  >  NZ_CP009273/1153982‑1154147

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