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%	72.6 / 24.0	51	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 (25/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

TCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGT  >  NC_000913/941574‑941855                                                                                                                                            |                                                                                                                                                 tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTTCt                                                                                                                                          >  1:107820/1‑148 (MQ=255)     cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCtt                                                                                                                                     >  1:19419/1‑149 (MQ=255)       cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                 <  2:199190/151‑1 (MQ=255)        aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTAACtttga                                                                                                                                  >  2:625915/1‑146 (MQ=255)              cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCCATCTTTaaaaag                                                                                                                              >  1:290509/1‑143 (MQ=255)                  aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg                                                                                                                      <  1:1138370/151‑1 (MQ=255)                  aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg                                                                                                                      <  1:94993/151‑1 (MQ=255)                   ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                     <  2:878057/151‑1 (MQ=255)                       tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                  <  1:892466/150‑1 (MQ=255)                        gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                 <  2:1021419/150‑1 (MQ=255)                        gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                 <  2:459544/150‑1 (MQ=255)                           cTATACCTACGGCGGCTGGGTAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTATTTTTTTAACAAGCACGCCCGAAAACTa                                                                                                               >  1:575948/1‑149 (MQ=255)                               aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATCTTTAACAGGCACGCCCGAAAACTAATGccc                                                                                                         >  2:24437/1‑151 (MQ=255)                                ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCCGAAAAATAATGccc                                                                                                         >  2:755824/1‑150 (MQ=255)                                ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAAAAAGCACGCCAGAAATATAATAcccc                                                                                                        >  1:268311/1‑150 (MQ=255)                                          ctggagaGATGGCGAGAGGCGGTCAGATATCGAATACAGTGGGATGGTGGTGCGGTTTTGTAATATCCCTGGCGAAACTGGAAGGAGAGGCGCTTGGGTGACGTACTATCTTGAACAGGCACGCCAGAAGGGTAACGCCCGCATT‑atca                                                                                               <  1:476530/143‑1 (MQ=255)                                            gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑atcatcat                                                                                           <  2:1030224/151‑1 (MQ=255)                                                 gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATc                                                                                      <  1:1040615/151‑1 (MQ=255)                                                   tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTATCTTGAACAGGCAAGCCAGAAATCTAATGCCCGCCCG‑AACATCACCGATcc                                                                                     >  2:607992/1‑150 (MQ=255)                                                       aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCt                                                                                <  2:937481/151‑1 (MQ=255)                                                         cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCtat                                                                              <  1:405805/151‑1 (MQ=255)                                                              cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCCCATG‑ATCATCATCGATCCGCGCTATAcc                                                                           >  2:592470/1‑149 (MQ=255)                                                               gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ac                                                                       <  2:157661/151‑1 (MQ=255)                                                                   gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCg                                                                   <  1:1079908/151‑1 (MQ=255)                                                                     tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTGACAGGCACGCCCGAAAACTAATGCCCGCCAG‑ATCATCATCGATCCGCGCTATAACC‑ACAcccg                                                                  >  2:226094/1‑148 (MQ=255)                                                                         gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCAAGCCAGAAATCTAATGCCCGCCTG‑ATAATCATCGATCCCCCCTATACCc‑acaaccggggc                                                              >  2:126586/1‑142 (MQ=255)                                                                          aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCgg                                                            <  1:827237/151‑1 (MQ=255)                                                                              cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTTTCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCTggggg                                                         >  2:975525/1‑148 (MQ=255)                                                                                    gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGAt                                                  <  2:590905/151‑1 (MQ=255)                                                                                     cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGGGCCGGGggcgga                                                     >  2:121130/1‑145 (MQ=255)                                                                                       ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAg                                               <  1:627717/151‑1 (MQ=255)                                                                                       ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCCGGAAACTAATGCCCGCCTG‑AATATCATCGATCCGCGCTATACCG‑CAACCGGTGCCggg                                                           >  2:1019386/1‑139 (MQ=255)                                                                                           gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGa                                           <  1:83384/151‑1 (MQ=255)                                                                                             gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTGAAAAAGCACGCCCGAAATATAAAGGCCGCCTG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGGGCCGGGCGCGAAGATGAgggg                                            >  1:1788/1‑148 (MQ=255)                                                                                                     gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                 <  1:1101043/151‑1 (MQ=255)                                                                                                       aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAAAAGGCACGCCAGAAAACTAAAGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCCATTCGTc                               >  2:832825/1‑151 (MQ=255)                                                                                                          aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                            <  1:835193/151‑1 (MQ=255)                                                                                                          aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAAGCACGCCAGAAATCTAATGGCCGCCGG‑ATTATTAACGATCCCCGCCATACCG‑ACACCGGTGCCGGGCGCGGAGAAGAGTGGGTCCCCTTTTGTCCgg                            >  2:532884/1‑151 (MQ=255)                                                                                                          aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCATG‑ATCATCATCCATCCCCGCTATACCG‑ACACCCGGGCCGGGGGCGAAGATGAGTGGATCCCCATTCGTCCgg                            >  2:624286/1‑151 (MQ=255)                                                                                                             ccTGGCGAAACGCGAATGAGTGGCGGCCGGGGGACCTTTTATCTCGAACAAGTACGCCTGAAATGTTCTGATCGCACGCACCATC‑TCCCTCCGCGCTGTTCCGCCCACCGGCACCGGG‑GCGAAGATGAGGGGATgccccatcgtc                               >  1:601739/1‑137 (MQ=255)                                                                                                                ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                      <  2:583185/151‑1 (MQ=255)                                                                                                                ggCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAAGCACGCCAGAAATCTAATGCCCCCCAT‑ATTATCATCGATCCCCGCCATACCG‑CCACCGGTGGCGGGCGGGAAGATGAGTGGATCCCCATTCGTcc                              >  2:506570/1‑143 (MQ=255)                                                                                                                      aCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTTGTCCGGGGACAGATGcc                 >  2:745358/1‑150 (MQ=255)                                                                                                                        gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa              <  2:822027/151‑1 (MQ=255)                                                                                                                        gcgAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGCCg                >  2:401795/1‑149 (MQ=255)                                                                                                                          gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc             <  1:992560/150‑1 (MQ=255)                                                                                                                           aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg           <  1:24435/151‑1 (MQ=255)                                                                                                                                  ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg    <  2:463275/151‑1 (MQ=255)                                                                                                                                   gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg   >  2:9514/1‑151 (MQ=255)                                                                                                                                    cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt  >  1:94615/1‑151 (MQ=255)                                                                                                                                    cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt  <  2:19421/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                 TCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATG‑ATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGT  >  NC_000913/941574‑941855

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