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

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

GACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCGCG  >  NC_000913/3567115‑3567385                                                                                                                                         |                                                                                                                                       gaCCTTTCTGGCTGGTCAGACGGCTCACCACCGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTaa                                                                                                                                      <  2:39468/139‑1 (MQ=255)     tttCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGt                                                                                                                                  <  2:740950/139‑1 (MQ=255)           gctggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTcc                                                                                                                            >  2:426349/1‑139 (MQ=255)           gctggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTcc                                                                                                                            >  2:181866/1‑139 (MQ=255)              ggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGtc                                                                                                                         >  2:670044/1‑139 (MQ=255)                tCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGtctc                                                                                                                       >  2:384334/1‑139 (MQ=255)                tCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGtctc                                                                                                                       >  1:551484/1‑139 (MQ=255)                                       aGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCa                                                                                                >  1:745476/1‑139 (MQ=255)                                                   tCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTa                                                                                    <  1:426349/139‑1 (MQ=255)                                                         aCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCgg                                                                              >  2:71155/1‑139 (MQ=255)                                                         aCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCgg                                                                              >  1:724926/1‑139 (MQ=255)                                                             tAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTATGACCGTCGCTGGACTCCAGATCTTCTCGTCCACGCCGTTCAGTACGCCGGAAAg                                                                          <  1:327376/139‑1 (MQ=255)                                                                     aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCtt                                                                  >  2:22648/1‑139 (MQ=255)                                                                           aTTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCTGTTCAGTACGCCGGAAAGACGCCCTTCACGGt                                                            <  1:496326/139‑1 (MQ=255)                                                                                    tGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtt                                                                                                          <  1:688309/84‑1 (MQ=255)                                                                                    tGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtt                                                                                                          >  2:688309/1‑84 (MQ=255)                                                                                      gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGc                                                 >  1:592678/1‑139 (MQ=255)                                                                                           tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                            <  2:479592/139‑1 (MQ=255)                                                                                           tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                            <  2:336390/139‑1 (MQ=255)                                                                                             ttttCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGAc                                          <  1:127528/139‑1 (MQ=255)                                                                                                             ccAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGc                          >  2:634582/1‑139 (MQ=255)                                                                                                                 cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCa                                                                <  2:796818/97‑1 (MQ=255)                                                                                                                 cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCa                                                                >  1:796818/1‑97 (MQ=255)                                                                                                                 cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       <  1:776537/106‑1 (MQ=255)                                                                                                                 cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       >  2:776537/1‑106 (MQ=255)                                                                                                                 cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAAc                      >  2:284006/1‑139 (MQ=255)                                                                                                                   tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg                    <  2:379450/139‑1 (MQ=255)                                                                                                                      cgcgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCgg                 <  1:618088/139‑1 (MQ=255)                                                                                                                          ggTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCg             >  2:243193/1‑139 (MQ=255)                                                                                                                            tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGt           >  1:473802/1‑139 (MQ=255)                                                                                                                                     aGGCTAATAGTAAGTCCGTCTCTGTACTCCAGATTTTCTCGTCCACGCCGTACAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGCTCGGTGATCTcgcg  >  2:489446/1‑139 (MQ=255)                                                                                                                                     aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCAGGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcgcg  >  1:508572/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       GACCTTTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCGCG  >  NC_000913/3567115‑3567385

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