Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG45.2% 60.3 / 15.6 43V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (2/21);  new base G (19/0);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 9.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCG  >  NC_000913/941578‑941861
                                                                                                                                       |                                                                                                                                                    
cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                       <  1:736153/151‑1 (MQ=255)
cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                       <  2:362565/151‑1 (MQ=255)
   aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGaa                                                                                                                                      >  2:506704/1‑149 (MQ=255)
          ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCtttgacaagcgagc                                                                                                                             >  1:787955/1‑139 (MQ=255)
           gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGACTTACTATTTTGAAAAAGCAAGcc                                                                                                                            >  2:622777/1‑151 (MQ=255)
              ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATTTTTAAAAGGCACGcc                                                                                                                            >  1:468909/1‑148 (MQ=255)
               gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGCACGcc                                                                                                                            >  2:144650/1‑147 (MQ=255)
                 ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                       <  2:658512/150‑1 (MQ=255)
                  tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                      <  2:220132/150‑1 (MQ=255)
                    aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTTAACCAGGACGcccgga                                                                                                                        >  1:466071/1‑142 (MQ=255)
                     aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                  <  2:108395/151‑1 (MQ=255)
                        ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGTAAAAGCAAGGCAGAAATATaa                                                                                                                  >  1:124275/1‑148 (MQ=255)
                           ccTGCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACTCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATg                                                                                                                <  1:829252/147‑1 (MQ=255)
                                 gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGat                                                                                                      <  1:808667/151‑1 (MQ=255)
                                     cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGGAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAAGGGCGGCTTTGTTTGTTTCACCGGCACGACGATAATTTAATGTCCGCATgcgc                                                                                                     >  2:515708/1‑145 (MQ=255)
                                      tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATTTTGAAAAAGGAAGCCCGAAAACTAAAGCCCGcc                                                                                                          >  1:45536/1‑141 (MQ=255)
                                                gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                       <  2:109660/151‑1 (MQ=255)
                                                       cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                <  2:466110/151‑1 (MQ=255)
                                                              gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCCATACCCACAcc                                                                          >  1:105299/1‑150 (MQ=255)
                                                               atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg                                                                        >  2:870365/1‑151 (MQ=255)
                                                                  tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGGACGCCAGGAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCGGTg                                                                      >  2:388779/1‑150 (MQ=255)
                                                                        aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                               <  2:676061/151‑1 (MQ=255)
                                                                         cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGAAAAAGCCACCCGAAAATATAAAAGCCGCATGATCAAACTCGATCCGCGCCATACCgcaacccgggcccgg                                                                 >  1:665522/1‑134 (MQ=255)
                                                                          aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                             <  2:60082/151‑1 (MQ=255)
                                                                             aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGaa                                                           >  1:345941/1‑150 (MQ=255)
                                                                                 tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAAGCACGCCCGAAAACTAATGCCCGCATGATTCTCCTCGATCCGCGCTATACCGAACCCGGGGCCGGGCGCGGAGAt                                                        >  1:814131/1‑149 (MQ=255)
                                                                                  ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                      <  1:515666/150‑1 (MQ=255)
                                                                                       tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                            <  1:591486/123‑1 (MQ=255)
                                                                                       tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATTTCTTTGAAAAGCACGCCAGAAATCTAAAGCCCGCATGATTATCATCGATCCCCGCCATACCCACACCGGGGCCGGGCGGGGAGATGAggg                                                   >  1:481711/1‑146 (MQ=255)
                                                                                           gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCGTCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                             <  1:152371/150‑1 (MQ=255)
                                                                                           gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                            <  1:459364/151‑1 (MQ=255)
                                                                                           gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATTATCGATCCGCGCTATACCGACAACGGTGCCGGGGGCGAAGATGAGGGGATccc                                             >  1:92617/1‑150 (MQ=255)
                                                                                               ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                         <  2:566175/150‑1 (MQ=255)
                                                                                                    taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATTTTTAAAATGCCCACCAGAAATCTAATGCCCCCATGATCATCATTGATCCGCGCTATACCCTGACCGGGGCCGGGCTCGAAGATGTGTGGACCCCTATTCGccc                                    >  2:196941/1‑150 (MQ=255)
                                                                                                     aaCCCTGGCGAAACGCGAATGAGTGGCGGTAGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                   <  2:761419/150‑1 (MQ=255)
                                                                                                          tGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTTTTATCTTTAACAGGCACGCCCGAAAAATAAAGCCCGCGTGATCATCATCGGTCCCCGCTATACCGACCCCGGTGCCGGGGGCGCAGGAGGGGGGGTACCTATTCGTCCGGGGAc                              >  1:198497/1‑150 (MQ=255)
                                                                                                            gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGa                           <  2:819838/151‑1 (MQ=255)
                                                                                                            gCGAAACGCGAATGAGTGGCGGTGGGGGGACCTTATATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCATTATCATCATAGACCCGCGCCATATCGACACCCGGGCCGGGCGGGGAGAAGAGGGGAACCCTATTTTTCCgg                                  >  2:34916/1‑144 (MQ=255)
                                                                                                             cGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACACGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCCATACCCACACCCGGGGCGGGCGCGAAGATGAGGGGGTCCCCATTCGTCCggg                                 >  1:138440/1‑144 (MQ=255)
                                                                                                                             ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg          <  2:757597/151‑1 (MQ=255)
                                                                                                                                 gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt      <  2:858478/151‑1 (MQ=255)
                                                                                                                                   ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg    >  2:23199/1‑151 (MQ=255)
                                                                                                                                     ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg  <  2:521507/151‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
CGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCG  >  NC_000913/941578‑941861

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑