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

GGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCG  >  NC_000913/3567124‑3567383                                                                                                                                |                                                                                                                                     ggctggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTc                                                                                                                           <  1:1060106/139‑1 (MQ=255)             ggCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAG‑CCCATTGCGATTAGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTc                                                                                                              <  2:857795/139‑1 (MQ=255)                         cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTtctc                                                                                                   >  2:926780/1‑139 (MQ=255)                              aGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCa                                                                                              <  2:923268/139‑1 (MQ=255)                                     cTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGtt                                                                                       <  1:1414004/139‑1 (MQ=255)                                                      aGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGtt                                                                                       <  1:116305/122‑1 (MQ=255)                                                      aGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGtt                                                                                       >  2:116305/1‑122 (MQ=255)                                                      aGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGAc                                                                      >  1:1513832/1‑139 (MQ=255)                                                      aGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGAc                                                                      <  2:1236408/139‑1 (MQ=255)                                                            aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCg                                                                                         >  2:490779/1‑114 (MQ=255)                                                            aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCg                                                                                         <  1:490779/114‑1 (MQ=255)                                                            aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGaaa                                                                         <  1:159615/130‑1 (MQ=255)                                                            aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGaaa                                                                         >  2:159615/1‑130 (MQ=255)                                                            aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCtt                                                                >  2:75570/1‑139 (MQ=255)                                                             ttGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTc                                                               >  2:969005/1‑139 (MQ=255)                                                                cGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACg                                                            >  2:754638/1‑139 (MQ=255)                                                                       tAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACGACAGATTTTCTCGTCCACGCCGTTGATTACGCCGGAAAGAAGCCCTTCACGGTGACGc                                                     >  1:161289/1‑139 (MQ=255)                                                                             gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACt                                                                                                               <  2:493635/75‑1 (MQ=255)                                                                             gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACt                                                                                                               >  1:493635/1‑75 (MQ=255)                                                                             gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGc                                               >  2:1009236/1‑139 (MQ=255)                                                                             gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGc                                               >  2:928731/1‑139 (MQ=255)                                                                                  tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          <  2:703849/139‑1 (MQ=255)                                                                                              tATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                            <  1:262299/77‑1 (MQ=255)                                                                                              tATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                            >  2:262299/1‑77 (MQ=255)                                                                                                          tATCGCGGGTGTAATGCGAGGCTAACAGTAAGTCCGTCTCTGGAc                                                                                                                >  1:106894/1‑45 (MQ=255)                                                                                                          tATCGCGGGTGTAATGCGAGGCTAACAGTAAGTCCGTCTCTGGAc                                                                                                                <  2:106894/45‑1 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                            <  1:526248/65‑1 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                            >  2:526248/1‑65 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          >  2:677581/1‑115 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          <  1:677581/115‑1 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATTCCGTAGGCAAACTg                  >  1:921342/1‑139 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg                  >  1:1245902/1‑139 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACCCCGGAAAg                                                                        >  2:1075252/1‑85 (MQ=255)                                                                                                          tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTACAGTACCCCGGAAAg                                                                        <  1:1075252/85‑1 (MQ=255)                                                                                                                   tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGt         <  2:1513832/139‑1 (MQ=255)                                                                                                                       aCGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGAtc     >  1:817095/1‑139 (MQ=255)                                                                                                                       aCGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGAtc     >  1:340155/1‑139 (MQ=255)                                                                                                                          cgAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcg  <  1:1009236/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                     GGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCG  >  NC_000913/3567124‑3567383

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