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

CGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTT  >  NC_000913/354149‑354423                                                                                                                                            |                                                                                                                                        cGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgcgc                                                                                                                              <  2:304753/151‑1 (MQ=255)    gACGCTGACGGTTATCACTTTATTCTCGGGCGCGCTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTGGGGGGGGAAGAGGGGCTgg                                                                                                                           >  1:378822/1‑150 (MQ=255)       gCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGATgggccg                                                                                                                            >  1:369503/1‑145 (MQ=255)         tGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGGGGAAGATGGGCCGAAAggg                                                                                                                      >  2:149386/1‑151 (MQ=255)          gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCAGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAAATGGGCt                                                                                                                             >  1:318520/1‑143 (MQ=255)            cGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGTGCTGGAAggg                                                                                                                      >  2:157945/1‑148 (MQ=255)               ttATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAGGAGCGGGGAAAggggcgggg                                                                                                                >  2:210294/1‑145 (MQ=255)                       ttATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAggtggcggtgg                                                                                                         <  2:463689/150‑1 (MQ=255)                         aTTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGAAAGGGGAggggggggggggg                                                                                                       >  2:43616/1‑148 (MQ=255)                           tctcGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGGGAGAGGCGCTGGAAGGGGGGGGGGGGGGGGCGctt                                                                                                    >  1:77526/1‑151 (MQ=255)                               gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCa                                                                                                <  1:100951/151‑1 (MQ=255)                                    cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCg                                                                                           <  2:800046/151‑1 (MQ=255)                                    cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGACAGGGCAGGTGGGGGGGGGGTTTGTGTTTcc                                                                                            >  1:395515/1‑150 (MQ=255)                                           gATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAagaga                                                                                    <  1:68409/151‑1 (MQ=255)                                            aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAAGGCAGGGGGCGGGGGCGGTTGTCATTTCGAAagag                                                                                     >  2:360769/1‑149 (MQ=255)                                            aTGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGGGGGGGGAAAAATGCGCTTGAAGGGGAGGGGGGGGGGGGGGTTGTCAATCCGgaa                                                                                        >  1:733610/1‑146 (MQ=255)                                                 aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGGGGTGGCGGTTGTCCTTTCGaaa                                                                                        >  2:481952/1‑141 (MQ=255)                                                    aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGt                                                                            <  1:661908/150‑1 (MQ=255)                                                       gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGGGCCGGAAGGGGAGGGGGGGGGGGGGGTTGGCATTTCGAAAGAGGGGGACAGGCTgg                                                                        >  1:740483/1‑151 (MQ=255)                                                         tGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTTTCATTCCGAAAGAGAGCGACAGTCTgg                                                                        >  2:758223/1‑149 (MQ=255)                                                               gtaccGGCTGGGGACGCGTGAGATGGAAGAGAGTATCTCCAGCCATCCGTGCGTTGCCGAAGTTGCGGTGGGTGTGGTGAAAGATGCGCTGAAAGTGCAAGTTGAGGTCGCGTTTGTCATTCCGAAAGAGAGCGACAGTCGGGAAGGTCg                                                                 <  1:703957/146‑1 (MQ=255)                                                                              gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                 <  2:315806/151‑1 (MQ=255)                                                                                  gagaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAag                                             <  2:107221/151‑1 (MQ=255)                                                                                     aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAgaga                                          <  2:601532/151‑1 (MQ=255)                                                                                     aTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTTGGGTGAAAGATGCGCTGGAAGGGCAGGGGGGGGGGGCGGTTGGCATTCCGAAAGAGAGGGACAGGCTGGGAGATCGGGGTGTGGCGGACTCGCaa                                              >  1:492627/1‑147 (MQ=255)                                                                                      ttGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTGTAAGGGCAGGGGGGGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGGCTGGAAGATCGGGGTGGGGCGCACTCGCAaga                                            >  1:495997/1‑148 (MQ=255)                                                                                        gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAg                                        <  2:517822/150‑1 (MQ=255)                                                                                           gagagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGa                                    <  2:553739/151‑1 (MQ=255)                                                                                                     ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGgcg                            <  2:407378/149‑1 (MQ=255)                                                                                                      cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCTGAAAGGGCAGGGGGCGGGGGCGGTTTTCATTCCGAAAGAGAGCGACAGGCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGGTTATGGCGCtg                         >  2:495799/1‑151 (MQ=255)                                                                                                          cATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                     <  2:337332/151‑1 (MQ=255)                                                                                                            tCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAAATGGGCTGGAAGGGGAGGGGGCGGGGGCGGTTGTCCTTCCGAAAGAGAGCGACAGTCTGGGAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGGGCcgggg                      >  1:749360/1‑146 (MQ=255)                                                                                                                      gCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCTTAAAGGGCAGGTGGCGGGGGCGGTTGGCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGGGGTGGGGGGCACTCGCAAGAGAAGGGGGTTATGGCGCTGGTGGACAGCCAGAtt          >  1:525929/1‑150 (MQ=255)                                                                                                                            gTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACt   <  1:801857/151‑1 (MQ=255)                                                                                                                             tGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACtt  >  1:190116/1‑151 (MQ=255)                                                                                                                                            |                                                                                                                                        CGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCTGAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTT  >  NC_000913/354149‑354423

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