Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,549,792	(C)6→7	coding (279/2085 nt)	malQ ←	4‑alpha‑glucanotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,549,786	1	.	C	100.0%	162.9 / NA	46	G95G (GGC→GGG)	malQ	4‑alpha‑glucanotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (20/26);  total (20/26)

TTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACT  >  NC_000913/3549647‑3549912                                                                                                                                             |                                                                                                                               tttGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTg                                                                                                                                 <  1:704656/140‑1 (MQ=255)      ttATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑ccccc                                                                                                                           <  1:99848/140‑1 (MQ=255)       tATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑cccccc                                                                                                                          <  2:317085/140‑1 (MQ=255)       tATTCAGCAACGCATGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑cccccc                                                                                                                          >  1:840196/1‑140 (MQ=255)          tCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑cccccccgt                                                                                                                       >  1:570950/1‑137 (MQ=255)               aaCGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACAt                                                                                                                   <  2:688149/140‑1 (MQ=255)               aaCGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACAt                                                                                                                   <  2:570950/140‑1 (MQ=255)                 cGCCTGCGGTTCGTAACAGCGTTTCGTGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATgg                                                                                                                 <  1:735311/140‑1 (MQ=255)                   ccTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACAt                                                                                                                   >  1:28632/1‑136 (MQ=255)                   ccTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACAt                                                                                                                   <  2:28632/136‑1 (MQ=255)                   ccTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGcc                                                                                                               >  2:481227/1‑140 (MQ=255)                   ccTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑ccccccc                                                                                                                         <  2:1272934/129‑2 (MQ=255)                   ccTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑ccccccc                                                                                                                         >  1:1272934/1‑128 (MQ=255)                        ggTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGGTAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGt                                                                                                          >  1:1166310/1‑140 (MQ=255)                         gTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTa                                                                                                         <  1:43642/140‑1 (MQ=255)                            cGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTg                                                                                                      >  2:552652/1‑140 (MQ=255)                                  aGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTcc                                                                                                <  1:552652/140‑1 (MQ=255)                                    cGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCctt                                                                                              >  1:437606/1‑140 (MQ=255)                                      tttCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCcttct                                                                                            >  1:576475/1‑140 (MQ=255)                                      tttCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCcttct                                                                                            >  1:219301/1‑140 (MQ=255)                                                     ctcccGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTAAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                                                             <  2:87096/138‑1 (MQ=255)                                                     cACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                                                             <  2:1153373/140‑1 (MQ=255)                                                        ccGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagccag                                                                          <  2:576475/140‑1 (MQ=255)                                                          ggCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTca                                                                                  <  1:1097967/130‑1 (MQ=255)                                                          ggCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTca                                                                                  >  2:1097967/1‑130 (MQ=255)                                                            cAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                                                             >  2:1314489/1‑133 (MQ=255)                                                            cAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                                                             <  1:1314489/133‑1 (MQ=255)                                                               tGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATAtt                                                                   >  2:1314070/1‑140 (MQ=255)                                                                    gcgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccg                                                              <  1:97931/140‑1 (MQ=255)                                                                     cgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccgc                                                             <  1:481227/140‑1 (MQ=255)                                                                       cTGATCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                                                           <  1:1096149/140‑1 (MQ=255)                                                                        tGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGc                                                          <  2:437606/140‑1 (MQ=255)                                                                                  gggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccaccac                                                >  1:720684/1‑140 (MQ=255)                                                                                    gTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                                                           <  2:610724/127‑1 (MQ=255)                                                                                    gTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                                                           >  1:610724/1‑127 (MQ=255)                                                                                         tgtCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATCGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGc                                         >  2:202042/1‑140 (MQ=255)                                                                                             aGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAttt                                     >  2:99560/1‑140 (MQ=255)                                                                                               cGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAttttt                                   >  2:1197476/1‑140 (MQ=255)                                                                                                     aTAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCg                             <  1:1314070/140‑1 (MQ=255)                                                                                                     aTAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCg                             <  2:840196/140‑1 (MQ=255)                                                                                                      tAACCTTCCGGCAGCTTCGTCGGTAGATCGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGc                            <  1:202042/140‑1 (MQ=255)                                                                                                        aCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTg                          >  2:1251892/1‑140 (MQ=255)                                                                                                           ttCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGta                       <  2:871839/140‑1 (MQ=255)                                                                                                              cGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATaa                    >  1:248393/1‑140 (MQ=255)                                                                                                               ggCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCg                             >  1:332598/1‑130 (MQ=255)                                                                                                               ggCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCg                             <  2:332598/130‑1 (MQ=255)                                                                                                               ggCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATaaa                   <  2:248393/140‑1 (MQ=255)                                                                                                                    cTTCGTCGGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATg              <  2:1163744/140‑1 (MQ=255)                                                                                                                         tcgGTAGATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACAtt         <  1:1251892/140‑1 (MQ=255)                                                                                                                               gATTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGAc   <  1:1306664/140‑1 (MQ=255)                                                                                                                                aTTGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACt  >  2:158055/1‑140 (MQ=255)                                                                                                                                  tGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGa    >  2:1041435/1‑136 (MQ=255)                                                                                                                                  tGAACGCTTTGCCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGa    <  1:1041435/136‑1 (MQ=255)                                                                                                                                             |                                                                                                                               TTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACT  >  NC_000913/3549647‑3549912

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