Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,548,915	T→C	T386A (ACC→GCC)	malQ ←	4‑alpha‑glucanotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,548,915	0	T	C	100.0%	94.2 / NA	30	T386A (ACC→GCC)	malQ	4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/12);  total (18/12)

GCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGA  >  NC_000913/3548791‑3549049                                                                                                                              |                                                                                                                                       gCACGCGCGGTGATGATATGCGGGTCCATTGGCGGGAATCCCCCGTTCTTCCCCAACGGCGCGGGGATAATCCGG‑GGCGCCCCCAGCCTGGCTTTCAAGGAATATAAGTCACGGCAAAACCAGGCTTCCGCCCCACCGTc                                                                                                                         >  2:376212/1‑140 (MQ=255)      cgcgGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTa                                                                                                                    >  2:453627/1‑140 (MQ=255)             atgatATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAAc                                                                                                             <  2:497389/140‑1 (MQ=255)                atatGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑TTCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGc                                                                                                          <  2:906856/140‑1 (MQ=255)                atatGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGc                                                                                                          <  1:984510/140‑1 (MQ=255)                atatGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGc                                                                                                          <  2:974461/140‑1 (MQ=255)                  atGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCa                                                                                                        <  2:224899/140‑1 (MQ=255)                     cGGGTCCATTGGCGGTTATTCCCTGTTCTGCCCCAACGGGCCGAGGAT‑AAACGGCGGCGCGCCCACCGATGCCTTTACGCAAAATATTTCACGGTCCCACCAAGCATCCGCCCCCCCTTCCCCTACGCCAACCGCCCGAt                                                                                                     >  2:374036/1‑140 (MQ=255)                       ggTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCa                                                                                                   <  1:398215/140‑1 (MQ=255)                       ggTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCa                                                                                                   <  1:376212/140‑1 (MQ=255)                       ggTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGGCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCa                                                                                                   >  2:123354/1‑140 (MQ=255)                          ccATTGGCGGTAATCCCCCGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGCTGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCCACCGCCAGATCAcca                                                                                                >  2:391677/1‑138 (MQ=255)                            aTTGGCGGCAATCCCCAGTTCTCCCCCAACGGGCCGAGGGT‑ATGCGGGGGCGCGGCAAACGAGGCTTTCCGGGAATATAAATAACGGGCACACCAGGCTTTCGCCCCACCTTTCCCTACGCCAACCGCCAGATCACGAtt                                                                                              >  2:930487/1‑139 (MQ=255)                             ttGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATAc                                                                                             >  1:752589/1‑140 (MQ=255)                                        ccccAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCgg                                                                                  <  2:662206/140‑1 (MQ=255)                                                 tGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATa                                                                         >  1:986686/1‑140 (MQ=255)                                                   ccccAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCCCCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGc                                                                       <  2:565149/140‑1 (MQ=255)                                                   ccccAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGc                                                                       >  2:1009750/1‑140 (MQ=255)                                                   ccccAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGc                                                                       >  1:436978/1‑140 (MQ=255)                                                        aCGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCTATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTgg                                                                  >  2:718200/1‑140 (MQ=255)                                                            gCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTa                                                              >  2:761527/1‑140 (MQ=255)                                                                       tCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCcagca                                                    >  2:832189/1‑140 (MQ=255)                                                                               cgcgcCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCaa                                            >  2:539129/1‑140 (MQ=255)                                                                                 cgcCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAAc                                          >  1:249405/1‑140 (MQ=255)                                                                                 cgcCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAAc                                          >  1:927368/1‑140 (MQ=255)                                                                                     aaCCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGc                                      <  2:986686/140‑1 (MQ=255)                                                                                           tGCTTTCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAag                                >  1:609692/1‑140 (MQ=255)                                                                                               ttCAGGCAATATAGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAagccag                            <  1:1009750/140‑1 (MQ=255)                                                                                                         ataGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAAc                  >  2:32439/1‑140 (MQ=255)                                                                                                         ataGTTCACGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAAc                  <  1:832189/140‑1 (MQ=255)                                                                                                                aCGGTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGAt           >  1:522083/1‑140 (MQ=255)                                                                                                                  ggTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATaa         >  2:188785/1‑140 (MQ=255)                                                                                                                  ggTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATaa         >  1:909006/1‑140 (MQ=255)                                                                                                                  ggTCACACCAGGCTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATaa         >  1:296544/1‑140 (MQ=255)                                                                                                                         ccAGGCTTCCGCCCCACCTCCCCCAACGCCACCCGCCAGATCCCGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGa  <  2:377171/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                       GCACGCGCGGTGATGATATGCGGGTCCATTGGCGGTAATCCCCAGTTCTGCCCCAACGGGCCGAGGAT‑ATCCGGCGGCGCGCCAACCGATGCTTTCAGGCAATATAGTTCACGGTCACACCAGGTTTCCGCCCCACCTTCCGCTACGCCAACCGCCAGATCACGATACAAGCCAATCGGCATTTCATAGCCCTGGCTTATCTCCCAGCAGGCGGCAAACTGGCTGTAAGCCAGCCACTGCAACCAGAGATAAAAATCGA  >  NC_000913/3548791‑3549049

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