Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	266,358	0	T	G	46.0%	59.9 / 13.0	51	T66P (ACC→CCC)	yafY	lipoprotein, inner membrane; degP regulator; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (7/20);  new base G (23/0);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 2.40e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266211‑266506                                                                                                                                                    |                                                                                                                                                      cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTggg                                                                                                                                                    <  2:587143/151‑1 (MQ=255) cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTgg                                                                                                                                                     <  2:118430/150‑1 (MQ=255)  ctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGTAGTAAAAACTGCCAgggggg                                                                                                                                                    >  1:738153/1‑150 (MQ=255)       ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGa                                                                                                                                              <  1:898716/151‑1 (MQ=255)       ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGTAAAAACTGCCAGGGGGGGCCgg                                                                                                                                               >  1:282930/1‑150 (MQ=255)          aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCGCAAGGGGGGGCCgggaa                                                                                                                                            >  2:661445/1‑147 (MQ=255)            tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCGCAAGGGGGGGCCGGGTaa                                                                                                                                           >  1:60002/1‑149 (MQ=255)              cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGCCgggaaa                                                                                                                                           >  1:642272/1‑143 (MQ=255)                      tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGCCGGATAAGTAGTCAA‑cc                                                                                                                                >  1:765942/1‑149 (MQ=255)                       gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCtgt                                                                                                                             <  2:901701/151‑1 (MQ=255)                          aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAACCCGTAAGGGGGGGCCCCAGAAAGAGTCAA‑CCagt                                                                                                                             >  1:640080/1‑148 (MQ=255)                               cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGTAAAAACTGCAAGGTGGGACCGGATAAATAGTCAA‑ACCGTGACCAtt                                                                                                                      >  1:709216/1‑150 (MQ=255)                                aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGGCGGGGAAATAGGCAA‑ACTTTGACCCtt                                                                                                                      >  1:277404/1‑149 (MQ=255)                                 aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGGAAAAACCGCAAAGGGGGGCCGGGAAAAAAGTCAA‑ACTGTGAACcttctt                                                                                                                   >  1:497391/1‑151 (MQ=255)                                   acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGCCGGATAAAAGGGCAA‑CCCGTGAGCcttctt                                                                                                                   >  2:952383/1‑149 (MQ=255)                                     acaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACCTAAAAACTGCAAAGGGGGACCGGGTAAGGAGTCAA‑CCTGTGACCAttctt                                                                                                                   >  1:941243/1‑147 (MQ=255)                                      caGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTg                                                                                                              <  2:427958/151‑1 (MQ=255)                                       aGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGa                                                                                                             <  2:228675/151‑1 (MQ=255)                                       aGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCACCAAGTGGGGCCGGGAAAGGAGTCAA‑ACTCTGACCCttcttctccg                                                                                                              >  1:503464/1‑146 (MQ=255)                                          gcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGAcat                                                                                                          <  2:228902/151‑1 (MQ=255)                                               catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGCCGGGTAAAGAAGCAA‑CCTGTGACCATTCTTTCCCGACATCATcc                                                                                                     >  1:895246/1‑151 (MQ=255)                                               catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACTTCAAAGTGGGACCCGCTAAGGAGTGAA‑ACTCGGACCCCTCTTCGCTgccggca                                                                                                        >  1:800777/1‑142 (MQ=255)                                                     tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGGCCGGATAAGTAGGCAA‑CCCGTGGCCAATCTTCGCTGACATCATCCttgtgg                                                                                               >  2:214949/1‑151 (MQ=255)                                                     tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCCAGGGGGGGCCGGGGAAATAGTCAA‑CCCGTGACCCTTCTTCGCCGACATCATCCtggtgg                                                                                               >  2:82366/1‑151 (MQ=255)                                                                     tCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCCAGGGGGGACCGGGGAAATAGGCAA‑ACTGTGACCCTTCTTCGCTGACATCCTCCTCGTGGCTGGCCACAt                                                                                     >  1:869434/1‑145 (MQ=255)                                                                        tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTATAAACTACAAAGGGGGGCCGGGAAAGAAGTCAA‑CCTGTGAACATTCTTTGCTGACATCATCCTGGTGGGTGGCCACAACAATGaaa                                                                             >  2:668636/1‑150 (MQ=255)                                                                         gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAAGGGGGAGCGGGGAAAGAGTCAA‑CCCGTGACCATTCTTCGCTGGCAAAATCCCGGTGGGGGACCACATCACCGaaa                                                                             >  2:103998/1‑149 (MQ=255)                                                                           atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgcg                                                                         <  2:642082/151‑1 (MQ=255)                                                                            tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGTGGGGCCGGGGAAGTAGTCAACCCTGT‑ACCCTTTTTCGCCTACACCCTCCTGGTGGCCGACCACCACACCTAAACGCgg                                                                        >  1:80526/1‑151 (MQ=255)                                                                             atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCgg                                                                        <  1:952180/150‑1 (MQ=255)                                                                             atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGt                                                                       <  2:632281/151‑1 (MQ=255)                                                                               ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAAGGGGGACCCGATAAGTAGTCAA‑ACCGTGACCATTCTTCGCCGACATCATCCCGGTGGCTGACCACATCACTGAAACGCGgggg                                                                     >  1:186327/1‑151 (MQ=255)                                                                                cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg                                                                    <  1:845811/151‑1 (MQ=255)                                                                                 aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGACCGGGTAAGTAGTCAA‑CCTGTGAGCATTCTTTGCTGACATCATCCTGGTGGGTGGCCACCTCACTGAAACGCGgggggg                                                                   >  1:620145/1‑151 (MQ=255)                                                                                     ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa                                                                <  2:582264/150‑1 (MQ=255)                                                                                        agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACTGCCAGGGGGGGCCGGATAAGTAGGCAA‑CCTGTGACCATTCTTTGCTGACATCATCCTGGGGGGTGACCACATCACTGAAAACCGggggggggggg                                                              >  1:65818/1‑144 (MQ=255)                                                                                          tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCCGCCAGGGGGGGCCGGATAAGAAGGCAA‑CCCGTGACCCTTCTTCGCTGACATCATCCTGGGGGCTGACCACATCACTGAAACGCGggggggggggggcc                                                           >  1:115992/1‑142 (MQ=255)                                                                                              aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGCAAAAACCGCAAGGAGGGGCCGGAAAAGGAGGCAA‑ACCGTGGCCCTTCTTTGCCGACCTCAACCCCGTGGCTGGGCCCATCACTGAAAAGCGggggggggggggccct                                                         >  2:373045/1‑138 (MQ=255)                                                                                                  gCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAtt                                                    <  1:93895/149‑1 (MQ=255)                                                                                                     aTCCTCGCAGCCATCGGGCGACTGGCGGGTATGTAAAAACTGCGAGGGGGGGGCGGATAAATAGGCAA‑ACCGTGGCCCTTCTTCTCCGTCATCATCCCGGGGGCTTCCCACAACACCGGAACCGGGGGGGGGGGGGGCCTTtaa                                                      >  1:14270/1‑144 (MQ=255)                                                                                                        cccgcgGCCATCGGGCGCCTGGCCGGTACGTAAAGACTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACGATTCTTCGCGGACATCATCCTGATGTCTGCCCACAGCACTGAAGCGGGGTGGGGTGTCGTCCTTATAtttc                                                  <  2:937917/143‑2 (MQ=255)                                                                                                        cTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAAGTGGGGCCGGATAAATAAGCAA‑ACTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGGCCACACCACTGAAAAGCGGTGGGGGGAGGGCCTTAAATTTTCTg                                               >  2:452678/1‑148 (MQ=18)                                                                                                          cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGc                                          <  2:59989/151‑1 (MQ=255)                                                                                                            cAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGAGAAATAGTCAA‑CCTGTGACCCTTCTTCGCTGACATCATCCTGGTGGCGGAACCCATCACTGAAACGCGGGGGGGGGGGGGCCTTAAAATTACTGAttggc                                         >  2:343919/1‑146 (MQ=18)                                                                                                               ccATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGGTAAGTAGTCAA‑CCTGTGACCATTTTTCGCTGACATCATCCTGGGGGGTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCTTTAATTTACTGATGGCCtt                                       >  1:325227/1‑149 (MQ=255)                                                                                                                 aTCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAgg                                   <  1:587088/151‑1 (MQ=255)                                                                                                                      gCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc                              <  2:577647/151‑1 (MQ=255)                                                                                                                                ggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa                    <  2:123408/151‑1 (MQ=255)                                                                                                                                    cGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc                <  2:277332/151‑1 (MQ=255)                                                                                                                                           ccTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt         >  1:541324/1‑151 (MQ=255)                                                                                                                                               cAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATgaga     <  1:967423/151‑1 (MQ=255)                                                                                                                                                   gTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  1:802178/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266211‑266506

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