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

GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354141‑354436                                                                                                                                                    |                                                                                                                                                     ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTg                                                                                                                                                     <  1:555489/149‑1 (MQ=255) ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTg                                                                                                                                                     <  2:159625/149‑1 (MQ=255) ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGa                                                                                                                                                    <  2:566734/150‑1 (MQ=255)   ggCATCCGTGACGCTGACGGTTATCCCTTTATTCTCGCGCGCGCTGACGATGTCATTAACGTTGCCGGGCTTCGGCTTTGGACGCGCGAGATTGACGAGAGTATCTCCATCCATGCGGGCGTTGCCGAACCGGCGGTGGTTGCGGTGAAAg                                                                                                                                                 <  1:26394/151‑1 (MQ=255)   ggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAg                                                                                                                                                 <  1:664657/151‑1 (MQ=255)    gCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                <  1:471160/151‑1 (MQ=255)           tGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATgcg                                                                                                                                            >  2:408103/1‑148 (MQ=255)             aCGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcgc                                                                                                                                           >  2:366241/1‑147 (MQ=255)             aCGCGGACGGTTATCACTTTATATTCGGGCTCACTGACGATTTAAGTAACGTTGCCGGGCAGCGGCTGGGGTCGCGCGAGATTGAAGAGAGTATCTCCAGCCTTCCGGGGGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGGTGaa                                                                                                                                       <  1:639793/151‑1 (MQ=255)                  gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAggg                                                                                                                                   <  2:688015/150‑1 (MQ=255)                      gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAgg                                                                                                                               >  2:200295/1‑150 (MQ=255)                         aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAggggg                                                                                                                            >  2:240957/1‑150 (MQ=255)                          tCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGCAggtggcg                                                                                                                          >  1:8793/1‑151 (MQ=255)                           cACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcg                                                                                                                          <  1:463638/150‑1 (MQ=255)                           cACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGAAAGGGGGgggggcgg                                                                                                                         >  2:424406/1‑151 (MQ=255)                               ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtg                                                                                                                       <  1:441841/149‑1 (MQ=255)                               ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCCCCGGAAAGGCAggggg                                                                                                                            >  1:72743/1‑144 (MQ=255)                                tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtggc                                                                                                                     <  2:374012/150‑1 (MQ=255)                                tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGAAGCGCCGAAAggggcggggggggggg                                                                                                                      >  1:216286/1‑136 (MQ=255)                                 aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCCTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAAGTGCGCTCACAGGGAagg‑gggggtggg                                                                                                                     >  1:434072/1‑147 (MQ=255)                                  ttCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTAACCCCAGCCATGCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATTCGCTGCACGGTGAGgggggggttgcggt                                                                                                                  >  1:335751/1‑149 (MQ=255)                                   tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGtt                                                                                                                  <  2:263075/150‑1 (MQ=255)                                   tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGtt                                                                                                                  <  2:511975/150‑1 (MQ=255)                                   tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGACGCGCTGTAAGGGGAGGGGGGGGGGGGGt                                                                                                                   >  1:498304/1‑149 (MQ=255)                                         gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGCATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCAtt                                                                                                           <  2:476612/151‑1 (MQ=255)                                            cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAAGCGCTGAAAGGGGCGGGGGGGGGGGGGGTTGTCATTcc                                                                                                         >  1:606526/1‑150 (MQ=255)                                                gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaaa                                                                                                     <  1:501962/150‑1 (MQ=255)                                                gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaaa                                                                                                     <  1:466673/150‑1 (MQ=255)                                                gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                                                                                    <  2:683686/151‑1 (MQ=255)                                                   gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagaga                                                                                                 <  2:515598/151‑1 (MQ=255)                                                   gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGaaa                                                                                                     >  2:505836/1‑147 (MQ=255)                                                    aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGGAaaag                                                                                                  >  2:173562/1‑149 (MQ=255)                                                              cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAAATGGGGTGGAAGGGCGGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGGGACAGTCt                                                                                       >  2:468612/1‑150 (MQ=255)                                                                ttGCCGGGCATCGGCTGGGGACGCGCGCGATTGACGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCTGTTGGGGGGGAAGAAGCGCCGGAAGGGCAGGTGGGGGGGGGGGTTGTTATTCCGCAAGAGAGCGACa                                                                                           >  1:85123/1‑144 (MQ=255)                                                                ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGTAAAGTGTGCTCGAAGGGGACGCGGGGGGGGCGGTTGTTATTCAGAAAGAGAGCGatcggctgg                                                                                     >  2:387240/1‑142 (MQ=255)                                                                                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgc                                                                     <  1:645149/150‑1 (MQ=255)                                                                                    aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACt                                                                <  2:561693/151‑1 (MQ=255)                                                                                     cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGAAAGGGGCGGGGGGGGGGGCGGTTGTCCTTTCGGAAGAGAGCGACAGGCTGGGAGAACGTGAGGTGGGGCAc                                                                 >  1:487800/1‑149 (MQ=255)                                                                                          gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGGGCACTCGc                                                             >  2:251523/1‑148 (MQ=255)                                                                                               tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGATAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGaa                                                      <  2:661065/150‑1 (MQ=255)                                                                                               tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGGGCTGGAAGGGGAGGGGGGGGGGGCGTTTGTCATTTCGGAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGGGCACTCGCAagag                                                        >  2:250410/1‑148 (MQ=255)                                                                                                gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAg                                                     <  2:505675/150‑1 (MQ=255)                                                                                                            tcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGAAAGGGCAGGTGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCCCTCGCAAGAGAAGGCGATTATgg                                           >  2:188939/1‑148 (MQ=255)                                                                                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACa                               <  1:37592/151‑1 (MQ=255)                                                                                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGGGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggggg                                  >  2:89802/1‑148 (MQ=255)                                                                                                                       gggCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACa                               >  2:212752/1‑149 (MQ=255)                                                                                                                               ccGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTgg                     >  2:111489/1‑151 (MQ=255)                                                                                                                                                 gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgc   <  1:150333/151‑1 (MQ=255)                                                                                                                                                 gggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgc   <  1:274451/151‑1 (MQ=255)                                                                                                                                                  ggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcg  >  2:416726/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354141‑354436

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