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

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941569‑941860                                                                                                                                                  |                                                                                                                                                    actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACt                                                                                                                                                 <  2:819418/149‑1 (MQ=255) actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTa                                                                                                                                               <  1:357271/151‑1 (MQ=255)  ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCtt                                                                                                                                                >  2:840835/1‑149 (MQ=255)        cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCtt                                                                                                                                        <  2:792231/151‑1 (MQ=255)          cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGATTTCTATCtt                                                                                                                                        >  2:714035/1‑149 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTACTAtttt                                                                                                                                        >  2:176964/1‑148 (MQ=255)            cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAAc                                                                                                                                    <  1:396241/151‑1 (MQ=255)              aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAg                                                                                                                                  <  2:643079/151‑1 (MQ=255)                ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCTTACTATCTTTaaaa                                                                                                                                   >  1:855175/1‑146 (MQ=255)                ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGAGGGGGGAGTTAATATCTTTaaaaag                                                                                                                                 >  1:89944/1‑146 (MQ=255)                       aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGCTACTATCTTTaaaa                                                                                                                                   >  1:823964/1‑139 (MQ=255)                             gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                    <  2:937992/150‑1 (MQ=255)                              aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                  <  1:854882/151‑1 (MQ=255)                                cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATTTTTAACAGGCACGCCAGAAAACTaa                                                                                                                 >  2:288453/1‑150 (MQ=255)                                  ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                              <  2:238348/151‑1 (MQ=255)                                   taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGc                                                                                                              <  1:685949/150‑1 (MQ=255)                                    aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc                                                                                                             <  2:2506/150‑1 (MQ=255)                                      cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGGTTACTATTTTGAAAACGCACCCCCCAAAACTaaaggcc                                                                                                            >  1:739624/1‑144 (MQ=255)                                        aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                         <  1:967498/150‑1 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGAGGC‑GGTGGGGGGATCTTCTATCTTGAACAAGCACACCACAAAACTaaaggcc                                                                                                            >  2:873212/1‑141 (MQ=255)                                          ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGACTTTCTATCTTGAACACGCACGCCCAAAAACTAATGCCCGCa                                                                                                         >  2:631443/1‑148 (MQ=255)                                                  ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAGCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                              <  1:179483/151‑1 (MQ=255)                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                           <  1:669222/151‑1 (MQ=255)                                                      gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                          <  2:172681/151‑1 (MQ=255)                                                        tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAttc                                                                                         >  2:638340/1‑148 (MQ=255)                                                         ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGGCTTACTATCTTTAACAGGCCCGCCAGAAATCTAATTCCCCCATGatca                                                                                                   >  1:68555/1‑139 (MQ=255)                                                          gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                      <  1:882328/151‑1 (MQ=255)                                                            aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTTCTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                     >  2:188437/1‑150 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCTTATTATCTTTAACAGGGACGCCAGAAAACTAATGACCCCATGATCATCATCTATCCGCGCt                                                                                    >  1:768027/1‑150 (MQ=255)                                                              cAGCCCGTCGGCTATCGAAAACAGATACCTGTTGGTGCTGTTTGGCAATAACCCAGGCCAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTTTTTTGAACAAGACAGCCGGAAATTTCATGCCCGCACGaacaccatcctttcgc                                                                                       >  1:865660/1‑139 (MQ=255)                                                                              gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCg                                                                  <  2:3867/151‑1 (MQ=255)                                                                                 aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                               <  2:521011/151‑1 (MQ=255)                                                                                    aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                            <  1:564214/151‑1 (MQ=255)                                                                                                 tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTTAC‑AACTTTAACAAGCCCCCCAAAAATCCAAAGCCCGCATGATCATCCTCGAGTCGCGCTTTACCGACACCCGGGCCGGGCGCGGAGAAGAGAgg                                                 >  2:309036/1‑149 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                              <  1:439561/151‑1 (MQ=255)                                                                                                   cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAAACTAAATCCCGCAAGAACATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGTGGGTcc                                             >  2:273303/1‑151 (MQ=255)                                                                                                   cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTcc                                             >  2:970547/1‑151 (MQ=255)                                                                                                       ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                          <  1:676881/150‑1 (MQ=255)                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                       <  2:187223/151‑1 (MQ=255)                                                                                                          gTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                      <  2:759906/151‑1 (MQ=255)                                                                                                           taataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGGCGGGGGCGAAGATGAGTGGGTCCCTAtt                                        >  2:217919/1‑148 (MQ=255)                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                                 <  1:713968/151‑1 (MQ=255)                                                                                                                aCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGGGGAACCCTATTCGTCCggg                                >  2:880926/1‑151 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                           <  1:120705/151‑1 (MQ=255)                                                                                                                         aaaCGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCCTGATCATCCTCGATCCGCGCTATACCGACACCCGTGCCGGGGGCGAAGATGAGTGGGTCCCTCTTCGTCCGGGGACAGa                          >  1:703698/1‑148 (MQ=255)                                                                                                                             gcgAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATTCGCGCTATACCCACACCGGTGCCGGGCGCGAAGATGAGGGGGACCCTATTCGTCCGGGGACAGATGCCGc                    >  1:772593/1‑150 (MQ=255)                                                                                                                                             tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg    >  2:350430/1‑151 (MQ=255)                                                                                                                                               gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  >  1:558709/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941569‑941860

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