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

CCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941575‑941858                                                                                                                                           |                                                                                                                                                   cggcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCGGGTCAGATGGCAACAGCTCGTCGGATATCGAAATCAGGAACCTGGTAGTGGTGTGTGGGAATAACCCTGGCGAAACGCGAATGAGTAGCGGTGGGGTgtttt                                                                                                                                               <  2:976847/142‑5 (MQ=255) ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                        <  2:265290/151‑1 (MQ=255)   gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTg                                                                                                                                      <  2:80788/151‑1 (MQ=255)   gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCtt                                                                                                                                       >  1:556895/1‑150 (MQ=255)   gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGACGGCAACAGCCCGTCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCCGGCGAAACGCGAATGAGTGGCGGTTGGGTgtttttctctt                                                                                                                                         >  2:533440/1‑138 (MQ=255)    cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTTTCtt                                                                                                                                       >  2:868600/1‑149 (MQ=255)           tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATTTTTAACAAGc                                                                                                                               >  2:1045321/1‑150 (MQ=255)             cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTGAACAAGCAc                                                                                                                             >  2:299686/1‑150 (MQ=255)                 aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAg                                                                                                                        <  1:394570/151‑1 (MQ=255)                 aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGAACAGGCACGcc                                                                                                                          >  1:94251/1‑149 (MQ=255)                      tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                    <  2:684501/150‑1 (MQ=255)                      tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATATTGGACAAGCACGCCCGaa                                                                                                                      >  2:518887/1‑148 (MQ=255)                          cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTATTATCTTTAAAAAGCGCGCCCGAAATCTaa                                                                                                                >  2:1027725/1‑150 (MQ=255)                             taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTTAAAAAGCAAGCCCGAAAACTaa                                                                                                                >  2:1047031/1‑147 (MQ=255)                               ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg                                                                                                          <  1:376915/151‑1 (MQ=255)                                  aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                       <  2:941105/151‑1 (MQ=255)                                   cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                       <  2:93955/150‑1 (MQ=255)                                     gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCCGAAAACTAATTccccc                                                                                                         >  1:455107/1‑144 (MQ=255)                                                 aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATcc                                                                                       <  2:1021002/151‑1 (MQ=255)                                                    gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATTTTGAACAGGCACGCCAGAAACCTAATGCCCGCATGATCAT‑CATCGATCcgcg                                                                                    <  2:690949/151‑1 (MQ=255)                                                         aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCtat                                                                                >  1:972276/1‑150 (MQ=255)                                                             cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAGGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGa                                                                           <  2:10651/151‑1 (MQ=255)                                                             cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAAACTAATGCCCCCATGATCAT‑CATCGATCCGCGCTATAcc                                                                             >  1:509973/1‑149 (MQ=255)                                                                     atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCAAGATCAT‑CATCGATCCGCGCTATACCGACACCgg                                                                     >  1:842948/1‑149 (MQ=255)                                                                      tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCgg                                                                     >  1:634591/1‑148 (MQ=255)                                                                       cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCTTACTATCTTGAATCAGCACCCCACAAAATTAATGGCCCGAAGACCAT‑TATCGAACCGCCCCATAGCCACAcccggg                                                                   >  2:104948/1‑145 (MQ=255)                                                                        gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGTCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCg                                                                <  2:753514/151‑1 (MQ=255)                                                                         aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCgg                                                               <  1:866125/151‑1 (MQ=255)                                                                         aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCgg                                                               <  2:969131/151‑1 (MQ=255)                                                                          aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGGCCGCCAGAAAACTAATGCCCCCATGATCAT‑CATCGATCCGCGCCATAACCACAACGGGGCCggg                                                              >  1:892001/1‑151 (MQ=255)                                                                           aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                             <  2:100697/151‑1 (MQ=255)                                                                              aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCggg                                                              >  2:793401/1‑147 (MQ=255)                                                                               gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGa                                                         <  1:486785/151‑1 (MQ=255)                                                                                 aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                        <  2:982299/150‑1 (MQ=255)                                                                                 aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                        <  2:167902/150‑1 (MQ=255)                                                                                 aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGGCAGAAAACTAATGCCCCCATGAACAT‑CATCGATCCGCGCTATACCGACACCGGGGCCGGGcgcg                                                          >  1:793118/1‑148 (MQ=255)                                                                                  aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAG‑CATCGATCCGCGCTATACCGGCACCGGTGCAGGGCGCGAAGa                                                      <  2:740636/151‑1 (MQ=255)                                                                                    cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCCGAAATCTAATGCCCGCATGATCAT‑AATCGATCCGCGCTATAACGACACCCGGGCCGGGGGCGAAGATg                                                    >  2:783993/1‑151 (MQ=255)                                                                                      ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                   <  1:795905/150‑1 (MQ=255)                                                                                      ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg                                                  <  1:388357/151‑1 (MQ=255)                                                                                       gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCAAGAACAT‑CATCGATCCCCGCTATACCGACACCCGGGCCGGGcgcgga                                                        >  1:625218/1‑142 (MQ=255)                                                                                       gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGAGGGGGGGGGGGTTAATATCTTGAAAAGGGAGGCCAGAAATATAAAGCCCGCATGGTCAT‑CAACGGTCCGCGCGATACCGACACCGGGGCCGGGCGGGGAGATGAg                                                  >  2:242045/1‑150 (MQ=255)                                                                                               gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                         <  2:894259/151‑1 (MQ=255)                                                                                                   ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                      <  1:880508/150‑1 (MQ=255)                                                                                                      aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                  <  1:438179/151‑1 (MQ=255)                                                                                                       ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCCTACTATTTTTAACAAGAACGCCAGAAAACTTATGGCCGCATGATCAT‑CATCGATCCGCGCCATACCGACACCCGTGGCGGGCGCGGAGATGAGGGGATCCCTCTTCGTcc                                 >  2:634761/1‑151 (MQ=255)                                                                                                            ccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAATTTCTATCTTGAAAAGGCAAGCCAGAAAAATAAAGCCCGCATGAGTATCCATGGA‑CCGCCCTATACCCACAACGGGGCCGGGCGCGAAGATGAGTGGCTCCCTAATCGTCCGGGa                             >  2:738570/1‑150 (MQ=255)                                                                                                              tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAAATTAATGCCCGCCTGATCCT‑CATCGATCCCCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGGGGATCCCCCTTCGTCCgggga                            >  2:595963/1‑147 (MQ=255)                                                                                                              tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGGGGATCCCCATTCGTCCGGGAAc                           >  2:73559/1‑150 (MQ=255)                                                                                                                 cGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAAGCCCCCCAGAAAACTAATGCCCGCACGATCAT‑CATCCATCCGCGCTCTACCCACACCGGGGCCGGGCGCGAAGAAGAGGGGATCCCCCTTCTTCCGGGGAc                           >  2:825754/1‑147 (MQ=255)                                                                                                                  gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                      <  1:235768/151‑1 (MQ=255)                                                                                                                    aaCGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGAACAGGCACGCCAGAAAACTAAAGCCCGCATGATCAT‑CATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGcc                    >  1:603234/1‑151 (MQ=255)                                                                                                                    aaCGCGAATGAGTGGCGGAGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCATGATCAA‑CATCCGTCCCCGCTATAACGCCACCGGGGGCGGGCGCGAAGATGAGTGGATCCCTATTCCTCCGGGGAAAGATGcc                    >  1:222576/1‑151 (MQ=255)                                                                                                                        cgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc                <  2:787772/151‑1 (MQ=255)                                                                                                                        cgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCAc                <  2:1127825/151‑1 (MQ=255)                                                                                                                                     gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt   >  1:130205/1‑151 (MQ=255)                                                                                                                                       ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg  <  1:754807/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                   CCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAT‑CATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTG  >  NC_000913/941575‑941858

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