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%	61.8 / 23.6	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

GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941565‑941860                                                                                                                                                     |                                                                                                                                                    gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                   <  2:473530/151‑1 (MQ=255)      ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGACtt                                                                                                                                                >  1:1028732/1‑149 (MQ=255)           ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTTCTa                                                                                                                                            >  1:1000390/1‑148 (MQ=255)            cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGTGTTACTAt                                                                                                                                           >  2:222911/1‑148 (MQ=255)              cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                        <  2:969631/149‑1 (MQ=255)                cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa                                                                                                                                     <  1:462814/150‑1 (MQ=255)                   aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                 <  1:333936/151‑1 (MQ=255)                      gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                               <  2:495425/150‑1 (MQ=255)                      gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAGTTTCAATTTTTAACAgg                                                                                                                                 >  2:670389/1‑148 (MQ=255)                          aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATTTTTAAAAAACACGcc                                                                                                                           >  2:276971/1‑150 (MQ=255)                          aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCTTACTATATTTAAAAgggaag                                                                                                                             >  2:243951/1‑144 (MQ=255)                             gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGCACGcc                                                                                                                           >  1:609388/1‑147 (MQ=255)                             gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGcc                                                                                                                           >  1:82521/1‑147 (MQ=255)                              tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTTAAAAGGCACGCCAGaa                                                                                                                       >  1:1062258/1‑150 (MQ=255)                                    cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGGACACGCACGCCAGAAAACTaa                                                                                                                 >  2:915790/1‑150 (MQ=255)                                      ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                              <  2:492871/151‑1 (MQ=255)                                       taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAAGCACGCCCGAAATCTAATTcc                                                                                                             >  1:874453/1‑151 (MQ=255)                                                  gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTGGTACTATTATTAACAAGCACGCCCGAAAACCTATGCCCGCACGagcat                                                                                                  >  2:93324/1‑151 (MQ=255)                                                      ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                              <  1:279288/151‑1 (MQ=255)                                                      ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGatcatca                                                                                                >  1:53213/1‑149 (MQ=255)                                                         aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                                  >  1:854772/1‑144 (MQ=255)                                                            tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTAATATTTTGAAAAGGCACGCCAGAAATCTAAAGCCCGCATGATCATCATCGATcc                                                                                         >  1:646093/1‑150 (MQ=255)                                                            tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTTCTATCTTGAACAGGCACGCCAGAAAACTAATGGCCGCCTGATCCTCATCGa                                                                                            >  2:534518/1‑147 (MQ=255)                                                                    gCCCGGCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATATCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCCCTCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCtattc                                                                                <  2:916258/151‑3 (MQ=255)                                                                     cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAATTAATATCTTGAACAGGCACGCCCGAAACCTAATGCCCGCCTGATCCACATCCATACCCGGTATAcc                                                                               >  2:469157/1‑151 (MQ=255)                                                                         tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCCCATTATCATCATCGATCCGCGCTATACCg                                                                              >  1:337408/1‑148 (MQ=255)                                                                           ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                         <  1:49837/151‑1 (MQ=255)                                                                           ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                         <  1:709219/151‑1 (MQ=255)                                                                            gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                        <  2:708713/151‑1 (MQ=255)                                                                                  gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATTCGCGCTATCCCGACACCGGTGcc                                                                   >  2:596277/1‑150 (MQ=255)                                                                                     aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                               <  1:898259/151‑1 (MQ=255)                                                                                        aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAAGATAATCCTCGATCCGCGCTATACCGACACCGGGGCCGGGggcg                                                            >  1:305126/1‑151 (MQ=255)                                                                                          tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTATTTTTTTTAACAAGCACGCCAGAAATCAAATGCCCGCATGAATAACATCGGTTCGCGCCATAACGACAACGGGGGCGGgcgc                                                             >  2:957808/1‑148 (MQ=255)                                                                                           aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                          <  1:18401/150‑1 (MQ=255)                                                                                             gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                       <  2:383307/151‑1 (MQ=255)                                                                                                 gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCCATACCCACACCCGGGCCGGGCGCGGAGATGAg                                                    >  1:642570/1‑150 (MQ=255)                                                                                                      gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                               <  1:202689/150‑1 (MQ=255)                                                                                                         gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                           <  2:305050/151‑1 (MQ=255)                                                                                                              gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                      <  2:967224/151‑1 (MQ=255)                                                                                                               taataaCCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACTTACTATCTTTAAAACGGCCGCCAGAAATCTAATGCCCGCCTGATCATCATCGAACCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGGGGGTACCTTTTCGt                                     >  1:707225/1‑151 (MQ=255)                                                                                                                  taaCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTAATAACTTCAAAAGAAACCGCAGAAAACTAATGGCCGCACGATTCTCCTTGATCCGCGGCCCTCCGACACCCGTGCCGGGCCGGAAGAAGAGTGGATGCCTACTccgcc                                   >  2:172857/1‑146 (MQ=255)                                                                                                                   aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGGGCGAAGATGTGTGGATCCCCATTTGTCCgg                                 >  1:792133/1‑151 (MQ=255)                                                                                                                    aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                <  1:1007568/151‑1 (MQ=255)                                                                                                                       cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCCgggga                              >  1:986384/1‑148 (MQ=255)                                                                                                                         ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                             <  1:936770/149‑1 (MQ=255)                                                                                                                           cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTTTCTTGAAAAGGCACGCCCGAAATCTAATTCCCGCCTGATCATCATCGATCCCCGCCATACCCACACCGGGGGCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCgggga                              >  1:318720/1‑144 (MQ=255)                                                                                                                             aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc                       <  2:609182/151‑1 (MQ=255)                                                                                                                              aaCGCGAATGAGTGGCGGTCGGATGACTTACTCTCTTGAACAGGCCCGCCAGAAATCTAATGCCCGCATGACGATCATGGATCCACGTTATGCCGTCCCCGGTGCCGGGAGCGAATATGAGTGGATCGTTATTTGTCCgggg                               <  2:685561/142‑2 (MQ=255)                                                                                                                               aCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCCGAAATCTAATGCCCCCCTGATCATCATCGATCCGCGCTATACCCCCACCGGGGCCGGGCGCGGAGATGAGTGGAACCCTATTCGTCCGGGGACAGATGcc                      >  1:967762/1‑150 (MQ=255)                                                                                                                                   gAATGAGTGGCGGGGGGGTGAATTATTATCTTGAACAAGCACACCAGAAAACTAAAGCCCGCATGATTATCATCGATCCGCGCTATACCGAAACCGGCGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCggggaacgaaggcgc                    >  2:348073/1‑136 (MQ=255)                                                                                                                                    aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg                <  2:378574/151‑1 (MQ=255)                                                                                                                                      tGAGTGGCGGTTGGGGGTCTTACTATCTTTAACAAGCACGCCAGAAAACTAATGGCCGCCCGATCATCATCGATCCGCCCTATACCGACACCCGTGCCGGGCGCGAAGATGAGTGGGGCCCCATTCGTCCCGGAACAGATGCCGCACTgg               >  2:691066/1‑150 (MQ=255)                                                                                                                                       gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt             <  1:754905/151‑1 (MQ=255)                                                                                                                                                  gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  >  2:81723/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941565‑941860

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