Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273937,1690TG58.8% 1.2 / 18.8 17V252G (GTG→GGG) dmsAdimethylsulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (10/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 5.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT  >  NZ_CP009273/937081‑937236
                                                                                        |                                                                   
ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                    <  1:331113/90‑1 (MQ=255)
   aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                 <  2:80965/90‑1 (MQ=255)
      aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                              <  2:631813/90‑1 (MQ=255)
                                      aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACCCCCGAAAACTAAAGc                              >  2:355499/1‑90 (MQ=255)
                                      aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATTTTTAAAAAGCACGCCAGAAAAATAATGc                              >  1:312128/1‑90 (MQ=255)
                                         gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCAGAAATCTAATGcccc                           >  1:681400/1‑89 (MQ=255)
                                           tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                         <  2:554261/90‑1 (MQ=255)
                                             tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCCCATg                       >  2:243181/1‑90 (MQ=255)
                                              ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTTTTATCTTTTAAAAGCACGCCAGAAAACTAATGCCCCCATGa                      >  1:687634/1‑90 (MQ=255)
                                                      aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCCCATGatcatcatc              >  2:410899/1‑90 (MQ=255)
                                                          cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc          <  2:544027/90‑1 (MQ=255)
                                                          cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc          <  2:611428/90‑1 (MQ=255)
                                                           tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAAAAAACTAAAGCCCCCATGATCATCCCCGATcc         >  2:307150/1‑90 (MQ=255)
                                                           tGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTTAACAGGCACCCCCGAAAACTAAATCCCGCCGGATCATCATCTATcc         >  1:259802/1‑90 (MQ=255)
                                                                  aCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCCCGCtat  >  2:542234/1‑90 (MQ=255)
                                                                          gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                <  1:35517/52‑1 (MQ=255)
                                                                          gAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                >  2:35517/1‑52 (MQ=255)
                                                                                        |                                                                   
GGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT  >  NZ_CP009273/937081‑937236

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