Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,086,4940TG64.7% 0.9 / 12.5 17Y102S (TAC→TCC) rhaML‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (11/0);  total (13/4)
Fisher's exact test for biased strand distribution p-value = 6.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAA  >  NZ_CP009273/4086421‑4086583
                                                                         |                                                                                         
gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAActtttttccg                                                                           >  2:573718/1‑86 (MQ=255)
gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAg                                                                           >  2:548550/1‑90 (MQ=255)
 tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGc                                                                          >  1:490541/1‑90 (MQ=255)
 tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAAAATCTTTCAGc                                                                          >  2:454619/1‑90 (MQ=255)
        tttCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACCTTTTTCCGCTCGCCAc                                                                   >  2:152825/1‑90 (MQ=255)
                  aaaTTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCt                                                         <  1:481990/90‑1 (MQ=255)
                                 tgCCGGAGGCATCTCTGGCACATTGGGCAATTCCGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCCCCGGGCTGTTATCCGGGTTAGc                                          <  2:577484/90‑1 (MQ=255)
                                     ggAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAACACCTCTTTCCGGCCGCGACCCCCCCGGCCGTTATACGGGGTAGCGGGc                                      >  1:577484/1‑90 (MQ=255)
                                       aGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACAATTCTTTCAGCTCGCCCCTCACCCGGCCGTTTTTCGGGGTAGCGGGCAt                                    >  2:231087/1‑90 (MQ=255)
                                       aGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACCTCTTTCCGCTCGCCACTCCCCCGGCTGGTTTCCCGGGTAGCGGGGAt                                    >  2:167904/1‑90 (MQ=255)
                                            gctctGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACAt                               <  1:288032/89‑1 (MQ=255)
                                             ctctGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATc                              <  1:454619/90‑1 (MQ=255)
                                                    acaTTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGGCata                       >  1:102765/1‑90 (MQ=255)
                                                       ttGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCACCCGGCTGTTTTCCGGGGTAGCGGGGATAAAAACGGGCatatac                    >  2:288032/1‑89 (MQ=255)
                                                         gggCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCCCCACTCCCCCGGCTGTTTTTCGGGGTAGCGGGCATAACATCGGGCATATAttt                  >  2:608976/1‑90 (MQ=255)
                                                               ttACGGCAGGGAAAAAACTTCTTTCCGCTCGCCACTCACCCGGCTGTTTTCCGGGGTAGCGGGCATAACATCGGGCATATATTTtcacca            >  2:355569/1‑90 (MQ=255)
                                                                   ggCAGGGAAAAAACTTCTTTCAGCTCGCTACTCACCCGGCTGTTTTTCGGGGTAGCGGGCATAACATCGGGCATATATTTCCCCCAACGt        >  2:620222/1‑90 (MQ=255)
                                                                         tAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCaa  <  2:3628/90‑1 (MQ=255)
                                                                         tAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCaa  <  1:675005/90‑1 (MQ=255)
                                                                         tAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCaa  <  1:573718/90‑1 (MQ=255)
                                                                         |                                                                                         
GTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAA  >  NZ_CP009273/4086421‑4086583

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