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

GTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCG  >  NC_000913/3567129‑3567383                                                                                                                           |                                                                                                                                     gTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACTTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGtct                                                                                                                      >  2:451349/1‑139 (MQ=255)         gCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTc                                                                                                              <  1:1268115/139‑1 (MQ=255)                 cTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      <  2:73601/139‑1 (MQ=255)                 cTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      <  2:340409/139‑1 (MQ=255)                    cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTctct                                                                                                                    >  1:329211/1‑122 (MQ=255)                    cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTctct                                                                                                                    <  2:329211/122‑1 (MQ=255)                    cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTcc                                                                                                             <  2:1239739/129‑1 (MQ=255)                    cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTcc                                                                                                             >  1:1239739/1‑129 (MQ=255)                                cTTTATCGACAATCTTCAGCCAAATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGCCTACCAGTATGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGtt                                                                                       <  2:1316577/139‑1 (MQ=255)                                           aCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCgg                                                                            >  2:246610/1‑139 (MQ=255)                                                      cATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCt                                                                 >  1:1273684/1‑139 (MQ=255)                                                       aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGt                                                                                                 <  2:97058/106‑1 (MQ=255)                                                       aTTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGt                                                                                                 >  1:97058/1‑106 (MQ=255)                                                                 gTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAg                                                                                    <  2:886783/109‑1 (MQ=255)                                                                 gTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAg                                                                                    >  1:886783/1‑109 (MQ=255)                                                                  tAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                     >  2:524640/1‑139 (MQ=255)                                                                     cTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGt                                                  <  1:528793/139‑1 (MQ=255)                                                                        gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGc                                               >  1:498461/1‑139 (MQ=255)                                                                        gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGc                                               >  1:85935/1‑139 (MQ=255)                                                                             tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          <  2:498461/139‑1 (MQ=255)                                                                             tATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          <  1:984625/139‑1 (MQ=255)                                                                                        ttATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATAc                               <  1:451349/139‑1 (MQ=255)                                                                                                cAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCa                       >  2:668849/1‑139 (MQ=255)                                                                                                   cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      <  1:1085524/57‑1 (MQ=255)                                                                                                   cGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      >  2:1085524/1‑57 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCATGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg                  >  2:706851/1‑139 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGc                                                                                           <  1:1026821/66‑1 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGc                                                                                           >  2:1026821/1‑66 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGc                                                                    >  1:333354/1‑89 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGc                                                                    <  2:333354/89‑1 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          >  2:1165002/1‑115 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          <  1:1165002/115‑1 (MQ=255)                                                                                                     tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg                  >  1:728454/1‑139 (MQ=255)                                                                                                       tCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      >  1:1339386/1‑53 (MQ=255)                                                                                                       tCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                      <  2:1339386/53‑1 (MQ=255)                                                                                                       tCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACg                                                                     >  2:1227821/1‑86 (MQ=255)                                                                                                       tCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACg                                                                     <  1:1227821/86‑1 (MQ=255)                                                                                                          cgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAg                                                                                                           >  2:282286/1‑45 (MQ=255)                                                                                                          cgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAg                                                                                                           <  1:282286/45‑1 (MQ=255)                                                                                                              tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          >  1:1174736/1‑106 (MQ=255)                                                                                                              tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                          <  2:1174736/106‑1 (MQ=255)                                                                                                               gtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGAc                                        >  1:880534/1‑107 (MQ=255)                                                                                                               gtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGAc                                        <  2:880534/107‑1 (MQ=255)                                                                                                                tAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGt                                                          <  2:1172765/88‑1 (MQ=255)                                                                                                                tAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGt                                                          >  1:1172765/1‑88 (MQ=255)                                                                                                                 aaCGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCATGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGAt      <  1:706851/139‑1 (MQ=255)                                                                                                                     cgAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCa                                  >  2:906410/1‑107 (MQ=255)                                                                                                                     cgAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCa                                  <  1:906410/107‑1 (MQ=255)                                                                                                                     cgAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcg  <  1:524640/139‑1 (MQ=255)                                                                                                                       aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACg                                                      <  1:922560/85‑1 (MQ=255)                                                                                                                       aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACg                                                      >  2:922560/1‑85 (MQ=255)                                                                                                                           |                                                                                                                                     GTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCG  >  NC_000913/3567129‑3567383

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