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

TACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941568‑941862                                                                                                                                                  |                                                                                                                                                       tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACtt                                                                                                                                                   >  1:636126/1‑151 (MQ=255)  actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                    <  1:946269/149‑1 (MQ=255)  actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                    <  2:140375/149‑1 (MQ=255)     cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                                                                               <  2:1014767/151‑1 (MQ=255)        ccGCGCAAATTGCGGAAGGTTTGAACTATTCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                            <  1:917034/151‑1 (MQ=255)          gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTg                                                                                                                                          <  1:373883/151‑1 (MQ=255)          gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTg                                                                                                                                          <  2:974009/151‑1 (MQ=255)                   gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                  <  2:840018/150‑1 (MQ=255)                      gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGcc                                                                                                                              >  1:133180/1‑151 (MQ=255)                           tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                          >  1:211899/1‑150 (MQ=255)                            ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaaa                                                                                                                         <  2:189640/150‑1 (MQ=255)                            ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGaaa                                                                                                                         >  1:134153/1‑150 (MQ=255)                              gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTGAACAGGGCCGCCAGgaa                                                                                                                         >  2:99757/1‑148 (MQ=255)                                        tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTTCTTTCCTTAACAAGCCCGCCCGAAAACTaaatacc                                                                                                               >  2:1035465/1‑143 (MQ=255)                                             cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAAGCCCGCATGat                                                                                                        >  1:450807/1‑150 (MQ=255)                                                      aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGACTTACTATCTTGAAAAGGGACGCCCGGAATATAATGCCCGCCCGATTATCATAGa                                                                                               >  2:768580/1‑150 (MQ=255)                                                      aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAAGCCCGCACGatcat                                                                                                     >  2:917817/1‑144 (MQ=255)                                                       gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATATAATGCCCGCCTGATCATCATCGAt                                                                                              >  2:851055/1‑150 (MQ=255)                                                          ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                           <  2:145476/150‑1 (MQ=255)                                                             aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                       <  1:738470/151‑1 (MQ=255)                                                              aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCt                                                                                       >  1:747593/1‑150 (MQ=255)                                                                aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                    <  1:613337/151‑1 (MQ=255)                                                                   ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCCCGCTATAcc                                                                                  >  1:336644/1‑150 (MQ=255)                                                                            atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGt                                                                        <  2:110544/150‑1 (MQ=255)                                                                            atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGt                                                                        <  1:1033042/150‑1 (MQ=255)                                                                              cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGcc                                                                     <  1:396921/151‑1 (MQ=255)                                                                              cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACA‑CCGGTGcc                                                                     >  1:293767/1‑151 (MQ=255)                                                                                  aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGgc                                                                 <  1:1059063/151‑1 (MQ=255)                                                                                      gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGa                                                             <  2:173000/151‑1 (MQ=255)                                                                                               tcgTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTGGACAGGGAAGCCAGAAATCTAAGTCCCGCATGATATTCATCGATCCGCGCTATACCAACA‑CCGGGGGCGGGCGGGAAGATGa                                                       >  1:182453/3‑148 (MQ=255)                                                                                                 gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGa                                                  <  1:979651/151‑1 (MQ=255)                                                                                                 gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGGGCCGGGCGCGAAGATGAGTgg                                                   >  1:197359/1‑150 (MQ=255)                                                                                                  tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAgggg                                                   >  1:895799/1‑149 (MQ=255)                                                                                                     tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAAGGACGCCAGAAAACTAATGCCCCCCTGATTATTATCGATCCCCGCTATAACCACACCCGGGGCCGGCGGCG‑AGAAGGGTGGATccc                                              >  2:682867/1‑151 (MQ=255)                                                                                                        ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                            <  1:375885/150‑1 (MQ=255)                                                                                                         tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  1:1006585/150‑1 (MQ=255)                                                                                                         tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                           <  1:357052/150‑1 (MQ=255)                                                                                                          ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCATGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                         <  1:807301/151‑1 (MQ=255)                                                                                                          ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                          <  1:299178/150‑1 (MQ=255)                                                                                                           gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                        <  1:734401/151‑1 (MQ=255)                                                                                                           gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCACAAATCTAATGCCCCCATGATCATCATCGATCCGCCCCATACCGACA‑CCGGGGCCGGGGGCGCAGATGAGTGGGTCCCTAtt                                          >  1:990836/1‑149 (MQ=255)                                                                                                           gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCCTTAACAGGCACGCCAGAAAATTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACC‑CCCGTGCCGGGCGCGGAGATGAGTGGATTCCTAtt                                          >  2:1018115/1‑149 (MQ=255)                                                                                                            taataaCCCGGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTGCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                         <  1:48320/149‑1 (MQ=255)                                                                                                                aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                                   <  2:17210/151‑1 (MQ=255)                                                                                                                  cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                  <  2:732869/150‑1 (MQ=255)                                                                                                                          aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTCTACCGACA‑CCGGGGCCGGGCGCGAAGATGAGTGGGACCCTATTCGTCCGGGAACAGATg                          >  2:494664/1‑150 (MQ=255)                                                                                                                               cgAATGAGTGGCGGGGGGGGGACTTTCTATCTTTAACAGGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGATCCCCGCCATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGGCGCAc                    >  2:906508/1‑151 (MQ=255)                                                                                                                                   tGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGCCGCACTgg                 >  1:637041/1‑150 (MQ=255)                                                                                                                                     aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  1:76380/151‑1 (MQ=255)                                                                                                                                          cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         <  2:188245/151‑1 (MQ=255)                                                                                                                                           ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc        <  1:403809/151‑1 (MQ=255)                                                                                                                                                 gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  2:1005194/1‑151 (MQ=255)                                                                                                                                                 gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  2:822570/1‑151 (MQ=255)                                                                                                                                                  tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:908526/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                       TACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACA‑CCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941568‑941862

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