Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,116,406	C→T	M272I (ATG→ATA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,406	0	C	T	100.0%	78.5 / NA	27	M272I (ATG→ATA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/13);  total (14/13)

TCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTA  >  NC_000913/4116266‑4116538                                                                                                                                             |                                                                                                                                     tCGCCCACCCACGGAATGGATGCGCCTGCCATAAACCCCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAagca                                                                                                                              <  1:118687/149‑1 (MQ=255)     gcAGCCACTGAATGGATGTGCCTGCCATAAACACCGCACCTTCCAACGCTGAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTGTTCTGATTTCACCGCTTTCTCGCCGGTGTTCATCAGTATAAAGCAGCCa                                                                                                                          <  2:141171/149‑1 (MQ=255)           aCTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGGGTTCATCAGTATAAAGCAGCCAGTGccg                                                                                                                    >  1:128211/1‑148 (MQ=255)                  gtaTGCGCCTCTCATACACCCCGCGCCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCATTGTTCATCAGTATAAAACAGCCAGTGCCATAGgtgt                                                                                                             <  2:403968/147‑1 (MQ=255)                   gacgcgcCTGCCACAAACCCCGCCCCTTCCATCGCATTGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGtt                                                                                                            <  1:156587/146‑1 (MQ=255)                         ccggCCATAAACCCCGCACCTTCCCACGCATCGTTCACCTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGc                                                                                                      <  2:51829/146‑1 (MQ=255)                         ccTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGc                                                                                                      <  2:273469/149‑1 (MQ=255)                              cATACACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATcc                                                                                                 <  1:82423/149‑1 (MQ=255)                                tAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTCCATCCGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCct                                                                                               >  2:22719/1‑149 (MQ=255)                                 aaaCACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCctt                                                                                              <  1:93017/149‑1 (MQ=255)                                  aaCACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCGTCCcttc                                                                                             >  1:5748/1‑149 (MQ=255)                                      cggCACGTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATGGGTTTTCTTCGCCATCCCTTCTTTc                                                                                         >  2:119171/3‑149 (MQ=255)                                      ccGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCCGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTc                                                                                         >  1:22741/1‑149 (MQ=255)                                      ccGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTc                                                                                         >  1:127155/1‑149 (MQ=255)                                                cAACGCATAGTTCACTTCGCCAGTCGGGCCGCAAGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACt                                                                               >  1:326762/1‑149 (MQ=255)                                                 aaCGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTg                                                                              <  2:417489/149‑1 (MQ=255)                                                                    cAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTAACGCACAACTGACCAAACAGCGCGGCctgc                                                           >  1:169726/1‑149 (MQ=255)                                                                       tCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTgg                                                        >  2:381491/1‑149 (MQ=255)                                                                            ccGCAGGCGATGGTGGTCAGCAGGAAGTTGTCTAATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAACCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCCCGCACAACTGACCAAACAGCGCGGCCTGCTGGTcccc                                                   >  2:431828/1‑149 (MQ=255)                                                                             cGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCg                                                  >  1:4662/1‑149 (MQ=255)                                                                                      tcgtggtCCGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCg                                         <  1:14247/147‑1 (MQ=255)                                                                                                   gCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATAc                            >  2:78132/1‑149 (MQ=255)                                                                                                        tttCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTg                       >  1:357121/1‑149 (MQ=255)                                                                                                                ttCGCCGCTTTGTCGCCAGGGTTCAACAGTATAAGGCAGCCAGTGCCATAGGTGTTCTTCGGCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGGGCCGCCttt               <  2:426233/149‑1 (MQ=255)                                                                                                                 tCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTg              >  1:8302/1‑149 (MQ=255)                                                                                                                        tttCTCCCCGGTGTCCTTCGGTATAAAGCAGCCAGGGCCATAGGTGTTCTTCGCCACCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa       <  2:127155/149‑1 (MQ=255)                                                                                                                             cGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa  <  2:369748/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                     TCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTA  >  NC_000913/4116266‑4116538

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