Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,254,4140AC52.0% 12.0 / 31.9 25V455G (GTG→GGG) tdcE2‑ketobutyrate formate‑lyase/pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (12/0);  new base C (0/13);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.92e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCA  >  NZ_CP009273/3254330‑3254495
                                                                                     |                                                                                 
cTGCCCATCTTTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCccg                                                                              <  1:35521/90‑1 (MQ=255)
        ctttttGTCGTAACCCACCCCGCCTCCCACCGGCGGTGCTTTTTTCGCCCCGACCTAAATCTTCAGCTTCCCCCCCCCCCCGCCGTTAAt                                                                      <  1:319674/87‑1 (MQ=255)
            ttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                  >  2:326772/1‑90 (MQ=255)
            ttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                  >  1:285737/1‑90 (MQ=255)
            ttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                  >  1:359096/1‑90 (MQ=255)
                    aTCCAGCACGTCGTCCATCAGCGGTGCTGGTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCTCCGTTAATTGCGTATagcag                                                          >  2:39139/1‑90 (MQ=255)
                        aGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                      >  2:281224/1‑90 (MQ=255)
                         gCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                     >  2:148792/1‑90 (MQ=255)
                             gtcgtcCACCAGGGGTGCTTTTTTCGCCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                 <  2:359096/90‑1 (MQ=255)
                              tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                >  2:281631/1‑90 (MQ=255)
                              tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGGAGAGCAGCGTTTTCGcc                                                >  1:261346/1‑90 (MQ=255)
                                        gCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcgc                                      >  1:308746/1‑90 (MQ=255)
                                          ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                    >  2:174647/1‑90 (MQ=255)
                                          ggTGCTTTTTCGGCCCCGCCCTAAATCTTCAGCTTCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                    <  1:117237/90‑1 (MQ=255)
                                                 tttcCGGCCCGCCCTAACTTTCCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                             <  2:285737/90‑1 (MQ=37)
                                                     cGGCCCGCCTTAAATCTTCAGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                         <  2:90146/90‑1 (MQ=37)
                                                     cGGCCCGACCTAAATCTTCAGCTTTTTGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                         <  2:58998/90‑1 (MQ=38)
                                                      ggCCCGCCCTAATTCTTCGGCTCCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAAc                        <  1:334945/90‑1 (MQ=255)
                                                           cGACCTGATTCTTCAGTTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                    <  2:304609/90‑1 (MQ=37)
                                                                  attCTTCAGTTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt             <  1:291250/88‑1 (MQ=25)
                                                                  aaTCTTCACCTTCCCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt             <  1:204691/90‑1 (MQ=25)
                                                                       tCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa        >  1:215839/1‑90 (MQ=255)
                                                                         tccTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc      <  1:269703/88‑1 (MQ=255)
                                                                         cgCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTTCATTTGCTTACCAATc      <  1:47177/89‑1 (MQ=255)
                                                                             tctcGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCa  >  1:249655/1‑90 (MQ=255)
                                                                                     |                                                                                 
CTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCA  >  NZ_CP009273/3254330‑3254495

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