Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,323,624	G→A	S149S (AGC→AGT)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,323,624	0	G	A	100.0%	73.3 / NA	25	S149S (AGC→AGT)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/11);  total (14/11)

GCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATAGAGAAGAATTTA  >  NC_000913/3323500‑3323738                                                                                                                             |                                                                                                                   gcgcGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCCACCCGCGCCATAAATGATACGACTGGCttt                                                                                                            >  1:196018/1‑133 (MQ=255)                     ccGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCaa                                                                                        <  1:348058/133‑1 (MQ=255)                      cGTTTGCACAATCCACCACAATCTGCAGCTCATTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAAc                                                                                       <  2:196018/133‑1 (MQ=255)                        tttGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGc                                                                                     <  2:310464/133‑1 (MQ=255)                          tGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAg                                                                                   >  1:585436/1‑133 (MQ=255)                                    accacAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCt                                                                         >  2:64151/1‑133 (MQ=255)                                           tCTTCAGTTCATTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCAt                                                                  >  1:618451/1‑133 (MQ=255)                                           tCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCAt                                                                  >  1:570301/1‑133 (MQ=255)                                             ttCAGTTCATTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCAttt                                                                <  2:618451/133‑1 (MQ=255)                                                        gAGGCTAAGTACGTTGGGGAACGTGGCTGTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATg                                                     <  2:222848/133‑1 (MQ=255)                                                                     ttcgGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCg                                        >  1:621440/1‑133 (MQ=255)                                                                     ttcgGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCg                                        >  2:520848/1‑133 (MQ=255)                                                                           aaCGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCg                                  >  2:37110/1‑133 (MQ=255)                                                                             cGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGc                                >  2:510091/1‑133 (MQ=255)                                                                               tGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAg                              >  2:520498/1‑133 (MQ=255)                                                                                ggCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGt                             <  1:510091/133‑1 (MQ=255)                                                                                    ttGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGTTGGCCTCTTCTGCCCCATCCGGCAGTTTgg                         >  1:215154/1‑133 (MQ=255)                                                                                     tGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGt                        <  1:517473/133‑1 (MQ=255)                                                                                     tGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGt                        <  2:555100/133‑1 (MQ=255)                                                                                          aaCTCGATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGt                   <  2:585436/133‑1 (MQ=255)                                                                                               gATATAGCGACCCGCGGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATa              <  2:570301/133‑1 (MQ=255)                                                                                                          ccGCGGCACCAACGATACTACTTGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATAGAGAAGAAttt   <  1:202579/133‑1 (MQ=255)                                                                                                           cgcgGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATAGAGAAGAATTTa  >  1:612795/1‑133 (MQ=255)                                                                                                           cgcgGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATAGAGAAGAATTTa  >  1:582679/1‑133 (MQ=255)                                                                                                           cgcgGCATCAACGATACGACTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATAGAGAAGAATTTa  >  1:456816/1‑133 (MQ=255)                                                                                                                             |                                                                                                                   GCGCGATGTGATAAGTCGCACCGTTTGCACAATCCACCACAATCTTCAGTTCACTGAGGCTAAGTTCGTTCGGGAACGTGGCTTTGCAAAACTCGATATAGCGACCCGCGGCATCAACGATACGGCTGGCTTTACCCAGTTCTGCCGAATCAACGCAGCTGATCTCCTTTTCCATTTCCGCTTCGATGGCCTCTTCTACCGCATCCGGCAGTTTGGTGCCGTCGATAGAGAAGAATTTA  >  NC_000913/3323500‑3323738

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