Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,556,8410TG57.9% 3.7 / 16.6 19V490G (GTG→GGG) glpDglycerol‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base G (11/0);  total (15/4)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACA  >  NZ_CP009273/3556753‑3556930
                                                                                        |                                                                                         
cgccgcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtg                                                                                          <  2:510541/90‑1 (MQ=255)
          acgacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc                                                                                <  1:612252/90‑1 (MQ=255)
          acgacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACgc                                                                                >  1:456562/1‑90 (MQ=255)
              cgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagca                                                                            >  1:369748/1‑90 (MQ=255)
              cgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagca                                                                            <  2:247138/90‑1 (MQ=255)
                    gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGtt                                                                      >  2:568744/1‑90 (MQ=255)
                          tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGTTATCGCt                                                                >  2:320088/1‑90 (MQ=255)
                          tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGGTATCGCt                                                                >  2:16202/1‑90 (MQ=255)
                            gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATCGCTgg                                                              >  2:497645/1‑90 (MQ=255)
                              acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAGGGTCTCGCAGGcg                                                            >  2:112787/1‑90 (MQ=255)
                                aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCGGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGcgtc                                                          >  1:192325/1‑90 (MQ=255)
                                  aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGcgtcgt                                                        <  2:360070/90‑1 (MQ=255)
                                           tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAAt                                                    <  1:38698/85‑1 (MQ=255)
                                           tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTAGAGAATACGGAGCAGAGGTTTTCGCTGGAGTCCTAAAATaaag                                               >  1:446445/1‑88 (MQ=255)
                                           tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAAt                                                    >  2:38698/1‑85 (MQ=255)
                                                        ggATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGGGATCAgg                                  >  2:229927/1‑90 (MQ=255)
                                                        ggATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAGGTTATCGCTGGGGGCGTAAATTAACGTAAGGTGGTCAgg                                  >  2:569757/1‑90 (MQ=255)
                                                                   ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa                       >  2:57802/1‑90 (MQ=255)
                                                                      gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGGGATCAGGTCAGATTTCAATCt                    >  2:144943/1‑90 (MQ=255)
                                                                      gtgtGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGGTATCGCTGGGGGCGTAAATTAACGGAAGGTGATCAGGTCAGGTTTCAATCt                    >  2:403638/1‑90 (MQ=255)
                                                                                      ggtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatga    >  1:367788/1‑90 (MQ=255)
                                                                                        tggAGTATGCGCAGCAGAGGTTATCGCTGGCGTCGAAAATTAACGTAAGGTGATAAGGTCAGATTTCAATCTGGCCTGAGACTGATGACa  <  2:370724/90‑1 (MQ=255)
                                                                                        |                                                                                         
CGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACA  >  NZ_CP009273/3556753‑3556930

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