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

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862                                                                                                                                                        |                                                                                                                                                      gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGt                                                                                                                                                       <  1:587481/151‑1 (MQ=255)   gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGa                                                                                                                                                     <  1:1114818/151‑1 (MQ=255)           ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACATAcca                                                                                                                                              >  2:639536/1‑148 (MQ=255)             ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTAt                                                                                                                                             <  2:346752/149‑1 (MQ=255)              cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCtt                                                                                                                                          <  1:1173447/151‑1 (MQ=255)               gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTg                                                                                                                                         <  1:564885/151‑1 (MQ=255)                cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGa                                                                                                                                        <  2:536562/151‑1 (MQ=255)                      ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGGTTACCATTTTGaaaaa                                                                                                                                    >  1:1247533/1‑146 (MQ=255)                         cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                               <  2:269367/151‑1 (MQ=255)                            aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTCAGCTGGTAGGCCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGCGGGGTGAGGTAATATcatttagaggaatcg                                                                                                                              >  1:31165/1‑135 (MQ=255)                              ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCGTACTATCTTTAAAAgg                                                                                                                                   >  2:464132/1‑142 (MQ=255)                                  tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                        <  1:230140/149‑1 (MQ=255)                                     aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCCGaaa                                                                                                                        >  2:161946/1‑146 (MQ=255)                                       tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGAGTTACTATCTTTAACAGGGACCCCAGCAAACTaa                                                                                                                   >  1:1107983/1‑149 (MQ=255)                                           ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGACTT‑CTATCTTTAAAAGGCACGCCACAAAACTAATGccc                                                                                                              >  1:117892/1‑150 (MQ=255)                                              aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                          <  2:907476/151‑1 (MQ=255)                                                ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAAAAAGGACGCCAGAAAACTAATGCCCGc                                                                                                            >  2:433322/1‑147 (MQ=255)                                                      tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACTATCTTGAACAAGCACGCCAGAAAAATAATGGCCGCCTGatgatc                                                                                                   >  1:552750/1‑150 (MQ=255)                                                       gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                                 <  1:1099701/151‑1 (MQ=255)                                                       gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                  <  1:759652/150‑1 (MQ=255)                                                        ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACTCGAATGAGTGGC‑GGTGGGGGGGCCTACTATTTTTAACACGGACGCCAGAAATCTAAAGGCCCCATGatcat                                                                                                    >  2:403250/1‑147 (MQ=255)                                                        ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGTGGGGGGCTTACTATTTTCAACAAGCACCGGAGAAATATAAAAACCGCACGatcatcat                                                                                                 >  2:1048391/1‑150 (MQ=255)                                                         gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                               <  1:815423/151‑1 (MQ=255)                                                           aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                             <  2:584470/151‑1 (MQ=255)                                                                       cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 <  1:122505/151‑1 (MQ=255)                                                                         cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGa                                                                               <  2:380180/151‑1 (MQ=255)                                                                         cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTAATATCTTGAAAAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 >  1:1045389/1‑149 (MQ=255)                                                                              gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                          <  2:815236/151‑1 (MQ=255)                                                                                 atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                        <  1:870089/150‑1 (MQ=255)                                                                                    gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGGGcc                                                                     >  1:543577/1‑150 (MQ=255)                                                                                      aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg                                                                  <  1:363498/151‑1 (MQ=255)                                                                                      aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTTTTATCTTGAAAAAGGACGCCAGAAATCTAATTCCCGCCAGATCATCCACGATCCGCGCTATACCCACACCCGGGGCggg                                                                  >  1:1020470/1‑151 (MQ=255)                                                                                       aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTCCTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCACGATTATCATCGATCCGCGCCATACCGACCCCCGTGCCggg                                                                  >  2:1158235/1‑150 (MQ=255)                                                                                                    tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTTCTTTCTTTTAAAAGCACGCCCGAAAACTAATAGCCGCCCGGTCAACAACGATCCGCGCTATACCCCCACCGGGGCCGGGGGCGGAGATGAGTg                                                    >  1:342939/1‑151 (MQ=255)                                                                                                     aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTgg                                                   <  1:223533/151‑1 (MQ=255)                                                                                                      gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGTgg                                                   >  2:353666/1‑150 (MQ=255)                                                                                                        gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                <  2:77155/151‑1 (MQ=255)                                                                                                                  aataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                      <  2:489928/151‑1 (MQ=255)                                                                                                                   ataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACCTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCAAGGTCATCATCGATACCGGCTATACCGAAACCCGTGCCGGGGGGGAAGGAGGGGGGATGCCTCTTc                                         >  1:644129/1‑147 (MQ=255)                                                                                                                     aaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                                   <  2:411414/151‑1 (MQ=255)                                                                                                                     aaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGCGCGCAGATGAGTGGATCCCTATTTCTCCgg                                   >  2:1238619/1‑151 (MQ=255)                                                                                                                     aaCCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCTTACTACCTTGAACCGGCACGCCAGAAATATAATGCCCGCACGGTTATCATCGGTCCGCGCATTACCCACACCGGGGCCCGGGGCGGAGAAGGGAGGAACCCTATTCGGCCg                                    >  1:357097/1‑150 (MQ=255)                                                                                                                       cccTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGGCTTTCTATCTTGAAAAGGCACGCCCGAAATCTCATGCCCGCATGATCATCATCGATCCCCGCCAAACCCACACCGGTGCCCGGCGCGGAAATGAGGGGATCCCCATTCGTCCcgg                                  >  2:935693/1‑150 (MQ=255)                                                                                                                       cccTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                                 <  1:797501/151‑1 (MQ=255)                                                                                                                        ccTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGaa                                <  2:552470/151‑1 (MQ=255)                                                                                                                         atGGCGAAACGCGAATGAGTGGC‑GGTGGGGCGGGTTACTTTCTTTGACAGGGACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCCTTACCGACAACAGGGGCGGGGGCGGAGATGGGGGGGGCCCTCTTCGTACggg                                  >  1:1036115/2‑148 (MQ=255)                                                                                                                           ggCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                              <  2:112668/150‑1 (MQ=255)                                                                                                                                aaCGCGAATGAGTGGC‑GGTGGGGGGAACTACTATCTTTAAAAAGGACCCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACCACACCGGGGCCGGGGGCGGAAATGAGTGGGTCCCCATTCGTCCCGGGACAGATGcc                        >  1:818700/1‑151 (MQ=255)                                                                                                                                   gcgAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                     <  1:144815/151‑1 (MQ=255)                                                                                                                                   gcgAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAACAAGCACGCCCGAAAACTAAAGCCCGCCTGATCATCATCGCACCCCGCTATAACGACACCGGTGCCGGGGGCGAAGATGAGTGGATCCCCATTCGTCCCGGGACAGa                            >  2:1217218/1‑144 (MQ=255)                                                                                                                                      aaTGAGTGGC‑GGTGGGGGGGCTTACTATCTTGGACAGGCACCCCAAAAATCTAATGCCCGCATGATCATCATCGATCCCCGCCATAACGCCACCGGGGCCCGGCGCGGAGAAGAGGGGATACCTCTTTCTCCCGGGACAGATGcccc                      >  1:1114270/1‑145 (MQ=255)                                                                                                                                           gTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGCACTGGtt               >  2:391598/1‑149 (MQ=255)                                                                                                                                                     gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc    <  2:940240/151‑1 (MQ=255)                                                                                                                                                        tGACTTACTATCTTGAACAGTCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGTCGt  >  2:550236/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  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 < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑