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

GGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354142‑354434                                                                                                                                                   |                                                                                                                                                      gggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGaaa                                                                                                                                                   <  1:429972/151‑1 (MQ=255)       ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcg                                                                                                                                             <  2:177699/151‑1 (MQ=255)       ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                                                                              <  2:333372/150‑1 (MQ=255)          tGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGGGGGGGGTAAGagg                                                                                                                                               >  1:921516/1‑144 (MQ=255)                      ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTG‑GAAA‑GGCAgggg                                                                                                                            >  1:625238/1‑149 (MQ=255)                       tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAggtg                                                                                                                            <  2:935597/150‑1 (MQ=255)                        aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGTCCTG‑AAAG‑GGCAggggg                                                                                                                           >  2:202188/1‑150 (MQ=255)                          cccTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAggtggcg                                                                                                                         <  1:69716/148‑1 (MQ=255)                          cACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAggtggcgg                                                                                                                        <  2:413398/151‑1 (MQ=255)                          cACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAggtggcg                                                                                                                         <  2:91641/150‑1 (MQ=255)                           aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAggtggcggt                                                                                                                       <  2:731098/151‑1 (MQ=255)                               tATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAggtggcggtggc                                                                                                                    <  2:680100/150‑1 (MQ=255)                                   ctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTG‑AAAG‑GGCAGGGGGGGGGGGGGttt                                                                                                                >  2:621406/1‑150 (MQ=255)                                    tcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGAGGGGGCG‑GAAG‑GGCAgggggcggggggggtt                                                                                                                >  1:907513/1‑144 (MQ=255)                                     cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGttt                                                                                                                <  2:1021818/148‑1 (MQ=255)                                     cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAGATGCGCTG‑AAAG‑GGGAGGTGGGGGGGGGGGTTGt                                                                                                              >  2:977052/1‑150 (MQ=255)                                       ggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCAt                                                                                                           <  1:552676/151‑1 (MQ=255)                                          gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTcc                                                                                                        <  1:175785/151‑1 (MQ=255)                                                aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAaga                                                                                                  <  2:175744/151‑1 (MQ=255)                                                 cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                                                                                 <  1:180863/151‑1 (MQ=255)                                                  gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCACCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTT‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagaga                                                                                                <  1:283983/151‑1 (MQ=255)                                                       aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTTGGGGGCAAGATCCGCTG‑GAAG‑GGGAGGGGGCGGGGGGGGGTTTCCTTCCGAAacagggag                                                                                             >  2:533951/1‑145 (MQ=255)                                                           aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGTTGGGGGGAAAGATGCGCGG‑AAAG‑GGCAGGAGGCGGTGGCGGTTGTCATTCCCGAAAAGA‑GCGACAg                                                                                        >  2:875673/1‑149 (MQ=255)                                                            aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCt                                                                                     <  2:115238/151‑1 (MQ=255)                                                            aCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTG‑AAAG‑GGCAGGGGGCGGGGGGGTTTGTCATTCCGAAAGAGA‑GCg                                                                                            >  2:429442/1‑144 (MQ=255)                                                                tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTG‑AAAG‑GGGAGGGGGGGGGGGGGGTTGTCATTTCGAAAAAGA‑GCGACCGTc                                                                                      >  1:823776/1‑146 (MQ=255)                                                                     ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAAAAATGCGCTG‑AAAG‑GGGCGGGGGGGGGGGGGTTTGTTATTCCGCAAAAGA‑GGGACAGTCTGGgaga                                                                               >  1:969677/1‑148 (MQ=255)                                                                     ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGCGGCGGCCGGTGGGGGGAAAAATGGGCTGAAAAG‑GGCA‑GTGGGGGCGGCGCTTGGTATCTCGCAagaga‑g                                                                                              >  1:968356/1‑133 (MQ=255)                                                                           ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTG‑AAAG‑GGGCGGGGGGGGGGGGGTTTGTCATTTCGCAAGAGA‑GCGGCAGTCTGGAAGATCGTgg                                                                         >  1:794647/1‑147 (MQ=255)                                                                            gCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGGTGGGGGGGAAGATGCGCTG‑GAAA‑GGGAGGTGGGGGGGGGGTTTGTTATTCCGAAAAAGA‑GCGACAGTCTGGGAGATCgtggtg                                                                       >  2:12132/1‑149 (MQ=255)                                                                                       gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCa                                                          <  1:925972/151‑1 (MQ=255)                                                                                        cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCaa                                                         <  2:841317/151‑1 (MQ=255)                                                                                        cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGTGAAAGAAGCGCAG‑AAAG‑GGGAGGTGGCGGCGGGGGTTGTTGTTGCGAAAGAGA‑GCGGCAATCTGGAAGATCGTGGTGTGGCGGACTCGCa                                                          >  1:879166/1‑150 (MQ=255)                                                                                         gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                                        <  2:558543/151‑1 (MQ=255)                                                                                              tGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTG‑AAAG‑GGCAGGCGGCGGGGGCGTTTGTCATTCCGCAAGAGA‑GCGACCGTCTGGAAGACCGTGATGTGGGGCACTCGCCAgaga                                                     >  1:162173/1‑149 (MQ=255)                                                                                                aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGc                                                 <  2:835342/151‑1 (MQ=255)                                                                                                aaGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAAAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGaa                                                    >  2:915145/1‑148 (MQ=255)                                                                                                 agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTG‑GAAG‑GGCAGGTGGCGGTGGGGTTTTTCATTCCGAAAGAGA‑GAGACAGTGTGGGAAAACGTGGTGTGGCGCACCCGCAAGAGaa                                                    >  2:735419/1‑147 (MQ=255)                                                                                                 agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGCGAAAAAAGCCCAG‑GAAGCGGCGGGGTGCGG‑GGCGGTTGTGCTTCCCCAAGAAAGGCTGCAGGCT‑GAAGATTGTGAGGTGGTGCACACGCAagag                                                      >  1:978542/1‑145 (MQ=255)                                                                                                   agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAt                                              <  2:574653/151‑1 (MQ=255)                                                                                                   agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGGGGGGGGGGGGTTTGTCATTTCGAAAGAGA‑GCGGCAGTCCGGGAGATCGTGATGTGGCGCACCCGCAAGAGAAGGCg                                                >  1:506054/1‑149 (MQ=255)                                                                                                   agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGTAAAATGGGCTG‑GAAG‑GGGAGGGGGCGGGGGCGGTTGTCATTTCGAAAGAGA‑GCGACAGTCTTGAAGATCGGGAAGTGGGGGCCTCGCAAGAGaa                                                    >  1:356207/1‑145 (MQ=255)                                                                                                      gTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTG‑AAAG‑GGCCGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAtt                                             >  2:318248/1‑149 (MQ=255)                                                                                                          ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCTG‑AAAG‑GGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGA‑GCGAAAGTCTGGAAGATCGTGATGTGGGGCACTCGCCAGAGAAGGGGATTATgg                                         >  2:408113/1‑149 (MQ=255)                                                                                                            ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcgc                                      <  2:272583/150‑1 (MQ=255)                                                                                                                ccATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtg                                 <  2:351475/151‑1 (MQ=255)                                                                                                                      gggCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTG‑AAAG‑GGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGGGATGTGGCGCACTCGCAAGAGAAGGCGGTTATGGCGCTGGTGGACAg                            >  2:162225/1‑150 (MQ=255)                                                                                                                          gTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAg                        <  1:742319/150‑1 (MQ=255)                                                                                                                                     ggCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt            <  1:155222/151‑1 (MQ=255)                                                                                                                                              tGGGGCGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCACGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGCTTGGCAACTTTGGCCGCCCg   >  1:1030395/1‑151 (MQ=255)                                                                                                                                               ggggTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg  <  1:991923/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTG‑AAAG‑GGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGA‑GCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354142‑354434

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