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

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862
                                                                                                                                                 |                                                                                                                                                         
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACtt                                                                                                                                                      <  2:249901/150‑1 (MQ=255)
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                       <  1:238760/149‑1 (MQ=255)
   cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                                   <  2:1091895/150‑1 (MQ=255)
   cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                                   <  1:345511/150‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTGATAACCCTGGCGAAACGGGAATAAGGCGCGGTGGGGGGaaaaa                                                                                                                                                     >  1:121672/1‑141 (MQ=255)
        gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGGTTACTATCtt                                                                                                                                             >  2:294059/1‑150 (MQ=255)
           cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATATCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                         <  1:530638/151‑1 (MQ=255)
           cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                         <  1:429878/151‑1 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                      <  2:754894/151‑1 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                    <  2:179822/150‑1 (MQ=255)
                  cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                                  <  2:100016/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                               <  1:767160/151‑1 (MQ=255)
                            gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAATTACTATACTTAAACAGCACGCCCGaa                                                                                                                            >  1:1088462/1‑147 (MQ=255)
                             aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGCAAACAGTAAGCTGGTAGTGGTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGggtgggtggtattagacattg                                                                                                                                             >  1:61771/1‑117 (MQ=255)
                                 ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                                                                                                     <  1:409830/149‑1 (MQ=255)
                                 ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                                   <  2:966376/151‑1 (MQ=255)
                                    ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg                                                                                                                <  1:1131964/151‑1 (MQ=255)
                                     cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTATTATTTTGAAAAAGGACGCCCGAAAACTAATGgcc                                                                                                                 >  2:401641/1‑149 (MQ=255)
                                      tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTAATTTCTTTTAAAAGCACACCCGAAAACTaa                                                                                                                      >  1:698067/1‑143 (MQ=255)
                                        cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                            <  2:698171/151‑1 (MQ=255)
                                          gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGat                                                                                                          <  2:1186804/151‑1 (MQ=255)
                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                        <  2:1106947/151‑1 (MQ=255)
                                              cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGGTTAATATCTTATAAAAGAACGCCCGAAATATAATGCCCGccggaacat                                                                                                       >  1:449090/1‑142 (MQ=255)
                                               tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGGGGTTCCCATATTTTACAGGCACGCCCAAAAACTAATTTCCGCCTGaac                                                                                                         >  1:955794/1‑145 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                    <  1:773554/150‑1 (MQ=255)
                                                   cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTAATATCTTCGACACGGACCCCAGAAATCCAATGCCCGcccgaccat                                                                                                       >  1:478067/1‑137 (MQ=255)
                                                       tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCAAGATCATCATCGATcc                                                                                             >  1:244048/1‑150 (MQ=255)
                                                         gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                           <  2:878658/150‑1 (MQ=255)
                                                           aaCAGCCATTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                        <  1:987517/151‑1 (MQ=255)
                                                             cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                       <  1:804802/150‑1 (MQ=255)
                                                              aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAGGCACGCCAGAAAACTAATGGCCGCATGATCATCATCGATCCGCGCtct                                                                                      >  1:1042876/1‑148 (MQ=255)
                                                                       gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAAAATAATGCCCGCATGAGCATCATTGGTTCCCGCCATACCGacaac                                                                             >  2:276698/1‑148 (MQ=255)
                                                                            cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                      <  2:155534/151‑1 (MQ=255)
                                                                            cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTTTCTTTAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCCTCGATCCCCGCCATACCGACACCGGTGcc                                                                      >  2:1230681/1‑151 (MQ=255)
                                                                                aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTTCTCTTTTGAACAAGCACGCCAGAAAAATAAATCCCCCCTGAGCATCACCCACCCGCCGTATACCCCCACCCGGGGCCgg                                                                    >  1:1189029/1‑150 (MQ=255)
                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCCTCATCCATCCCCCCTATACCGACACCGGTGCCgg                                                                    >  2:1180169/1‑147 (MQ=255)
                                                                                       aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCCCATGATCATCATCGATCCGCACTATACCCACACCCGGGCCGGgcg                                                                 >  2:700142/1‑145 (MQ=255)
                                                                                         cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGACTTCTATCTTTAACAAGCACACCCGCAATTTTATAGCCGCATGGTCATTAACCGATCACGCCATACCCACAACGGCGGCGGGCGCGCAGa                                                           >  2:935939/1‑149 (MQ=255)
                                                                                           ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGAATTTATCACTTGAAAAGGCACGCCAGAAAAATAAATCCCCCCCGATCCTCATCGATTCGCGCTATACCGACACCGGGGGCggggg                                                                 >  2:636901/1‑139 (MQ=255)
                                                                                              aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGGACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATAACGACACCCGGGCCGGGCGCGGGGATGAgggg                                                    >  1:773500/1‑151 (MQ=255)
                                                                                                 gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACCTACTATATTGAAAACGCGCGCCCCCAATCTAATGGCCCCCCGATGATCATCGACGCCCGCTTTACCCTGACCAGTGCCGGGGGCGgagagatgagg                                                    >  2:462430/1‑141 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                              <  2:984027/150‑1 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCCTGATCATTCTCCATACGCGCTCTAACCCCACCGGTGTCGGGGGCGGAGATGAGTgggaccc                                               >  1:1096768/1‑145 (MQ=255)
                                                                                                       tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  2:1096851/150‑1 (MQ=255)
                                                                                                             taaCCCTGGCGAAACGCGAATGAGTGGCGGTAGGGGGGACTACTTTTTTTAAAAAGCAAGCCCGGAATCTAAATCCCCCATGGGCATCATCGATCCCCGGCATAACGCCACCCCGGCCGGGCGCCGAGATTAGAgggtacccactgttcc                                      >  1:205134/1‑136 (MQ=255)
                                                                                                               aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCCATACCGACACCGGTGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCCggg                                  >  2:80979/1‑151 (MQ=255)
                                                                                                                cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                                 <  1:664128/151‑1 (MQ=255)
                                                                                                                   tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGGACGCCAGAAAAATAATGCCCGCATGATCATTATCGATCCGCGCTATACCCACACCGGTGCCGGGGGCGGAGAAGAGGGGATCCCCATTTCTCCgggga                                >  1:1076977/1‑147 (MQ=255)
                                                                                                                      cGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCAGAAATCTAAAGCCCGCATGATCCTCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTCCcgg                                  >  2:430199/1‑144 (MQ=255)
                                                                                                                       ggaaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGACCGCCCGAAAACTAATGACCGCCTGATCATCATCGGTCCCCGCTATACCCACACCGGCGCCGGGGGCGAAGATGAGGGGGTCCCTCTTCGTCCGGGGACAGATg                          >  1:203110/3‑151 (MQ=255)
                                                                                                                       gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTTTCTTTAAAAAGCACGCCAGAAAACTAATGGCCGCCTGATCATCATCGCTCCGCGCCATACCCACACCGGGGCCCGGCGCGGAGATGAGTGGATCCCTATTCGTCCCGGGACAGAt                           >  1:557735/1‑150 (MQ=255)
                                                                                                                             cgAATGAGTGGCGGGGGGGGGGCTTTATAACTTTAAAAGGCACGCCAGAAAAATAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAAATGGGGGGCATCCCATTCTTCCGGGGAAAGATGCCGc                      >  2:684775/1‑149 (MQ=255)
                                                                                                                                        cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         >  1:456644/1‑151 (MQ=255)
                                                                                                                                        cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         >  1:1245598/1‑151 (MQ=255)
                                                                                                                                             ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg     <  1:891984/151‑1 (MQ=255)
                                                                                                                                                 tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:930289/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                         
ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862

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