Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,615,733	0	T	A	46.4%	52.7 / 32.7	57	Y50F (TAT→TTT)	yfgO	putative UPF0118 family inner membrane permease
Reads supporting (aligned to +/- strand):  ref base T (12/18);  new base A (26/0);  total (39/18)
Fisher's exact test for biased strand distribution p-value = 4.15e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCCTGTTGCCAGGCGATAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCA  >  NC_000913/2615593‑2615880                                                                                                                                             |                                                                                                                                                    tGCCCTGTTGCCAGGCGATAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGcac                                                                                                                                           <  2:272696/151‑1 (MQ=255)      tGTTGCCAGGCGATAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAg                                                                                                                                      <  2:782820/151‑1 (MQ=255)      tGTTGCCAGGCGATAGGCAGTACCACGAACGCCATCAGTAGCAATAGACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAg                                                                                                                                      <  1:338942/151‑1 (MQ=255)               gCGATAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGGGCAGCCAATAGATTGCAGACGCCCGGTTGGCCATTCCAGCAAAAAGGCCAGCACAAAAGCCACCAGCa                                                                                                                             >  2:185448/1‑151 (MQ=255)               gCGATAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGGGCAACCAAAAGAATGCAGACGCCCGGGTGGGCATTCCAGCAAAAAGGCCAGCACAAAAACCACCAGCa                                                                                                                             >  2:539259/1‑151 (MQ=255)                  aTAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGGGCAGCCAATAAATTTCAGACGCCCGGTTGGCCATTCCAGCAAAAAGGCCAACACAAAAACCACCCGCAACg                                                                                                                          >  2:796749/1‑151 (MQ=255)                   tAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACgg                                                                                                                         <  1:657818/151‑1 (MQ=255)                   aaGGCAGTCCCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGCCGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACgg                                                                                                                         <  2:794698/150‑1 (MQ=255)                     ggCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAAAAGAATGCAAACGCACGGTTGGCCATTTCAGCAAAAAGGCCAGCACAATAGCCACCCGCAACgggg                                                                                                                       >  2:696502/1‑149 (MQ=255)                      gCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGCGCAGCCCATAAATTGTAGACGCACGGTTGGCCATTCCAGCCAAAAGGCCAGCACAAAAGCCACCCGCAACgggg                                                                                                                       >  1:626451/1‑148 (MQ=255)                       cAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCCGACGCACGGTTGGCCATTCCAGCAAAAAGGCCAGCACAAAAGCCACCAGCAACggggga                                                                                                                     >  2:1051778/1‑147 (MQ=255)                         gTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCCGCCAATAGATTGCAAACGCACGGGTGGCCATTCCCGCAAACAGGCCCCCACAATAGCCACCAGCAACAGagcaag                                                                                                                   >  2:626992/1‑151 (MQ=255)                         gTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCCGCAAAAAGATTGCAGACGCACGGTTGGCCCTTCCAGCAAAAAGGCCAGCACAATTGCCACCACCCACGTagcaag                                                                                                                   >  1:553125/1‑151 (MQ=255)                         gTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAAAAGGCCCGCACAAAAACCACCAGCAACGGagcaag                                                                                                                   >  2:232461/1‑151 (MQ=255)                         gTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACCGTTGGCCCTTCCAGCCAAAAGGCCAGCACAATAACCCCCCGCAACGGagcaaa                                                                                                                   >  1:1014524/1‑150 (MQ=255)                         gTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAAACGCACGGTTGGCCATTCCCGCAAAAAGGCCAGCACAATAACCACCAGCAACGGagcaag                                                                                                                   >  2:320079/1‑151 (MQ=255)                          tACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGagcaagc                                                                                                                  <  1:791625/151‑1 (MQ=255)                          tACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGagcaagc                                                                                                                  <  1:951185/151‑1 (MQ=255)                                       aTCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAAAAGATTGCAAACGCACGGTTGGCCATTCCAGCAAAAAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGa                                                                                                     >  2:499022/1‑151 (MQ=255)                                       aTCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAACCAATAAATTGCAAACCCACGGTTGGCCCTTCCAGCAAAAAGGCCAGCACAAAAGCCCCCAGCAACGGAGCAAGCAGGCCACTAAAGa                                                                                                     >  2:372629/1‑151 (MQ=255)                                         cAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGGGCAGCCAATAGATTGCAGACGCACGGGTGGCCATTCCAGCAAAAAGGCCCGCACAATAACCCCCAGCAACGGAGCAAGCAGGCCCCTAaaaaaa                                                                                                   >  2:762502/1‑151 (MQ=255)                                         cAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTTCAAACCCACCGTTGGCCCTTCCCGCAAATAAGCCCGGACCAAAGCCACCAACAACGGAGCAAACACGCCCCCAAAGaaa                                                                                                   >  2:55593/1‑151 (MQ=255)                                           gTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGGGCAACCAATAGATTGCAGACCCACGGTTGGCCATTCCCGCAAAAAGGCCAGCACAAAAGCCACCAGCAAAGGGGCCAGCAGGCCCCTAAAGaaaaa                                                                                                 >  2:540106/1‑151 (MQ=255)                                           gTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAAAAGATTTCAGACGCACGGTTGGCCATTCCAGCAAAAAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCCGGCCACTAAAGaaaaa                                                                                                 >  1:608984/1‑151 (MQ=255)                                              gCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAAATTGCAGAAGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGaa                                                                                              >  1:264839/1‑150 (MQ=255)                                              gCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAAATTGCAAACGCACGGTTGGCCATTCCCGCAAAAAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAAAt                                                                                              >  1:284725/1‑151 (MQ=255)                                              gCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAAAAGATTGCAGACGCACGGGTGGCCCTTCCACCAAAAAAGCTAGCACAAAAACCACCAGCAACGGGGGAAACCGGGCACTCaaaaaaaaaaa                                                                                              >  2:704172/1‑148 (MQ=255)                                              gCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAACCAAAAAATTGCCCAAGCACGGTTGGCCATTCCCCCAAAAAAGCCCGCACAAAAGCCACCACCAAAGGAGCAAGCAGGCCACTAAAGaaaaaaaa                                                                                              >  2:482723/1‑148 (MQ=255)                                               cAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGGGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCCGCAAAAAGGCCAGCACAAAAGCCACCCGCAACGGAGCAAACAGGCCACTAaaaaaaaaaaaa                                                                                             >  2:61871/1‑147 (MQ=255)                                               cAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGAt                                                                                              <  1:964644/150‑1 (MQ=255)                                               cAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGAt                                                                                              <  1:681147/150‑1 (MQ=255)                                               cAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAAAAGATTTCAGACGCACGGGTGGCCATTCCAGCAAAAAGGCCAGCACAAAAACCACCAGAAACGGGGCAAGCAGGCCAACAaaaaaaaaaaaa                                                                                             >  2:441786/1‑147 (MQ=255)                                                   ataCCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGGAGACAACAGATTGCAAAAGCACGGCTGGCCCTTCCCGCAAATAGGCCCGCACAATAACCACCCGCAACGGAGCAAGCAGGCCCCCAAAAAAAAAGATAATg                                                                                          >  2:519769/1‑150 (MQ=255)                                                   ataCCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGAATTCAGACGCACGGTTGGCCATTCCAGCAAAAAGGCCAGCACAAAAACCCCCAGCAACGGAGCCAGCAGGCCACTAAAGAAAAAAAAAATg                                                                                          >  2:116865/1‑150 (MQ=255)                                                        gacgaAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAAc                                                                                                                           <  2:876269/112‑1 (MQ=255)                                                        gacgaAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAAAAAAATGCAGACGCACGGGTGGCCATTCCAGCAAAAAGGCCAGCACAAAAACCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAAATAAAGCCaaaa                                                                                    >  2:829464/1‑150 (MQ=255)                                                         ccgaAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATc                                                                                   <  1:619953/150‑1 (MQ=255)                                                          cgaAACCCACCAATACAATCGACGTCGCCCAGCAGCGGGAGCAGCCAATAGAATGCAGACGCACGGTGGGCCATTCAAGCCAATAGGCAAGCACAATAGCCACCCGCAACGGAGTAAGCAGGCCACTAAAGAAAAAGATAATGCCAACTcc                                                                                  <  1:314068/151‑1 (MQ=255)                                                                ccaccaATACAATCGACGTCGCCCAGCGGCGGGGGCAGCCAAAAAATTGCAGACCCACGGTTGGCCCTTCCAGCAAAAAGGGCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACCAAAAAAAAAAAAAAAGCCAAAACCCGcca                                                                             >  2:1014919/1‑148 (MQ=255)                                                                ccaccaATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGGTGGCCATTCCAGCAAAAAGGCCAGCACAATAACCACCAGCAACGGAGCAAACAGGCCACTAAAAAAAAAGATAATGCCAAATCCCGCaa                                                                             >  2:436556/1‑150 (MQ=255)                                                                 caccaATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAAATTGTAGACGCACGGGTGGCCATTTCAGCAAAAAGGCCAGCACAATAGCCACCAGCAACGGAGCCAGCAGGCCACTAAAGAAAAAGATAATGACAAATCCGGCaa                                                                             >  1:87157/1‑149 (MQ=255)                                                                      aTACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAAACGCACGGTTGGCCATTCCAGCAAAAAGGCCAGCCCAATAGCCACCAGCAACGGAGCAAACAGGCCACTAAAGAAAAAAATAATGCCAAATCCGGCAACTaaaaa                                                                      >  2:337630/1‑150 (MQ=255)                                                                         cAATCGACGTCGCCCAGCGGCGGGGGCGGCCAATAGATTGCAGACGCACGGGTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGGAACGGAGCAAGCAGGCCACAAAAGAAAAAGAAAAAGACCAAACCGGGAAATAaaaaaa                                                                    >  2:464888/1‑150 (MQ=255)                                                                          aaTCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAAcc                                                                  <  1:278247/151‑1 (MQ=255)                                                                          aaTCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAAcc                                                                  <  1:449860/151‑1 (MQ=255)                                                                                                   cagcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCg                                          >  2:234729/1‑150 (MQ=255)                                                                                                    agcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCg                                          >  2:701287/1‑149 (MQ=255)                                                                                                     gcCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGACAAAGATAATGCCAAATCCGGCAACTAAAATNNCCAGCAAGGCAATCGCTTCCGGGTCGCTa                                       >  1:842044/1‑151 (MQ=255)                                                                                                         aTAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTaaaa                                    >  2:881473/1‑150 (MQ=255)                                                                                                                 cAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCTCTTCCGGGTCGCTAAAACggcggcca                           >  2:335066/1‑149 (MQ=255)                                                                                                                     cGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATAcc                       <  2:95115/151‑1 (MQ=255)                                                                                                                         cGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTg                   <  1:16007/151‑1 (MQ=255)                                                                                                                         cGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATAACATTg                   >  2:228895/1‑151 (MQ=255)                                                                                                                                ccATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAAc            <  1:319982/151‑1 (MQ=255)                                                                                                                                         cAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAg    <  1:1051397/150‑1 (MQ=255)                                                                                                                                          aaaTAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCa  >  2:487132/1‑151 (MQ=255)                                                                                                                                          aaaTAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCa  >  2:957840/1‑151 (MQ=255)                                                                                                                                             |                                                                                                                                                    TGCCCTGTTGCCAGGCGATAGGCAGTACCACGAACGCCATCAGTAGCAATATACCGACGAAAACCACCAATACAATCGACGTCGCCCAGCGGCGGGAGCAGCCAATAGATTGCAGACGCACGGTTGGCCATTCCAGCAAATAGGCCAGCACAATAGCCACCAGCAACGGAGCAAGCAGGCCACTAAAGAAAAAGATAATGCCAAATCCGGCAACTAAAATAACCAGCAAGGCAATCGCTTCCGGGTCGCTAAAACGGCGGCGATACCATTGCATCAACATTTCGAGCA  >  NC_000913/2615593‑2615880

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