Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,572,091	T→C	D471G (GAT→GGT)	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,091	0	T	C	100.0%	207.3 / NA	60	D471G (GAT→GGT)	glgB	1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/28);  total (32/28)

TAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAAT  >  NC_000913/3571959‑3572225                                                                                                                                     |                                                                                                                                       tAGTTGTGGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCgg                                                                                                                                 >  2:1188465/1‑140 (MQ=255)  aGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGa                                                                                                                                >  2:37493/1‑140 (MQ=255)         gagaATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGaaa                                                                                                                         >  1:13849/1‑140 (MQ=255)         gagaATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGaaa                                                                                                                         >  2:379026/1‑140 (MQ=255)         gagaATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGaaa                                                                                                                         <  2:872023/140‑1 (MQ=255)           gaATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACgaga                                                                                                                           <  1:576999/136‑1 (MQ=255)           gaATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACgaga                                                                                                                           >  2:576999/1‑136 (MQ=255)           gaATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACg                                                                                                                       >  1:220191/1‑140 (MQ=255)            aaTCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGc                                                                                                                      >  1:1478283/1‑140 (MQ=255)              tCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCa                                                                                                                    <  2:920248/140‑1 (MQ=255)              tCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCa                                                                                                                    >  1:820015/1‑140 (MQ=255)              tCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCa                                                                                                                    <  2:1290253/140‑1 (MQ=255)                 cGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCa                                                                                                                    <  2:1407920/137‑1 (MQ=255)                 cGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCa                                                                                                                    >  1:1407920/1‑137 (MQ=255)                     ggTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAAt                                                                                                             <  2:319990/140‑1 (MQ=255)                        cAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCgg                                                                                                          >  2:1457891/1‑140 (MQ=255)                          gTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTa                                                                                                        >  1:876649/1‑140 (MQ=255)                           tttATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAg                                                                                                       <  2:1116972/140‑1 (MQ=255)                            ttATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGa                                                                                                      <  2:320061/140‑1 (MQ=255)                                 gTGATGATCCTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGActcct                                                                                                 <  2:220191/140‑1 (MQ=255)                                 gTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGActcct                                                                                                 <  1:379026/140‑1 (MQ=255)                                         aCTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCAtt                                                                                         <  2:820015/140‑1 (MQ=255)                                              cGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcac                                                                                    <  1:968742/140‑1 (MQ=255)                                                aTAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcaccg                                                                                  <  1:662179/140‑1 (MQ=255)                                                aTAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAAGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcaccg                                                                                  <  1:164391/140‑1 (MQ=255)                                                  aaaCCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcaccgca                                                                                >  2:896210/1‑140 (MQ=255)                                                  aaaCCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcaccgca                                                                                >  1:817078/1‑140 (MQ=255)                                                       gggTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGCGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGa                                                                           <  2:1413046/140‑1 (MQ=255)                                                       gggTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGa                                                                           >  2:905897/1‑140 (MQ=255)                                                         gTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGCTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGaaa                                                                         >  1:349469/1‑140 (MQ=255)                                                              gCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGc                                                                    >  2:1126822/1‑140 (MQ=255)                                                               cTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcaccgcaccg                                                                             >  1:618914/1‑130 (MQ=255)                                                               cTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTcaccgcaccg                                                                             <  2:618914/130‑1 (MQ=255)                                                                ttCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTc                                                                  <  1:702250/140‑1 (MQ=255)                                                                   aTGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCAcc                                                               <  1:1188465/140‑1 (MQ=255)                                                                        gTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGaa                                                          >  2:948417/1‑140 (MQ=255)                                                                         tCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAAt                                                         <  1:37493/140‑1 (MQ=255)                                                                         tCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAAt                                                         <  2:817078/140‑1 (MQ=255)                                                                         tCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAAt                                                         >  1:1079608/1‑140 (MQ=255)                                                                                  tCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGaa                                                          >  1:576888/1‑130 (MQ=255)                                                                                  tCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGaa                                                          <  2:576888/130‑1 (MQ=255)                                                                                       cATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTAc                                           >  2:265752/1‑140 (MQ=255)                                                                                            aGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGc                                         <  1:1416208/137‑1 (MQ=255)                                                                                            aGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGc                                         >  2:1416208/1‑137 (MQ=255)                                                                                            aGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAg                                      >  2:707171/1‑140 (MQ=255)                                                                                                  ggTTCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCa                                <  1:265752/140‑1 (MQ=255)                                                                                                    ttCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCTCCGCACCGGAACCCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAAt                              <  2:261853/140‑1 (MQ=255)                                                                                                    ttCCACTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAAt                              <  2:349469/140‑1 (MQ=255)                                                                                                         cTTGTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCtt                         <  1:905897/140‑1 (MQ=255)                                                                                                            gTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCaa                      >  1:303345/1‑140 (MQ=255)                                                                                                            gTACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCaa                      >  1:772161/1‑140 (MQ=255)                                                                                                             tACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAg                     >  1:741657/1‑140 (MQ=255)                                                                                                             tACCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAg                     >  2:362153/1‑140 (MQ=255)                                                                                                              aCCAGAAGCCCAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGCCTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGa                    >  2:163124/1‑140 (MQ=255)                                                                                                                       ccAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCgcgc           >  2:1315083/1‑140 (MQ=255)                                                                                                                        cAGACCGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGcgcc          >  1:486018/1‑140 (MQ=255)                                                                                                                          ggccGCCCATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGcgcccg        <  2:211872/138‑1 (MQ=255)                                                                                                                               cccATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCaaa   >  1:881453/1‑140 (MQ=255)                                                                                                                               cccATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCaaa   >  2:272986/1‑140 (MQ=255)                                                                                                                                ccATACCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGTCTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAAt  <  1:163124/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATGATACTGACGATAAACCGGGTCGAGCTTCATGTAGTCCAGGGTGTCATGCATCCAGCCGAGGTTCCACTTGTACCAGAAGCCCAGACCGCCCATATCCTGCGGACGAGAAACGCCAGGGAAATCGGTAGACTCCTCAGCCATTGTCACCGCACCGGAAACCTGCTCACCAAGAATACGGTTGGTATTACGCAAGAATTCAATCGCTTCAAGATTCTCGCGCCCGCCAAAT  >  NC_000913/3571959‑3572225

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