Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,774,2150TG56.2% 6.4 / 13.3 18G343G (GGT→GGGlldDquinone‑dependent L‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (9/0);  total (14/4)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATT  >  NZ_CP009273/3774136‑3774296
                                                                               |                                                                                 
gTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGAGAAGCCAAcc                                                                         >  2:110460/1‑90 (MQ=255)
               gTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGgt                                                                                   <  1:28479/65‑1 (MQ=255)
               gTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTa                                                                    <  2:249259/80‑1 (MQ=255)
               gTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTa                                                                    >  1:249259/1‑80 (MQ=255)
               gTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCg                                                          >  2:111196/1‑90 (MQ=255)
                   cGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGaaaa                                                      <  2:222527/90‑1 (MQ=255)
                         cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGAt                                                >  2:102972/1‑90 (MQ=255)
                         cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGCTAACCTGCTAAAACTGGTCGAAAAAAAGAt                                                >  2:230234/1‑90 (MQ=255)
                            aCTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGaa                                             <  2:217939/90‑1 (MQ=255)
                               gctgGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGAAGAAAgg                                          >  2:155759/1‑89 (MQ=255)
                                    gTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGGTCGAAAAAAAGATGAAAGTGGGGa                                     >  2:226600/1‑90 (MQ=255)
                                    gTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCCGCTAAATCTGATCGAAAAAAAGATGAAAAGGGCGa                                     >  1:287558/1‑90 (MQ=255)
                                    gTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGa                                     >  2:21354/1‑90 (MQ=255)
                                                 gTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGTTCGAAAAAGAGATGAAAAGGGCGATGACGCTGACTgg                        >  1:367608/1‑90 (MQ=255)
                                                 gTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCCAAATCTGTAAGAAAAAGAGATGAAAGTGGCGCAGAGGctggctgg                        >  2:566720/1‑90 (MQ=255)
                                                   aTGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGgcg                      >  1:292039/1‑90 (MQ=255)
                                                   aTGCGCTGGCAACAGCGGGCCAGGCGGGCGTAGCTAACCTGCTAAATCTGATCGAAAAAGAAATGAAAATTGCGATGAAGCTGGCTGgcg                      >  1:286072/1‑90 (MQ=255)
                                                                 gCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATTTGTTCGAAAAAGAGATGAAAATGGGGATGACGCCGAATGGGGCGAAAACGATCAGc        >  1:442528/1‑90 (MQ=255)
                                                                       cAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAAtt  <  2:286072/90‑1 (MQ=255)
                                                                               |                                                                                 
GTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATT  >  NZ_CP009273/3774136‑3774296

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