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

GTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATT  >  NZ_CP009273/3254332‑3254497
                                                                                  |                                                                                   
gTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                              >  1:72317/1‑90 (MQ=255)
      cccTTTGTCGTACCCCCCCCCGTCGCCCACCGGCGGTGTTTTTTTCGCCCCGACTTGATTCTTCACCTTCCCGCCCCCCCCGCCGTTAAt                                                                        <  1:766480/88‑1 (MQ=255)
      cccTTTGTCGTAACCCAGCACGTCTCCCACCAGGGGGGCTTTTTTCGGCCCGCCCTAAATTTTCATCTTCCCGCCCCCCCCGCCGTTAAt                                                                        <  1:574514/88‑1 (MQ=255)
          ttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                    >  2:802806/1‑90 (MQ=255)
                  acccAGCGCGTCGTCCATCGGGGGTGCTTTTTTCGGCCCGCCCTAATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGagcag                                                            <  2:535597/88‑1 (MQ=255)
                  acccAGCCCGTCGTCCATCGGGGGTGTTGTTTTCGCCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGagcag                                                            <  2:766982/88‑1 (MQ=255)
                  aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagcag                                                            >  1:535597/1‑90 (MQ=255)
                      aGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                        >  2:370625/1‑90 (MQ=255)
                           gtcgtcCATCAGCGGTGTTTTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                   <  2:186131/90‑1 (MQ=255)
                           gtcgtcCATCAGCGGTGTTGTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                   <  2:72317/90‑1 (MQ=255)
                           gtcgtcCATCAGCGGTGCTTTTTTCGGCCCGCCCTAAATTTTCAGCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                   <  2:100832/90‑1 (MQ=255)
                           gtcgccCACCGGGGGTGTTTTTTCCGCCCCCACTTAATTCTTCAGCTCCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                   <  1:662971/90‑1 (MQ=255)
                            tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                  >  1:186131/1‑90 (MQ=255)
                            tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                  >  2:39978/1‑90 (MQ=255)
                            tcgtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                  >  2:312554/1‑90 (MQ=255)
                                        ggTGTTTTTTTCGCCCCGCCCTAAATTTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                      <  2:88634/90‑1 (MQ=255)
                                        ggTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                      >  1:182353/1‑90 (MQ=255)
                                          tGCTGTTTTCGCCCCGACCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc                                    <  2:18867/90‑1 (MQ=255)
                                              gTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                >  2:36288/1‑90 (MQ=255)
                                               ttttCGGCCCGTCCAGATTCTTCCCCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                               <  1:370625/90‑1 (MQ=25)
                                                        cGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                      >  2:778866/1‑90 (MQ=255)
                                                               aaTCTTCTGCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt               <  1:443242/90‑1 (MQ=37)
                                                               aaTCTTCAGCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt               <  2:88471/90‑1 (MQ=37)
                                                                      aGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc        <  2:47326/90‑1 (MQ=255)
                                                                          tctcGTCCACACCGCCGTCAATTACGTAGAGCTGCGTGTTCGCCGGGTTAGCGCGTACCCCAAGGAACTGCATTTGCTTACCAACCActa    >  1:343679/1‑88 (MQ=255)
                                                                            tcGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt  <  2:486043/90‑1 (MQ=255)
                                                                                  |                                                                                   
GTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATT  >  NZ_CP009273/3254332‑3254497

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