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

TTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAA  >  NC_000913/3566469‑3566730                                                                                                                                    |                                                                                                                                   ttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCa                                                                                                                             <  2:1428085/139‑1 (MQ=255) ttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCa                                                                                                                             <  1:319008/139‑1 (MQ=255)     tGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCg                                                                                                                         >  2:1464531/1‑139 (MQ=255)         gacgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAg                                                                                                                     >  1:425560/1‑139 (MQ=255)                       ttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                          >  2:233571/1‑120 (MQ=255)                       ttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGc                                                                                                                          <  1:233571/120‑1 (MQ=255)                                         gCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTg                                                                                     >  2:1325551/1‑139 (MQ=255)                                         gCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTg                                                                                     <  2:425560/139‑1 (MQ=255)                                          cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                    <  2:595961/139‑1 (MQ=255)                                          cTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGc                                                                                    <  1:1464531/139‑1 (MQ=255)                                              tAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGc                                                                                >  2:1572472/1‑139 (MQ=255)                                                 gCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTaa                                                                             <  1:1066653/139‑1 (MQ=255)                                                   gATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTaaaa                                                                           >  1:1457131/1‑139 (MQ=255)                                                    aTAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAAt                                                                          >  1:1492612/1‑139 (MQ=255)                                                    aTAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATA‑GCGTTTCCTGAAAATTATTTCa                                                                                                                             <  2:839411/87‑1 (MQ=255)                                                    aTAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATA‑GCGTTTCCTGAAAATTATTTCa                                                                                                                             >  1:839411/1‑87 (MQ=255)                                                       gagTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                       <  2:1343093/139‑1 (MQ=255)                                                        agTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCAc                                                                                                           >  1:1556403/1‑102 (MQ=255)                                                        agTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCAc                                                                                                           <  2:1556403/102‑1 (MQ=255)                                                                               aGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACaa                                               >  1:1150150/1‑139 (MQ=255)                                                                                 cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                          <  2:498366/94‑1 (MQ=255)                                                                                 cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                          >  1:498366/1‑94 (MQ=255)                                                                                 cGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAAc                                             >  2:1137807/1‑139 (MQ=255)                                                                                     ggCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGcc                                         >  2:35235/1‑139 (MQ=255)                                                                                         aTGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAg                                     <  2:1150150/139‑1 (MQ=255)                                                                                                tGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGa                              >  1:323840/1‑139 (MQ=255)                                                                                                    aaCGGAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACggg                          >  2:1300779/1‑139 (MQ=255)                                                                                                        gAGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                      <  2:1492612/139‑1 (MQ=255)                                                                                                         aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCTATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAg                                <  2:981314/128‑1 (MQ=255)                                                                                                         aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                            >  1:844734/1‑100 (MQ=255)                                                                                                         aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGcc                                                            <  2:844734/100‑1 (MQ=255)                                                                                                         aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAg                                >  1:981314/1‑128 (MQ=255)                                                                                                         aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                      <  2:157266/138‑1 (MQ=255)                                                                                                         aGCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACca                      >  1:157266/1‑138 (MQ=255)                                                                                                          gCATTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcaca                    >  2:739739/1‑139 (MQ=255)                                                                                                            aTTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                     >  1:902449/1‑136 (MQ=255)                                                                                                            aTTCATATAGGCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACcac                     <  2:902449/136‑1 (MQ=255)                                                                                                                     ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGa                                                                                                    >  1:403828/1‑48 (MQ=255)                                                                                                                     ggCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGa                                                                                                    <  2:403828/48‑1 (MQ=255)                                                                                                                      gCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGa                                                                                         >  1:603204/1‑58 (MQ=255)                                                                                                                      gCGTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGa                                                                                         <  2:603204/58‑1 (MQ=255)                                                                                                                        gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                       <  1:1184943/74‑1 (MQ=255)                                                                                                                        gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCa                                                                       >  2:1184943/1‑74 (MQ=255)                                                                                                                        gTTTCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGt      >  2:57251/1‑139 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                          <  2:968103/52‑1 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCcgcg                                                                                          >  1:968103/1‑52 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                        <  2:406545/118‑1 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGAc                        >  1:406545/1‑118 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa   >  1:3286/1‑139 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa   >  1:1417465/1‑139 (MQ=255)                                                                                                                           tCCTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGa   >  2:963316/1‑139 (MQ=255)                                                                                                                            ccTGAAAATTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa  <  1:1572472/139‑1 (MQ=255)                                                                                                                                aaaaTTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa  <  2:806766/135‑1 (MQ=255)                                                                                                                                aaaaTTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGaa  >  1:806766/1‑135 (MQ=255)                                                                                                                                    |                                                                                                                                   TTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGCCGCGACCTGCCAGCTAAAATCCATTGCCATAGCCTGACGTTGCACAAACCGCCACAGTGAAGGACGGGACCACAGTACAAAAGCACGTCGAA  >  NC_000913/3566469‑3566730

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