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

TGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC  >  NZ_CP009273/3738291‑3738465
                                                                                        |                                                                                      
tGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATgcgggcg                                                                                       >  2:242091/1‑90 (MQ=255)
       aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAACATGAATGCGGGGGGCGTggc                                                                                >  2:617245/1‑88 (MQ=255)
         cGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCaa                                                                              <  2:221163/90‑1 (MQ=255)
           cTGGTGAACGGAGCCGATAGCTTCTCACTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCATGAATGCGGGGGGCTTGTCaaaa                                                                            >  1:306843/1‑90 (MQ=255)
              gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGc                                                                         <  1:631837/90‑1 (MQ=255)
                gAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAAAATGAATGGGGGGGGCTTGTCAAAACGCAt                                                                       >  1:233350/1‑90 (MQ=255)
                        ccGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGAcc                                                               <  1:242091/90‑1 (MQ=255)
                          aaTAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTTTCAAAAAGCATTTTTGACCtt                                                             >  1:526812/2‑89 (MQ=255)
                               cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGAt                                                        >  1:82790/1‑90 (MQ=255)
                               cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAAAATGAATGCGGGGGGCCTTTCCAAACGCATTGTTGGCCCGCCCAt                                                        >  2:734058/1‑90 (MQ=255)
                                ttCTCACTGCTGGCGATTCCGTTTTTTGTGCGGGGGGGTGAAATCATGAATGCGGGGGGTTTGTCAAAACGCATTGTTGCCCTGCCGATg                                                       <  2:306843/90‑1 (MQ=255)
                                 tctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGa                                                      >  1:636297/1‑90 (MQ=255)
                                    cccTGCTGGCGTTTCCTTTCTTTGTGCTGGGGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAAc                                                   <  2:297197/90‑1 (MQ=255)
                                     ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACt                                                  <  1:601701/90‑1 (MQ=255)
                                        gctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAAAGCATTGTTGACCTGCCGGTGAAACtggt                                               >  1:289045/1‑90 (MQ=255)
                                         ctgGCGTTTCCGTTTTTTGTGCTGGGGGGTGAAATCATGAAGGCGGGCGGTTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                              <  2:75485/90‑1 (MQ=255)
                                         ctgGCGATTCCTTTTTTTGTGCGGGGGGGGAAAATCATGAAGGGGGGGGGTTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACtggtg                                              <  2:546060/90‑1 (MQ=255)
                                          tgGCGTTTCCGTTTTTTGTGCTGGGGGGTGAAATCATGAATGCGGGCGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACtggtgg                                             <  2:263115/90‑1 (MQ=255)
                                              gATTCCTTTTTTTGTGCTGGGGGGTAAAATCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCa                                         <  2:233350/90‑1 (MQ=255)
                                                ttCCGTTCTTTGTGCTGGCGGGTGAAATCATGAAGGGGGGGGGTTTGTCAAACCGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATa                                       <  2:540667/90‑1 (MQ=255)
                                                            tgCTGGGGGGTGAAATCATGAAGGGGGGGGGTTTGTAAAAACGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                           <  2:803298/90‑1 (MQ=255)
                                                            tgCTGGGGGGTGAAAACAAAAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                           <  1:780724/90‑1 (MQ=255)
                                                            tgCTGGCGGGTGAAATCATGAATGCGGGCGGTTTGTCAAAACCCTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                           <  1:668493/90‑1 (MQ=255)
                                                            tgCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                           <  1:314508/90‑1 (MQ=255)
                                                                             aTGAAGGGGGGGGGTTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTg          <  2:800192/90‑1 (MQ=255)
                                                                                     gggcggCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc  >  2:31615/1‑90 (MQ=255)
                                                                                        |                                                                                      
TGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGC  >  NZ_CP009273/3738291‑3738465

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