Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,547,986	C→T	*695* (TAG→TAA)	malQ ←	4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,547,986	0	C	T	100.0%	208.9 / NA	64	*695* (TAG→TAA)	malQ	4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/27);  total (37/27)

TCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCC  >  NC_000913/3547851‑3548120                                                                                                                                        |                                                                                                                                       tCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTAc                                                                                                                                     >  2:761017/1‑139 (MQ=255)    gctgctgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtc                                                                                                                                  <  1:1001509/139‑1 (MQ=255)    gctgctgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtc                                                                                                                                  >  2:1001509/1‑139 (MQ=255)         tgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                             <  1:901824/139‑1 (MQ=255)           aaTATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCg                                                                                                                           >  2:888581/1‑139 (MQ=255)                  tGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                    >  1:469171/1‑139 (MQ=255)                         aTAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTgcgc                                                                                                             >  2:64528/1‑139 (MQ=255)                             cacTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTc                                                                                                         <  1:761017/139‑1 (MQ=255)                               gtgactgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGAGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTg                                                                                                       >  1:807338/2‑139 (MQ=255)                                tgactgCCGGATGCGGCGTGACCGCCTTATCCAGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGt                                                                                                      >  2:781905/1‑139 (MQ=255)                                  actgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTcc                                                                                                    >  1:633370/1‑139 (MQ=255)                                            gcggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCa                                                                                          >  2:1224621/1‑139 (MQ=255)                                            gcggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCa                                                                                          <  2:469171/139‑1 (MQ=255)                                            gcggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCa                                                                                          <  1:750241/139‑1 (MQ=255)                                                cgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                   >  1:202736/1‑126 (MQ=255)                                                cgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                   <  2:202736/126‑1 (MQ=255)                                                 gTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAAc                                                                                     >  2:368932/1‑139 (MQ=255)                                                      cGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGtt                                                                                >  1:3951/1‑139 (MQ=255)                                                         cTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCAc                                                                             >  2:760485/1‑139 (MQ=255)                                                          ttATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                            >  1:664252/1‑139 (MQ=255)                                                          ttATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                            <  1:781905/139‑1 (MQ=255)                                                                   cTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCg                                                                   >  1:158993/1‑139 (MQ=255)                                                                       gATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTc                                                                               >  2:1013537/1‑123 (MQ=255)                                                                       gATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTc                                                                               <  1:1013537/123‑1 (MQ=255)                                                                            gggATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATc                                                          <  1:368932/139‑1 (MQ=255)                                                                                 gAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                            >  2:361518/1‑116 (MQ=255)                                                                                 gAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                            <  1:361518/116‑1 (MQ=255)                                                                                    ttAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAg                                                  <  1:888581/139‑1 (MQ=255)                                                                                      agtagGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCg                                                <  2:633370/139‑1 (MQ=255)                                                                                          ggCCGGATAAGACGCGTCAAGCATCGCATCCGACATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGc                                            >  2:1126700/1‑139 (MQ=255)                                                                                           gCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCTCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCaa                                                   <  2:896314/131‑1 (MQ=255)                                                                                           gCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCaa                                                   >  1:896314/1‑131 (MQ=255)                                                                                           gCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCg                                           >  1:1264253/1‑139 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTg                                                                                                                        >  1:71941/1‑61 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTg                                                                                                                        <  2:71941/61‑1 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                          >  1:361468/1‑139 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                          >  2:618851/1‑139 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                          >  2:118478/1‑139 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                          >  1:123975/1‑139 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                          >  1:141098/1‑139 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAAc                                                                                     <  2:393204/95‑1 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAAc                                                                                     >  1:393204/1‑95 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGa                                         >  2:77543/1‑139 (MQ=255)                                                                                                aTAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGc                                                                                        <  2:89686/89‑1 (MQ=255)                                                                                                aTAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGc                                                                                        >  1:89686/1‑89 (MQ=255)                                                                                                  aaGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCt                                    >  1:114265/1‑139 (MQ=255)                                                                                                     aCGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTg                                             >  2:1226231/1‑127 (MQ=255)                                                                                                     aCGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTg                                             <  1:1226231/127‑1 (MQ=255)                                                                                                         gTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                   <  1:618406/69‑1 (MQ=255)                                                                                                         gTCAAGCATCGCATCCGGCATCAACCGCACTTTACTACTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                   >  2:618406/1‑69 (MQ=255)                                                                                                          tCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGc                                                                 <  1:1305282/102‑1 (MQ=255)                                                                                                          tCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGc                                                                 >  2:1305282/1‑102 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                           <  2:884254/139‑1 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                           <  2:807338/139‑1 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                           <  1:1224621/139‑1 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                           <  2:114265/139‑1 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                           <  2:141098/139‑1 (MQ=255)                                                                                                              gcatcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGt                        >  2:650835/1‑139 (MQ=255)                                                                                                                  cgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGttttt                    <  1:760485/139‑1 (MQ=255)                                                                                                                         ggCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTgg             >  2:50210/1‑139 (MQ=255)                                                                                                                              cAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACt        <  2:3951/139‑1 (MQ=255)                                                                                                                               aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGt            >  2:725660/1‑134 (MQ=255)                                                                                                                               aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGt            <  1:725660/134‑1 (MQ=255)                                                                                                                                    cACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGCTTTTATACTGGTAACTGGTGcc  >  2:806157/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       TCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCC  >  NC_000913/3547851‑3548120

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