Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	54.0%	35.5 / 52.0	38	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (6/11);  new base G (20/0);  total (27/11)
Fisher's exact test for biased strand distribution p-value = 1.45e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-05

CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354145‑354437                                                                                                                                                |                                                                                                                                                      cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                                                <  1:323723/151‑1 (MQ=255)         aCGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTg                                                                                                                                          >  2:539340/1‑149 (MQ=255)           gCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAg                                                                                                                                      <  2:169075/151‑1 (MQ=255)            cTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAgg                                                                                                                                     >  2:332170/1‑151 (MQ=255)              gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAgg                                                                                                                                     <  2:178137/149‑1 (MQ=255)              gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAgg                                                                                                                                     <  2:98991/149‑1 (MQ=255)              gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAAAGGGCCGGAAggg                                                                                                                                    >  1:2105/1‑150 (MQ=255)                    tATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtg                                                                                                                              <  2:258934/150‑1 (MQ=255)                     aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggc                                                                                                                            <  2:523983/151‑1 (MQ=255)                           ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggggg                                                                                                                       >  1:439223/1‑150 (MQ=255)                             aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAAAGCGCTGAAAGGGGAgggggcggggggg                                                                                                                     >  2:533763/1‑148 (MQ=255)                                  cGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTc                                                                                                               <  1:409718/151‑1 (MQ=255)                                   gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGGGCTGAAAGGGCAgggggcggggggggtt                                                                                                                  >  2:465092/1‑142 (MQ=255)                                       gcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCCGGAAGGGCAGGTGGCGGGGGGGGTTGTCAtt                                                                                                            >  1:343434/1‑149 (MQ=255)                                              cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagag                                                                                                   <  2:141192/151‑1 (MQ=255)                                              cGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGTGGCGTTTGTCATTCCGaaa                                                                                                      >  1:514589/1‑148 (MQ=255)                                                  gTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGc                                                                                                <  2:66371/150‑1 (MQ=255)                                                   tAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATGCGCTGAAAGGGGGGGGGGGGGGGGGGTTTGTTATTCCGAAAGAGAGCg                                                                                               >  1:400624/1‑150 (MQ=255)                                                        aaCGTTGCCGGACATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGGAGGGGGCGGGGGGGGTTGTCATTTCGAAAGAGAGCg                                                                                               >  1:548750/1‑145 (MQ=255)                                                          cGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAGGCGCTGAAAGGGGAGGGGGGGGGGGGGGTTGTCATTTCGAAAGAGAGCGACAGTCTg                                                                                       >  1:231559/1‑151 (MQ=255)                                                           gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                      <  1:11990/151‑1 (MQ=255)                                                           gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                      >  2:278070/1‑151 (MQ=255)                                                           gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGTCTgg                                                                                      >  2:133740/1‑151 (MQ=255)                                                                cGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGGAAAGGCAGGGGGGGGGGGCGTTTTTTATTTCGAAAAAGAGCGACAGTCTgg                                                                                      >  2:421179/1‑146 (MQ=255)                                                                        ggCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTg                                                                             >  1:648774/1‑147 (MQ=255)                                                                               gACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGGAAGGGCAGGTGGCGGTGGGGGTTGTCATTTCGAAAGGGAGCGACAGTCTGGAAGATCGTGGTGTGGCGCAc                                                                  >  1:123452/1‑151 (MQ=255)                                                                                aCGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGCGAAAGAGGCGCTGAAAAGGCAGGGGGCGGTGGGGTTTGTGATTCCGAAAGAGAGGGACAGCGTGGAAGATCGGGGTGTGgcg                                                                     >  2:591942/1‑147 (MQ=255)                                                                                    gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGAAGGGCTGAAAGGGCAGGGGGCGGTGGGGTTTGTCCTTCCGAAAGAGAGCGGCAGTCTGGAAGATCGGGATGTGGCGCACTCGc                                                              >  2:224235/1‑150 (MQ=255)                                                                                      gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCCGAAAGGGGAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCAccccc                                                              >  1:539700/1‑144 (MQ=255)                                                                                                       ctcCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGTTGGGGGGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGGGGTTGTCATTCCGGAAGAGAGGGACAGTCTGGAAGATCGTGGTGGGGCGCACTCGCAAGAGAAGGCGATTATTgc                                           >  1:233067/1‑150 (MQ=255)                                                                                                        tcCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGGAGGGGGCGGGGGCGTTTGTCCTTCCCAAAGAGAGGGACAGGCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGGGATTTTTgcgc                                         >  1:523087/1‑151 (MQ=255)                                                                                                             ccATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGGTGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtgggg                                    >  2:39298/1‑149 (MQ=255)                                                                                                             ccATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGAAGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGGGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtg                                    >  2:222608/1‑151 (MQ=255)                                                                                                              cATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGTGGAAGATGCGCTGAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGGGCACTCGCAAGAGAAGGCGATTATGGCGCtggggg                                   >  1:507582/1‑151 (MQ=255)                                                                                                                       gTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGcc                             >  2:648931/1‑148 (MQ=255)                                                                                                                                   gCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg              <  2:208822/151‑1 (MQ=255)                                                                                                                                     ggtggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgc            >  2:623787/1‑151 (MQ=255)                                                                                                                                      gtggtTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc           >  2:561636/1‑151 (MQ=255)                                                                                                                                                tGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  >  2:584191/1‑150 (MQ=255)                                                                                                                                                tGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcgc  >  1:429539/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCGC  >  NC_000913/354145‑354437

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