Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,572,871	Δ1 bp	coding (632/2187 nt)	glgB ←	1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,572,869	0	T	.	100.0%	215.4 / NA	52	coding (634/2187 nt)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (24/28);  total (24/28)

TGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTA  >  NC_000913/3572734‑3572996                                                                                                                                        |                                                                                                                                tgtgGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCAttgcg                                                                                                                             <  1:580435/140‑4 (MQ=255)    ggCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCtt                                                                                                                         >  2:577727/1‑140 (MQ=255)       gACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCaa                                                                                                                      <  1:577727/140‑1 (MQ=255)           ccAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGc                                                                                                                  >  2:595029/1‑140 (MQ=255)                aaCCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAg                                                                                                             <  1:435095/140‑1 (MQ=255)                   ccAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGt                                                                                                          <  2:370502/140‑1 (MQ=255)                   ccAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGt                                                                                                          <  1:772146/140‑1 (MQ=255)                   ccAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGt                                                                                                          <  1:538553/140‑1 (MQ=255)                            cTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCa                                                                                                 >  1:244239/1‑140 (MQ=255)                             ttCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAg                                                                                                >  2:322851/1‑140 (MQ=255)                                       gagaTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTg                                                                                      >  1:453867/1‑140 (MQ=255)                                       gagaTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTg                                                                                      >  1:71449/1‑140 (MQ=255)                                       gagaTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTg                                                                                      >  1:324083/1‑140 (MQ=255)                                         gaTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGcc                                                                                    >  2:625412/1‑140 (MQ=255)                                          aTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGACGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCa                                                                                   <  1:513284/140‑1 (MQ=255)                                          aTTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGACGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCa                                                                                   >  2:513284/1‑140 (MQ=255)                                           ttGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCAt                                                                                  <  1:230804/140‑1 (MQ=255)                                           ttGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCAt                                                                                  <  1:192552/140‑1 (MQ=255)                                              gcgcATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCAt                                                                                  >  1:81980/1‑137 (MQ=255)                                              gcgcATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCAt                                                                                  <  2:81980/137‑1 (MQ=255)                                              gcgcATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTgg                                                                               >  2:493386/1‑140 (MQ=255)                                               cgcATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGc                                                                              >  2:113701/1‑140 (MQ=255)                                               cgcATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGc                                                                              <  2:472578/140‑1 (MQ=255)                                                 cATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCat                                                                            <  1:578276/140‑1 (MQ=255)                                                   tCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCa                                                                                   >  2:263957/1‑131 (MQ=255)                                                   tCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCa                                                                                   <  1:263957/131‑1 (MQ=255)                                                   tCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCatca                                                                          <  1:568108/140‑1 (MQ=255)                                                   tCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCatca                                                                          >  2:588813/1‑140 (MQ=255)                                                     aaaCTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCatcaat                                                                        <  1:292206/140‑1 (MQ=255)                                                          gATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCAtct                                                                   <  1:595029/140‑1 (MQ=255)                                                             tCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGt                                                                <  1:690196/140‑1 (MQ=255)                                                             tCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGt                                                                <  2:437261/140‑1 (MQ=255)                                                                cTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTAtt                                                             >  1:764377/1‑140 (MQ=255)                                                                       gcgcTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATagag                                                      <  2:453867/140‑1 (MQ=255)                                                                                       aCAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGcc                                      >  2:548420/1‑140 (MQ=255)                                                                                        cAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGccc                                     <  1:113701/140‑1 (MQ=255)                                                                                        cAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCACA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGccc                                     <  2:407232/140‑1 (MQ=255)                                                                                           ccTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAg                                  >  1:336043/1‑140 (MQ=255)                                                                                            cTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAgg                                 >  2:419750/1‑140 (MQ=255)                                                                                              tttCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGa                               >  1:148255/1‑140 (MQ=255)                                                                                                tCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATa                             >  1:362133/1‑140 (MQ=255)                                                                                                 ccGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATaa                            <  1:540063/140‑1 (MQ=255)                                                                                                 ccGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATaa                            <  2:71449/140‑1 (MQ=255)                                                                                                  cGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATaaa                           <  2:148255/140‑1 (MQ=255)                                                                                                  cGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATaaa                           <  1:322851/140‑1 (MQ=255)                                                                                                                   gagaCGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGc          <  2:336043/140‑1 (MQ=255)                                                                                                                   gagaCGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGc          <  1:588813/140‑1 (MQ=255)                                                                                                                    agaCGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCCTTGGCATCAATCCTCTCGTATTTATAGAGCTTACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGct         <  2:171220/140‑1 (MQ=255)                                                                                                                    agaCGCGGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGct         <  1:469931/140‑1 (MQ=255)                                                                                                                       cgcgGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGctctc      >  2:250452/1‑140 (MQ=255)                                                                                                                        gcgGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGctctct     >  1:75485/1‑140 (MQ=255)                                                                                                                        gcgGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGctctct     <  1:642717/140‑1 (MQ=255)                                                                                                                         cgGTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCtt    >  2:49079/1‑140 (MQ=255)                                                                                                                           gTTTCCGGGCGCA‑TTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTa  >  1:540273/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                TGTGGCGACGCCAGGAACCCAGGTGAACTTCATAAATAGAGATTGGCGCATCAAACTGATTCGCTTTTTTGCGCTCTTCAGTCTGTACAACCTTTTCCGGCAGCCCGCAAATAAGAGACGCGGTTTCCGGGCGCATTTGCGCTTCAAAGGCATAAGGGTCGGACTTCAGACGCAAGTTGCCATTGGCATCAATCATCTCGTATTTATAGAGCTGACCGTTATGCGCCCCAGGGATAAACAGTTCCCAGATGCCGCTCTCTTTA  >  NC_000913/3572734‑3572996

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