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

GGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGAT  >  NC_000913/3572‑3865                                                                                                                                                     |                                                                                                                                                    ggTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTAc                                                                                                                                                   <  2:460094/151‑1 (MQ=255)       tccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGACTgg                                                                                                                                                >  1:650340/1‑148 (MQ=255)        ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaaa                                                                                                                                             <  2:106372/150‑1 (MQ=255)          ggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGACTGGgaaa                                                                                                                                            >  1:489600/1‑149 (MQ=255)                ggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGaaaaaaaaa                                                                                                                                       >  1:338246/1‑143 (MQ=255)                           gTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt                                                                                                                         <  1:810273/151‑1 (MQ=255)                            ttCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAt                                                                                                                         <  2:56722/150‑1 (MQ=255)                              cGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCt                                                                                                                       <  2:152914/150‑1 (MQ=255)                                  ctGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGaa                                                                                                                    <  2:851643/149‑1 (MQ=255)                                    gtgtgACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGAt                                                                                                                <  2:389195/151‑1 (MQ=255)                                       tgACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCGCGAGGACAGGGAAAATAAAAAAAAATCTACAAtaataaa                                                                                                                   >  2:526249/1‑139 (MQ=255)                                       tgACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGAATGGGAAAAAAAATGAAAAACTACAAACTTaaaaaaca                                                                                                              >  1:875194/1‑145 (MQ=255)                                              ccGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGAATGGGAAAATAAAAAAAAACCTACAAAATGAAAAATCACAAc                                                                                                          >  2:346180/1‑147 (MQ=255)                                                ggAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGAGGGGCCAGGACTTGGAAAATAAAATAAACTCCACAAACTGGAAAATCACaa                                                                                                           >  2:202992/1‑144 (MQ=255)                                                 gAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGGGCGAGGAATGGGAAAATAAAAGAAAATCTACAAAATGAAAGAAAACAACGaggagg                                                                                                   >  2:803973/1‑151 (MQ=255)                                                    aCCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCa                                                                                                <  2:278724/151‑1 (MQ=255)                                                        cccAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCtt                                                                                            <  2:557032/151‑1 (MQ=255)                                                            gcgcGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTgcg                                                                                        <  1:458200/151‑1 (MQ=255)                                                                gTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGACGAGGACTGGGAAAATAAATGAAAATCTACAATCTGAAAGATCACAACAAGGAGGGCAGCtttt                                                                                          >  1:182336/1‑144 (MQ=255)                                                                      gACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑aa                                                                            <  1:198025/151‑1 (MQ=255)                                                                            ttGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACCAGGAATGGGAAAATAAAAAAAAATCTTCAATCTTAAAGAACAAAACGAGCAGGGGAGCCTTTCGCCAG‑CCGT‑AAACCCgg                                                                      >  1:767561/1‑151 (MQ=255)                                                                              gggTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAggg                                                                     <  2:752241/150‑1 (MQ=255)                                                                              gggTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACCGGGACTTGAAAAAAAAATGAAAAACTACAATCTTAAAAAACCCAACGAGCAGGGCCGCTTTTCGCAAA‑CCCC‑AAcccccggg                                                                    >  2:650858/1‑146 (MQ=255)                                                                                gTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGtt                                                                  <  2:658044/151‑1 (MQ=255)                                                                                  aaGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTgg                                                                <  2:144855/151‑1 (MQ=255)                                                                                   aGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGAGGACGAGGAACGGGAAAATAAAAAAAACAATTCAAACTGAAAAAAACAAACGAGCAGGTCAGCATTTCGCAAG‑CCCTAAACCCGGGGGTTgg                                                                >  2:483624/1‑151 (MQ=255)                                                                                    gAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGGGCTGGGAAAAAAAATGAAAAACTACAATCTGAAAGATCACAAACAGCAGGGCAGCTTTGCGCAAGCCCGA‑AACCAAGGGGTgggg                                                               >  2:119093/1‑151 (MQ=255)                                                                                       cTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAACTAAAAGAAAATCTACAATCTGAAAAATCACAAAGAACAGGGCAGCTTTTCGCAAG‑CCGT‑AAACCAGGGGGTggggca                                                            >  1:875531/1‑147 (MQ=255)                                                                                        tACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGACGAGGACTGGGAAAATAAAAAAAAATCTACAAACTGAAAGAACACAAACAGCAGGGCAGCCTTTCGCCAG‑CCGC‑AACCCCGGGGGTGGGGaaa                                                           >  2:247271/1‑150 (MQ=255)                                                                                        tACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGAAAAATAAATGAAAATCTACAATCTGAAAGATCACAACGAGCAGGGCAGCTTTTCGCAAG‑CCCT‑AAACCAGGGGGTGGGGaaa                                                           >  1:167762/1‑150 (MQ=255)                                                                                         aCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCaaa                                                           <  1:75802/149‑1 (MQ=255)                                                                                         aCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCaaa                                                           <  2:249762/149‑1 (MQ=255)                                                                                              cAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAgg                                                    <  2:456736/151‑1 (MQ=255)                                                                                              cAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAACTAAAAGAAAATCTACAATCTGAAAGATCACAAAGAGCAGGTCAGCTTTGCGCAAA‑CCGT‑AAACCAGGGGGTGGGGAAAAAtaa                                                      >  1:408314/1‑147 (MQ=255)                                                                                               aaaaTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCCGGTACTGGAAAAATAAAAGAAAATCCACAATCTTAATGATCACAACGAGATGGTGAGCTTTGTGCCAC‑GCCT‑AATCCATGGGGTTGGCAAAAATCAagg                                                   >  1:510574/1‑151 (MQ=255)                                                                                                    cAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGAATGGAAAAATAAATGAAAATCTACAATCTGAAAGATCACAAACAGCAGGTCAGCTTTTCGCAAG‑CCGT‑AACCCCGGGGGTGGGGAAAAAACAgggg                                                  >  1:791555/1‑147 (MQ=255)                                                                                                      ggAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGttt                                            <  1:871840/151‑1 (MQ=255)                                                                                                              ttGTTCATATTTGCCGGCTGGATACGGCGGGCGCACCAGGACTGGAAAACTAAAAGAAACTCTACAATCTGAAAGATCACAACGAGCAGGGCAGCTTTGTGCAAG‑CCGT‑AACCCAGGGGGTGGGCGAAAAAAAGGGGCTGGTTTTTTCGc                                     >  1:137905/1‑150 (MQ=255)                                                                                                                     tatTTGCCGGCTGGATACGGCGGGGGCCAGAGGCATTGGAAACAAAATGAAAATCTACAATATGAAAGAATACCACGAGGAGGGCAGACTTGCGCAAGACCCA‑AACCCAGGGG‑GGGGGAAAAATAAAAGGGTGGTTTTTTTGCTCGACCt                              >  2:535672/1‑150 (MQ=255)                                                                                                                      atTTGCCGGCTGGATACGGCGGGCGCACGCGGACTGGGAAAAAAAATAAAACCCTACAAACTGAAAGAACACAACGAGCAGGGCAGCTTTTCGCAAG‑CCCT‑AACCCAGGGGGTGGGCAAAAAAAAGGGGGGGtttttttcccccccccccc                            >  2:85873/1‑138 (MQ=255)                                                                                                                         tGCCGGCTGGATACGGCGGGCGCCCGAGGACTTGGAAAAAAAATTAAACTCTACCAACTTAAAGATCCCAACGCGCAGGGCAGCTTTGCGCAAG‑CCCT‑CACCCCCGGGGTGGGCCAAAATAAGGGGGGGgttttttctccccccccgcc                           >  2:564872/1‑135 (MQ=255)                                                                                                                                     aCGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctg              <  1:162171/150‑1 (MQ=255)                                                                                                                                         cgggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgactg          >  1:494831/1‑150 (MQ=255)                                                                                                                                          gggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgactga         >  1:402768/1‑150 (MQ=255)                                                                                                                                           ggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGaa        >  1:691603/1‑150 (MQ=255)                                                                                                                                               cACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtg    >  1:264727/1‑150 (MQ=255)                                                                                                                                                aCGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgat  >  2:602661/1‑151 (MQ=255)                                                                                                                                                aCGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgat  <  2:523262/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAG‑CCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGAT  >  NC_000913/3572‑3865

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