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

CTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354139‑354434                                                                                                                                                      |                                                                                                                                                   cTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGt                                                                                                                                                    <  2:648676/150‑1 (MQ=255)     ggCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTTGTTGGGGTGAAAg                                                                                                                                               <  2:480705/151‑1 (MQ=255)       cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGAt                                                                                                                                             <  2:554933/151‑1 (MQ=255)        aTCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGa                                                                                                                                              >  1:40244/1‑149 (MQ=255)          ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                                                                           <  1:741891/150‑1 (MQ=255)          ccGTGACGCTGACGGTTATCACGTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                                                                           <  2:170459/150‑1 (MQ=255)           cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAggaggg                                                                                                                                          >  1:609062/1‑145 (MQ=255)                    gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAAGGCGCTTAAAggg                                                                                                                                 >  1:709432/1‑150 (MQ=255)                            tCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggc                                                                                                                         <  1:373515/150‑1 (MQ=255)                            tCACGTTATTCTCGGGCGCACTGACGATGTAATTAACGGTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtgg                                                                                                                          <  2:414375/149‑1 (MQ=255)                              aCTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAgggggggg                                                                                                                       >  1:1361/1‑150 (MQ=255)                                   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                                                                  <  2:632771/151‑1 (MQ=255)                                                                                  ggggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGgcg                                                                  <  1:644434/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)                                                                                                                                                      |                                                                                                                                                   CTGGGGCATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGG  >  NC_000913/354139‑354434

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