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%	175.2 / NA	54	*695* (TAG→TAA)	malQ	4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/25);  total (29/25)

TCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCA  >  NC_000913/3547851‑3548121                                                                                                                                        |                                                                                                                                        tCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTAc                                                                                                                                      <  1:428465/139‑1 (MQ=255) tCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTAc                                                                                                                                      >  1:862165/1‑139 (MQ=255)                   gATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGct                                                                                                                    <  1:918951/139‑1 (MQ=255)                              actgactgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCt                                                                                                         >  2:261032/1‑139 (MQ=255)                                  actgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTcc                                                                                                     >  1:877795/1‑139 (MQ=255)                                   ctgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                    >  1:1038514/1‑139 (MQ=255)                                   ctgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                    >  1:206958/1‑139 (MQ=255)                                   ctgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                    >  2:89187/1‑139 (MQ=255)                                     gCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCTTCGCATCCGGCATCATCCGCACTTTATTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCaaa                                                                                                  <  1:672120/139‑1 (MQ=255)                                     gCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCaaa                                                                                                  <  1:261032/139‑1 (MQ=255)                                                cgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCaa                                                                                       >  1:160779/1‑139 (MQ=255)                                                                ggCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcatca                                                                       <  2:45424/139‑1 (MQ=255)                                                                  ccTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCatcatc                                                                     >  1:58177/1‑139 (MQ=255)                                                                         ttCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAAc                                                              <  2:973934/139‑1 (MQ=255)                                                                           cGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACAt                                                            >  2:460327/1‑139 (MQ=255)                                                                            gggATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATc                                                           <  2:58177/139‑1 (MQ=255)                                                                                tGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTAc                                                                                                                                      >  1:567894/1‑60 (MQ=255)                                                                                tGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTAc                                                                                                                                      <  2:567894/60‑1 (MQ=255)                                                                                   aTTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCaa                                                    >  1:168545/1‑139 (MQ=255)                                                                                        tagGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACa                                                             <  2:494565/125‑1 (MQ=255)                                                                                        tagGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACa                                                             >  1:494565/1‑125 (MQ=255)                                                                                         agGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                              <  2:320823/59‑1 (MQ=255)                                                                                         agGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                              >  1:320823/1‑59 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                     >  2:811778/1‑65 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                     <  1:811778/65‑1 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACAt                                                            <  1:872105/122‑1 (MQ=255)                                                                                            ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACAt                                                            >  2:872105/1‑122 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                              <  1:585972/55‑1 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGc                                                                  <  1:971428/115‑1 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGc                                                                  >  2:971428/1‑115 (MQ=255)                                                                                             cGGATAAGACGCGTCAAGAATCGCATCCGGCATCAACCGCACTTTActtcttctt                                                                                                                              >  2:585972/1‑55 (MQ=25)                                                                                              ggATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcat                                                                         <  1:802934/107‑1 (MQ=255)                                                                                              ggATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcat                                                                         >  2:802934/1‑107 (MQ=255)                                                                                               gATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttc                                                                                                                                <  1:651351/51‑1 (MQ=255)                                                                                               gATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActtcttc                                                                                                                                >  2:651351/1‑51 (MQ=255)                                                                                               gATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCa                                                                                           <  1:1011465/88‑1 (MQ=255)                                                                                               gATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCa                                                                                           >  2:1011465/1‑88 (MQ=255)                                                                                               gATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGaaa                                        >  1:487826/1‑139 (MQ=255)                                                                                                    gACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTg                                   >  2:721296/1‑139 (MQ=255)                                                                                                         gTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCg                              >  2:869034/1‑139 (MQ=255)                                                                                                         gTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCg                              >  2:357552/1‑139 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                            <  2:393381/139‑1 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                            <  2:1038514/139‑1 (MQ=255)                                                                                                           cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                            <  2:862165/139‑1 (MQ=255)                                                                                                               catcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGtt                        >  2:468398/1‑139 (MQ=255)                                                                                                                atcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGttt                       <  1:89187/139‑1 (MQ=255)                                                                                                                 tcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                           >  2:1024624/1‑54 (MQ=255)                                                                                                                 tcgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                           <  1:1024624/54‑1 (MQ=255)                                                                                                                  cgcatcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGttttt                     <  1:468398/139‑1 (MQ=255)                                                                                                                          gCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGt             >  2:375825/1‑139 (MQ=255)                                                                                                                               aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                            <  2:846104/119‑1 (MQ=255)                                                                                                                               aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                            >  1:846104/1‑119 (MQ=255)                                                                                                                                   gCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGc    <  2:451889/139‑1 (MQ=255)                                                                                                                                     aCTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCa  >  2:664216/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        TCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGCCA  >  NC_000913/3547851‑3548121

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