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%	114.6 / NA	36	L261L (TTG→TTA)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/19);  total (17/19)

GGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGT  >  NC_000913/3567124‑3567376                                                                                                                                |                                                                                                                              ggctggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTc                                                                                                                    <  1:790901/139‑1 (MQ=255)       tCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGtctc                                                                                                              >  2:85318/1‑139 (MQ=255)                tCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCa                                                                                                     >  2:1081355/1‑139 (MQ=255)                      cTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                               <  1:85318/139‑1 (MQ=255)                         cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTtctc                                                                                            >  2:126688/1‑139 (MQ=255)                                            ttCAACCTTAAGCCCCATTGCGATTTGTAACTGCCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACg                                                                         <  1:776541/138‑1 (MQ=255)                                             tCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                     <  1:319617/126‑1 (MQ=255)                                             tCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                     >  2:169651/1‑126 (MQ=255)                                             tCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                     <  1:169651/126‑1 (MQ=255)                                             tCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                     >  2:319617/1‑126 (MQ=255)                                                       gCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACg                                                              <  1:1081355/139‑1 (MQ=255)                                                           cATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCt                                                          >  1:252429/1‑139 (MQ=255)                                                                     tGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACg                                                              >  1:440739/1‑125 (MQ=255)                                                                     tGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACg                                                              <  2:440739/125‑1 (MQ=255)                                                                     tGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGAc                                                <  2:252429/139‑1 (MQ=255)                                                                     tGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTGTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGAc                                                <  1:35855/139‑1 (MQ=255)                                                                            ggCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                              <  1:121705/134‑1 (MQ=255)                                                                            ggCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                              >  2:121705/1‑134 (MQ=255)                                                                              cgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGt                                                                                          <  2:497647/88‑1 (MQ=255)                                                                              cgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGt                                                                                          >  1:497647/1‑88 (MQ=255)                                                                                  tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACTCGAGGCTAACAGTAAGTCCGTCTCTGTACTCCAGATTTTCTCGTCCACACCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGTCGCTGTTGCAACAg                                   <  1:676700/139‑1 (MQ=255)                                                                                  tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                   <  2:144561/139‑1 (MQ=255)                                                                                        ccGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTc                                                                                                                    <  2:247484/52‑1 (MQ=255)                                                                                        ccGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTc                                                                                                                    >  1:247484/1‑52 (MQ=255)                                                                                                        cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGt                                                                            >  1:799472/1‑76 (MQ=255)                                                                                                        cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATAc                        >  1:320648/1‑128 (MQ=255)                                                                                                        cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATAc                        <  2:320648/128‑1 (MQ=255)                                                                                                        cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCCGGACTCCAGATTTTCTCGTCCACGCCGTTCAGt                                                                            <  2:799472/76‑1 (MQ=255)                                                                                                             cgcgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                              <  2:951071/101‑1 (MQ=255)                                                                                                             cgcgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                              >  1:951071/1‑101 (MQ=255)                                                                                                             cgcgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCgg        >  2:1088469/1‑139 (MQ=255)                                                                                                                   tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGt  >  1:582781/1‑139 (MQ=255)                                                                                                                     tAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg           <  2:606429/128‑1 (MQ=255)                                                                                                                     tAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg           >  1:606429/1‑128 (MQ=255)                                                                                                                            aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTa                                                                           >  2:532075/1‑57 (MQ=255)                                                                                                                            aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTa                                                                           <  1:532075/57‑1 (MQ=255)                                                                                                                                |                                                                                                                              GGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGT  >  NC_000913/3567124‑3567376

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