Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,151	G→A	P217S (CCT→TCT)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,151	0	G	A	97.0%	93.6 / NA	33	P217S (CCT→TCT)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (18/14);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCC  >  NC_000913/4306911‑4307392                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gcgcCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACATATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTcc                                                                                                                                                                                                                                         <  1:168595/251‑1 (MQ=255)     ccatatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGCCTCTGGCCGAAAGAGGCGTTTTCCGCAt                                                                                                                                                                                                                                     >  2:217197/3‑251 (MQ=255)      aaTATTGCTGCACCCCGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGTCGTTTTCCGCATg                                                                                                                                                                                                                                    <  1:207420/251‑1 (MQ=255)       atatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATgg                                                                                                                                                                                                                                   >  1:145750/1‑251 (MQ=255)              tgcACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGCCGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGCTCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGaca                                                                                                                                                                                                                             >  1:246920/1‑248 (MQ=255)                  cccTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATc                                                                                                                                                                                                                        <  2:257550/251‑1 (MQ=255)                   ccTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGaaa                                                                                                                                                                                                                             >  2:218431/1‑245 (MQ=255)                           ccATTAACGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATAACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATgg                                                                                                                                                                                                               >  1:202516/1‑251 (MQ=255)                             aTTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATgggg                                                                                                                                                                                                             >  2:76752/1‑250 (MQ=255)                               tAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGTCTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTa                                                                                                                                                                                                            <  2:168146/250‑1 (MQ=255)                                                   gcAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTa                                                                                                                                                                                       <  2:14119/251‑1 (MQ=255)                                                      aTAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCa                                                                                                                                                                                    <  2:105703/251‑1 (MQ=255)                                                       tAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCAACGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTAcc                                                                                                                                                                                     >  1:170868/1‑249 (MQ=255)                                                            cgcACGACCTGATGCGGCAGTGGGCGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCt                                                                                                                                                                              <  2:246920/251‑1 (MQ=255)                                                                              aGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGTGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTTGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAAGTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCa                                                                                                                                                            <  1:38617/251‑1 (MQ=255)                                                                                  gACGGAGCAGGTACTTTTGGGGCATGGCCACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCCCCGCATCGGGATCGAACAGATTAATGATGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGACAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGACATGGGTAATTTCGGGGATGTTGTTCGTAACAGCGTCTTAGCGCCATTCCAGAGCAATTg                                                                                                                                                        >  2:98050/1‑251 (MQ=255)                                                                                   aCGGCGCAGGTACGTTTGGGTCAGGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGTGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAAGTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTgg                                                                                                                                                       <  2:101153/251‑1 (MQ=255)                                                                                     ggCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGACAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGtt                                                                                                                                                     >  1:262700/1‑251 (MQ=255)                                                                                              aCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGACAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCAc                                                                                                                                            >  2:2518/1‑251 (MQ=255)                                                                                                       tCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTg                                                                                                                                   >  1:134824/1‑251 (MQ=255)                                                                                                                                gAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCCGTGTTGGGTCatat                                                                                                          >  2:76182/1‑251 (MQ=255)                                                                                                                                     cGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTccc                                                                                                     <  1:49180/251‑1 (MQ=255)                                                                                                                                                             ccAGGATCACCGCATCGTGATCGAACAGGTTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGAATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGACGTGTTGGGTCATATCTCCCCGGGGGGTATGACCCAGTTCCCCt                                                                             >  2:164204/1‑251 (MQ=255)                                                                                                                                                                    cACCTCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACGCTGGACGAAAGAGGCGTTTGCCGCAGGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTTTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCATTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCacac                                                                      <  1:263415/251‑1 (MQ=255)                                                                                                                                                                             gggATCGAACAGGTTAATGCTGGTGGCACTGGCCCGAGCCGCGTTTTCACGCCGACTCTGGACGGAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCCGCGTCTTAACGCCATTCCAGAGCACTGGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGGGTTGGGACATCTCTACCAGGGGGATATGGACCAGTTCGCCGGCCACACCGGGTGCAc                                                             >  2:149352/1‑251 (MQ=255)                                                                                                                                                                                    aaCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTccc                                                      >  1:147784/1‑250 (MQ=255)                                                                                                                                                                                        gATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGCCTCGGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGACCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATGGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGGTCGCCTGCCACACCGTGTGCACCCGTCCACGgc                                                  <  1:217197/251‑1 (MQ=255)                                                                                                                                                                                                   gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAg                                                                                    <  2:127418/206‑1 (MQ=255)                                                                                                                                                                                                   gtggCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAg                                                                                    >  1:127418/1‑206 (MQ=255)                                                                                                                                                                                                        aaTGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTAGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                                   >  1:25712/1‑250 (MQ=255)                                                                                                                                                                                                                              aGCAGACTCTGGCCGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTCCCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc            <  1:76182/251‑1 (MQ=255)                                                                                                                                                                                                                                      ctGGACGAAAGGGGCGTTTTCCGAATGGACAAAAAGACCGCGCAATGTGTAATTTCGGGGCTTTTGTTTGTACCAACGGCTTAGCGCCATCCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACCCGCGCTGTGTTGGGTCATATCCCCCAGGGGGATATGACTAAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa      <  1:86655/249‑1 (MQ=255)                                                                                                                                                                                                                                        ggACGAAAGAGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGCTAGgcc  >  2:80616/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCACCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCC  >  NC_000913/4306911‑4307392

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