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

CTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941570‑941857                                                                                                                                                |                                                                                                                                                 ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTa                                                                                                                                            <  2:180882/150‑1 (MQ=255)    ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                             >  1:344034/1‑146 (MQ=255)            aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATTTTTaa                                                                                                                                  >  1:41208/1‑149 (MQ=255)             aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaaa                                                                                                                                 >  1:45339/1‑148 (MQ=255)               ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAGCTACTATCTTTAACa                                                                                                                                >  2:644953/1‑148 (MQ=255)                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                           <  1:590588/151‑1 (MQ=255)                    gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGGGGCGGTGGGGGGGATTATTATCTTGTACAGGCAcac                                                                                                                         >  2:269437/1‑148 (MQ=255)                        gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                    <  1:264197/151‑1 (MQ=255)                            gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                                <  1:446339/151‑1 (MQ=255)                            gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                                <  2:431306/151‑1 (MQ=255)                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                      <  2:517410/150‑1 (MQ=255)                                           cggcTGGGCAGATTGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGTCGAAACTCGAAGGAGTGTCGGTGCGGTGACTTTCTATCTTGGACAGGCTCGCCAGAAATCTAATGCCCGCATGttc                                                                                                 <  2:524452/151‑1 (MQ=255)                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTTTCTTTAACAGGGACGCCCGAAAACTAATGCCCGCACGatc                                                                                                 >  2:499039/1‑150 (MQ=255)                                             gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                               <  2:407966/151‑1 (MQ=255)                                               tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                               >  2:7545/1‑149 (MQ=255)                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGGTTACTATCTTGAAAAGGCACGCCAGAAAAATAATGGCCGCATGgtca                                                                                                >  1:6802/1‑147 (MQ=255)                                                        ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                     <  2:101720/150‑1 (MQ=255)                                                         gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTTTCTTGAACAAGGACGCCAGAAATCTAATTGCCGCCTGTTCATCAACGATc                                                                                       >  1:423826/1‑147 (MQ=255)                                                               gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                             <  1:341138/151‑1 (MQ=255)                                                                 ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                                                                           <  2:300013/151‑1 (MQ=255)                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGGGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAAACTAAAGCCCGCCTGAACATCATCCATCCCCGCTATCCCGaca                                                                        >  2:623016/1‑150 (MQ=255)                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGAACAAGCACGCCCGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                              >  1:638772/1‑151 (MQ=255)                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                              >  1:140241/1‑151 (MQ=255)                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGGTTTTTATATTTAAAAAGCACGCCCGAAAACTAATGCCCCCCCGATTATCATCGATCCGCGCCAAACCCACAcccggggcgg                                                              >  2:60200/1‑143 (MQ=255)                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATTCTGAAAAGGCCCGCCCGAAAAATAATGCCCCCATGATCATAATCGAAGCGCGCCCTACCgccccccggggcggg                                                             >  2:205210/1‑134 (MQ=255)                                                                                      aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCCCGCCATACCCACACCGGGGCCggggg                                                           >  2:660196/1‑144 (MQ=255)                                                                                       aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTTCTATCTTGAACACGCACGCCAGAAATCTAATGGCCCCAAGAATATCCTCGATCCGCGCTATACCGACACCGGTTCCCGGcgcgga                                                       >  2:677349/1‑147 (MQ=255)                                                                                         cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGACTTACTATCTTTAAAAGGCACCCCAGAAATATAATGGCCGCATGATCATCATCGATTCCCGCCATTCCGCCACCGGTGCCGGGCGCGaa                                                       >  2:457322/1‑147 (MQ=255)                                                                                                         gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCCTTAACAGGGACCCCCGAAATCTAATGCCCCCCTGATCATCAACGAACCGCGCTCTACCGACACCGGGGCCGGGCGCGGAAAGGAGTGGAACCCTAt                                      >  2:538309/1‑148 (MQ=255)                                                                                                          taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGGGGATCCCTATTCGt                                  >  1:73608/1‑151 (MQ=255)                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTTCTATCTTTAACAGGCCCGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                >  1:109142/1‑151 (MQ=255)                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTATTATCTTGAACAGGCACGCCAGAAATTTAAAGCCCGCCTGATCATCATCGATCCGCGCTATACCCACACCGGGGGCGGGCGCGAAGATGAGTGGGTCCCTATTCGTcc                                >  2:188772/1‑151 (MQ=255)                                                                                                              aaCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGACTTACTATCTTTAACAAGCCCGCCAGAAATCTAATGCCCGCATGATCATCATCGGTTCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTcccg                              >  2:569115/1‑149 (MQ=255)                                                                                                                 ccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGTaa                           <  2:659439/151‑1 (MQ=255)                                                                                                                      cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGa                       <  1:287618/150‑1 (MQ=255)                                                                                                                        aaaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGAAAGATg                     >  1:345835/1‑150 (MQ=255)                                                                                                                         aaCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGTGGATCCCTATTCGTCCGGGGAAAGATGcc                   >  1:310791/1‑151 (MQ=255)                                                                                                                                aTGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGGACGCCCGAAAACTAATGGCCGCCTGATCATCCTTGCTCCGCGCTATACCGACACCCGGGCCGGGGGCGAAGATGAGTGGATCCCCATTCGTCCGGGAACAGATGCCGCACTgg            >  2:417799/1‑151 (MQ=255)                                                                                                                                  gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt          <  2:343996/151‑1 (MQ=255)                                                                                                                                      ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc       <  1:504512/150‑1 (MQ=255)                                                                                                                                       gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg     <  1:352073/151‑1 (MQ=255)                                                                                                                                          gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  <  1:590467/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 CTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941570‑941857

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