Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	43.4%	63.6 / 12.8	54	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (12/18);  new base G (23/0);  total (36/18)
Fisher's exact test for biased strand distribution p-value = 1.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508511‑3508806                                                                                                                                                    |                                                                                                                                                      tCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCGC‑CGcaaa                                                                                                                                                    >  2:251687/1‑150 (MQ=255)    aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTc                                                                                                                                                <  2:937991/151‑1 (MQ=255)     cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGgaaaa                                                                                                                                                   >  1:745612/1‑147 (MQ=255)        tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTcaca                                                                                                                                             <  1:114832/150‑1 (MQ=255)             ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGGGTCTCCTGCAAACCGGAGGGAGGAACGAAAAGAGCGACCATTTTGTGGCGGTGGACCGC‑GGGAAAGATTacacaca                                                                                                                                         >  2:116842/1‑148 (MQ=255)              tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGt                                                                                                                                       <  1:1010603/150‑1 (MQ=255)               aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGC‑GGGCAAATTCacacac                                                                                                                                          >  2:55518/1‑146 (MQ=255)                  gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGGGC‑GGGAAAAATTacacacactaaacc                                                                                                                                  >  2:467139/1‑143 (MQ=255)                       agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCcaca                                                                                                                             <  1:43603/151‑1 (MQ=255)                       agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCcaca                                                                                                                             <  1:463260/151‑1 (MQ=255)                        gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAAACACGGTAGccccccc                                                                                                                              >  2:874300/1‑147 (MQ=255)                        gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATAACAAACAGGAGccccc                                                                                                                                >  2:960253/1‑147 (MQ=255)                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑GGGAAAAATTAAACACCGGAGccccccc                                                                                                                              >  1:1056637/1‑145 (MQ=255)                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTACACCCGGGAGccccccc                                                                                                                              >  2:945395/1‑145 (MQ=255)                          gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGGAGccccccc                                                                                                                              >  1:1009953/1‑145 (MQ=255)                           aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACCGGAGCCCCcac                                                                                                                              >  1:104677/1‑146 (MQ=255)                             cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAAAAACGGGAAccccccc                                                                                                                              >  2:789488/1‑142 (MQ=255)                               atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCCC‑GGGAAAAATCACAAACCGTATGCCCCACAACAGaaa                                                                                                                      >  1:83420/1‑150 (MQ=255)                                  tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACCCCGGAGCCCCCCCCCCAGaaa                                                                                                                      >  2:1063841/1‑147 (MQ=255)                                   gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCGC‑GGGAAAAATTAAACACGGTAGACCCCCCATCAAAAAACggg                                                                                                                 >  2:464584/1‑149 (MQ=255)                                         gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAAACACGGGAGCCCCCCCATCAGAAAGCgggggggg                                                                                                            >  1:707721/1‑146 (MQ=255)                                            cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAACAATGCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                        <  2:177394/151‑1 (MQ=255)                                               gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACAcc                                                                                                      <  1:63450/150‑1 (MQ=255)                                                    aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTTACAAACGGTAACCCCCCCCTCCGAAAGCGTGGGGGGcaca                                                                                                        >  2:192567/1‑143 (MQ=255)                                                    aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                <  2:105232/151‑1 (MQ=255)                                                      ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCg                                                                                              <  2:620/151‑1 (MQ=255)                                                          tAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTACACACGGTAGCCCCCCCAACAGAAAACGGGGGGGGGACCCCGAGCCGCGTAg                                                                                           >  1:56192/1‑150 (MQ=255)                                                              ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATACAAACACGGGAGCCCCCCCAACCGAAAGAGGGGGGCGGCCACCCAGCTGCGGCGCGatc                                                                                      >  2:244961/1‑151 (MQ=255)                                                              ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                      <  1:26475/151‑1 (MQ=255)                                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGTGCGCACACCCAGCTGCGTAGCGATCATccc                                                                                 >  1:797032/1‑151 (MQ=255)                                                                   ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCCCGCGGGAAAAATACA‑AAACGGCAACCCCCCCAACCGGCAACCGGGGGGGCACACCCAACTTCGTGGCGATCATccc                                                                                 >  2:82609/1‑151 (MQ=255)                                                                        aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCAAACACGGTAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCCAGCTGCGGAGCGATCATcccct                                                                               >  1:470381/1‑146 (MQ=255)                                                                                    gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                                                                <  2:789378/151‑1 (MQ=255)                                                                                            tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACAAACCGGAGCCCCCACATCAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTGCTTATAAAACAGCCCTCGcacc                                                        >  2:39461/1‑151 (MQ=255)                                                                                            tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACACAGTAGCCCCCACATCAGAAAACGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGcacc                                                        >  1:543560/1‑151 (MQ=255)                                                                                               aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGAAAAATTTACACACGGGAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCCCGCTGCGGAGCCGTCATCCCCTGCTGATAAAACAACACCGGgaccac                                                      >  1:816936/1‑150 (MQ=255)                                                                                                 aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTCAAACACGGGAGCCCCCACATCAGAAAGCGTGGGGGGCACACCGAGCTTCGTAGCGATCATCCCCTTCTGGTAAAACAGCCCTGGCACCCCTTc                                                   >  2:788746/1‑151 (MQ=255)                                                                                                       gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACCTCGCggg                                              >  2:669718/1‑150 (MQ=255)                                                                                                        ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATTAAAAAAGGGAGCCACCACAACAGAAAGAGTGGTGGGGAAACCCAGATGGGTGGCGATTATCCCTTTCTGGTAAAACAGCACTGGGCACACTTCGCggg                                              >  2:432532/1‑149 (MQ=255)                                                                                                             aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATTACAAACGGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACAACCAGCTGCGTAGGGAACATCCCCTGCTGGTAAACCAGCACTGGCCCCCCTTCGCGGGGAtgggg                                       >  2:897228/1‑149 (MQ=255)                                                                                                              aagaaaagaGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAATTTACACACGGTAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCCGCTGCGGAGCGCTCATCCCCTGTTGCTAAACCCGCACTCGCCCCACCTCGCGGGTGtgg                                         >  2:833105/1‑148 (MQ=255)                                                                                                               agaaaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACGGGAGCCCCCCCATCCGAAAGGGGGGGGCGCCCACCCAGACGCGGAGCGATCATCCCTTGCTGTTAAAACAGCACAGGCACCCCTTCGCGGGGAtggggg                                      >  2:612258/1‑150 (MQ=255)                                                                                                                   aagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGCAAAAATCACACACGGTAGCCCCCCCCTCAGAAAGCGGGGGGGGGCCCCCCAGCTGCGTAGCGATCATCCCCTGCTGGTAAAACAACACTGGGCCCACCTCGCGGGGAtgg                                         >  2:42709/1‑143 (MQ=255)                                                                                                                         gACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACTTCGCGGGGATGGTGGCTGTGACCaa                            >  1:60670/1‑150 (MQ=255)                                                                                                                          aCCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                          <  1:58987/151‑1 (MQ=255)                                                                                                                                  gtggcggtggcgCGC‑GTAAAAATTCACCCACGGTAGCCCCCTCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                  <  1:306859/151‑1 (MQ=255)                                                                                                                                  gtggcggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                   <  2:358396/150‑1 (MQ=255)                                                                                                                                   tggcggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                  <  1:117003/150‑1 (MQ=255)                                                                                                                                    ggcggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                >  2:377541/1‑151 (MQ=255)                                                                                                                                      cggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAACGCGTGGTGCGCACACCGAGCTGCGTAGCGATCCTCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGCCGGATCGTTg              >  2:174717/1‑151 (MQ=255)                                                                                                                                             gcgc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt       <  2:910670/151‑1 (MQ=255)                                                                                                                                             gcgc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt       >  2:171440/1‑151 (MQ=255)                                                                                                                                              cgc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc      >  2:617990/1‑151 (MQ=255)                                                                                                                                               gc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc      <  2:122809/150‑1 (MQ=255)                                                                                                                                                c‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc     >  2:581672/1‑150 (MQ=255)                                                                                                                                                   gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  <  1:464827/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508511‑3508806

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