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

TTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCGCG  >  NC_000913/3567120‑3567385                                                                                                                                    |                                                                                                                                       ttCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTa                                                                                                                                 >  1:649617/1‑139 (MQ=255)         ggTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGtc                                                                                                                         >  1:647549/1‑139 (MQ=255)                    tCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTg                                                                                                                     >  1:714114/1‑132 (MQ=255)                    tCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTg                                                                                                                     <  2:714114/132‑1 (MQ=255)                    tCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCa                                                                                                              >  2:134515/1‑139 (MQ=255)                          cTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGAtttt                                                                                                        <  2:1192278/139‑1 (MQ=255)                           tGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTtc                                                                                                       >  2:151746/1‑139 (MQ=255)                             cAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTtctc                                                                                                     >  1:242123/1‑139 (MQ=255)                                  aGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCa                                                                                                >  2:1186766/1‑139 (MQ=255)                                  aGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCa                                                                                                >  2:384204/1‑139 (MQ=255)                                    cGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                              <  2:647549/139‑1 (MQ=255)                                    cGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                              <  1:151746/139‑1 (MQ=255)                                       cACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCg                                                                                           <  1:134515/139‑1 (MQ=255)                                            tATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAg                                                                                      >  1:1314392/1‑139 (MQ=255)                                                   aaCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCg                                                                               <  1:1186766/139‑1 (MQ=255)                                                    aCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTg                                                                                                                     >  1:1236193/1‑100 (MQ=255)                                                    aCCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTg                                                                                                                     <  2:1236193/100‑1 (MQ=255)                                                           gCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACg                                                                       <  1:168421/139‑1 (MQ=255)                                                               cATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCt                                                                   >  2:692867/1‑139 (MQ=255)                                                                  tGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCa                                                                <  1:692867/139‑1 (MQ=255)                                                                        ttGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGa                                                          <  2:1314392/139‑1 (MQ=255)                                                                         tGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGAc                                                         >  1:432096/1‑139 (MQ=255)                                                                          gTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACg                                                        <  2:1394320/139‑1 (MQ=255)                                                                                ggCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCTAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       >  1:929089/1‑134 (MQ=255)                                                                                ggCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       <  2:929089/134‑1 (MQ=255)                                                                                 gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       <  2:998283/133‑1 (MQ=255)                                                                                 gcgcTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       >  1:998283/1‑133 (MQ=255)                                                                                    cTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAAc                                              <  2:649617/139‑1 (MQ=255)                                                                                    cTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAAc                                              <  1:791472/139‑1 (MQ=255)                                                                                                          aaCGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTc                                                                                                                             <  1:1412888/38‑1 (MQ=255)                                                                                                          aaCGTATCGCGGGTGTAACGCGAGGCTAACAGTAAGTc                                                                                                                             >  2:1412888/1‑38 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGATTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCaa                                               <  2:1015036/112‑1 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACt                                                                                                                 >  1:552021/1‑46 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACt                                                                                                                 <  2:552021/46‑1 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGAc                                                         <  2:1240766/102‑1 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGAc                                                         >  1:1240766/1‑102 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCaa                                               >  1:1015036/1‑112 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg                    >  1:1160374/1‑139 (MQ=255)                                                                                                              tATCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTg                    >  1:612587/1‑139 (MQ=255)                                                                                                               aTCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACg                                                              >  1:1468640/1‑96 (MQ=255)                                                                                                               aTCGCGGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACg                                                              <  2:1468640/96‑1 (MQ=255)                                                                                                                 cgcgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGc                                                                                 <  1:1230572/75‑1 (MQ=255)                                                                                                                 cgcgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGc                                                                                 >  2:1230572/1‑75 (MQ=255)                                                                                                                   cgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTa                                                                                    <  2:776772/70‑1 (MQ=255)                                                                                                                   cgGGTGTAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTa                                                                                    >  1:776772/1‑70 (MQ=255)                                                                                                                      gtgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCgg            >  1:208203/1‑139 (MQ=255)                                                                                                                       tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                              <  1:1126740/56‑1 (MQ=255)                                                                                                                       tgtAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACg                                                                                              >  2:1126740/1‑56 (MQ=255)                                                                                                                         tAACGCGAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGa         >  2:769978/1‑139 (MQ=255)                                                                                                                              cgAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcg    <  2:1160374/139‑1 (MQ=255)                                                                                                                               gAGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcgc   <  2:242123/139‑1 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       <  2:1406567/86‑1 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGc                                                       >  1:1406567/1‑86 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCTACAg                                            >  2:115970/1‑97 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCTACAg                                            <  1:115970/97‑1 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                            <  2:524891/97‑1 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAg                                            >  1:524891/1‑97 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcgcg  >  2:1120521/1‑139 (MQ=255)                                                                                                                                aGGCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTcgcg  >  2:1142748/1‑139 (MQ=255)                                                                                                                                 ggCTAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGAc                                                         >  1:370085/1‑83 (MQ=255)                                                                                                                                    |                                                                                                                                       TTCTGGCTGGTCAGACGGCTCACCACTGCAAAAAGCGGCACTTTATCGTCAACCTTAAGCCCCATTGCGATTTGTAACTGGCGCTTATTTTCCGCTTTATCTTCCAACGTATCGCGGGTGTAACGCGAGGCCAACAGTAAGTCCGTCTCTGGACTCCAGATTTTCTCGTCCACGCCGTTCAGTACGCCGGAAAGACGCCCTTCACGGTGACGCTGTTGCAACAGACCTTCCATACCGTAGGCAAACTGCGGTTCGGTGATCTCGCG  >  NC_000913/3567120‑3567385

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