Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,740	0	C	A	42.4%	53.5 / 10.6	33	L3I (CTC→ATC)	thrC	L‑threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (11/8);  new base A (14/0);  total (25/8)
Fisher's exact test for biased strand distribution p-value = 1.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGAT  >  NC_000913/3598‑3883                                                                                                                                               |                                                                                                                                                 gTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACaa                                                                                                                                           <  2:734458/150‑1 (MQ=255)    cGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCt                                                                                                                                        <  2:194615/150‑1 (MQ=255)              gACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATcac                                                                                                                              <  2:803878/150‑1 (MQ=255)              gACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATcac                                                                                                                              <  2:649523/150‑1 (MQ=255)                    ccGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGCGGAATGGGAAAAAAAATGAAAATCTACAAACTTAAAGATCaaaa                                                                                                                            >  1:315337/1‑146 (MQ=255)                          aCCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCCCGGGGACTGGGAAAATAAAAGAAAATCTACAATCTGAAAGAACACaa                                                                                                                            >  1:655747/1‑140 (MQ=255)                              cccAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGCGGACTGGAAAAATAAAAGAAAATCTACAAACTGAAAGAAAACAACGAGCAGGGGAGCtt                                                                                                             >  2:243747/1‑151 (MQ=255)                               ccAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACGAGGAATGGGAAAATAAATGAAAATCTACCATCTGAAAGATCACAACCAGCCGGGCAGCgtt                                                                                                            >  1:58909/1‑151 (MQ=255)                               ccAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCCAGGACCGGGAAAAAAAATTAAAAACTACAATATGAAAAAACACAAAAAGCAGGGCAGCttt                                                                                                            >  2:35685/1‑151 (MQ=255)                                 aGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAAATAAATGAAAATCTACAAACTGAAAGATCACAAACAGCAGGGCAGCtttt                                                                                                           >  1:170130/1‑149 (MQ=255)                                     cgTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCaa                                                                                                      <  2:413674/151‑1 (MQ=255)                                             aCTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑c                                                                                             <  1:677157/151‑1 (MQ=255)                                               tGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGCGGACTGGGAAAATAAATGAAAATCTACAAACTTGAAGATAACAAAGAGCAGGGCAGCTTTTCGCAAACCGCAA‑ccc                                                                                           >  2:572303/1‑151 (MQ=255)                                                           aaCTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGCGACCTGGAAAACTAAATGAAAATCTAAAAAATCGAAGATCACAAACAGCCGGGCAGCTTTTCGCCAGCCGTAA‑ACCAGGGGGTg                                                                                   >  1:331223/1‑147 (MQ=255)                                                               aCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGTGCACGAGGAATGGGAAAAAAAAAAAAAAACTACAAAATTAAAAATAACCACAAGGAGGGCAGACTTTTGCAAGCCCGAA‑ACCACGGGGTggggacaa                                                                            >  2:579297/1‑145 (MQ=255)                                                                  tGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGGACTGGGGAACAAAAAAAAAATATCCAATTTGAAAAATAACAACAAGAAGGGCAGACTTTTGCCAGACGTAA‑CCCCGGGGGTggggcca                                                                             >  2:104141/1‑142 (MQ=255)                                                                    cAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGTAAAATAAATTAAAATCTACAAACTTAAAGAACACAACAAGCAGGGCCGCTTTTTGCAAACCGCAACCCCAGGGGTGGGGGaaaa                                                                            >  2:545719/1‑146 (MQ=255)                                                                     aaaaTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAAATAAAAGAAACTCTACAAACTGAAAGAAAACAACGAGCACGGCAGCGTTTTGCAAGCCCTAA‑ACCAGGGGGTGGGGAAAAAAAAggg                                                                     >  1:384259/1‑151 (MQ=255)                                                                           aGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACGAGGACTTGAGAAATAAATGAAACTCTCCAATCTGAAAGATCACAACGAGCCGGTCAGCCTTTCGCCAGCCGGAA‑ACCAAGGGGTGGGGGAAAAAAAGGGGGTg                                                                 >  2:16435/1‑149 (MQ=255)                                                                           aGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGtt                                                               <  2:230164/151‑1 (MQ=255)                                                                           aGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGGATGGGAAAAAAAATGAAAATCTACAATCTGAAAGAACACAAACAGCAGGGCAGCTTTGTGCAAACCGCAA‑CCCAGGGGGTGGGGGAAAAAAAGGGGGTg                                                                 >  2:263413/1‑149 (MQ=255)                                                                              aaGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGGAATGGGAAAAAAAATAAAAAACTACAAACTGAAAGAAAACAACGAGCAGGGCAGCTTTTCGCAAGACCTAA‑ACCCCGGGGTGGGCAAAAAAAAAGGGGTGgtttt                                                            >  2:84420/1‑151 (MQ=255)                                                                                   tttGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTTGGAAAATAAATAAAAATCTACAAACTGAAAAAAAACAAAAAGCAGGGCAGCATTTTGCAAACCGGAA‑CCCCAGGGGTGGGGAAAAAAAAGGGGGTGttttttt                                                          >  2:615083/1‑147 (MQ=255)                                                                                    ttGTTCATATTTGCCGGCTGGATACGGCGGGGGCACGAGAACTGGGAAACTAAAATAAACTCCACAAACTTAAAAAACAAAAAAAGCAGGTCAGCTTTTTGCAAACCGTAA‑CCCAAGGGGTGGGGAAAAAACAGGGGGTGTTTTTTCTGc                                                       >  1:525602/1‑150 (MQ=255)                                                                                                   ggCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAAtt                                       <  1:716462/151‑1 (MQ=255)                                                                                                           aCGGCGGGCGCACGAGTACTGGAAAAATAAATGAAACTCTACAATCTGAAAGATCACAACGAGCCGGGCAGCTTTGCGCAAACCGTAA‑CCCAGGGGGTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGGATTCAGcc                                  >  1:115302/1‑148 (MQ=255)                                                                                                                     cACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtg                      >  2:695544/1‑150 (MQ=255)                                                                                                                      aCGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGGAATtgat                    >  1:640755/1‑151 (MQ=255)                                                                                                                       cGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgatg                   >  2:83203/1‑151 (MQ=255)                                                                                                                         aGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgatg                   >  1:745673/1‑149 (MQ=255)                                                                                                                          gTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgatg                   >  1:328125/1‑148 (MQ=255)                                                                                                                                gAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTg           >  1:197950/1‑150 (MQ=255)                                                                                                                                        aaTGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGAt  >  2:462910/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                 GTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAA‑CCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGAT  >  NC_000913/3598‑3883

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