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%	76.4 / 17.3	51	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/26);  new base G (22/0);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 3.11e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354144‑354436                                                                                                                                                  |                                                                                                                                                     gCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGa                                                                                                                                                <  2:706018/151‑1 (MQ=255)  cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGAt                                                                                                                                               <  1:301763/151‑1 (MQ=255)       gTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTGTGGGGGGaggaagcg                                                                                                                                             >  1:450187/1‑142 (MQ=255)       gTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCt                                                                                                                                          <  2:441513/151‑1 (MQ=255)       gTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCt                                                                                                                                          <  2:247602/151‑1 (MQ=255)            gCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGG‑TGGGGGGAAAGATGCGCTGaa                                                                                                                                       >  2:369918/1‑149 (MQ=255)             cTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGG‑TTGGGGGGAAGAGGGGCTGGAAgg                                                                                                                                    >  1:424322/1‑151 (MQ=255)              tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTTGA‑TGGGGGGGAAGATGCGCTGAAAggg                                                                                                                                   >  2:492939/1‑151 (MQ=255)              tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGAAAGATGGGCTGAAAggg                                                                                                                                   >  2:778692/1‑151 (MQ=255)               gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGGAAGAAGGGCTGGAAgggg                                                                                                                                  >  1:143055/1‑150 (MQ=255)                aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGG‑TGGGGGGGAAGATGGGCTGGAAggg                                                                                                                                   >  1:246313/1‑149 (MQ=255)                 cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAg                                                                                                                                <  1:505296/151‑1 (MQ=255)                    ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGGAAAATGGGCTGGAAGGGCggggg                                                                                                                             >  1:63141/1‑149 (MQ=255)                      aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAggtgg                                                                                                                            <  2:299502/150‑1 (MQ=255)                       tCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGGGGAAAATGCGCTGAAAAGGCAggaggcg                                                                                                                          >  2:455262/1‑151 (MQ=255)                              aTTCTCGGGCGCACTGACGATGTAATTACCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtggc                                                                                                                     <  1:396989/149‑1 (MQ=255)                                tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGG‑TTGGGGGGAAGAAGCGCCGGAAGGGGAGggggcgggggcggt                                                                                                                  >  1:174021/1‑148 (MQ=255)                                  tcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGGAAGATGCGGCGGAAGGGGAggggggggggggggtt                                                                                                                 >  1:303909/1‑144 (MQ=255)                                   cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGt                                                                                                               >  2:108104/1‑150 (MQ=255)                                      gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGAAAGAAGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTTTCAtt                                                                                                           >  2:506987/1‑151 (MQ=255)                                       cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGGGGAAGATGCGCTGGAAGGGGGGGGGGCGGGGGGGGTTGTCAtt                                                                                                           >  1:857633/1‑150 (MQ=255)                                       cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCAtt                                                                                                           >  2:512393/1‑150 (MQ=255)                                         cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCg                                                                                                        <  1:611918/151‑1 (MQ=255)                                             gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                                                                                    <  1:856773/151‑1 (MQ=255)                                               cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                                                                                  <  1:579541/151‑1 (MQ=255)                                                  tGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGAAAGAAGGGCTGGAAGGGGAGGGGGCGGGGGGGGTTGTCATTCCGAAagaga                                                                                                 >  2:193423/1‑149 (MQ=255)                                                     aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTTCGAAagagaggg                                                                                              >  2:524292/1‑147 (MQ=255)                                                         aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTc                                                                                        <  2:408121/151‑1 (MQ=255)                                                         aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGGTTTTCATTCCGAAAGAGAGCGACa                                                                                           >  2:557555/1‑148 (MQ=255)                                                               gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGat                                                                                   <  1:755393/150‑2 (MQ=255)                                                                 cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGGAAGATGCGCCGGAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                     >  1:601250/1‑146 (MQ=255)                                                                 cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTTCGAAAAAGAGCGACAGTCTgg                                                                                     >  2:159462/1‑146 (MQ=255)                                                                      aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTg                                                                            <  1:198221/150‑1 (MQ=255)                                                                      aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGa                                                                           <  1:522575/151‑1 (MQ=255)                                                                        cGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGTCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAt                                                                          <  1:310103/150‑1 (MQ=255)                                                                               ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa                                                                  <  2:246337/151‑1 (MQ=255)                                                                                gACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGGGAAAGAAGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGAAGATCGGGATGTGGCGCa                                                                  >  2:93874/1‑150 (MQ=255)                                                                                  cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                               <  2:487087/151‑1 (MQ=255)                                                                                   gcgcgAGATTTAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                              <  1:367002/151‑1 (MQ=255)                                                                                       gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                                          <  1:108081/151‑1 (MQ=255)                                                                                       gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGG‑TGGGGTGGAAAATGCGCTGAAAGGGCAGGTGGGGGGGGGGGTTGTCATTCCGGAAAAGAGGGACAGGCTGGGAGATCGGGATGTGGCGCCCTCGcca                                                           >  1:627902/1‑148 (MQ=255)                                                                                       gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGCGGCGGCGGG‑TGGGGCGAAACAACCGCCCACAGGGCCCGGGGCGGCGGCCCTTGTCCTTCCCAACCAGAGCGCCAGACCGGAAGATCCTGCTGCgggcccccccc                                                             >  1:616143/1‑139 (MQ=255)                                                                                             gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAg                                                     <  2:524628/150‑1 (MQ=255)                                                                                              aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGc                                                   <  1:677680/151‑1 (MQ=255)                                                                                               agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCg                                                  <  1:54013/151‑1 (MQ=255)                                                                                                  gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAtt                                               <  2:499518/151‑1 (MQ=255)                                                                                                          ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgc                                          <  2:774579/148‑1 (MQ=255)                                                                                                                 tCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGGGGTTGTTATTCCGAAAGAGAGCGACAGTCTTGGAGATCGTGGTGTGGGGCACTCGCAAGAGAAGGGGATTATGGGGCtggggg                                  >  2:413982/1‑149 (MQ=255)                                                                                                                                       gtggt‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGGCATTCCGAAAGAGAGCGACAGTCTGGAAGAGCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc          <  1:365828/151‑1 (MQ=255)                                                                                                                                             tGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCg     >  1:726978/1‑151 (MQ=255)                                                                                                                                                ggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcg  <  1:709451/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGT‑TGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354144‑354436

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