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

GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑354437                                                                                                                                                    |                                                                                                                                                       ggggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGa                                                                                                                                                      <  1:27835/150‑1 (MQ=255)  gggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGaaa                                                                                                                                                    <  1:372733/151‑1 (MQ=255)     cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                                                 <  2:356407/151‑1 (MQ=255)      aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATg                                                                                                                                                <  1:743190/151‑1 (MQ=255)      aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAGGTg                                                                                                                                                >  1:248965/1‑151 (MQ=255)                 tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCT‑GGAAggg                                                                                                                                    >  1:203961/1‑151 (MQ=255)                       ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAAATGTGCT‑GAAAGGGGAggtg                                                                                                                              >  2:516146/1‑151 (MQ=255)                        tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGCAggtgg                                                                                                                             >  2:752641/1‑151 (MQ=255)                         aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAggtgg                                                                                                                             <  1:44322/150‑1 (MQ=255)                             cTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGCCgggggcgggg                                                                                                                        >  2:90534/1‑149 (MQ=255)                              tttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCC‑GGAAggggggggggggg                                                                                                                          >  1:558168/1‑138 (MQ=255)                                tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGGAgggggcggggggg                                                                                                                     >  2:931179/1‑149 (MQ=255)                                      cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GGAAGGGGGGGGGGGGGGGGCGgtt                                                                                                                  >  1:866622/1‑148 (MQ=255)                                       gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAAAGGGCT‑GGAAGGGCAggggggggggggggtt                                                                                                                  >  1:22024/1‑142 (MQ=255)                                           gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTcc                                                                                                          <  1:219208/151‑1 (MQ=255)                                              cTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAAGGGAGGTGGCGGGGGCGGTTGTCATTTCGa                                                                                                        >  2:314567/1‑150 (MQ=255)                                                 aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGaa                                                                                                       <  1:431921/148‑1 (MQ=255)                                                             aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCt                                                                                        <  2:31351/151‑1 (MQ=255)                                                             aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGGAGGGGGCGGGGGGGTTTGTCCTTTCGAAAAAGAGCGGAAGTCt                                                                                        >  1:965381/1‑151 (MQ=255)                                                               gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTg                                                                                       <  2:354823/150‑1 (MQ=255)                                                                ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTGGGGGGGGAAGAAGCGCTGGAAAGGGC‑GGGGGGGGGGGGGATTGGCATTTCGGAagagcgggagcggaccgg                                                                                     >  1:234858/1‑137 (MQ=255)                                                                 tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa                                                                                     <  2:409943/150‑1 (MQ=255)                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGaa                                                                                    <  1:610810/150‑1 (MQ=255)                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                   <  2:808001/151‑1 (MQ=255)                                                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAAAAGAGCT‑GAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGGGCGACAGTCTgg                                                                                      >  2:818730/1‑148 (MQ=255)                                                                        cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGGAGGTGGCGGGGGGGGTTTTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCggg                                                                             >  2:782351/1‑149 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGGAGGGGGGGGGGGGGGTTGTTATTTCGGAAGAGAGCGACCGTCTGGAAGATCg                                                                               >  2:70003/1‑148 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCT‑GGAAAGGGAGGGGGGGGGGGGGGTTGTTATTCCGAAAGAGAGGGACAGGCTGGGAGGTc                                                                                >  1:855569/1‑147 (MQ=255)                                                                             gCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCT‑GAAAGGGGCGGGGGGGGGGGGGGTTTTCATTTCCAAAAAGAGCGACAGGCTGGAAGGTCGTGAtg                                                                          >  1:554865/1‑149 (MQ=255)                                                                                    aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGCGCAc                                                                  >  2:260692/1‑150 (MQ=255)                                                                                     cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                                <  2:279023/151‑1 (MQ=255)                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCCGAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGGGCACTCGc                                                              >  1:602529/1‑149 (MQ=255)                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCC‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGGGCACTccc                                                              >  1:389955/1‑147 (MQ=255)                                                                                         cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTTGAAGATCGTGATGTGGCGCACTc                                                                >  2:688523/1‑147 (MQ=255)                                                                                          gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGGAGGTGGCGGGGGCGTTTGGCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATTTGGCGCACTCGc                                                              >  2:119816/1‑148 (MQ=255)                                                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAaga                                                          <  2:902200/151‑1 (MQ=255)                                                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTTCGAAAGAGAGCGACAGGCTGGGAGATCGTGATGTGGCGCCCTc                                                                >  1:160050/1‑145 (MQ=255)                                                                                             aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCC‑GGAAGGGCAGGGGGCGGGGGGGTTTGTCATTCCGGAAGAGAGCGACAGTCTGGAAGATCGGGGTGTGGCGCACTCGcca                                                            >  1:961356/1‑145 (MQ=255)                                                                                                         aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATg                                             <  2:736242/150‑1 (MQ=255)                                                                                                           ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCT‑GGAAAGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAAATCGTGGTGTGGCGCACTCGCAAGAGAAGGGGATTATGgcg                                          >  2:336972/1‑151 (MQ=255)                                                                                                              cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtg                                       <  2:380212/151‑1 (MQ=255)                                                                                                              cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtg                                       <  2:334444/151‑1 (MQ=255)                                                                                                                 ccATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGCAGGGGGGGGGGGCGGTTGTCCTTCCGGAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCCCTCGCAAGAGAAGGCGGTTATGGCGCtgggg                                    >  1:495411/1‑149 (MQ=255)                                                                                                                  cATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                                   >  1:769472/1‑151 (MQ=255)                                                                                                                     ccGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGAc                                 <  2:129892/150‑1 (MQ=255)                                                                                                                      cGGGCGTTGCCGAAGTGGCGGAGGGTGGGGTGGAAAATGCGCT‑GAAAAGGGAGGGGGCGGTGGGGGTTGTCGTTTCGAAAAAAAACGACAGTCTGGGAGAACGGGGTGTGGGGGACTCGCGAGAGAAAGCGAATATGGCGCTGGGGGAc                                 >  1:96040/1‑149 (MQ=255)                                                                                                                                  aaGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCa                    >  1:137012/1‑150 (MQ=255)                                                                                                                                  aaGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCa                    >  2:791016/1‑150 (MQ=255)                                                                                                                                       gCGGTGGTTGGGGTGAAAGATGCTCT‑CAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGCGGCCCACTCGCATGAGAAGGCGATTATGGCCCTGGTGGACAGCCTGATCGGCAACtttc              <  1:549333/151‑2 (MQ=255)                                                                                                                                       gCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg              <  2:635342/151‑1 (MQ=255)                                                                                                                                        cGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTgg             <  1:90530/151‑1 (MQ=255)                                                                                                                                            ggtTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccgc         >  1:121810/1‑151 (MQ=255)                                                                                                                                               tGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCg      >  2:167061/1‑151 (MQ=255)                                                                                                                                                   gTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  <  2:22024/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑354437

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