Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,356,1520CA40.9% 28.5 / 10.2 22A321S (GCC→TCC) yfaYnicotinamide mononucleotide deamidase‑related protein YfaY
Reads supporting (aligned to +/- strand):  ref base C (10/3);  new base A (0/9);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 4.62e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGAAAGTGCCGTCTGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACGCCAGCAA  >  NZ_CP009273/2356073‑2356239
                                                                                |                                                                                       
gCGAAAGTGCCGTCTGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTg                                                                               >  2:406798/1‑90 (MQ=255)
gCGAAAGTGCCGTCTGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTg                                                                               >  2:157604/1‑90 (MQ=255)
    aaGTGCCGTCTGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGccc                                                                           >  1:366370/1‑90 (MQ=255)
        cccGTCGGGGGGGGCTGGCGAAAAGTTGAGATCCTCTTTCCCAAAACCAAAACCAGCCAGTGCCAGCCCGGCAAAATGGTTGGCccgccg                                                                       <  2:212839/89‑1 (MQ=255)
         ccGTCGGGGGTGGCTAGCGAAAAGTGGAGATGCTCTTTCTCAAACCCAAAAACAGCCAGTGCCACCCCGGAAAAATGGTTGGCCCGCCGt                                                                      <  2:290216/90‑1 (MQ=255)
          cGTCTGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGtt                                                                     >  1:212839/1‑90 (MQ=255)
           gTCTGGCGTGGCTAGCGAAAAGTTGAGATGCTCGTTCTCAAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGTTc                                                                    <  1:683200/90‑1 (MQ=255)
           gTCGGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGAAAAATGGTTGGCCCGCCGTTc                                                                    <  1:34864/90‑1 (MQ=255)
           gTCGGGCGTGGCTAGCGAAAAGTTGAAAACCTCTTCCCCAAACCCGAAAACAGCCAGTGCCAGCCCGGAAAAATGGTGGGCCCGCCGTTc                                                                    <  2:617647/90‑1 (MQ=255)
              gggggTGGCGGGCGAAAAGTGGAGTTTCTCGTTCCCAAACCCCAAACCAGCAAGTGCCAGCCGGAAAAAATGGTTGCCCCGCCGTTCTGt                                                                 <  2:649502/86‑1 (MQ=255)
                                        gCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAggg                                       >  1:467245/1‑90 (MQ=255)
                                          tCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGtt                                     >  1:546312/1‑90 (MQ=255)
                                                      cccAAAACACCCGGTGCCACCCCGGAAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGaa                         <  1:350144/90‑1 (MQ=255)
                                                             aCAGCCGGGGCCAGCCCGGAAAAAGGGTTGCCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAAcc                   <  2:93260/90‑1 (MQ=255)
                                                                 ccAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACtt               >  1:144576/1‑90 (MQ=255)
                                                                       ccAGCCCGGAAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACg         <  1:108444/90‑1 (MQ=255)
                                                                         aGCCCGGAAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACGcc       <  1:105318/90‑1 (MQ=255)
                                                                          gCCCGGCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACGCCa      >  2:8263/1‑90 (MQ=255)
                                                                             cGGAAAAAGGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACGCCAGCa   <  1:154974/90‑1 (MQ=255)
                                                                              ggCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGa                          <  1:62078/66‑1 (MQ=255)
                                                                              ggCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGa                          >  2:62078/1‑66 (MQ=255)
                                                                              ggCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACGCCAGCaa  >  2:593706/1‑90 (MQ=255)
                                                                                |                                                                                       
GCGAAAGTGCCGTCTGGCGTGGCTAGCGCAAAGTTGAGATGCTCGTTCTCGAAACCCGAAACAGCCAGTGCCAGCCCGGCAAAATGGTTGGCCCGCCGTTCTGTAATCCAGTGCGCAGTTTGCGCCAGGGTTTCCTCCTGTGAAGGAACCACTTCACACGCCAGCAA  >  NZ_CP009273/2356073‑2356239

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