Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	40.0%	76.6 / 18.8	50	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (2/28);  new base G (20/0);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 4.90e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862                                                                                                                                                        |                                                                                                                                                      gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGt                                                                                                                                                       <  1:1003495/151‑1 (MQ=255)       actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGACGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCGGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACtt                                                                                                                                                  <  2:1257327/150‑1 (MQ=255)             ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCt                                                                                                                                           <  2:686282/151‑1 (MQ=255)             ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTACTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATc                                                                                                                                            <  2:541815/150‑1 (MQ=255)             acgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGAGGCTGGGCAGGTGGCAACCGCCCGTGGGGAATCGAAAAAAGTAGGCTGGTAGTGCTGTCTGGTAATAACCCTGGCGAAGCACGAATGAGTGGCGGTTGGGGTgggtagt                                                                                                                                                >  2:135429/2‑141 (MQ=255)              cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCtt                                                                                                                                          <  2:1155778/151‑1 (MQ=255)              cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCtt                                                                                                                                          <  1:669051/151‑1 (MQ=255)                             aGGCTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGGAATAACCCCGGCGAAACGCGAATGAGTGGCGG‑CGGGGGGATGTACATTCTTGaaaaagcta                                                                                                                                >  2:398240/1‑139 (MQ=255)                              ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                          <  2:479166/151‑1 (MQ=255)                                tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGCTTACTATCTTGAAAAGGGACGCCCGaa                                                                                                                         >  1:258531/1‑150 (MQ=255)                                      cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGAGTTACTATCTTTAAAAAGCACGCCCGAAAtat                                                                                                                     >  1:907670/1‑146 (MQ=255)                                         taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGGTTTCTATTTTTTAAAAGGACGCCCGAAAACTaaatcc                                                                                                               >  2:47396/1‑147 (MQ=255)                                          aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTCCATTTTTTAAAAGCCCGCCAGaaaaataa                                                                                                                   >  2:974686/1‑141 (MQ=255)                                               cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                         <  2:1131556/151‑1 (MQ=255)                                                  cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                      <  1:1048900/151‑1 (MQ=255)                                                   ggcTGGGCAGATGGCAACAGCCCGTCGGATGTCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                     <  1:695330/151‑1 (MQ=255)                                                   ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                     <  1:483392/151‑1 (MQ=255)                                                    gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                                    <  1:650322/151‑1 (MQ=255)                                                         gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                               <  2:749564/151‑1 (MQ=255)                                                           aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                             <  1:516973/151‑1 (MQ=255)                                                           aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTTATATCTTTGAAAGGCACGCCAGAAAACCAATGCCCCCATGATCCTCTTCGa                                                                                              >  1:430107/1‑150 (MQ=255)                                                           aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGGCCGCCTGATCATCATCGa                                                                                              >  1:380468/1‑150 (MQ=255)                                                                gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                        <  1:1066170/151‑1 (MQ=255)                                                                 cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGGTTACCATTTTGAACAGGCAATCCACGAAAATAAAATCCGCATGGTCAACATCGATTcgcg                                                                                        >  2:990092/1‑150 (MQ=255)                                                                     aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                   <  2:1096941/151‑1 (MQ=255)                                                                      gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGGTTATTATCTTTAAAAAGGCCCCCAGAAAATTAATACCCGCCCAATAATAATCGAACCccgcctta                                                                                   >  2:125335/1‑146 (MQ=255)                                                                        ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                                                                                <  1:1232539/151‑1 (MQ=255)                                                                          gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑ATGGGGGGCTTTCTATCCTTAAAAGGCACCGCAAAAAACTAAATCCCCCCGGATCATCATTGATCCGCCCGATACCAac                                                                              >  1:340456/1‑151 (MQ=255)                                                                                 atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTATCTTTAACAAGCCCGCCCGAAAACTAATGCCCGCAAGATCATCATCGATCCCCGCCATACCCacac                                                                            >  2:1076283/1‑146 (MQ=255)                                                                                   cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGCTTTCTATCTTTAACAAGCCCGCCAAAAAACTAAAGCCCGCCCGATCACCATCGATTCGCCCTATACCCACAcc                                                                           >  2:1261908/1‑145 (MQ=255)                                                                                       aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTTCTATCTTTAACAGGGACGCCAGAAAACTAATGCCCCCCTGATCATCATCGATACGCGCTATACCGACACCGGGGCCggg                                                                  >  1:781507/1‑150 (MQ=255)                                                                                       aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                                 <  2:868994/151‑1 (MQ=255)                                                                                              aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                          <  1:433720/151‑1 (MQ=255)                                                                                                  ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                       <  1:831801/150‑1 (MQ=255)                                                                                                   gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                     <  1:249755/151‑1 (MQ=255)                                                                                                   gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGGCTTTCTATCTTTAACAAGCACGCCCGAAAACTAATGCCCGCATGATCATCAACGAACCGCGCCATAACCACACCCGGGGCGGGGGCGGAGATggg                                                      >  2:1008551/1‑148 (MQ=255)                                                                                                     aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGTgg                                                   >  2:1132431/1‑151 (MQ=255)                                                                                                         cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGATTTATATTTTGAAAAAGGACACCAGAAAACTAAATCCCGCACGAACATCAATCATCCGCGCTATACCGACACCGGGGCCGGGCGGGaaagagaggtgggcc                                               >  2:462445/1‑138 (MQ=255)                                                                                                         cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑GGGGGGGAAATACTCTATTTAAAAGGCCCGCACCAAAACTAAAATCCGCAAGAATCTCATTGCTCCCCCCAATACCCACACACGGGCCGGGGGCGgagagttggggatcc                                               >  1:1135807/1‑140 (MQ=255)                                                                                                           gtctGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCGCACCGGTGCCCGGGGCGAAGATGAGTGGGTccc                                              >  1:766717/4‑150 (MQ=255)                                                                                                           gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATAACCACACCGGGGCCGGGCGCGAAGAAGAGTGGGTccc                                              >  1:339371/1‑150 (MQ=255)                                                                                                              tGGTAATAACCCTGGCGAAACGCGAATGAGTGGTGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  2:325975/150‑1 (MQ=255)                                                                                                               ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                          <  2:461946/150‑1 (MQ=255)                                                                                                                 taataaCCCTGGCGAAACGCGAATGAGTGGCGG‑CAGGGGGAGCTACTATCTTGAACACGCGCGGCCGAAAAATAATGCCCCCGCGATCATAATCGAACCCCGCTTAACCGGAACACGGGCCGGGGCCGCAGATGAGAGGGCCCCTATTc                                         >  2:183399/1‑149 (MQ=255)                                                                                                                  aataaCCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                      <  1:33245/151‑1 (MQ=255)                                                                                                                   ataaCCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                     <  2:75741/151‑1 (MQ=255)                                                                                                                     aaCCCTGGCGAAACGCTAATGAGTGGCGG‑TGGGGTGACTTACTGTCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCATAGTCGTcctg                                   <  2:8035/151‑3 (MQ=255)                                                                                                                      aCCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                  <  2:1116551/151‑1 (MQ=255)                                                                                                                        ccTGGCGAAACGCGAATGAGTGGCGG‑TGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCCATTCGTCCggg                                  >  1:127569/1‑149 (MQ=255)                                                                                                                                aaCGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                        <  2:364787/151‑1 (MQ=255)                                                                                                                                   gcgAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                     <  1:540645/151‑1 (MQ=255)                                                                                                                                                        tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:676602/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGG‑TGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862

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: ‑