Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,318	0	T	G	53.8%	19.0 / 33.1	26	V538G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base G (14/0);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 1.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-03

ATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCACGTCTGGTTTGTCTCGCAATTG  >  NC_000913/354173‑354460                                                                                                                                                  |                                                                                                                                               aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCCGGAAGGGGcgggggcgggggggg                                                                                                                                           >  2:455455/1‑148 (MQ=255)          cgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCCGAAAGGGCAgggggggggggggggtgtgatt                                                                                                                                   >  2:441626/1‑139 (MQ=255)           gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAGGGCCGGGGGCGGGGGCGGTTGTCATTcc                                                                                                                                 >  1:632915/1‑151 (MQ=255)                gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGAAGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTTCGaaa                                                                                                                             >  1:361197/1‑150 (MQ=255)                  cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGAGTGAAGAGAGTATCTCCAGCCTTCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGACAGATGCGCTGAAATGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                                                                                                          <  2:591327/151‑1 (MQ=255)                       tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGGGCTGCAAGGGCAGGGGGCGGGGGCGTTTGGGATTTCGAAaaaga                                                                                                                         >  2:149732/1‑147 (MQ=255)                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCTGGAAGGGCAGGGGGGGGGGGGGGTTGTCATTCCGAAAGAGAGAGACAGTCTgg                                                                                                             >  1:705514/1‑148 (MQ=255)                                  gCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGGAAGGGCAGGTGGCGGGGGCGCTTGTCATTTCCAAAGAGAGCGagccgctgg                                                                                                             >  1:275939/1‑140 (MQ=255)                                    cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGAt                                                                                                        <  2:693179/151‑1 (MQ=255)                                      ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCg                                                                                                      <  2:90426/151‑1 (MQ=255)                                               tGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGAGACAGGCTGGAAGAGCGGGATGTgg                                                                                              >  2:835028/1‑150 (MQ=255)                                                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgcg                                                                                            <  1:644434/151‑1 (MQ=255)                                                ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgcg                                                                                            <  2:632771/151‑1 (MQ=255)                                                      gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAAATGCGCTGAAAGGGCAGGGGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGGCAGTCTGGAAGATCGTGATGTGGCGCACTc                                                                                       >  1:780115/1‑150 (MQ=255)                                                       cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                                                     <  1:198360/151‑1 (MQ=255)                                                           agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGAGCTGGAAGGGCAGGGGGCGGGGGGGGTTGTTATTTCGAAAGAAGGCGACAGTCTGGGAGATAGGGGTGTGGGGCACTCGCaa                                                                                   >  2:560170/1‑149 (MQ=255)                                                              ttGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGACAGGGCAGGGGGGGGGGGGGTTTGGCGTTCCGAAAGGGAGGGCCAGGCTGGAAGATAGTGGTGTGGGGCACTCGCCag                                                                                  >  2:26192/1‑147 (MQ=255)                                                                                      cGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGGAAGGGCAGGGGGGGGGGGGGGTTGTTATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGGTGGGGCGCCCTCGCAAGAGAAGGCGGTTATGGCGCtggggg                                                          >  2:836664/1‑147 (MQ=255)                                                                                       gggCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGGGGGCa                                                       >  1:493358/1‑149 (MQ=255)                                                                                          cGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGAAGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCCTTTCGAAAGAGAGCGACAGGCTGGAAGATCGGGATGTGGCGGACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGAAAGcc                                                    >  2:743285/1‑149 (MQ=255)                                                                                               acGAAGTGGCGGGGGGTGGGGGGAAAGATGCGCCAAAAGTGCAGGGGGGGGGGGCGGGTGTTTTTTCGCAAAAGAGGGACAGTCTGGAAGATCGTGATGTGGCGCCCTCGCAAGAGGAGGGGGATTTGGTGGTGGGGGGCCGCCAGATTgg                                             >  2:49071/2‑151 (MQ=255)                                                                                                  aaGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGGGGCGGTGGCGGTTGTCATTTCGAAAGAGAGCGACAGGCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGAACGCCAGATTGGCa                                           >  2:297492/1‑150 (MQ=255)                                                                                                     tGGCGGAGGTTGGGGAGAAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCCAGCGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACtt                                       >  2:538976/1‑151 (MQ=255)                                                                                                          gtggtTGGGGTGAAAGAGGGGCTGAAAGGGAAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc                                  <  2:175555/151‑1 (MQ=255)                                                                                                                ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCgg                            <  1:248525/151‑1 (MQ=255)                                                                                                                ggggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCg                             <  2:840617/150‑1 (MQ=255)                                                                                                                                 gAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCACGTCTGGTTTg             <  2:179229/149‑1 (MQ=255)                                                                                                                                   aaGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCACGTCTGGTTTGtctc         >  1:455135/1‑151 (MQ=255)                                                                                                                                          ggtggcggtggcGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCACGTCTGGTTTGTCTCGCAATTg  <  2:248462/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               ATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCACGTCTGGTTTGTCTCGCAATTG  >  NC_000913/354173‑354460

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: ‑