Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,141,2410CA48.0% 31.0 / 13.0 25N29K (AAC→AAA) yegDmolecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (2/11);  new base A (12/0);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 2.02e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATG  >  NZ_CP009273/2141167‑2141320
                                                                           |                                                                               
tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGctgctg                                                                  <  1:835855/90‑1 (MQ=255)
tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGctgctg                                                                  <  1:768416/90‑1 (MQ=255)
tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGctgctg                                                                  <  1:603318/90‑1 (MQ=255)
  tGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGcc                                                                <  1:192291/90‑1 (MQ=255)
       aCGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCC‑TTTTCTAAAAAAAGAAAAAGACAGCACGCCGCTGCCTTCaa                                                           >  2:192291/1‑90 (MQ=255)
                    cTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCC‑TTTTCTAAAAATGGGAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGGGCGcc                                              >  1:639675/1‑90 (MQ=255)
                    cTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCA‑TTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGcc                                              >  1:862022/1‑90 (MQ=255)
                    cTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCA‑TTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGcc                                              >  2:248922/1‑90 (MQ=255)
                    cTGTTCAGTGGCGGTCATGCGGGACGGGAAAACCCA‑TTTTCTAAAAAAGGAAAAAGACAGCACCCCGCTGCCCTCAATGCTTTTCGCGcc                                              >  1:212163/1‑90 (MQ=255)
                     tGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCa                                             <  1:734061/90‑1 (MQ=255)
                     tGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCa                                             <  2:639675/90‑1 (MQ=255)
                      gTTCAGTGGCGGTCTTGCGTGGGGGTTAAACGCAGTTTTTTAAAAAT‑GAAAAAGAAAGCACCCTCCTGCCTTCGATTCTTTTCTCGccca                                            >  2:858062/1‑88 (MQ=255)
                      gTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCC‑TTTTCTAAAAAAAGAAAAAGACAACACCCCGCTGCCCTCAAATCTTTTCGCGccca                                            >  2:421092/1‑88 (MQ=255)
                      gTTCAGTGGCGGTCATGCGTGACGGGAAACCGCA‑TTTGCTAAAAAATGAAAAAGACAGCACGCCGCTGCCTTCAATGCTTTGCGCGCCaa                                            >  1:432613/1‑90 (MQ=255)
                      gTTCAGTGGCGGTCATGCGTGACGGGAAACCGCA‑TTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCaa                                            >  1:228116/1‑90 (MQ=255)
                                ggTCATGCGTGGCGGGAAAACGCC‑TTTTCTAAAAAAGGAAAAAGACAGCACGCCGGTGCCTTCAATGCTTTGCGCGGCAACGCGTGAAGc                                  >  2:775353/1‑90 (MQ=255)
                                ggTCATGCGTGACGGTAAACCCCC‑TTTTCTAAAAATGGAAAAAGACAACACGCCGCCGCCTTCAATGCTTTTCGCGCCAACCCGTGAAGc                                  >  2:345199/1‑90 (MQ=255)
                                ggTCATGCGTGACGGGAAACCGCC‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAAAGCGTGAAGc                                  >  1:124420/1‑90 (MQ=255)
                                ggTCATGCGTGACGGGAAAACGCA‑TTTTCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATTCTTTTCGCGCCAACGCGTGAAGc                                  >  2:88282/1‑90 (MQ=255)
                                    aTGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTa                              <  1:870444/90‑1 (MQ=255)
                                           aCGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAAt                       <  1:88282/90‑1 (MQ=255)
                                           aCGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAAt                       <  1:581800/90‑1 (MQ=255)
                                            cGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATg                      <  1:775353/90‑1 (MQ=255)
                                                              tAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatg     <  2:536478/90‑1 (MQ=255)
                                                                 aaaTGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgatg  >  1:207807/1‑90 (MQ=255)
                                                                           |                                                                               
TTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCA‑TTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATG  >  NZ_CP009273/2141167‑2141320

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