Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	46.9%	76.2 / 11.4	49	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (3/23);  new base G (23/0);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 5.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCC  >  NC_000913/354138‑354431                                                                                                                                                       |                                                                                                                                                aCTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGt                                                                                                                                                 <  2:708815/151‑1 (MQ=255)          tCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATg                                                                                                                                         >  2:559959/1‑150 (MQ=255)          tCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGATg                                                                                                                                         >  2:504675/1‑150 (MQ=255)                   cTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAgg                                                                                                                               >  2:335235/1‑151 (MQ=255)                            aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggc                                                                                                                      <  2:974810/151‑1 (MQ=255)                              cACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGGAAggggcgggggggg                                                                                                                    >  2:946190/1‑141 (MQ=255)                                  ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtg                                                                                                                  <  2:293501/149‑1 (MQ=255)                                  ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGTTGGGGGGGAAGATGCGCCGAAAGGGCAgggggcggggg                                                                                                                 >  2:525123/1‑150 (MQ=255)                                    aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCTggtagcggttgc                                                                                                                <  2:705942/149‑1 (MQ=255)                                      tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGGAggggggggggggggtt                                                                                                            >  1:319680/1‑146 (MQ=255)                                         cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGt                                                                                                          <  1:309845/150‑1 (MQ=255)                                                aCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAAGGGCAGGGGGGGGGGGCGGTTGTCATTCCg                                                                                                   >  2:99199/1‑150 (MQ=255)                                                   gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaaa                                                                                                <  1:773805/150‑1 (MQ=255)                                                   gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGaa                                                                                                 >  1:704563/1‑149 (MQ=255)                                                       aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagagag                                                                                           <  2:705366/151‑1 (MQ=255)                                                           aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGa                                                                                        <  2:76932/150‑1 (MQ=255)                                                                aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAATGCGCTGGAAGGGGAGGGGGGGGGGGGGGTTGTCATTTCGAAAGAGGGCg                                                                                         >  2:609737/1‑144 (MQ=255)                                                                aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGTGGGGGTTGTCATTCCGAAAGAGAGGGACAGTc                                                                                   >  1:275706/1‑150 (MQ=255)                                                                 cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGGCCTTTCGAAAGAGAGCGACAGt                                                                                    >  2:54212/1‑148 (MQ=255)                                                                  gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAGATGCGCTGCAAGGGCAGGGGGCGGGGGGGTTTGTTATTCCGAAAGAGAGCGACAGTCTgg                                                                                >  2:224405/1‑151 (MQ=255)                                                                     gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGaa                                                                              <  1:582809/150‑1 (MQ=255)                                                                     gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                             <  2:505467/151‑1 (MQ=255)                                                                      ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGa                                                                            <  1:528238/151‑1 (MQ=255)                                                                       cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                >  1:253720/1‑146 (MQ=255)                                                                        gggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATc                                                                          <  1:155768/151‑1 (MQ=255)                                                                               ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGGGTTTGTTATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTggt                                                                     >  1:162518/1‑147 (MQ=255)                                                                               ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGGAGGGGAGGTGGGGGTGGGGGTTGTCATTCCGAAAGAGAGGGACAGTCTGGgagaac                                                                          >  2:175510/1‑142 (MQ=255)                                                                                gCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGTTTGTCATTCCGGAAGAGAGCGACAGTCTGGAAGATCGTGAtg                                                                    >  2:990057/1‑149 (MQ=255)                                                                                          cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                        <  2:121959/151‑1 (MQ=255)                                                                                           gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCACTc                                                          >  2:679989/1‑148 (MQ=255)                                                                                            cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGTGATGTGGCGCACTCGCa                                                       >  1:596015/1‑150 (MQ=255)                                                                                               gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAaga                                                    <  2:690585/150‑1 (MQ=255)                                                                                               gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGAAGATCGTGATGTGGCGCACTCGCaa                                                      >  2:152358/1‑148 (MQ=255)                                                                                                  tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGaa                                                 <  2:199797/150‑1 (MQ=255)                                                                                                   gAAGAGAGTATCTCCAGCCATCCGGGCGTTGGCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAg                                                <  2:864262/150‑1 (MQ=255)                                                                                                    aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGc                                              <  2:703219/151‑1 (MQ=255)                                                                                                    aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAgg                                               >  2:809637/1‑150 (MQ=255)                                                                                                    aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAAGGGCAGGGGGCGGGGGGGGTTGTCATTTCGAAAGAGAGGGACAGTCTGGAAGATCGGGATGTGGCGCACTCGCCAgaga                                                  >  2:952510/1‑147 (MQ=255)                                                                                                      gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGa                                            <  2:941177/151‑1 (MQ=255)                                                                                                        gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGTTTGTCATTTCGAAAGAGAGCGACAGGCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGGtt                                          >  1:622307/1‑151 (MQ=255)                                                                                                             tctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGGAAGGGGAGGGGGCGGGGGCGTTTGTCATTTCGAAAGAGAGGGAAAGTCTGGAAGATCGTGATGTGGGGCACTCGCAAGAGAAGGCGATTATgg                                      >  2:571300/1‑150 (MQ=255)                                                                                                                ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgc                                     <  2:802397/148‑1 (MQ=255)                                                                                                                    ccATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtg                              <  1:585373/151‑1 (MQ=255)                                                                                                                        ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACa                          <  1:18864/151‑1 (MQ=255)                                                                                                                            gCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCa                      <  2:552051/151‑1 (MQ=255)                                                                                                                                   cGAAGTGGCGGTGGTTGAGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGc               <  2:828441/151‑1 (MQ=255)                                                                                                                                    gAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTGTCATTCCGGAAGAGAGGGACAATCTGGGAGAACGGGATGTGGCGCACTCGCAAGAGAAAGCGATTATGGGGCTGGTGGACAGCCAGATTgg                >  2:954065/1‑149 (MQ=255)                                                                                                                                         ggCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt         <  1:285748/151‑1 (MQ=255)                                                                                                                                          gCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg        >  2:301966/1‑151 (MQ=255)                                                                                                                                                gtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgcc  >  1:262392/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                ACTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCC  >  NC_000913/354138‑354431

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