Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273937,1690TG46.0% 60.8 / 23.5 50V252G (GTG→GGG) dmsAdimethylsulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (1/26);  new base G (23/0);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.22e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGC  >  NZ_CP009273/937080‑937257
                                                                                         |                                                                                        
tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                          <  2:873071/90‑1 (MQ=255)
tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                          <  2:1211581/90‑1 (MQ=255)
tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                          <  1:1773447/90‑1 (MQ=255)
 ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                         <  2:1531828/90‑1 (MQ=255)
  gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                        <  2:70665/90‑1 (MQ=255)
  gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                        <  1:1774280/90‑1 (MQ=255)
   cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAc                                                                                       >  1:252531/1‑90 (MQ=255)
       aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATACCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                   <  2:1542738/90‑1 (MQ=255)
       aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                   <  1:1606819/90‑1 (MQ=255)
        gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACt                                                                                  >  1:325331/1‑90 (MQ=255)
           cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                               >  1:1349125/1‑90 (MQ=255)
                    tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTGCTATCTTTaaaaag                                                                      >  2:2049708/1‑86 (MQ=255)
                     cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                     <  1:1799550/90‑1 (MQ=255)
                     cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                     <  1:1366849/90‑1 (MQ=255)
                          aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                <  2:433665/90‑1 (MQ=255)
                            aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg                                                              <  2:1349125/90‑1 (MQ=255)
                             gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTGAACAGGCACGCCAGa                                                             >  1:2051089/1‑90 (MQ=255)
                                    ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTTAAAAAGGACGCCCGAAAACTaa                                                      >  2:254140/1‑90 (MQ=255)
                                     gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                   <  2:1996292/76‑1 (MQ=255)
                                     gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCAc                                                                   >  1:1996292/1‑76 (MQ=255)
                                     gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTTAAAAGGCACGCCCGAAAACTAAt                                                     >  2:569910/1‑90 (MQ=255)
                                       aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                 <  2:422066/60‑1 (MQ=255)
                                       aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTa                                                                                 >  1:422066/1‑60 (MQ=255)
                                             gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                             <  2:1525550/90‑1 (MQ=255)
                                              tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTGACAGGGACCCCCGAAAACTAATTCCCCCATg                                            >  1:978954/1‑90 (MQ=255)
                                              tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTTTTATCTTTAACAAGGGCGCCCGAAAATTAATGCCCCCATg                                            >  2:1437084/1‑90 (MQ=255)
                                               ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                           <  1:569910/90‑1 (MQ=255)
                                                 ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                         <  1:1890230/90‑1 (MQ=255)
                                                       aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                   <  1:1249459/90‑1 (MQ=255)
                                                       aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCATGatcatcatc                                   >  1:2196993/1‑90 (MQ=255)
                                                       aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                   >  1:565066/1‑90 (MQ=255)
                                                       aaCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGatcatcatc                                   >  1:2045595/1‑90 (MQ=255)
                                                           cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                               <  1:854282/90‑1 (MQ=255)
                                                           cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                               <  1:532126/90‑1 (MQ=255)
                                                           cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                               <  1:1737621/90‑1 (MQ=255)
                                                           cTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                               <  2:2177725/90‑1 (MQ=255)
                                                           cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCAGAAAACTAATTCCCGCCTGATCATCATTGATc                               >  2:854282/1‑90 (MQ=255)
                                                           cTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTGAAAAGGCACGCCAGAAAACTAAAGCCCGCATGATCATCATCGATc                               >  1:1928367/1‑90 (MQ=255)
                                                            tGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                 <  1:1547084/39‑1 (MQ=255)
                                                            tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCAGAAAACTAATGGCCCCATGATCATCATCGATcc                              >  2:763896/1‑90 (MQ=255)
                                                            tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATcc                              >  2:302834/1‑90 (MQ=255)
                                                            tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACCCCAGAAAATTAATGCCCCCATGATCATCatccatcc                              >  1:1332692/1‑90 (MQ=255)
                                                            tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTa                                                                                 >  2:1547084/1‑39 (MQ=255)
                                                            tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc                              >  1:769951/1‑90 (MQ=255)
                                                              gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                            <  1:1876452/90‑1 (MQ=255)
                                                                 aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                         <  1:2049708/90‑1 (MQ=255)
                                                                 aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                         <  1:2199580/90‑1 (MQ=255)
                                                                   aCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGGACCCCAGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCtat                       >  1:1879541/1‑90 (MQ=255)
                                                                       gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                   <  2:769951/90‑1 (MQ=255)
                                                                                  ggTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc        >  1:810168/1‑90 (MQ=255)
                                                                                    tGGGGGGACTTTCTATCTTTAACAGGGACGCCAGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCCATAACCACACCCGGGCCgg      >  2:567136/1‑90 (MQ=255)
                                                                                        gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc  <  2:537893/90‑1 (MQ=255)
                                                                                        gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc  <  1:407352/90‑1 (MQ=255)
                                                                                         |                                                                                        
TGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGC  >  NZ_CP009273/937080‑937257

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