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%	204.8 / NA	63	*478* (TAG→TAA)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/27);  total (36/27)

GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATA  >  NC_000913/3566484‑3566732                                                                                                                     |                                                                                                                                     gggTCGTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGCCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCt                                                                                                                <  2:913771/139‑1 (MQ=255)    tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCAc                                                                                                             >  1:591614/1‑139 (MQ=255)                           cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                      >  2:196599/1‑139 (MQ=255)                            ttGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     <  2:243342/139‑1 (MQ=255)                               tAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTTAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGc                                                                                  >  2:888014/1‑139 (MQ=255)                                aaGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGACGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCt                                                                                 >  2:592285/1‑139 (MQ=255)                                  gCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTaa                                                                               <  2:591614/139‑1 (MQ=255)                                       agagTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTCAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATcc                                                                          >  1:849275/1‑139 (MQ=255)                                         agTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGa                                                                                                                          >  1:1172512/1‑89 (MQ=255)                                         agTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGa                                                                                                                          <  2:1172512/89‑1 (MQ=255)                                         agTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCAt                                                                        >  2:1266495/1‑139 (MQ=255)                                                     cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     >  2:49977/1‑114 (MQ=255)                                                     cTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGcc                                                                                     <  1:49977/114‑1 (MQ=255)                                                       tCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                            <  2:397015/73‑1 (MQ=255)                                                       tCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                            >  1:397015/1‑73 (MQ=255)                                                                aGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                                  <  2:1467588/138‑1 (MQ=255)                                                                aGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                                  >  1:1467588/1‑138 (MQ=255)                                                                    tGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCg                                             <  1:196599/139‑1 (MQ=255)                                                                      ggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGcc                                           >  2:129890/1‑139 (MQ=255)                                                                       gCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCca                                          <  1:888014/139‑1 (MQ=255)                                                                          aTGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCa                                                                   <  2:1266223/111‑1 (MQ=255)                                                                          aTGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCa                                                                   >  1:1266223/1‑111 (MQ=255)                                                                            gAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                               <  2:1214800/129‑1 (MQ=255)                                                                            gAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                               >  1:1214800/1‑129 (MQ=255)                                                                                     aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                                  <  2:663264/117‑1 (MQ=255)                                                                                     aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGcaca                                                  >  1:663264/1‑117 (MQ=255)                                                                                     aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTg                                     <  1:800256/130‑1 (MQ=255)                                                                                     aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTg                                     >  2:800256/1‑130 (MQ=255)                                                                                     aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACggg                            >  2:904686/1‑139 (MQ=255)                                                                                       cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                              >  1:892982/1‑103 (MQ=255)                                                                                       cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                              <  2:892982/103‑1 (MQ=255)                                                                                       cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  1:904686/139‑1 (MQ=255)                                                                                       cGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          >  2:1357731/1‑139 (MQ=255)                                                                                       cGGAGCATTCATATAGGCGTTTCCAGAAAATTATTTCAAGAGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          >  2:897322/1‑139 (MQ=255)                                                                                          aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                       >  2:1305671/1‑139 (MQ=255)                                                                                          aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                       >  2:1054704/1‑139 (MQ=255)                                                                                          aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                        >  1:7953/1‑138 (MQ=255)                                                                                          aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                        <  2:7953/138‑1 (MQ=255)                                                                                            cATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgc                                                                                                <  2:1278614/64‑1 (MQ=255)                                                                                            cATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgc                                                                                                >  1:1278614/1‑64 (MQ=255)                                                                                            cATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAg                     >  2:468641/1‑139 (MQ=255)                                                                                             aTTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGt                    >  1:589418/1‑139 (MQ=255)                                                                                                      ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCa           >  1:226685/1‑139 (MQ=255)                                                                                                        cGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                               <  2:1418245/101‑1 (MQ=255)                                                                                                        cGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                               >  1:1418245/1‑101 (MQ=255)                                                                                                          tttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCca                                          <  2:202785/104‑1 (MQ=255)                                                                                                          tttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCca                                          >  1:202785/1‑104 (MQ=255)                                                                                                           ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCTCGACCTGCCAGCTAAAATCCATTGCCGTAGCCTGACGTTGCACAAACCGCCACAGTg                                     <  1:297639/108‑1 (MQ=255)                                                                                                           ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTg                                     >  2:297639/1‑108 (MQ=255)                                                                                                           ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGa                           <  2:30068/118‑1 (MQ=255)                                                                                                           ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGa                           >  1:30068/1‑118 (MQ=255)                                                                                                           ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          >  2:426043/1‑119 (MQ=255)                                                                                                           ttCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  1:426043/119‑1 (MQ=255)                                                                                                            tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          <  2:861828/118‑1 (MQ=255)                                                                                                            tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                          >  1:861828/1‑118 (MQ=255)                                                                                                            tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACg         <  2:351947/135‑1 (MQ=255)                                                                                                            tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACg         >  1:351947/1‑135 (MQ=255)                                                                                                            tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa     >  1:471159/1‑139 (MQ=255)                                                                                                            tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa     >  1:856372/1‑139 (MQ=255)                                                                                                             ccTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa    <  1:1357731/139‑1 (MQ=255)                                                                                                              cTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                            <  1:811863/50‑1 (MQ=255)                                                                                                              cTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                            >  2:811863/1‑50 (MQ=255)                                                                                                               tGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATa  <  2:856372/139‑1 (MQ=255)                                                                                                                     |                                                                                                                                     GGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAATA  >  NC_000913/3566484‑3566732

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