Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,774,2150TG56.8% 24.6 / 34.7 37G343G (GGT→GGGlldDquinone‑dependent L‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (8/8);  new base G (21/0);  total (29/8)
Fisher's exact test for biased strand distribution p-value = 3.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAA  >  NZ_CP009273/3774126‑3774302
                                                                                         |                                                                                       
gATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGgt                                                                                         <  2:22291/90‑1 (MQ=255)
gATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGgt                                                                                         <  2:449979/90‑1 (MQ=255)
gATGTCGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGgt                                                                                         <  2:536920/90‑1 (MQ=255)
     cGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGc                                                                                    >  1:188312/1‑90 (MQ=255)
     cGTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGc                                                                                    >  1:682896/1‑90 (MQ=255)
      gTGCGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCt                                                                                   >  2:752847/1‑90 (MQ=255)
         cGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAAc                                                                                <  2:743172/90‑1 (MQ=255)
         cGTATGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAAc                                                                                <  2:252658/90‑1 (MQ=255)
             tGATTGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGc                                                                            >  2:664955/1‑90 (MQ=255)
                 tGCGCTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTGGCTAACCTGCTaaa                                                                        >  1:347829/1‑90 (MQ=255)
                  gcgcTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTGGGTAAACTGCTAAAt                                                                       >  1:341016/1‑90 (MQ=255)
                  gcgcTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGCTAACCTGCTaaaa                                                                       >  2:436644/1‑89 (MQ=255)
                    gcTCGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCt                                                                     >  2:472049/1‑90 (MQ=255)
                       cGGTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGAt                                                                  >  2:27432/1‑90 (MQ=255)
                        ggTGCCGACACCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATc                                                                 <  2:5125/90‑1 (MQ=255)
                                 acCGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTGGCTAACCTGCTAAATCTGATCGAAAAaaag                                                        >  1:487458/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGAt                                                      >  2:346583/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTTATCGAAAAAgagag                                                      >  1:765394/1‑89 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGGTCGAAAAAGAGAt                                                      >  1:163410/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGAt                                                      >  1:312919/1‑90 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAAACTGCTAAATCTGATCGAAAAAaagag                                                      >  2:592297/1‑89 (MQ=255)
                                   cGTACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGCTAACCTTCTAAATCTTCTCGAAAAAGAGAt                                                      >  1:818554/1‑90 (MQ=255)
                                         gctgGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGt                                                >  2:40639/1‑90 (MQ=255)
                                         gctgGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAAACTGATCGAAAAAAAGATGAAAgg                                                >  1:255527/1‑89 (MQ=255)
                                            gggTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGc                                             >  2:565345/1‑90 (MQ=255)
                                              gTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCCGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGa                                           >  1:279434/1‑90 (MQ=255)
                                              gTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGGAGGTAACCTGCTAAATCTGATCGAAAAAGAGAGGAAAggggggg                                           >  1:491499/1‑87 (MQ=255)
                                               tCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGAt                                          <  1:772251/90‑1 (MQ=255)
                                                 gTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGa                                        >  2:151236/1‑90 (MQ=255)
                                                 gTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCCGCTAAATCTGATCGAAAAAAAGATGAAAGTGGGGATGa                                        >  2:589188/1‑90 (MQ=255)
                                                         ttGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACt                                <  1:346583/90‑1 (MQ=255)
                                                           gTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTgg                              >  2:801209/1‑90 (MQ=255)
                                                             aTGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGgcg                            >  2:54445/1‑90 (MQ=255)
                                                                  cTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGGGGGGGTGAAGCTGACTGGCGCGaaa                       >  1:432547/1‑90 (MQ=255)
                                                                           gCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGc              >  1:131204/1‑90 (MQ=255)
                                                                             gggCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGa            <  1:589187/90‑1 (MQ=255)
                                                                              ggCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGGGGGGGTGACGCTGACTGGCGCGGAATCGATCAGCGaa           >  1:789799/1‑90 (MQ=255)
                                                                              ggCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAAAGATGAAAGGGGGGATGACGCTGACTGGCGCGAAATCGATCAGCGaa           >  2:405826/1‑90 (MQ=255)
                                                                                     cGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGc    <  2:302972/90‑1 (MQ=255)
                                                                                       ggTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCaa  <  2:789800/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 < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: