Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,774,2150TG52.0% 22.3 / 31.8 25G343G (GGT→GGG) lldDquinone‑dependent L‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (3/9);  new base G (13/0);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAA  >  NZ_CP009273/3774126‑3774302
                                                                                         |                                                                                       
gATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGgt                                                                                         <  1:275175/90‑1 (MQ=255)
     cGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGc                                                                                    >  2:229140/1‑90 (MQ=255)
      gTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGGGTTTTTTTGTATGGGCGGGAAACAGCGGGCCAGGCGGGGGTGGCt                                                                                   >  2:75662/1‑90 (MQ=255)
         cGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAAc                                                                                <  2:86122/90‑1 (MQ=255)
         cGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAAc                                                                                <  1:778540/90‑1 (MQ=255)
                ttGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGTGCTAACCTGCTaa                                                                         >  1:397691/1‑90 (MQ=255)
                  gcgcTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAAt                                                                       >  1:871599/1‑90 (MQ=255)
                    gcTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGCTAAACTTCTAAAACt                                                                     >  1:525148/1‑90 (MQ=255)
                         gTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAAACTGCTAAATCTGATCg                                                                >  1:363130/1‑90 (MQ=255)
                            ccGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGaaa                                                             <  2:854793/90‑1 (MQ=255)
                             cGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGaaaa                                                            <  1:278649/90‑1 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGAt                                                      >  2:697005/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTGGCTAACCTGCTAAATCTGATCGAAAAAaagaa                                                      >  1:877514/1‑89 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTGGCTAACCTGCTAAATCTGATCGAAAAAGAGAt                                                      >  1:784932/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTTCTAAATCTGATCGAAAAAGAGAt                                                      >  1:452495/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGCTAACCTGCTAAATCTGGTCGAAAAAGAGAt                                                      >  1:480692/1‑90 (MQ=255)
                                            gggTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGc                                             >  2:113914/1‑90 (MQ=255)
                                                   gCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACg                                      <  2:60799/90‑1 (MQ=255)
                                                   gCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAAAGATGAAAGGGGCGATGACg                                      >  2:407047/1‑90 (MQ=255)
                                                         ttGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACt                                <  2:525148/90‑1 (MQ=255)
                                                         ttGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACt                                <  1:350178/90‑1 (MQ=255)
                                                                cgcTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAAAGATGAAAGTGGCGATGACGCTGACTGGCGCGa                         >  2:938807/1‑90 (MQ=255)
                                                                   tGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAAt                      <  1:229140/90‑1 (MQ=255)
                                                                           gCGGGCCAGGCGGGGGTAGCTAACCTGCTAAAACTGATCGAAAAAGAGAAGAAAGGGGGGAAGAGGCTGGCTGGCGCGAAATCGATTAGc              >  1:820875/1‑90 (MQ=255)
                                                                              ggCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGaa           >  2:596849/1‑90 (MQ=255)
                                                                                       ggTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCaa  <  2:824417/90‑1 (MQ=255)
                                                                                         |                                                                                       
GATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAA  >  NZ_CP009273/3774126‑3774302

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