Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,566,600	C→T	*478* (TAG→TAA)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,600	0	C	T	100.0%	200.1 / NA	62	*478* (TAG→TAA)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/30);  total (32/30)

TTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGA  >  NC_000913/3566472‑3566729                                                                                                                                 |                                                                                                                                  ttGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                         <  2:1267225/139‑1 (MQ=255)             gggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCt                                                                                                             <  1:1100855/139‑1 (MQ=255)                                    cAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                         <  2:1246986/104‑1 (MQ=255)                                    cAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                         >  1:1246986/1‑104 (MQ=255)                                    cAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                         <  1:500015/104‑1 (MQ=255)                                    cAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                         >  2:500015/1‑104 (MQ=255)                                                 aTAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAAt                                                                         >  1:410747/1‑139 (MQ=255)                                                         cTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGcc                                                                 >  2:577646/1‑139 (MQ=255)                                                               agCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                           <  2:410747/139‑1 (MQ=255)                                                                 cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGAAAAGCTCACGGTACGACTTCGACGCGACCTGCCAACTAAAATCCATTGCCATAGCCTg                                                         >  1:1590623/1‑139 (MQ=255)                                                                   tCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGAc                                                       >  1:517512/1‑139 (MQ=255)                                                                   tCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGAc                                                       >  2:1579905/1‑139 (MQ=255)                                                                         cTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGca                                                 <  1:149334/139‑1 (MQ=255)                                                                                 gggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                         >  1:672548/1‑139 (MQ=255)                                                                                  ggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCt                                                                                                             >  1:207322/1‑70 (MQ=255)                                                                                  ggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCt                                                                                                             <  2:207322/70‑1 (MQ=255)                                                                                   gCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCca                                       <  2:1590623/139‑1 (MQ=255)                                                                                        gAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTg                                  <  2:672548/139‑1 (MQ=255)                                                                                             tGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGa                             >  1:900156/1‑139 (MQ=255)                                                                                                 aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                  >  2:728073/1‑82 (MQ=255)                                                                                                 aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                  <  1:728073/82‑1 (MQ=255)                                                                                                 aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                   <  2:762899/81‑1 (MQ=255)                                                                                                 aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                   >  1:762899/1‑81 (MQ=255)                                                                                                   cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                           >  2:867208/1‑55 (MQ=255)                                                                                                   cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                           <  1:867208/55‑1 (MQ=255)                                                                                                     gAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                           <  1:1099914/101‑1 (MQ=255)                                                                                                     gAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                           >  2:1099914/1‑101 (MQ=255)                                                                                                     gAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                         >  2:988108/1‑119 (MQ=255)                                                                                                     gAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGc                                         <  1:988108/119‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGATAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                    <  2:900156/139‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                           <  1:1529499/52‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCa                                                                                                           >  2:1529499/1‑52 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                                                       >  2:222481/1‑72 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGAc                                                                                       <  1:222481/72‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                      >  2:257223/1‑89 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                      <  1:257223/89‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                     <  1:1209094/138‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                     <  1:1375106/138‑1 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                     >  2:1209094/1‑138 (MQ=255)                                                                                                      aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                     >  2:1375106/1‑138 (MQ=255)                                                                                                                taGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAg          >  2:1335797/1‑139 (MQ=255)                                                                                                                 aGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGc         <  1:1579905/139‑1 (MQ=255)                                                                                                                  ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCa        >  1:337846/1‑139 (MQ=255)                                                                                                                   gCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCca                                       >  1:428189/1‑107 (MQ=255)                                                                                                                   gCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCca                                       <  2:428189/107‑1 (MQ=255)                                                                                                                    cGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACg      <  1:577646/139‑1 (MQ=255)                                                                                                                     gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACt                                                                                                 >  1:921427/1‑47 (MQ=255)                                                                                                                     gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACt                                                                                                 <  2:921427/47‑1 (MQ=255)                                                                                                                       ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                      <  1:1075073/72‑1 (MQ=255)                                                                                                                       ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTg                                  >  1:226859/1‑108 (MQ=255)                                                                                                                       ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTg                                  <  2:226859/108‑1 (MQ=255)                                                                                                                       ttCCTGAAAATTATTTCAAACGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                      >  2:1075073/1‑72 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                   >  2:1103409/1‑58 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                   <  1:1103409/58‑1 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCa                                                                <  2:47107/77‑1 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCa                                                                >  1:47107/1‑77 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                               <  1:959849/94‑1 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                               >  2:959849/1‑94 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa  >  2:197785/1‑139 (MQ=255)                                                                                                                        tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa  >  2:1134516/1‑139 (MQ=255)                                                                                                                         ccTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACt                                                                                                 >  2:691244/1‑43 (MQ=255)                                                                                                                         ccTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACt                                                                                                 <  1:691244/43‑1 (MQ=255)                                                                                                                                 |                                                                                                                                  TTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGA  >  NC_000913/3566472‑3566729

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: ‑