Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,720	0	T	G	55.6%	27.8 / 41.9	28	V307G (GTA→GGA)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/12);  new base G (15/0);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 5.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-06

ATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTG  >  NC_000913/3575‑3857                                                                                                                                                  |                                                                                                                                           aTCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTg                                                                                                                                        <  1:576975/150‑1 (MQ=255)     ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACGAGGACTGGaaa                                                                                                                                    >  1:604847/1‑150 (MQ=255)          tccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGgaaaaaaa                                                                                                                               >  2:182431/1‑146 (MQ=255)              gcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGACTGGgaaaaaaaa                                                                                                                              >  2:516043/1‑142 (MQ=255)              gcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGaa                                                                                                                          <  2:400727/151‑1 (MQ=255)                 cGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCGCGAGGACTGGgaaaaaaaaaaaaa                                                                                                                         >  2:326406/1‑139 (MQ=255)                    ccTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTa                                                                                                                    <  2:146052/151‑1 (MQ=255)                                         aGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACGAGGAATGGGAAAATAAATAAAAATCTCAAATTTTAAAAATCACaa                                                                                                  >  1:349095/1‑148 (MQ=255)                                                  ccGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACAGGAAAAATAAAAAAAACACCACAATCTTAAAAAACACaaacaacagg                                                                                          >  2:104829/1‑139 (MQ=255)                                                        aGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTgc                                                                                <  1:564206/151‑1 (MQ=255)                                                          cgcgTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAACTAAAAAAAAAACTACAATCTGAAAAAAAACAAAGAGCAGGTCAGCGTTgcg                                                                               >  2:35458/1‑150 (MQ=255)                                                          cgcgTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGAAAACTAAAAGAAAATCTACAATCTGAAAAAACACAACAAGCAGGGCAGCTTTtcg                                                                               >  1:350162/1‑150 (MQ=255)                                                            cgTTGCTGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACCAGGACTGGAAAAAAAAATAAAAATCTTCAAACTGAAAAATAAAAAAAAGCAggggagattt                                                                                  >  1:82946/1‑137 (MQ=255)                                                                   gACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGAACTGGGAAAAAAAATAAAAATTTACAAACTGAAAGAACACAAAGAGCAGGGCAGCTTTGTGCAAACCgga                                                                      >  2:236711/1‑148 (MQ=255)                                                                                aGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTggg                                                       <  1:132686/151‑1 (MQ=255)                                                                                  aaCTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGc                                                      <  2:548711/150‑1 (MQ=255)                                                                                      aCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCCAGGTATTGAAAAATAAAAGAAAAACTTCCACCTGAAAATTTACAAAAAGGAGGAAAGGGTTTCGCAAGCCGGAACACCAGGGGTGGGGGaa                                                    >  2:338308/1‑149 (MQ=255)                                                                                          gCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAg                                             <  1:326413/151‑1 (MQ=255)                                                                                               aTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCGCGAGGAATGGGAAAAAAAAAAAAAATCTACAAAATTAAAAATCACAACGAACAAGGCAACTTTTCGGAAGCCCCAAA‑CCACGGGGGGGGGGAAAAAAAgggg                                          >  2:220191/1‑149 (MQ=255)                                                                                                   ggAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGttt                                    <  1:300883/151‑1 (MQ=255)                                                                                                   ggAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGGCAGCTTTTCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGttt                                    >  1:275709/1‑151 (MQ=255)                                                                                                    gAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGGACTGGGAAAAAAAATGAAAATCTAAAAAATGAAAAAACACAACGAGCAGGGCAGCTTTGCGCAAACCCTAAC‑CCAGGGGGTGGGCAAAAAACAGGGGGTGtttt                                   >  1:300318/1‑151 (MQ=255)                                                                                                    gAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGGAATTGGAAAAAAAATTAAAATCTACAAAATTAAAAAAAACAACGAGCAGGGGAACTTTTCGCAAACCCGAAC‑CCACGGGGTGGGGaaaaaaaaaggggggtttt                                   >  2:92988/1‑137 (MQ=255)                                                                                                         ttttGTTCATATTTGCCGGCTGGATACGGCGGGGGCCCGAGGACTTGGAAAAAAAAAAAAAATCTACAATCTGAAAAAAAACAACAAGCAGGGCAGCTTTTCGCAAAACCGAAA‑CCAGGGGGTGGGGAAAAAACAGGGGGGGttttttt                                >  2:632825/1‑148 (MQ=255)                                                                                                                 atatTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCt                      <  1:625885/151‑1 (MQ=255)                                                                                                                  tatTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCaa                                                    <  1:350962/120‑1 (MQ=255)                                                                                                                          ggCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAAtt             <  1:105190/151‑1 (MQ=255)                                                                                                                                     gcgggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgactg  <  2:145957/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                           ATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAAC‑CCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTG  >  NC_000913/3575‑3857

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