Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	52.1%	43.7 / 69.8	48	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (5/18);  new base G (25/0);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 4.60e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-05

GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354141‑354437                                                                                                                                                     |                                                                                                                                                       ggggCATCCGGGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGGCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTg                                                                                                                                                       <  2:408362/149‑1 (MQ=255)    gCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGa                                                                                                                                                  <  2:472646/151‑1 (MQ=255)     cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGAt                                                                                                                                                 <  1:196737/151‑1 (MQ=255)      aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGGAAGATg                                                                                                                                                >  2:952898/1‑151 (MQ=255)        ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATgcg                                                                                                                                              <  2:648799/151‑1 (MQ=255)          gggACGCTGACGGGTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATgcg                                                                                                                                              <  2:441559/147‑1 (MQ=255)                   aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCa                                                                                                                                  <  1:550482/151‑1 (MQ=255)                      gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGGAAGAAgcggcgcaaagggcgg                                                                                                                                 >  2:105810/1‑137 (MQ=255)                         aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAggtgg                                                                                                                             <  1:414189/150‑1 (MQ=255)                               ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAggtggcggtgg                                                                                                                       <  2:5012/150‑1 (MQ=255)                                tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑AGGTTGGGGGGGAAGATGCGCTGAAA‑GGGCCGgtggcgggggcg                                                                                                                     >  1:569897/1‑151 (MQ=255)                                 aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGGAggggggggtg                                                                                                                        >  1:914199/1‑147 (MQ=255)                                     tcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGCGGCGG‑GGGGTGGGGGGAAAGAAAGCCGGAAA‑GGGCAGGTGGCGGGGGCGGTTGt                                                                                                                >  2:82924/1‑151 (MQ=255)                                       gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCa                                                                                                              <  1:762780/151‑1 (MQ=255)                                       gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGAAAGATGCGCTGGAA‑GGGGCGGGGGGGGGGGGGttt                                                                                                                  >  2:333293/1‑147 (MQ=255)                                         gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑GGGGTGGGGGGGAAGATGCGCTGGAA‑GGGCAGGGGGGGGGGGCGttttttctt                                                                                                            >  2:159503/1‑145 (MQ=255)                                           gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTcc                                                                                                          <  1:306779/151‑1 (MQ=255)                                           gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGGGGGGGGGGGGGTTGTCAtt                                                                                                            >  1:198429/1‑149 (MQ=255)                                           gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGACGAAGTGGCGGCGGGTGGGGGGGACAAATCGGCGGAAA‑AGG‑AGGAGGCGGGGGCGTTTggtctt                                                                                                            >  2:989474/1‑144 (MQ=255)                                            cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGAAAGATGCGCTGGAA‑GGGCAGGGGGGGGTGGCGGTTGTCAtt                                                                                                            >  1:122197/1‑148 (MQ=255)                                             aCTGACGATGTAATTAACGATGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGGAAGAAGCGCTGGAA‑GGGCAGGTGGGGGGGGGGTTTGGCATTcc                                                                                                          >  1:401806/1‑149 (MQ=255)                                                 aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑AGGATGGGGGGAAAAAGGGGCCGAAA‑AGGAAgggggggggggggcttttgatttcg                                                                                                         >  1:341897/1‑132 (MQ=255)                                                     tGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGAAAGATGCGCTGAAA‑AGGGAGGGGGCGGTGGCGGTTGTTATTCCGAAaaagag                                                                                                 >  2:827374/1‑150 (MQ=255)                                                     tGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAagagag                                                                                                 >  1:531074/1‑150 (MQ=255)                                                           tAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAg                                                                                           <  1:633575/150‑1 (MQ=255)                                                             aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑GGGGTGGGGGGAAAGAAGCGCTGGAA‑GGGCAGGTGGGGGGGGCGGTTGTCATTTCCAAaaaga                                                                                                  >  2:100143/1‑141 (MQ=255)                                                              cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑GGGTTGGGGGGGAAGATGCGCTGGAA‑GGGCAGGGGGGGGGGGGGGTTTGCATTTCGGAag                                                                                                     >  2:138724/1‑137 (MQ=255)                                                                    cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                                 <  2:882759/151‑1 (MQ=255)                                                                        cATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGGGGCGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATc                                                                                >  1:629523/1‑148 (MQ=255)                                                                         aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCCCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑GGGGTGGGGGGAAAAAAGAGCTGAAAGGGGCGGGGGGGGGGGGCG‑GTGTCATTCCGAAAGAGAGTGAGAGTCTGGGGGATCg                                                                               >  2:792338/1‑148 (MQ=255)                                                                            ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAtgt                                                                         <  1:952947/151‑1 (MQ=255)                                                                               tGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTgg                                                                       >  1:50600/1‑150 (MQ=255)                                                                               tGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑GGGGTTGGGGGAAAGATGCGCTGAAA‑AGGGAGGGGGGGGGGGGGTTTTTCATTCCGAAAGAGAGCGGCAGGCTGGAAGATTGGGATGTgg                                                                       >  1:1012387/1‑150 (MQ=255)                                                                                  ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGGAAGAAGCGCTGAAA‑GGGGAGGGGGCGGGGGCGTTTGTCATTCCGAAAAAGAGCGACAGGCTGGGAGATCGTGATGTGgggc                                                                    >  1:912015/1‑150 (MQ=255)                                                                                  ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgcgc                                                                    >  1:6697/1‑150 (MQ=255)                                                                                  gaaCGCGCGAGAGTGAAGAGAGTATCTCCTGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATACGCTGAAA‑GGGCAGGTGGCGGTTGCCTATATCATTCCGAAAGAGAGCGACAGCCTGGAAAATCGTGCTGTGGCGCa                                                                   <  2:276812/149‑1 (MQ=255)                                                                                   gACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGGAAGAAGCGCTGAAA‑GGGCAGGGGGCGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTTGAAGATCGTGGTGTGgcg                                                                     >  1:883000/1‑148 (MQ=255)                                                                                              ttGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGGAAAATGCCCTTAAA‑GGGCAGGGGGCGGTGGCGGTTGGCCTTCCCGAAGAAAGCGACAGTCTTGAAGAGCGCGGTGTGGCGCACCCGcca                                                            >  2:63004/1‑144 (MQ=255)                                                                                                 aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGc                                                    <  1:27777/151‑1 (MQ=255)                                                                                                   gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGa                                                  <  1:544939/151‑1 (MQ=255)                                                                                                         aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATg                                             <  2:562646/150‑1 (MQ=255)                                                                                                          tctcCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGAAAGATGCGCTGAAA‑GGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATgg                                            >  1:185919/1‑150 (MQ=255)                                                                                                           ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGGGGAAAATGCGCTGAAA‑GGGGAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGAGGGTCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcg                                          >  1:119096/1‑151 (MQ=255)                                                                                                              cAGCCATCCGGGCGTTGCCGAAGTGGCGG‑CGGGTGGGGGGAAAGATGCGCCGCAA‑CGGCAGGTGGGGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGGCTTGGAGATGGCGGTGTGGCGCACTCGCAAGAGAAGGCGCTTATGgcgccg                                       >  1:133814/1‑149 (MQ=255)                                                                                                                      cGGGCGTTGCCGAAGTGGCGG‑TGGGTGGGGGGAAAGATGCGCTGAAA‑GGGGAGGGGGGGGGGGGGTTTGTCATTTCGAAAGAGAGGGACAGTCTGGAAGATCGGGGTGTGGGGCGCTCGCAAGAGAAGGGGATTATGGGGCtggggg                                   >  1:786553/1‑147 (MQ=255)                                                                                                                               ccGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTgg                      <  1:1067541/151‑1 (MQ=255)                                                                                                                                                 tgggTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg     >  2:142931/2‑151 (MQ=255)                                                                                                                                                  gggTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgc    >  1:914268/1‑151 (MQ=255)                                                                                                                                                     tGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  >  2:579103/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGG‑TGGTTGGGGTGAAAGATGCGCTGAAA‑GGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  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 < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑