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

CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCA  >  NC_000913/354145‑354438                                                                                                                                                |                                                                                                                                                       cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                  <  2:160657/150‑1 (MQ=255) cATCCGTGACGCTGACGGTTATCACGTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGa                                                                                                                                                  <  1:70050/150‑1 (MQ=255)                  gTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGGGCCGAAAGGGCgggg                                                                                                                                >  2:363086/1‑150 (MQ=255)                        aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGACGTGGCGGGGGGTGGGGGGGAAGACGGGCCGAAAggggcggggggggg                                                                                                                          >  2:621826/1‑140 (MQ=255)                        aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCCGGAAGGGGAggggggggg                                                                                                                          >  2:426537/1‑150 (MQ=255)                        aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGGTGCCGAAGGGGCGGTGGGTGGGGGGGAAAATGGGCCGGAAGGGGAggggggggg                                                                                                                          >  2:599693/1‑150 (MQ=255)                             aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTTGGGGGGAAGAAGGGCCGGAAGGGGCggggggggtggggg                                                                                                                     >  1:244020/1‑148 (MQ=255)                             aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGGGCTGAAAGGGCAgggggggggggggg                                                                                                                     >  2:71132/1‑148 (MQ=255)                                  cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAGGTGCGCCGAAAGGGGGGGTGGGGGGGGGGGTTGt                                                                                                                 >  2:413323/1‑150 (MQ=255)                                    ggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCa                                                                                                               <  1:84742/150‑1 (MQ=255)                                                aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagagag                                                                                                  <  1:426387/151‑1 (MQ=255)                                                  gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCTGAAAGGGGCGGGGGGGGGGGGGGTTGTCATTTCGAAagagaggg                                                                                                >  2:539980/1‑149 (MQ=255)                                                  gTAATTAACGTTGCCGGGCATCGGCTGGGGACCCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGCGGCGGAGGTGGGGGGGAAAGAAGCGCTGGAAGGGCCGGCGGGGGGGGGGGTTGTCATTTCGATagagggtg                                                                                                >  1:1859/1‑147 (MQ=255)                                                   tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGTGCCGGAAGGGGAGGGGGGGGGGGGGGGTGTGATTTCGACAGAGAGCgg                                                                                               >  1:15373/1‑150 (MQ=255)                                                   tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCTGGAAGGGGAGGGGGCGGTGGCGGGTGTCATTTCCGAAGAGAGGGa                                                                                               >  1:269044/1‑151 (MQ=255)                                                    aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGAc                                                                                              <  2:432231/151‑1 (MQ=255)                                                    aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGGGGCGGGGGGTGGGGGGAAAGAAGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGGCATTTCGAAAGGGGGCgg                                                                                               >  2:353473/1‑149 (MQ=255)                                                        accGTTGCCGGGCATCGGCTGGGGGCGCGCGGGATTGAAGAGAGGATCTCCAGCCATCCGGGGGTTGCCGGAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTc                                                                                          <  2:269129/149‑1 (MQ=255)                                                               ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                    <  2:270823/150‑1 (MQ=255)                                                               ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                    <  2:159859/150‑1 (MQ=255)                                                               ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAg                                                                                    <  1:241713/150‑1 (MQ=255)                                                                  ggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAAAGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTGATTTCGAAAGAGAGGGAAAGGTTGGAAGAACg                                                                                >  2:106539/1‑151 (MQ=255)                                                                      tCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGAt                                                                            <  2:533243/151‑1 (MQ=255)                                                                        ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAGGGGAGGGGGCGGGGGGGGTTGTTGTTCCGGAAGAGAGCGACAGTCTGGGAGATCggggtgg                                                                          >  2:455499/1‑145 (MQ=255)                                                                        ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGAAAAATGGGCTGAAAAGGGGGGGGGCGGGGGGGGTTGGCATTCCGGAAGAGGGCGACAGTCTGGGAGATGGTggggt                                                                          >  2:98549/1‑147 (MQ=255)                                                                                       agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGGAAGGGGAGGTGGCGGGGGCGTTTTTCATTCCGGAAGAGAGCGGCAGTCTGGAAGATCGTGATGTGGGGCACTCGCAaga                                                           >  1:337109/1‑151 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAagag                                                          <  2:109726/150‑1 (MQ=255)                                                                                         aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGAAGCCCTGAAAGGGCCGGGGGCGGGGGCGTTTGTCATTGCGAAAAAGAGCGACAGCCTGGAAAATTGGGATGTGGCGGACTCGCAAgaga                                                         >  2:424960/1‑151 (MQ=255)                                                                                                 gagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGCGCTGGAAGGGGAGGGGGCGGGGGGGGTTGTCATTTCGAAAGAGAGCGACAGTCTGGGAGATCGTGATGTGGCGCACTCGCAAGAGAAGGGGAtt                                                 >  1:188792/1‑151 (MQ=255)                                                                                                    tATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTAt                                               <  2:419808/150‑1 (MQ=255)                                                                                                       ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGTGGGGGGGGCGTTTGTGATTCCGAAAGAGAGGGACAGGCTGGAAGATCGTGAGGTGGGGGACTCGCAAGAGAAGGCGATTATgggg                                           >  2:63488/1‑149 (MQ=255)                                                                                                       ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTCATTTCGAAAGAGAGGGACAGGCTGGGAGAACGGGGTGTGGGGCACCCGCAAGAGAAGGGGATTATgggg                                           >  2:84907/1‑149 (MQ=255)                                                                                                               aTCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGGGGTTGTTATTCCGGAAGAGGGCGACAATCTGGAAGAACGTGATGGGGGGCACTCGCAAGAGGAGGGGGTTATTGGGCcggtggg                                   >  1:534568/1‑150 (MQ=255)                                                                                                                        ttGCCGAAGTGGCGGTGGTTTGGGGGAAAGATGCGCTGAAAGGGCAGGGGGGGGGGGCGTTTGTCATTCCGAAAGAGAGGGACAGGCTTGGAGATGGTGGGGTGGGGCACTCGCAAGAGAAGGGGATTATGGCGCTGGTGGACAGCCAGAt                          >  1:23589/1‑151 (MQ=255)                                                                                                                                  ggCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACttt                >  1:251969/1‑151 (MQ=255)                                                                                                                                              ggTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcg    >  2:351896/1‑151 (MQ=255)                                                                                                                                              ggTGAAAGATGCGCTGAAAGGGCAGGGGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcg    >  2:67357/1‑151 (MQ=255)                                                                                                                                                tGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCa  <  1:605765/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                       CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGCA  >  NC_000913/354145‑354438

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