Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,678,039	0	T	G	43.1%	82.8 / 11.3	51	N108T (AAC→ACC)	yphE	putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (5/24);  new base G (22/0);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 8.69e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677892‑2678187                                                                                                                                                     |                                                                                                                                                     aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCA‑GAGGgttt                                                                                                                                                   >  2:924087/1‑151 (MQ=255)  gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTc                                                                                                                                                  <  2:893658/151‑1 (MQ=255)  gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTc                                                                                                                                                  <  2:822034/151‑1 (MQ=255)      ctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGGACGGGGGA‑GAGGgtttttgcc                                                                                                                                              >  1:444964/1‑146 (MQ=255)       ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCA‑GAGGGTTTCcccc                                                                                                                                              >  2:1110478/1‑150 (MQ=255)        tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGCGGCCCCGACCGCGGCAGGGGGTTTTcccc                                                                                                                                               >  1:516441/1‑147 (MQ=255)               cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCA‑GAGGGTTTTCGcc                                                                                                                                              >  1:767014/1‑142 (MQ=255)                 ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGGTTTCCGCCACTGTCAGccc                                                                                                                                   >  2:703286/1‑151 (MQ=255)                 ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCC‑GAGGGTTTTCGCCACTGTCAGccc                                                                                                                                   >  2:363221/1‑151 (MQ=255)                     tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCtt                                                                                                                                 <  1:45805/149‑1 (MQ=255)                      aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTcc                                                                                                                               <  2:702941/150‑1 (MQ=255)                       gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTccac                                                                                                                             <  1:832289/151‑1 (MQ=255)                        cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTccacc                                                                                                                            <  1:647346/151‑1 (MQ=255)                            gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGGCCGAGGCA‑GAGGGTTTTCGCCCCTGTCAGCCCCTTCACCCGa                                                                                                                         >  1:868144/1‑150 (MQ=255)                              aaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACtt                                                                                                                      <  1:891269/151‑1 (MQ=255)                              aaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCC‑GAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGAAtt                                                                                                                      >  1:906317/1‑151 (MQ=255)                                cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCA‑GAGGGTTTTCGCCACTGTCAGCCCCTCCCCCCGACtt                                                                                                                      >  2:267529/1‑149 (MQ=255)                                cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTaa                                                                                                                    <  1:349444/151‑1 (MQ=255)                                   gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTc                                                                                                                 <  1:1007814/151‑1 (MQ=255)                                    ttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGTGCCCCCGACCCGGGGC‑GAGGGTTTGCGCCCCTGTCAGCCCGTccacccaaatttaacac                                                                                                                >  2:990983/1‑139 (MQ=255)                                      gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCt                                                                                                               <  1:396816/150‑1 (MQ=255)                                         cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCCTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGAt                                                                                                            <  2:320854/150‑1 (MQ=255)                                          aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCA‑GAGGGTTTTCCCCCCCGTCAGCCCTTCCCCCCAACTTAACTcc                                                                                                                >  2:40059/1‑145 (MQ=255)                                              cTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAcc                                                                                                       <  1:86012/150‑1 (MQ=255)                                               tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAAcc                                                                                                       >  2:392418/1‑149 (MQ=255)                                                          cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGTGCCACCGCCCGAGGCA‑GAGGGTTTTCCCCGCCTTCACCCCTTCCCCCAGACTTAACTCCTGATAAACCGCCCGCcccccc                                                                                           >  1:213040/1‑150 (MQ=255)                                                           ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCCAGGCA‑GGGGGTTTTCGCCACTGTCACCCCCTCCCCCACACCTAACTCCTTATAAAACACCCGCcccccc                                                                                           >  1:83207/1‑149 (MQ=255)                                                            cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGt                                                                                        <  2:44437/151‑1 (MQ=255)                                                             aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGCGGGA‑GAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCACACTTAAATCCTGGTAAAACCCGCGCAccccc                                                                                           >  1:656094/1‑147 (MQ=255)                                                                       aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGGCCCAGGCA‑GAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACTTTACTCCTGATAAAACCCGCGCACCCCCCGCTCCGCGGCAc                                                                              >  1:986713/1‑150 (MQ=255)                                                                         cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGc                                                                           <  2:1103960/151‑1 (MQ=255)                                                                          aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGCCCCCCGCCCCAGGAA‑GAAGGTTTTCGCCCCCGTCAGCCCCTTCACCAGACCTAAAACCCGATAAACCCCGCGCACCCCCCCTTCCGCGGCACGGCg                                                                          >  1:1081352/1‑151 (MQ=255)                                                                           aCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCTACCCAGGCA‑GAGGGTTTTCCCCACCGTCAGCCCCTCCCCCCGACTTAAATCCCGATAAACCGCGCGCCCCCCCCGTTCAGCGGCCCgg                                                                            >  1:640359/1‑148 (MQ=255)                                                                                gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                     <  1:687014/150‑1 (MQ=255)                                                                                      ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAg                                                               <  2:506400/150‑1 (MQ=255)                                                                                        tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGCGGCA‑GAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACTTTACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGGACCGCGGGTCAGCgga                                                                >  2:150612/1‑145 (MQ=255)                                                                                         gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc                                                           <  2:825540/151‑1 (MQ=255)                                                                                             cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCa                                                       <  2:631427/151‑1 (MQ=255)                                                                                               tttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCA‑GAGGGTTTCCCCCACCGTCAGCCCTTCCACCCGACCTCACTCCTGCTAAACCGCCCGCACCCCCCCTTCACCGGCACGGCGAGTCCGCGTAGCCTCGTCAcc                                                     >  2:300018/1‑151 (MQ=255)                                                                                                 tGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCCAGGCA‑GAGGGTTTTCCCCCCTGTCAGCCCCTCCACCAGACCTAAATCCCGATAAAACGCGCGCACCCCCCCTTCAGCGGCACGGCGAGGCAACGTAGCTTCGTCAACtt                                                   >  2:39091/1‑151 (MQ=255)                                                                                                   cAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGAACCCGGGA‑GAGGGTTTTCCCCCCCGTCCGCCCCTCCACCAGACCTAAATCCTGTTAAAAACCGCGCCCCCCCCACTCCGCGGCCCCGCGCGGCAGCGCAGCTTTCTCAccct                                                   >  1:785577/1‑147 (MQ=255)                                                                                                        aatcaatcATGCCGTTGCGGCGGGGCCCCCGACCGAGGCA‑GAGGGTTTTCGCCACCGTCAGCCCCTCCACCCGACTTAAATCCTGATAAAACCCGCGCCCCCCCCGTTCAGCGGCACGGCGGGTCAGCGGAGCTTCGTCACCCTCCAGTcg                                            >  2:1088278/1‑151 (MQ=255)                                                                                                           caatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgcg                                          <  1:602172/150‑1 (MQ=255)                                                                                                            aatcaTGCCGTTGCGGCGGGGCCACCGACCGCGGCA‑GAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCCGACTTAAATCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGGCAGCGGAGCTTCGTCACCCTCCACTcgcg                                          >  2:136597/1‑149 (MQ=255)                                                                                                                aTGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGc                                    <  2:647035/151‑1 (MQ=255)                                                                                                                       tGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                              <  2:755091/150‑1 (MQ=255)                                                                                                                              gggCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                        >  1:345236/1‑149 (MQ=255)                                                                                                                                  cACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc                  <  1:1108022/151‑1 (MQ=255)                                                                                                                                       aCCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATcc             >  1:217677/1‑151 (MQ=255)                                                                                                                                        ccGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg            >  2:748257/1‑151 (MQ=255)                                                                                                                                         cGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg           <  1:574399/151‑1 (MQ=255)                                                                                                                                              cA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGtt      <  1:894302/151‑1 (MQ=255)                                                                                                                                                   ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  >  2:918965/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCA‑GAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677892‑2678187

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