Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,903,129	(G)5→6	coding (1234/1479 nt)	ygcE →	putative sugar kinase YgcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,903,124	1	.	G	100.0%	59.6 / NA	18	L410R (CTG→CGG)	ygcE	putative sugar kinase YgcE
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (9/9);  total (9/9)

AATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCT‑GGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGAT  >  NC_000913/2902991‑2903223                                                                                                                                       |                                                                                                    aatgaaatgaATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCT‑gggggg                                                                                               >  1:295819/1‑139 (MQ=255)   tgaaatgaaTCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCT‑ggggggca                                                                                             <  2:233617/140‑4 (MQ=255)                          aTGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCt                                                                       >  1:112299/1‑140 (MQ=255)                            gtgATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGt                                                                             >  2:306334/1‑132 (MQ=255)                            gtgATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGt                                                                             <  1:306334/132‑1 (MQ=255)                                tCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCt                                                                 <  2:112299/140‑1 (MQ=255)                                   ttACTGGCGGCGGTTCGAACAGCGATCTGTTTATACAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCtata                                                              >  1:284504/1‑140 (MQ=255)                                         gcggcgGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGa                                                        >  1:187445/1‑140 (MQ=255)                                           ggcggTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGAtt                                                      >  1:99171/1‑140 (MQ=255)                                            gcggTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTa                                                     >  1:9813/1‑140 (MQ=255)                                                 tCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCaa                                                <  2:135496/140‑1 (MQ=255)                                                 tCGAACAGCGATCTGTTTATACAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCaa                                                <  2:284504/140‑1 (MQ=255)                                                    aaCAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACgg                                             <  2:316046/140‑1 (MQ=255)                                                          gATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCg                                       <  2:295819/140‑1 (MQ=255)                                                                        gATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTcgc                         >  2:96135/1‑140 (MQ=255)                                                                          tttttGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGt                       <  2:1684/140‑1 (MQ=255)                                                                            tttGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGa                     <  2:9813/140‑1 (MQ=255)                                                                                     gtgtTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                 >  2:113565/1‑135 (MQ=255)                                                                                     gtgtTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGa                 <  1:113565/135‑1 (MQ=255)                                                                                               ttCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCTGGGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGAt  >  1:22982/1‑140 (MQ=255)                                                                                                                                       |                                                                                                    AATGAAATGAATCACTTTGCGAAGCATGTGATCATTACTGGCGGCGGTTCGAACAGCGATCTGTTTATGCAGATTTTTGCCGACGTGTTCAACCTTCCGGCACGACGTAACGCCATTAACGGTTGTGCAAGCCT‑GGGGGCAGCGATTAATACAGCGGTAGGTCTGGGGCTATACCCGGATTACGCAACGGCTGTCGATAACATGGTTCGCGTGAAAGATATCTTTATACCGAT  >  NC_000913/2902991‑2903223

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