Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	51.7%	58.8 / 95.5	58	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (4/24);  new base G (30/0);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.60e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-10

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941569‑941859                                                                                                                                                 |                                                                                                                                                   actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                <  2:1097780/149‑1 (MQ=255)    cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACtt                                                                                                                                           <  1:156277/151‑2 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGTCTTACTATCt                                                                                                                                        >  2:267825/1‑151 (MQ=255)        cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                       <  1:227669/151‑1 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa                                                                                                                                    <  1:808682/151‑1 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCtt                                                                                                                                       >  2:661908/1‑148 (MQ=255)                    ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGc                                                                                                                           <  1:985223/151‑1 (MQ=255)                        ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                       <  1:427983/151‑1 (MQ=255)                             gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                   <  1:378362/150‑1 (MQ=255)                              aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                 <  2:126662/151‑1 (MQ=255)                              aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAAAAAGCACGCCCGaaaaat                                                                                                                  >  1:760690/1‑147 (MQ=255)                                     ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTGAACAAGGACGCCAGAAAACTAATGccc                                                                                                           >  2:1186652/1‑150 (MQ=255)                                        aCGGCGGCTGTGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGTCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                       <  2:181888/151‑1 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCCGAAAACTAATGCCCGCATg                                                                                                      >  2:482701/1‑151 (MQ=255)                                          ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCATg                                                                                                      >  2:321330/1‑150 (MQ=255)                                             ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGGACGCCAGAAAACTTATGCCCGc                                                                                                         >  1:718055/1‑144 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                              <  2:1183106/151‑1 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatc                                                                                                <  2:534491/149‑1 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAGAAATCTAATTCCCGCCTGatcatca                                                                                               >  1:817653/1‑150 (MQ=255)                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAAAAGGCCCGCCCGAAATATAATGCCCCCATGATCATCCTCGa                                                                                           >  2:1217474/1‑150 (MQ=255)                                                      gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAGCTAATGCCCGCATGATCATCATCGATc                                                                                         <  1:740749/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                    >  2:56487/1‑149 (MQ=255)                                                                gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                               <  1:653177/151‑1 (MQ=255)                                                                  ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAATCTAATTCCCGCATGATCATTCTCGATCCgcgccct                                                                                 >  1:1190551/1‑144 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCAGAAAAATAAAGCCCGCCTGATCATCATCGAGCCGCGCCTTACCGa                                                                            >  1:762404/1‑148 (MQ=255)                                                                              gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCCCCTGATCATCATCGATCCGCGCCATACCGACACCGGGGcc                                                                  >  2:205443/1‑150 (MQ=255)                                                                               aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                <  1:532989/151‑1 (MQ=255)                                                                                 aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                              <  1:447327/151‑1 (MQ=255)                                                                                  aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                             <  1:129485/151‑1 (MQ=255)                                                                                    aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                           <  2:681356/151‑1 (MQ=255)                                                                                        aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                       <  1:773543/151‑1 (MQ=255)                                                                                         gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCACGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGa                                                          >  2:244993/1‑147 (MQ=255)                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGGGCGGAGATGAgggg                                                >  2:415264/1‑150 (MQ=255)                                                                                                 tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTATTTTCTTTAACAGGCACGCCAAAAATCTAATGCCCGCATGAACAACATCCATTCGCGCTATACCGACACCCGGGTCGGGCGCGGAAAAGAGTgg                                                >  2:603362/1‑149 (MQ=255)                                                                                                 tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGCGCGAAGATGAGTgg                                                >  2:642727/1‑149 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCCGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGGGCGGAGATGAGTgg                                                >  2:602831/1‑148 (MQ=255)                                                                                                     gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATAAGTGGATcc                                            <  1:148016/149‑1 (MQ=255)                                                                                                       ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATAACGACACCCGGGCCGGGCGCGAAGATGAGTGGATccc                                           >  2:271251/1‑148 (MQ=255)                                                                                                        tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAATCTAATACCCCCCAGATCATCATCCATCCGCGCTATAACGACACCGGGGCCCGGCGCGGAGATGAGGGGGTCCCCAtt                                       >  1:16431/1‑151 (MQ=255)                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                       <  2:68335/150‑1 (MQ=255)                                                                                                           taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGTGGGTCCCTAtt                                       >  2:149041/1‑148 (MQ=255)                                                                                                             ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCCGGGCCGGGGGCGAAGATGAGGGGGGCCCCATTCGTcc                                  >  2:626885/1‑151 (MQ=255)                                                                                                              taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAATCTAAATCCCGCCTGATCATCATCGATCCCCGCTATAACGACACCCGGGCCGGGCGCGGAGATGAGGGGGTCCCCATTTGTcc                                  >  2:933988/1‑150 (MQ=255)                                                                                                              taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGCCACCGGTGCCGGGCGCGAAAATGAGGGGGTCCCCCTTCGTcc                                  >  1:329561/1‑150 (MQ=255)                                                                                                              taaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTTCTATCTTGAAAACGCCCGCCAGAAAAATAATGCCCCCATGATCATCATTGATCCGCCCTATACCCACAACGGTGCCCGGCGCGCAGATGAGGGGAGCGCTATTcct                                    >  1:1215105/1‑146 (MQ=255)                                                                                                                 cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAGGGACGCCAGAAAATTAATGCCCGCGTGATCATCCTCGATCCGCGCCATACCCACCCCCGGGCCGGGCGCGAAGATGAGGGGGACCCCATTCGCCCggg                               >  1:95070/1‑150 (MQ=255)                                                                                                                 cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCCATCCGCGCTATACCGACACCGGTGCCGGGGGCGGAGATGAGTGGATCCCCATTCGTCCggg                               >  1:420686/1‑150 (MQ=255)                                                                                                                 cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                               >  2:1221018/1‑150 (MQ=255)                                                                                                                  ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCCCTATACCGACACCGGGGCCGGGCGCGAAGATGAGGGGATCCCTCTTCGTCCggg                               >  1:163805/1‑149 (MQ=255)                                                                                                                    tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCgggga                             >  2:193798/1‑147 (MQ=255)                                                                                                                        gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                       <  1:973597/151‑1 (MQ=255)                                                                                                                        gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTTATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGGGCCGGGGGCGGAGATGAGTGGGTCCCCATTTGTCCCGGAACAGa                         >  1:1122881/1‑149 (MQ=255)                                                                                                                            cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATAAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                   <  2:809280/151‑1 (MQ=255)                                                                                                                                    aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa           <  1:785127/151‑1 (MQ=255)                                                                                                                                      tGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc         >  1:1004508/1‑151 (MQ=255)                                                                                                                                         cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt      >  2:985385/1‑151 (MQ=255)                                                                                                                                             ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg  >  1:1117845/1‑151 (MQ=255)                                                                                                                                             ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg   <  2:190282/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941569‑941859

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