Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273937,1710AG62.5% 0.5 / 11.9 16T253A (ACT→GCT) dmsAdimethylsulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base A (1/5);  new base G (10/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 1.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCC  >  NZ_CP009273/937087‑937251
                                                                                    |                                                                                
aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                             <  1:112772/90‑1 (MQ=255)
     gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCt                                                                        >  2:475971/1‑90 (MQ=255)
            atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTATTATCTTTAAAAg                                                                 >  1:255068/1‑90 (MQ=255)
                     aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg                                                        <  2:296305/90‑1 (MQ=255)
                       tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAAAAGGGACGCCAGaa                                                      >  2:321711/1‑90 (MQ=255)
                              gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTTTTTTTAAACAGCACCCCCGAAAACTAAt                                               >  1:348291/1‑90 (MQ=255)
                                   caggTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg                                          <  2:300320/87‑1 (MQ=255)
                                          ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                   <  2:195384/90‑1 (MQ=255)
                                           gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGGTTAATATCTTGAAAAGGCACGCCCGAAAACTAATGCCCGCGTGagca                                  >  1:357947/1‑90 (MQ=255)
                                                     tGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTGAACAGGGACGCCAGAAAACTAATTCCCGCATTAGCATCATCGATcc                        >  2:188991/1‑90 (MQ=255)
                                                     tGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTTCTATCTTTAAAAAGGACGCCCAAAAACTTATGCCCCCCTGATTATCATTGGTcc                        >  2:154636/1‑90 (MQ=255)
                                                          aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                   <  1:321711/90‑1 (MQ=255)
                                                            aCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCCCGCtat                 >  2:596100/1‑90 (MQ=255)
                                                            aCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGGACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCtct                 >  1:328499/1‑88 (MQ=255)
                                                                           ggTGGGGGGGCTTTCTATCTTTTACAAGCACGCCAGAAATATAATGCCCGCAAGATCATCATCGGTCCGCGCGATACCCACACCCGGGcc  >  1:281181/1‑90 (MQ=255)
                                                                           ggTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACACCCGGGcc  >  1:310366/1‑90 (MQ=255)
                                                                                    |                                                                                
AGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCC  >  NZ_CP009273/937087‑937251

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