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

CCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  NC_000913/266211‑266506                                                                                                                                                     |                                                                                                                                                     cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTggg                                                                                                                                                   <  2:261521/151‑1 (MQ=255)  ctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAA‑CTGCAAgggggg                                                                                                                                                   >  1:700114/1‑150 (MQ=255)          aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATaa                                                                                                                                          <  1:392000/151‑1 (MQ=255)            tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAA‑CTGCCAGGGGGGGCCgg                                                                                                                                              >  2:590823/1‑145 (MQ=255)             ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAagt                                                                                                                                        <  1:510947/150‑1 (MQ=255)               aGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAA‑CCGCCAGGGGGGGCCGGATaaaga                                                                                                                                       >  2:628068/1‑146 (MQ=255)                   tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                                  <  1:103127/150‑1 (MQ=255)                   tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCaa                                                                                                                                  <  1:410459/150‑1 (MQ=255)                    gATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAAc                                                                                                                                 <  2:277706/150‑1 (MQ=255)                      tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAAGGGGGACCGGGTAAAGAGTCAAcc                                                                                                                                >  1:718286/1‑149 (MQ=255)                       gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAC‑CTGCAAGGGGGGACCGGATAAATAGTCAACCtg                                                                                                                              >  1:27330/1‑150 (MQ=255)                             ggCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAA‑CTGCAAGGGGGGGCCGGGTAAGTAGTCAAACTGTGAcc                                                                                                                         >  2:794183/1‑149 (MQ=255)                                aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttct                                                                                                                    <  2:690285/151‑1 (MQ=255)                                 aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAA‑CTCCAAAGGGGGGACGGATAAGGAGTCAACCTGTGGACAttctt                                                                                                                   >  2:857680/1‑151 (MQ=255)                                   acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttcttc                                                                                                                  <  1:40163/150‑1 (MQ=255)                                   acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAA‑CCGCAACGGGGGACCGGGGAAAAAATCTACCTGTGagccatctt                                                                                                                   >  2:469889/1‑141 (MQ=255)                                    cacaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGc                                                                                                                <  1:104236/151‑1 (MQ=255)                                      caGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTg                                                                                                              <  2:785622/151‑1 (MQ=255)                                       aGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAA‑CCGCAAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGc                                                                                                                >  1:368286/1‑148 (MQ=21)                                          gcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAC‑CCGCAAGGGGGGGCCGGGGAAATAGTCAAACTGTGACCCTTTTTTGCTGAca                                                                                                           >  1:158409/1‑150 (MQ=255)                                                     tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAA‑CCGCCAGGTGGGGCCGGATAAGTAGTCAACCTGTGGCCATTCTTCGCTGACATCATCCtggtgg                                                                                               >  2:369234/1‑151 (MQ=21)                                                       gCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAA‑CCGCCAGGGGGGACCGGATAAATAGTCAACCTTTGACCATTCTTCGCTGACATCATCCcggtg                                                                                                >  2:792424/1‑148 (MQ=21)                                                           gACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAA‑CTGCAAGGGGGGGCCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGAcc                                                                                         >  1:99453/1‑151 (MQ=255)                                                            aCGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGGGGGACCGGATAAGTAGTCAACCCGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGAcc                                                                                         >  1:58379/1‑150 (MQ=255)                                                                  ttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAA‑CTCCAAAGGGGGGCCGGAGAAGTAGACAACCTGTGACCCATCTTCGCCGGAATCACCCTGGTGGCGGaccacaccac                                                                                  >  2:178074/1‑151 (MQ=255)                                                                   tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAA‑CTGCAAGGTGGGACCGGGTAAGGAGGCAACCTGTGACCATTCTTCGCTGACATCATCCCGGTGGCTGACcac                                                                                       >  2:74295/1‑145 (MQ=255)                                                                      ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa                                                                              <  1:27775/151‑1 (MQ=255)                                                                        tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGAAACAGAAAAC‑CACCCAGGGGGGACCGCAAAAAAAGTGAAACCGGGAGCAATACTCCCCGACAACAACCTGGGGGGTGACCCCATCACTGaaa                                                                             >  2:373555/1‑150 (MQ=255)                                                                              tCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGGGGGACCGGATAAGTAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGACCACATCACTGAAACGCGGTg                                                                      >  1:450780/1‑151 (MQ=255)                                                                                cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg                                                                    <  2:83313/151‑1 (MQ=255)                                                                                 aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGAATGGAAAACGCCGC‑AGGGGGGGCCGGGGAAGAAATCAAAGTTTTGCCCTTCCTCGCTGACATCAACCTGGTGGAGGACCCCCTCACTGTAACGCGgggggg                                                                   >  2:845591/1‑151 (MQ=255)                                                                                         gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGGGGGGCCGGATAAATAGTCAAACTGTGACCATTCTTCGCTGACATCATCCTTGGGGGTGACCCCATCACTGAAAAGCGGGGGGGGGAGGGcc                                                           >  1:489105/1‑151 (MQ=18)                                                                                              aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTTCAAGGGGGGACCGGGTAAGTAGTCCACCCGTGACCCTTTTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCtt                                                         >  1:363133/1‑148 (MQ=255)                                                                                                       ccTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAA‑CTGCCAGGGGGGGCCCGATAAAGAGTCAAACCGGGGCCATTTTTCGCTGACATCATCCTGGGGGCTGGCCCCATCACTGAAACGCGGGGGGGGGGGGGCCTTTAATTTACt                                                >  2:357798/1‑148 (MQ=255)                                                                                                          cGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATg                                            <  1:409230/149‑1 (MQ=255)                                                                                                           gCAGCCATCGGGCGACTGGCGGGTACGGAAAAA‑CTGCAAGGGGGGGCCGGGAAAAGAGTCAACCCGTGACCCCTTTTTTCCGAACCCCTCCCGGGGGGTGTCCACCTCACTCAAACCCGGTGGGGGGAGGGCCTTCAATTTAATGAt                                             >  2:845421/1‑147 (MQ=255)                                                                                                              gCCATCGGACGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTc                                      <  1:304880/151‑1 (MQ=255)                                                                                                               ccATCGGGCGACTGGCGGGGACGTAAAAA‑CCGCAAGGGGGGGCCGGATAAGTAGTCAACCTGTGACCCTTCTTCGCTGACATCATCCTGGGGGCTGACCACAACACTGAAAAGCGGGGGGGGGAGGGCCTTTAATTTACTGATGGCCTTc                                      >  2:718580/1‑150 (MQ=18)                                                                                                                cATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGGGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGGGAGGGCCTTAAATTTACTGATGGCCtt                                       >  1:791948/1‑148 (MQ=255)                                                                                                                           tGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTTAATTTACTGATGGCCTTCAGGGCCTCACtt                          >  2:681604/1‑150 (MQ=255)                                                                                                                              cGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt                      >  1:148849/1‑151 (MQ=255)                                                                                                                                   aCGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc                  <  1:614236/150‑1 (MQ=255)                                                                                                                                   aCGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc                  <  2:128126/150‑1 (MQ=255)                                                                                                                                         aaC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCg           <  1:814599/151‑1 (MQ=255)                                                                                                                                          aC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCg           >  1:71817/1‑150 (MQ=255)                                                                                                                                                  aGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGaa    >  2:806501/1‑150 (MQ=255)                                                                                                                                                   ggTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAc  >  2:617755/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAC‑CTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGATGAGAAAC  >  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: ‑