Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	45.5%	78.1 / 12.4	55	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (7/23);  new base G (25/0);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941565‑941857                                                                                                                                                     |                                                                                                                                                  gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                 <  2:322938/151‑1 (MQ=255)           ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGTCAGTTGCATTACTTACTATCt                                                                                                                                       <  1:84699/151‑1 (MQ=255)                cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaa                                                                                                                                   >  2:819121/1‑150 (MQ=255)                   aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTaaaa                                                                                                                                 >  2:916594/1‑147 (MQ=255)                    ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTTAACAgg                                                                                                                               >  2:734695/1‑150 (MQ=255)                      gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                            <  1:851334/151‑1 (MQ=255)                       cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTgtaaaagg                                                                                                                              >  1:398459/1‑141 (MQ=255)                       cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCTTACTATCTTGAACAGGCAc                                                                                                                            >  1:316269/1‑150 (MQ=255)                          aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATTTCTTTAACAGGCACGcc                                                                                                                         >  1:216788/1‑150 (MQ=255)                                 gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGCACGCCAGaaa                                                                                                                    >  1:599249/1‑148 (MQ=255)                                  aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCt                                                                                                                 >  2:99106/1‑150 (MQ=255)                                   aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                               <  1:85359/151‑1 (MQ=255)                                      ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                                                                                              <  1:557723/149‑1 (MQ=255)                                        aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGccc                                                                                                          <  2:406771/151‑1 (MQ=255)                                            aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                        <  2:623334/149‑1 (MQ=255)                                            aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                      <  2:177648/151‑1 (MQ=255)                                              ggcggcTGGGCAGATGGCAACAGCCCGTCGTATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                                                                                    <  1:213570/151‑1 (MQ=255)                                                 ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                 <  1:23170/151‑1 (MQ=255)                                                 ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                 <  1:612010/151‑1 (MQ=255)                                                     gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                             <  2:377456/151‑1 (MQ=255)                                                      ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAAGGCCGCCCAAAAACTAATTCCCCCATg                                                                                                     >  1:426535/1‑142 (MQ=255)                                                         aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCggt                                                                                         >  2:561045/1‑149 (MQ=255)                                                             ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                                     <  2:19505/151‑1 (MQ=255)                                                              gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                    <  1:275719/151‑1 (MQ=255)                                                                aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTTAACAAGGACGCCCGAAAACTAATGCCCGCATGATTATCCTCGAACCgcgc                                                                                   >  1:821482/1‑150 (MQ=255)                                                                  cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                  >  1:30667/1‑149 (MQ=255)                                                                  cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                   >  1:239603/1‑148 (MQ=255)                                                                    gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                              <  2:867303/151‑1 (MQ=255)                                                                                tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCCGTg                                                                  >  1:690080/1‑150 (MQ=255)                                                                                  gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTAATATCTTGAACAAGCACGGCCGAAAACTAATGCCCGCCAGATCATCATCGATTCGCGCTATACCGA‑cccccg                                                                    >  1:826952/1‑144 (MQ=255)                                                                                   aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCgg                                                              >  2:843913/1‑151 (MQ=255)                                                                                       cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGgcgc                                                          <  1:369893/151‑1 (MQ=255)                                                                                          tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGcgcg                                                         >  1:856991/1‑149 (MQ=255)                                                                                               tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCAGTGCCGGGCGCGAAGATg                                                   >  1:68835/1‑150 (MQ=255)                                                                                                ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCCGAAAAATAATGCCCGCCTGATCATCATCGATCCGCGCCATAACGG‑CACCGGTGCCGGGGGCGGAGATg                                                   >  1:80636/1‑149 (MQ=255)                                                                                                 gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAg                                                 >  2:428313/1‑150 (MQ=255)                                                                                                     tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGA‑CACCGGTGCCGGGGGCGAAGATGAGTgg                                              >  2:14136/1‑149 (MQ=255)                                                                                                      gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGAt                                            <  1:22942/150‑1 (MQ=255)                                                                                                      gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATc                                           <  2:63206/151‑1 (MQ=255)                                                                                                            tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                     >  2:874169/1‑151 (MQ=255)                                                                                                               taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                     >  2:826473/1‑148 (MQ=255)                                                                                                                  taaCCCTGGCGAAACGCGAATGAGTGGCGGAGGGGTGGCGTTCTACCATGCAAAGGCCCGCAACAATTCCCAGGCCAGCGTGATCATTATAGATCCGCGCAGTACCCATAACCAGGTCCGGGC‑CGAAGATCGGTGGACGCCTTCTCGTcc                                >  1:621496/1‑150 (MQ=255)                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                              <  1:480286/151‑1 (MQ=255)                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATTCCCGCTATACCGA‑CACCGGGGCCGGGGGCGAAGATGAGGGGGTCCCTATTCGTCCgg                              >  2:536059/1‑151 (MQ=255)                                                                                                                      ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                             >  2:903187/1‑149 (MQ=255)                                                                                                                        tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGA‑CACCGGGGCCGGGGGCGAAGATGAGTGGATCCCCATTCGTCCGGGGAc                          >  2:279535/1‑150 (MQ=255)                                                                                                                        tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCA‑CACCGGTGCCGGGGGCGAAGATGAGGGGATCCCTATTCGTCCGGGaa                           >  2:91304/1‑149 (MQ=255)                                                                                                                           cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                      <  2:341367/151‑1 (MQ=255)                                                                                                                           cGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCCATCCGCGCTATACCGA‑CACCGGTGCCGGGGGCGAAGATGAGGGGATCCCTATTCGTCCgggga                           >  2:182588/1‑144 (MQ=255)                                                                                                                                 gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGA‑CACCCGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGCCGc                 >  2:461023/1‑150 (MQ=255)                                                                                                                                  cgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc               <  1:706848/151‑1 (MQ=255)                                                                                                                                   gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc               >  2:39435/1‑150 (MQ=255)                                                                                                                                       gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt          <  2:570155/151‑1 (MQ=255)                                                                                                                                          tGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc       >  1:811761/1‑151 (MQ=255)                                                                                                                                               gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  <  2:426574/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGA‑CACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941565‑941857

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