Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	51.2%	28.6 / 54.5	41	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (5/15);  new base G (21/0);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 2.44e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-06

GCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354144‑354434                                                                                                                                                 |                                                                                                                                                     gCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAg                                                                                                                                                 <  2:1103785/150‑1 (MQ=255)    tCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGaagg                                                                                                                                             >  1:1093218/1‑148 (MQ=255)                      aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGA‑AAGGGCAgggggg                                                                                                                          >  1:833880/1‑150 (MQ=255)                         aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAAAATCGCTGG‑AAGGGCCggggggggg                                                                                                                       >  1:878886/1‑150 (MQ=255)                              aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGCCTTA‑AAGGGGAgggggggggggggg                                                                                                                  >  1:289626/1‑148 (MQ=255)                                   cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAAAGGGCTGG‑AAGGGGAGGGGGGGGGGGGGGTTggt                                                                                                             >  1:918666/1‑149 (MQ=255)                                              aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAag                                                                                                   <  2:1018861/150‑1 (MQ=255)                                               cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGTGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                                                                                 <  2:773296/151‑1 (MQ=255)                                                   gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGgtggcggtggcg                                                                                                                   <  1:395522/129‑1 (MQ=255)                                                   gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGGAAGATGCGCTGG‑AAAGGGAGGGGGGGGGGGGGGTTGTCATTTCGGAAGAGAGCg                                                                                             >  1:365722/1‑151 (MQ=255)                                                     aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGGGGGGGTTGTCATTTCGAAAGAGAGGGa                                                                                            >  1:987816/1‑150 (MQ=255)                                                     aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGA‑AAAGGGGGGGGGGGGGGGGGGTTGTCATTTCGAAagagagggg                                                                                            >  1:102337/1‑147 (MQ=255)                                                         aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTc                                                                                       <  2:92780/151‑1 (MQ=255)                                                             ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑GAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGa                                                                                   <  2:587743/151‑1 (MQ=255)                                                                ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                 <  1:642581/150‑1 (MQ=255)                                                                 cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                               <  2:86866/151‑1 (MQ=255)                                                                        cGGCTGGGGACGCGCGAGCTTGAAGAGAGTATCTCCAGCCATCCGGGCGCTGCCGAAGTGGCGGTGGTTCGGGGGGAAGAGAAGCTGATAAGGGCCGCGGGC‑GCGGCGTTTTGCCTTTCGAAAGAGAGGCACAGCCTCGGAGAACggg                                                                           >  2:197866/1‑146 (MQ=255)                                                                        cGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAtg                                                                        <  1:127415/151‑1 (MQ=255)                                                                         ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGAGGCGCTGA‑AAGGGCAGGGGGCGGGGGCGGTTGTTATTCCGGAAGAGAGGGACAGTCTGGGAGAACGGGaggg                                                                       >  1:238644/1‑148 (MQ=255)                                                                                        agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGATGCGCTGA‑AAGGGCCGGGGGGGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGGCTGGAAGATCGGGGGGTGGGG‑CACTCGCAaga                                                       >  1:743561/1‑151 (MQ=255)                                                                                        agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGA‑AAGGGGAGGTGGCGGGGGGGGTTGTCATTCCGGAAGAGAGCGGCAGGCTGGAAGATCGGGAGGGGGCG‑CACTCGCAaga                                                       >  2:440407/1‑151 (MQ=255)                                                                                        agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGA‑AAGGGGAGGGGGCGGGGGGGTTTTTCATTTCGGAAAAGAAGGACAGACTGGAAGATCGGGGTGTGGGG‑CACTCGCAaga                                                       >  2:1256905/1‑151 (MQ=255)                                                                                        agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGACGCGGTGG‑AAGGGGAGGGGGCGGCGGCGGTTGTTGTTCCGAAAGAAGGCGGCCGCCTGGCAGGTAGGGGTGGGGGG‑CGCTCTCAaga                                                       >  2:298542/1‑151 (MQ=255)                                                                                          aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAAATGCGCTGA‑AAGGGCCGGGGGCGGGGGGGGTTGTCATTCCGGAAGAGAGGGGCAGTCTTGGAGATCGTGGTGGGGGG‑CCCTCCCAAggga                                                     >  2:907296/1‑151 (MQ=255)                                                                                           ttGAAGAGAGTATCTCCAGCCACCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGaa                                                    <  2:816554/151‑1 (MQ=255)                                                                                              aaGAGAGTATCTCCAGTCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGCGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCCTGATGTGGCG‑CACTCGCAAGAGAAGGc                                                 <  2:657258/151‑1 (MQ=255)                                                                                               agagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGA‑AAGGGCAGGGGGGGGGGGGGGTTGTCATTTCGAAAGAGAGTGACAGTCTGGAAGATCGTGATGTGGGG‑CACCCGCAAGaggagggg                                                >  1:522595/1‑149 (MQ=255)                                                                                                  gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAAATGCGCCGG‑AAGGGGAGGGGGGGGGGGGGTTTGTCATTTCGGAAGAGAGCGAAAGGCTGGAAGATCGTGATGTGGGG‑CACCCCCAAGAGAAGGGGAtt                                             >  2:1186826/1‑151 (MQ=255)                                                                                                   agTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGAAGCGCTGA‑AAGGGCCGGTGGGGGGGGGGTTTGTTATTCCGAAAAAGAGCGAAAGGTTGGAAGATAGTGGTGGGGGGCCACTC‑CAAAAGAAAGCGAATa                                            >  1:779841/1‑151 (MQ=255)                                                                                                     tATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTAt                                           <  1:672845/150‑1 (MQ=255)                                                                                                      aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGA‑AAGGGCAGGGGGGGGGGGGGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGGTGTGGCG‑CACTCGCAAGAGAAGGGGGTTATgg                                         >  1:778564/1‑151 (MQ=255)                                                                                                        ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGgcg                                       <  1:635465/151‑1 (MQ=255)                                                                                                          ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGA‑AAAGGCGGGGGGCGGGGGGGTTTGGGAGTTCGAAAGAGAGGGGCAGTCCGGGAGAGCGGGATGTGGCG‑CACTCGCAAGAGAAGGCGGTTATGgcgcg                                     >  2:657526/1‑150 (MQ=255)                                                                                                             gCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGA‑AAGGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGGGGCAGTCTGGAAGATCGGGATGTGGCG‑CACTCGCAAGAGAAGGGGATTATGGGGCtggg                                  >  2:592147/1‑150 (MQ=255)                                                                                                               cATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGGGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                                <  1:635763/151‑1 (MQ=255)                                                                                                                      gCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGcc                          <  1:43406/150‑1 (MQ=255)                                                                                                                         ttGCCGAAGTGGCGGTGGTTGGGGGGAAAAATTCGCTGA‑AAGGGCAGGGGGCGGGGGCGGTTGTCTTTTCGAAAGAGAGGGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGGGGACAGGCAGAt                      >  2:240986/1‑151 (MQ=255)                                                                                                                                        tggtTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg       >  1:987770/1‑151 (MQ=255)                                                                                                                                         ggtTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgccg       >  2:280423/1‑150 (MQ=255)                                                                                                                                              gggTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg  >  1:1120498/1‑150 (MQ=255)                                                                                                                                              gggTGAAAGATGCGCTGA‑AAGGGCAGGGGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGGGGACAGCCAGATTGGCAACTTTGGCCGCCCgg  >  1:230197/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGA‑AAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCG‑CACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354144‑354434

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