Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,738,3790CG58.8% ‑0.5 / 11.8 17G72G (GGC→GGGyiaN2,3‑diketo‑L‑gulonate transporter large permease YiaN
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base G (5/5);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGC  >  NZ_CP009273/3738293‑3738462
                                                                                      |                                                                                   
gCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGggcggc                                                                                  <  2:36926/90‑1 (MQ=255)
   caAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTg                                                                               >  2:104926/1‑90 (MQ=255)
     aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCATGAATGCGGGGGGCTTGTc                                                                             >  2:7566/1‑90 (MQ=255)
            gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGc                                                                      <  2:58140/90‑1 (MQ=255)
            gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGc                                                                      <  2:397615/90‑1 (MQ=255)
                   gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGGAAACATGAATGAGGGGGGGTTTTCAAAACGCCttgttg                                                               >  1:78237/1‑90 (MQ=255)
                     gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTGTCAAAAAGCATTGTTGAc                                                             >  2:69646/1‑90 (MQ=255)
                     gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGCGGGGGGCTTGTCTAAACGCCttgttggc                                                             >  1:393894/1‑88 (MQ=255)
                                       cgggCGATTCCGTTCTTTGTGCGGGGGGGTAAAATCATAAATGGGGGGGGCTTGAAAAAACGCATTGTTGCCCTGCCGATGAAACtggtg                                           <  1:246826/88‑1 (MQ=255)
                                        tgGCGATTCCTTTTTTTGTGCGGGGGGGAGAAATCATGAATGCGGGGGGTTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACtggtgg                                          <  1:383821/90‑1 (MQ=255)
                                         ggCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATGCGGGCGGGTTGTCCAAACGCATTGTTGACCTGCCGATGGAACTGGTggg                                         >  1:68992/1‑90 (MQ=255)
                                         ggCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGGGGGCGGCTTGTCAAAAAGCATTGTTGAGCTGCCGGTGAAACTGGTggg                                         >  2:246826/1‑90 (MQ=255)
                                                     cTTTGTGCTGGGGGGTGAAATCATGAATGGGGGGGGCTTGCAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAgg                             <  2:402510/90‑1 (MQ=255)
                                                       ttGTGCTGGGGGGTAAAATCATGAATGCGGGGGGTTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCg                           <  1:130545/90‑1 (MQ=255)
                                                          tgCTGGGGGGTGAAATCATGAATGGGGGGGGTTTGTAAAAACGCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                        <  1:104926/90‑1 (MQ=255)
                                                                            tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGc      >  2:486861/1‑90 (MQ=255)
                                                                            tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGc      >  1:307098/1‑90 (MQ=255)
                                                                                tGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGc  <  2:166475/90‑1 (MQ=255)
                                                                                      |                                                                                   
GCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGC  >  NZ_CP009273/3738293‑3738462

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