Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,960	G→A	A280A (GCC→GCT)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,960	0	G	A	94.1%	87.1 / ‑1.8	34	A280A (GCC→GCT)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (13/19);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

ATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGC  >  NC_000913/4306723‑4307182                                                                                                                                                                                                                                              |                                                                                                                                                                                                                               aTTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCa                                                                                                                                                                                                                   <  1:354523/251‑1 (MQ=255)                tgtgCTGAAGAGGGAAAATTGCAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTGATGTGTGTATTATCTATATTTTACCTGTCGTTGGCACAGCGGGACATTGCGCTTCTGTTCCTACGTGTTACTCATTGGGCTTTGGCACAGAACTGGGGCAAAAACTGTTGATGCGCCTATATTGCTGCACCCTGAGCGCGATTAAAGTCAGCAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGc                                                                                                                                                                                                     <  2:110309/250‑1 (MQ=255)                        agagGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACgg                                                                                                                                                                                            >  1:141155/1‑251 (MQ=255)                                   tACAGTAACACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATCTATTGTACCTGTCGTTGGCAAAGCGGGACATTGAGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGATGATGAGACAATATTGCTGCACACTGAGCGACATTAAAGTAAGAGGAAGGGGCATCTATAAAGAGCACGACCTGACTCGGCAGTGGACgggacgggaaca                                                                                                                                                                                  >  2:338374/1‑240 (MQ=255)                                         aaaaCTATTATTGCAAAAAATAAGCTTAGACGTTTGTTGTGTTTATTATATATATTGTTCCTGTCGTTGGCACAGTGGGTCATTGCGCTACTGTGCCAACGTGTTACTCATGGGGCTTTATCACAGATCTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTAAGCGCCATTAAAGTCAGAAGATGATGCAGCAATAAAGCTCACGACCTGATGCGGCAGTGGTCGCCGCAGGTCCTTTTggtt                                                                                                                                                                           <  2:29118/251‑3 (MQ=255)                                                  aTTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATTTATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAAc                                                                                                                                                                  <  1:234136/251‑1 (MQ=255)                                                         aaaaTAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGAACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAACAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTACAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACATGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGtc                                                                                                                                                           >  1:29796/1‑251 (MQ=255)                                                                                tttttagtctctatatTGTATCTTACGTTGGCACAGCGGGACATTGCTCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGCACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGCTGCGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCCGGTACTTTTTGGTCATGGCAACGAGAGTCTCGAGTGGGAAGGAGGGCCTATc                                                                                                                                    >  1:351235/12‑251 (MQ=255)                                                                                    tattatATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGAGGGAAGGCGGGCATATCCATc                                                                                                                                >  1:231311/1‑251 (MQ=255)                                                                                       tatatataTTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCAc                                                                                                                              <  1:365278/250‑1 (MQ=255)                                                                                               ttGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATACCCATCACGCCACCGcc                                                                                                                     >  1:304367/1‑251 (MQ=255)                                                                                                         cGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAcc                                                                                                           <  1:310370/251‑1 (MQ=255)                                                                                                             ggcacggagggACATTGCGCTACTGTGCCAATGTTTTACTCATGTGGTTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCACCAATATTGCTGCACCCTGATCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGGCCTGATGCGTCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCGCAGGATCACCGCa                                                                                                        <  2:92010/242‑1 (MQ=255)                                                                                                              gCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACAGCCCAGGACCAACGCATc                                                                                                      >  1:257750/1‑251 (MQ=255)                                                                                                                    cGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAACCGTTGCTGCGCCACTATTGCTGCACCCTGATCTCCATTAAAGTCAGAAGATGAGGCAGCACTAAAGCGCACGACCTGCTGCTGCCGTTGACGGCCCAGGTACTTTTGGGTCATGGCCAAGCGCGTCTCGCGTTGGCACGCGGGCCTCTCCATCACCCCACCTCCCaagctcggcagagcgtcgtg                                                                                                >  2:35980/1‑232 (MQ=255)                                                                                                                    cGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCa                                                                                                                   <  1:35980/232‑1 (MQ=255)                                                                                                                       gccATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGTCCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGaa                                                                                             <  1:153115/249‑1 (MQ=255)                                                                                                                                      tgCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGACGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGc                                                                                                                                                                     <  1:99610/164‑1 (MQ=255)                                                                                                                                        ccAACGTTTTACTCATGGGGTTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCtggtgg                                                                            <  1:255434/251‑1 (MQ=255)                                                                                                                                                                   acaGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCa                                                 <  1:300118/251‑1 (MQ=255)                                                                                                                                                                     aGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCGATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGc                                                  >  1:113470/1‑248 (MQ=255)                                                                                                                                                                             ggCAAAAAACGTTTATGCGCCAATATTGCTGCGCCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTGCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCTCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCTCGTTTTCAAGCAGACTCTGGAc                                       <  2:151192/251‑1 (MQ=255)                                                                                                                                                                                aaaaaaCGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGa                                     >  1:204742/1‑250 (MQ=255)                                                                                                                                                                                 aaaatCGTTGATGCGCCAATATTGCTGCTCCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGTCAAGGGCTCGTGCCGCGTTTTCAAGCAGACTCTGGACGaaa                                   <  1:127740/251‑1 (MQ=255)                                                                                                                                                                                    aaCGGTGATGCGCCTATATTGCTGCACCCTGAGCGCTATTAAAGTCATAAGATGAGGCAGCAATAAAGCGCTCGACCGGATGCGTCAGTGGACGGCGCAGGTACGTTTGGGGCATGGCAACGATAGTCTCGCGTGGGAAGGCGTGCATATCCAGCACGCCCCCGCCCAGGATCACCGCATCGGGAGCGAACAGATTAAGGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAgaa                                <  2:204705/251‑3 (MQ=255)                                                                                                                                                                                          gATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGttt                          >  1:178846/1‑251 (MQ=255)                                                                                                                                                                                                 cAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCAt                   <  1:1257/251‑1 (MQ=255)                                                                                                                                                                                                   atatTGCTGCGCCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATAGCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATgg                 <  1:23022/251‑1 (MQ=255)                                                                                                                                                                                                     atttcttatctcTGAGCGCCATTAAAGTCAGAAGAGGTGGCATCAATAAAGCGCGCGCCCTGATGCGGCAGAGGGCGGCGCAGGCACTTTTGGGTCAGGGCAACGATAGTCGCGCGTGTGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGa                <  2:304367/239‑1 (MQ=255)                                                                                                                                                                                                     atTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCCGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTg                                          >  1:99467/1‑224 (MQ=255)                                                                                                                                                                                                     atTGCGGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCCGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTg                                          <  2:99467/224‑1 (MQ=255)                                                                                                                                                                                                       tgctgcACACTGAGCGCCATTAAAGTCAGAAGGTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCCGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCCTCGGGATCGAAAAGGTTAATGCTGGTGGCAAGGGCCAGTGCCGCGTTTTCAAACAGACTCTGGACGAAAGAGACGTTTTCCGCATGGACg              >  2:219250/1‑250 (MQ=255)                                                                                                                                                                                                        gctgcACCCTGATCGCCATTAAAGGCAGAATATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCATGTACATTTGGGTCATGGCAACGAGAGTCTCTCGTGGGAAGGCGGGCATAGCCATCACGCCACCGCCCAGGATCCCCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGGGTTGTCAATCAGACTCTGGACGAAAGAGGCGTTTTCCTCATGGACGaa            <  1:333161/251‑1 (MQ=255)                                                                                                                                                                                                                   aGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGAACGAACAGATTAATGCTGGTGGCAATGGCCAGTGCAGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCacgc  >  2:193770/1‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                               ATTTTATTAGGCATTTGTGCTGAAGAGGGAAAATTACAGTAAAACTATCATTGCAGAAAATAATATGAGTCGTTTTATGTGTTTATTATATATATTGTACCTGTCGTTGGCACAGCGGGACATTGCGCTACTGTGCCAACGTTTTACTCATGGGGCTTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGC  >  NC_000913/4306723‑4307182

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