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

CCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941572‑941860                                                                                                                                              |                                                                                                                                                      cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                                                                              <  2:40038/151‑1 (MQ=255)    ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                             <  2:41144/149‑1 (MQ=255)     cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                          <  1:284280/151‑1 (MQ=255)     cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCtt                                                                                                                                          >  1:609351/1‑151 (MQ=255)              tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCtt                                                                                                                                          >  1:450645/1‑142 (MQ=255)               gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTAATATCTTT‑AACAgg                                                                                                                                  >  2:356907/1‑148 (MQ=255)                  gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGcc                                                                                                                            <  1:124097/151‑1 (MQ=255)                    aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTT‑AACAAGGACGcc                                                                                                                            >  1:37191/1‑149 (MQ=255)                         tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTTTTTTT‑AAAAAGGACGCCCGaaa                                                                                                                       >  2:191824/1‑149 (MQ=255)                               ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAAt                                                                                                                <  1:309854/149‑1 (MQ=255)                               ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTG‑AAAAGGGACGCCAGAAAAC‑Taa                                                                                                                 >  1:596408/1‑148 (MQ=255)                                  ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCg                                                                                                           <  1:470887/151‑1 (MQ=255)                                  ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTATTATTTTG‑AAAAGGGACGcccaaaaaataaaatcc                                                                                                             >  2:400074/1‑135 (MQ=255)                                      cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTG‑AACAGGCACGCCAGAAAAC‑TAATGCCCGc                                                                                                          >  2:301370/1‑148 (MQ=255)                                      cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTAATATTTTGCACCGGGCAGCCCAGAAATT‑Taaacccgc                                                                                                           >  1:581847/1‑142 (MQ=255)                                       ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGa                                                                                                      <  1:632927/151‑1 (MQ=255)                                       ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTG‑AACAGGCACGCCAGAAAAC‑TAATGCCCGc                                                                                                          >  1:528098/1‑147 (MQ=255)                                              gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGatcatcat                                                                                               <  2:659377/151‑1 (MQ=255)                                                  aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGAGAAACGCGAATGAGTGGCGGAGGGGGGAAATATTATCCTG‑AAAAGGCACGCCAGAAAAC‑TAATGCCCGCATGatgaccaacg                                                                                             >  2:559604/1‑141 (MQ=255)                                                    aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATcc                                                                                         <  2:9236/151‑1 (MQ=255)                                                      ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCcgc                                                                                       <  1:489796/151‑1 (MQ=255)                                                      ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCcgc                                                                                       <  2:581826/151‑1 (MQ=255)                                                      ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCcg                                                                                        <  2:525610/150‑1 (MQ=255)                                                      ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCcg                                                                                        <  2:140858/150‑1 (MQ=255)                                                           cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGTTTTCTATCTTG‑AAACGGCAAGCCAGAAAAC‑TAAAGCCCGCATGATCATCCTCCATCcgc                                                                                       >  2:53532/1‑146 (MQ=255)                                                                    ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTG‑AAAAGGCACGCCAGAAAAC‑TAATGCCCGCAAGATCATCATCGAAACGCGCTATACCGacac                                                                          >  1:187384/1‑150 (MQ=255)                                                                     gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGGTTACTATCTTT‑AACAGGCACGCCAGAAATC‑TAAAGCCCGCACGATCATCATAGCTCCGCGCCATACCCACAcc                                                                         >  2:375400/1‑150 (MQ=255)                                                                     gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTT‑AACAGGCACGCCAGAAATC‑TAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACAcc                                                                         >  2:272613/1‑150 (MQ=255)                                                                     gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTT‑AAAAGGCACGCCCGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCCATACCGAcccc                                                                         >  1:91022/1‑150 (MQ=255)                                                                     gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                         >  1:4605/1‑150 (MQ=255)                                                                        atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                     <  1:272622/151‑1 (MQ=255)                                                                             aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTT‑AAAAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg                                                                >  1:394415/1‑151 (MQ=255)                                                                                cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                             <  1:85641/151‑1 (MQ=255)                                                                                    aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAg                                                         <  2:658445/151‑1 (MQ=255)                                                                                    aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTT‑AACCGGCCCCGCAGAAAAC‑TTATTCCCCCCTCATCATAATCGCTCCCCCCTATACCGAAACCTGCGGCggggg                                                              >  1:280098/1‑144 (MQ=255)                                                                                         ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑GACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                      <  2:277802/149‑1 (MQ=255)                                                                                         ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                     <  1:191831/150‑1 (MQ=255)                                                                                           tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTT‑TACCAGCACGCCCGAAAAC‑TAATGCCCCCCTGAACATCCTCCATCCGCCGTATACCGCCAACCGTGcc                                                                   >  2:37603/1‑134 (MQ=255)                                                                                           tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATTTTG‑AACAGGCGCGCCAGAAAACTTAATG‑CCGCCTGATCATCCTTGATCCGCGCGTTACCGAAACCCGTGCCGGGcgcg                                                            >  2:454457/1‑141 (MQ=255)                                                                                              tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGTGTCTTTATTACTTG‑AATAAGCAAACCCAAAAAG‑TAATGCCCGCATGAAACTCATTGAAGCGTGCTATAACGACACACGGGCCCGGGGCGaaaaaggggg                                                  >  2:271427/1‑140 (MQ=255)                                                                                                          ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                   <  1:77926/151‑1 (MQ=255)                                                                                                             aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTG‑AAAAGGCACGCCAGAAAAC‑TAATGCCCGCATGAACATCATCGATCCGCGCTATACCGACACCGGGGCCCGGCGCGGAAATGAGTGGGTCCCCCTTCGTCCggg                                >  2:220734/1‑151 (MQ=255)                                                                                                                    cGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCCTGATCATCCTCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                >  1:65/1‑144 (MQ=255)                                                                                                                     gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTT‑AACAAGCACGCCCGAAAAC‑TAATGGCCGCATGATCATCACCGATCCGCGCCATAACGACACCCGGGCCGGGGGCGAAAATGAGTGGGTCCCTATTCGTCCGGGGACa                            >  1:187328/1‑147 (MQ=255)                                                                                                                        aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                     <  2:558045/151‑1 (MQ=255)                                                                                                                                     gCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg        <  1:532765/151‑1 (MQ=255)                                                                                                                                          ggggTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg   >  2:409220/1‑151 (MQ=255)                                                                                                                                          ggggTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg    <  2:19179/150‑1 (MQ=255)                                                                                                                                           gggTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  >  1:384514/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                      CCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTG‑AACAGGCACGCCAGAAATC‑TAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941572‑941860

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