Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,498,4400CT60.0% ‑0.0 / 12.4 15E314K (GAA→AAA) yhfSaminotransferase class I/II‑fold pyridoxal phosphate‑dependent enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base T (9/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 2.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAAT  >  NZ_CP009273/3498353‑3498515
                                                                                       |                                                                           
tCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCat                                                                           <  2:239554/90‑1 (MQ=255)
         tGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGAct                                                                  <  1:86966/90‑1 (MQ=255)
         tGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGAct                                                                  <  2:257249/90‑1 (MQ=255)
            gAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTCActctc                                                               >  1:303711/1‑89 (MQ=255)
                          tCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTTACTCTGCCCCCCCCCGGGaa                                                 >  1:506763/1‑90 (MQ=255)
                          tCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTaa                                                 <  1:93948/90‑1 (MQ=255)
                                 gtgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTTATATTTTTACTCTGCACCCCCCCGGGaaagcaagg                                          >  2:344262/1‑90 (MQ=255)
                                  tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTTTATTTTTACTCTTCCCCCCCCCGGGaaaggaaggc                                         >  1:379394/1‑90 (MQ=255)
                                  tgGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTTACTCTGCACCCCCCCGGGaaagcaaagc                                         >  2:382613/1‑90 (MQ=255)
                                     aTTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTGACTCTGCACCCCCCCGGGAAGGCAAGGCAcc                                      >  2:62735/1‑90 (MQ=255)
                                       tCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTTATATTTTGACTCTGCACCCACCGGGGAAGGCAAGGCACcgc                                    >  2:137346/1‑90 (MQ=255)
                                                   aCGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCct                        <  1:137346/90‑1 (MQ=255)
                                                       tCCGGAAAGGCGATAAAAGAGCGGCGGGGTTTTATATTTTTACCTTGCACCCCCCCGGGGAGGCAAGGGCACGCGCTTTTTCGCCCCTTc                    >  1:456047/1‑90 (MQ=255)
                                                         cGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCa                  <  2:379394/90‑1 (MQ=255)
                                                                         gagCGGCGGGATTTTATATTTCGACTCTGCACCCCCCGGGGAAGGGAAGGCACCGCGCTTTTTCGCCTCTTCCAGCACTCTGGCGGGaaa  >  2:532609/1‑89 (MQ=255)
                                                                                       |                                                                           
TCGGGTTAATGCGAATCGCACAATGTTCTGATTGTGGATTCGCCTGGCGAAACGTTCCGGAAAGGCGATAAAAGAGCGGCGGGATTTCATATTTCGACTCTGCACCCACCGGGTAAGGCAAGGCACCGCGCTTTTGCGCCTCTTCCAGCACTCTGGCGGCAAT  >  NZ_CP009273/3498353‑3498515

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