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

GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTA  >  NC_000913/941565‑941850                                                                                                                                                     |                                                                                                                                          gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                         <  2:631152/151‑1 (MQ=255)    tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                      >  2:477076/1‑151 (MQ=255)        cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                                                                    <  2:146089/149‑1 (MQ=255)                 aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAtttt                                                                                                                              >  2:629246/1‑146 (MQ=255)                  aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAg                                                                                                                        <  1:213949/151‑1 (MQ=255)                       cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                   <  2:728609/151‑1 (MQ=255)                       cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTAtttt                                                                                                                              >  1:596353/1‑140 (MQ=255)                          aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTGTCTTTAAAAGGCACGcc                                                                                                                 >  1:422826/1‑150 (MQ=255)                          aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGcc                                                                                                                 >  2:10553/1‑150 (MQ=255)                                         ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                         <  1:72369/143‑1 (MQ=255)                                         ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg                                                                                                 <  2:332414/151‑1 (MQ=255)                                              ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCCTCCTATCTTGAACCGGCACACCAGAAATCTAATGctc                                                                                                  >  2:756163/1‑143 (MQ=255)                                                cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                          <  2:383683/151‑1 (MQ=255)                                                 ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTATTTCTTTAAAAGGCAAGCCCGAAATATAATGcccccca                                                                                              >  1:185088/1‑142 (MQ=255)                                                    tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTATTATCTTTAAACGGCACGCCAGCAAAATAATGTCCCCGTGattatc                                                                                       >  1:827657/1‑150 (MQ=255)                                                     gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTAATATCTTTAACAAGCACGCCAGAAATATAATGCCCGCACGatcctca                                                                                      >  1:39584/1‑150 (MQ=255)                                                             ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgc                                                                             <  1:596770/151‑1 (MQ=255)                                                               cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                           <  2:30003/151‑1 (MQ=255)                                                                      ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                                                                    <  2:393601/151‑1 (MQ=255)                                                                          cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGTTTACTATCTTGAAAAAGGACGCCAGAAAACTAATGCCCGCATGATTATCATTCATCCGCGCTATACCGccac                                                                >  1:161816/1‑151 (MQ=255)                                                                             atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg                                                             >  2:934003/1‑151 (MQ=255)                                                                             atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTAATATCTTTAAAAAGGACGCCAGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCCATACCGCCAcc                                                               >  1:676400/1‑149 (MQ=255)                                                                               atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATCCCGACACCGGTg                                                           <  2:936031/151‑1 (MQ=255)                                                                                tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATACCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGACCAGGCACGGCAGAAATCTAATGCCCGCATGATCATGATCGATCCGCGCTGTCCCGACACCGGTg                                                           <  1:26317/150‑1 (MQ=255)                                                                                     aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                     <  1:218657/151‑1 (MQ=255)                                                                                     aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGAATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                     <  2:378986/151‑1 (MQ=255)                                                                                       cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                   <  1:579652/151‑1 (MQ=255)                                                                                       cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCCGAAAACTAAATCCCGCCAGGTCATTCTTGAACCGCGCCATAACGACAACCGTGCCGGgcg                                                    >  1:835049/1‑150 (MQ=255)                                                                                           aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                               <  2:445121/151‑1 (MQ=255)                                                                                           aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                               <  1:702557/151‑1 (MQ=255)                                                                                                ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTCTCTTTAACCAGCACGCCAGAAATCTAATACCCGCCTGAGCATCATTCATCCGCGCTATAACGACACCGGGggcgggcgggga                                                >  2:417753/1‑141 (MQ=255)                                                                                                  tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGTg                                        >  2:724255/1‑151 (MQ=255)                                                                                                    gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                      <  1:161108/151‑1 (MQ=255)                                                                                                        tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCTATACCCACACCCGGGCCGGGCGCGAAGATGGGTGGATccc                                  >  1:829183/1‑151 (MQ=255)                                                                                                        tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAGAAAACTAATGCCCCCATGATCATCCTCGATCCGCGCTATACCGACACCCGGGGCGGGCGCGAAGATGAGTGGATccc                                  >  1:479956/1‑151 (MQ=255)                                                                                                         gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                 <  1:417748/151‑1 (MQ=255)                                                                                                         gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAAGCACGCCAGAAAACTAAAGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATccc                                  >  2:926570/1‑150 (MQ=255)                                                                                                          tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATccc                                  >  2:502564/1‑149 (MQ=255)                                                                                                            tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTAAAGCCCGCCTGATTATCATCCATCCGCGCTATAACCACACCCGTGCCGGGGGCGGAGAAGAGTGGATCCCTAtt                              >  2:940727/1‑151 (MQ=255)                                                                                                             ggTAATAACCCTGGCGAAGCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAACGCCCGCATGATCATCATCGATCCGCTCTATACCGCCACCGGTGCCGGGCGCGAAGATGAGTGGTTCCCTAtta                             <  1:259959/151‑2 (MQ=255)                                                                                                              gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTGAAAAGGCACGCCAGAAAACTAATGGCCGCATGATAATCATCGATCCCCGCTATACCGACACCCGGGCCCGGCGCGAAGATGAGTGGAACCCCAtt                              >  1:314971/1‑149 (MQ=255)                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                       >  2:808317/1‑151 (MQ=255)                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTGAAAAGGCACGCCACAAATCTAATGCCCCCCCGATAATCATTGATCCGCGCTATACCCACACCCGTGCCGGGCGCGAAGAAGAGTGGATCCCTATTCGTcc                         >  2:693811/1‑149 (MQ=255)                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCCCATGATCATTATCGATCCGCGCTATACCGACACCAGGGGCGGGGGCGAAGAAGAGTGGGACCCTATTCGTCCgg                       >  2:235401/1‑151 (MQ=255)                                                                                                                    aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCACGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGCGCGGAGATGAGTGGAACCCTATTCGTCCggg                      >  2:518849/1‑151 (MQ=255)                                                                                                                        tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGCGCGAAGATGAGTGGATCCCTTTTTGTCCcggga                    >  2:892006/1‑147 (MQ=255)                                                                                                                        tGGCGAAACGCGAATGAGTGGCGGCGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGCGCGAAGATGAGTGGATCCCTATTTGTCCGGGAAc                   >  1:118336/1‑150 (MQ=255)                                                                                                                            gAAACGCGAATGAGTGGCGGTGGGGTGGATTACTTTCTTTAAAAGGCCCGCCCGAAAACCAATGCCCGCATGATCATCATCGATCCGCGCCATAACCACACCCGCGCCGGGGGCGAAGATGAGGGGGTCCCCATTCCTCCGGGGACAGATg              >  1:73884/1‑151 (MQ=255)                                                                                                                            gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCTAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg              <  2:896625/151‑1 (MQ=255)                                                                                                                              aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc            <  2:52113/151‑1 (MQ=255)                                                                                                                              aaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACGCCAGAAAACTTATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGcc            >  1:504664/1‑151 (MQ=255)                                                                                                                                      tGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAGAAATCTAAAGCCCGCATGATAATCATCGATCCGCGCCATACCGACACCCGGGCCCGGGGCGGAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGCCGCCCTgg     >  1:287986/1‑150 (MQ=255)                                                                                                                                        aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa  <  1:767104/151‑1 (MQ=255)                                                                                                                                         gTGGCGGGGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATACCGACACCGGGGCCCGGCGCGAAGATGAGGGGGTCCCTATTCGTCCGGGGACAAATACCGCCCTGGtt   >  1:863387/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                          GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTA  >  NC_000913/941565‑941850

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