Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,567,257	T→C	L259L (TTA→TTG)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,567,257	0	T	C	100.0%	67.1 / NA	20	L259L (TTA→TTG)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/5);  total (15/5)

GGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCG  >  NC_000913/3567128‑3567383                                                                                                                                  |                                                                                                                               ggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAg                                                                                                                             >  2:353538/1‑133 (MQ=255) ggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAg                                                                                                                             >  1:589095/1‑133 (MQ=255)                     cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTg                                                                                                                   <  1:613671/123‑1 (MQ=255)                     cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTg                                                                                                                   >  2:613671/1‑123 (MQ=255)                     cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGAt                                                                                                         >  2:60215/1‑133 (MQ=255)                        aaaGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      <  1:375715/133‑1 (MQ=255)                          aGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTtct                                                                                                    >  1:434744/1‑133 (MQ=255)                                           aaCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGt                                                                                   >  2:566250/1‑133 (MQ=255)                                                       cATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGa                                                                       >  1:349618/1‑133 (MQ=255)                                                        aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGAc                                                                      >  2:100211/1‑133 (MQ=255)                                                        aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGAc                                                                      >  2:572659/1‑133 (MQ=255)                                                                                      gCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGAc                                        <  1:353538/133‑1 (MQ=255)                                                                                      gCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGAc                                        <  1:14487/133‑1 (MQ=255)                                                                                       cTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGAcc                                       >  1:528811/1‑133 (MQ=255)                                                                                           aTCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTcc                                   >  1:80034/1‑133 (MQ=255)                                                                                                                  aaCGCGAGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACATACCTTCCCTACTGTCGGCAAAGTGCGGTTc            >  1:293825/1‑133 (MQ=255)                                                                                                                        aGGCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGCAGGAAATCTGCGGCTCGGTGAt      >  1:221365/1‑133 (MQ=255)                                                                                                                         ggCCAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCGACAGACCTTCCAAACCGCAGGCAAAGTGCGGATCGGTGAtc     >  2:295151/1‑133 (MQ=255)                                                                                                                           ccAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGAtctc   >  1:2731/1‑133 (MQ=255)                                                                                                                            cAACAGCAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcg  <  2:349618/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                               GGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCG  >  NC_000913/3567128‑3567383

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