Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,498,4470CT46.0% 52.4 / 26.0 50S311S (TCG→TCA) yhfSaminotransferase class I/II‑fold pyridoxal phosphate‑dependent enzyme
Reads supporting (aligned to +/- strand):  ref base C (2/25);  new base T (22/1);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 9.07e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGA  >  NZ_CP009273/3498358‑3498532
                                                                                         |                                                                                     
ttAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTc                                                                                       <  2:615655/90‑1 (MQ=255)
   aTGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGAc                                                                                    >  1:1490586/1‑90 (MQ=255)
    tGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGAct                                                                                   <  2:1199517/90‑1 (MQ=255)
     gCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGActc                                                                                  <  2:209456/90‑1 (MQ=255)
     gCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGActc                                                                                  <  1:702044/90‑1 (MQ=255)
       gAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGActctt                                                                                >  1:945738/1‑89 (MQ=255)
       gAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGGTTTTTATTTTTActctc                                                                                >  2:60720/1‑89 (MQ=255)
        aaTCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGc                                                                               >  1:434300/1‑90 (MQ=255)
        aaTCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTCATATTTTGActcttc                                                                               >  1:734843/1‑88 (MQ=255)
        aaTCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGGTTCATATTTTTActctcc                                                                               >  2:1132827/1‑88 (MQ=255)
        aaTCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTTGACTCTGc                                                                               >  1:398090/1‑90 (MQ=255)
           cGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacc                                                                            <  2:497170/90‑1 (MQ=255)
            gCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTTACTCTCcaccc                                                                           >  2:1077678/1‑90 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  2:1490586/90‑1 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  2:1256087/90‑1 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  2:1012971/90‑1 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  1:472324/90‑1 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  1:112201/90‑1 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  1:660355/90‑1 (MQ=255)
              acaATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGcacccac                                                                         <  1:310457/90‑1 (MQ=255)
                   gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGt                                                                    <  2:447137/90‑1 (MQ=255)
                   gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGt                                                                    <  1:1315812/90‑1 (MQ=255)
                   gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGt                                                                    <  1:932666/90‑1 (MQ=255)
                   gTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGt                                                                    <  1:31502/90‑1 (MQ=255)
                     tCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaa                                                                  <  2:217876/90‑1 (MQ=255)
                       tGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTTACTCTGCACCCCCCCGGGaagg                                                                >  2:472324/1‑90 (MQ=255)
                           tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTGACTCTGCACCCCCCCGGGaaggcaag                                                            >  1:860682/1‑90 (MQ=255)
                           tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTCATATTTTGTCTCTGCACCCCCCCGGTaaggcaag                                                            >  2:151427/1‑90 (MQ=255)
                           tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTCATATTTTGACTCTGCACCCCCCCGGGaaggcaag                                                            >  2:475532/1‑90 (MQ=255)
                           tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGGTTTTTATTTTTACTCTTCACCCCCCCGGGaaagcaag                                                            >  2:717469/1‑90 (MQ=255)
                           tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTTATATTTTCACTCTGCACCCCCCCGGGaaggcaag                                                            >  1:1494539/1‑90 (MQ=255)
                           tgtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACCCTGCACCCCCCCGggaagggaag                                                            >  2:52145/1‑90 (MQ=255)
                             tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCACCCACCCGGGaaagcaaggc                                                          >  2:520594/1‑90 (MQ=255)
                             tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGGTTTATATTTTGACTCTTCCCCCCCCCGGGaaaggaaggg                                                          >  2:1072994/1‑89 (MQ=255)
                              ggATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCACCCCCCCGGGAaggcaaggca                                                         >  1:1210668/1‑90 (MQ=255)
                                  tCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTTTTTTTACTCTGCACCCCCCCGGGAAAGGAAGGGACcgc                                                     >  2:1398371/1‑90 (MQ=255)
                                  tCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTTATATTTTGACTCTGCACCCCCCGGGGAAGGGAAAGCACcgc                                                     >  1:1444758/1‑90 (MQ=255)
                                              aCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCCCCCCCCCGGGAAGGCAAAGCACCCCGCCTTTTCGccc                                         >  1:779938/1‑89 (MQ=255)
                                              aCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCct                                         <  2:859129/90‑1 (MQ=255)
                                              aCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCct                                         <  1:1072994/90‑1 (MQ=255)
                                                     ggAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAg                                  <  2:284568/90‑1 (MQ=255)
                                                      gAAAGGCGATGAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGc                                 <  2:1054978/90‑1 (MQ=255)
                                                      gAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGc                                 <  2:667514/90‑1 (MQ=255)
                                                                      gcggcgGGGTTTTATATTTTGACTCTGCCCCCCCCCGGGAAGGGAAAGGACCCCGCCTTTTCGCCCCTTTCCGCACTCTGGCGGCAAACg                 >  2:1136790/1‑90 (MQ=255)
                                                                        ggcggGATTTCATATTTTGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCct                                         <  2:291403/64‑1 (MQ=255)
                                                                        ggcggGATTTCATATTTTGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCct                                         >  1:291403/1‑64 (MQ=255)
                                                                            ggATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCtgat           <  1:1398371/90‑1 (MQ=255)
                                                                            ggATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCtgat           <  1:1393446/90‑1 (MQ=255)
                                                                              gttttATATTTTTACTCTGCACCCCCCCGGGAAGGCAAGGCACCCCGCTTTTTTGCCCCTTCCCGCCCCCTTGCGGCAATCGGCtggtga         >  2:1106535/2‑90 (MQ=255)
                                                                                     atTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGa  <  2:434300/90‑1 (MQ=255)
                                                                                     atTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGa  <  1:264733/90‑1 (MQ=255)
                                                                                         |                                                                                     
TTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAATCGGCTGATGAAACTCGA  >  NZ_CP009273/3498358‑3498532

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