Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,549,583	T→G	D163A (GAT→GCT)	malQ ←	4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,549,583	0	T	G	100.0%	181.9 / NA	52	D163A (GAT→GCT)	malQ	4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (20/32);  total (20/32)

TTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGG  >  NC_000913/3549442‑3549730                                                                                                                                              |                                                                                                                                                    ttCAGCCACCGGCGAGAAGGCGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGACCGCATTAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATa                                                                                                                                             <  1:799426/150‑1 (MQ=255)  tCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATAcc                                                                                                                                           >  2:315408/1‑151 (MQ=255)      ccAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAg                                                                                                                                       <  2:20920/151‑1 (MQ=255)      ccAACGGCGAGAAGACGGGCTGTATGGGCGGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCa                                                                                                                                        <  2:295404/150‑1 (MQ=255)         aCGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGttt                                                                                                                                    >  2:642919/1‑151 (MQ=255)           ggCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGttttt                                                                                                                                  >  2:864169/1‑151 (MQ=255)               agaagaCGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTcc                                                                                                                              >  2:767728/1‑151 (MQ=255)                 aagaGGGGCTGTATGGGCTGGCGTTCTCCGGTTTTGCCGGATAGAGCGAATGAATCGGGTTCAGGCCAAGGAGCGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGtttttt                                                                                                                                 <  2:540895/146‑1 (MQ=255)                 aagaCGGGCTGTATGGGCTGGCGCTCTCCGGTTTGGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGa                                                                                                                            <  2:85489/151‑1 (MQ=255)                 aagaCGGGCGGTATGGGCTGGTGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGa                                                                                                                            <  1:328031/151‑1 (MQ=255)                      gggCTGTAGGGGCTGGCGCTCCCCGGTTTTGCCGGATAGAGCCCATGAATCGGGTTCAGCCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCa                                                                                                                       <  2:523899/151‑1 (MQ=255)                        gCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGc                                                                                                                     <  2:208670/151‑1 (MQ=255)                          tGTATGGCCTGGCGCGCCCCGGATTTGCGGGATAGACCGCATGAATCGGGTTCAGGCCAATAAACGCCCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCg                                                                                                                    <  2:624134/150‑1 (MQ=255)                              tGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGtat                                                                                                                 <  2:408035/149‑1 (MQ=255)                               gggCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATaa                                                                                                               >  1:266427/1‑150 (MQ=255)                                ggCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAg                                                                                                             <  2:758361/151‑1 (MQ=255)                                       gcgccccGGAGTAGCCGGATAGGGAGCACGATTGGGATGCCGGCCAAGTAACGCCCCGCCACGTTTTGCCACATCCACCAGCATCGATTTGAGATCGCCAAAAGCCCCAATTCCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACg                                                                                                      <  1:529458/146‑1 (MQ=255)                                         tctcCGAATTTGCCGGATAGGGCGCATAAATCGGGTTCAGACCAATAAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACg                                                                                                      <  1:875559/149‑1 (MQ=255)                                          ctcCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGc                                                                                                     <  2:760794/149‑1 (MQ=255)                                               gATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCAcc                                                                                              >  1:630991/1‑151 (MQ=255)                                                 tttGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCAcccc                                                                                            >  1:440114/1‑151 (MQ=255)                                                  ttGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCAcccc                                                                                            >  2:297997/1‑150 (MQ=255)                                                     ccGGAGAGGGCGAATGAATCGGGTGCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAg                                                                                        <  1:526221/151‑1 (MQ=255)                                                           agagCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTg                                                                                  >  1:85821/1‑151 (MQ=255)                                                            gagCGAATGAAGCGGGTTCAGGCAAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGt                                                                                 <  1:859237/151‑1 (MQ=255)                                                                 cATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTAt                                                                             <  1:693620/150‑1 (MQ=255)                                                                          ggTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGc                                                                   <  2:662816/151‑1 (MQ=255)                                                                             tCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTg                                                                >  1:527189/1‑151 (MQ=255)                                                                              cAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGc                                                               <  2:413056/151‑1 (MQ=255)                                                                                ggCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCg                                                              <  2:731062/150‑1 (MQ=255)                                                                                   cAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTc                                                          >  2:713643/1‑151 (MQ=255)                                                                                    aaTGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCg                                                         >  2:785385/1‑151 (MQ=255)                                                                                     aTGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGt                                                        <  1:769326/151‑1 (MQ=255)                                                                                              ccGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGtt                                               <  2:631011/151‑1 (MQ=255)                                                                                               cGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGttt                                              <  2:930922/151‑1 (MQ=255)                                                                                                  cACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCgg                                           <  2:864064/151‑1 (MQ=255)                                                                                                     gTTTTGCCACATCCACCAGCATCGCTTTGATATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGc                                        <  1:653374/151‑1 (MQ=255)                                                                                                     gTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGc                                        >  2:710750/1‑151 (MQ=255)                                                                                                           ccACATCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCGCCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAAt                                  <  2:349482/151‑1 (MQ=255)                                                                                                               aTCCACCAGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCAcc                              >  1:386596/1‑151 (MQ=255)                                                                                                                 ccaccaGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCAccc                             <  2:37969/150‑1 (MQ=255)                                                                                                                 ccaccaGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCg                            >  2:414344/1‑151 (MQ=255)                                                                                                                    ccaGCATCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCa                         <  1:7964/151‑1 (MQ=255)                                                                                                                         aTCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATgcg                    <  2:70641/151‑1 (MQ=255)                                                                                                                         aTCGCTTTGAGATCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATgcg                    >  2:274923/1‑151 (MQ=255)                                                                                                                                    aTCGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGgtcg         >  1:896766/1‑151 (MQ=255)                                                                                                                                      cGCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGgtcgtc       <  1:713261/151‑1 (MQ=255)                                                                                                                                       gCCAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTcc      >  1:901910/1‑151 (MQ=255)                                                                                                                                         cAAAAGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTg    <  1:775580/151‑1 (MQ=255)                                                                                                                                          aaaaGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTgg   >  2:785298/1‑151 (MQ=255)                                                                                                                                          aaaaGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTgg   >  2:396287/1‑151 (MQ=255)                                                                                                                                           aaaGCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTggg  <  1:676129/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    TTCAGCCAACGGCGAGAAGACGGGCTGTATGGGCTGGCGCTCTCCGGATTTGCCGGATAGAGCGCATGAATCGGGTTCAGGCCAATGAACGACCCGCCACGTTTTGCCACATCCACCAGCATCGCTTTGAGATCGCCAAAATCCCCAATACCCCAGTTTTTTTCCGATCGCAGCGTATAAAGCTGAACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGG  >  NC_000913/3549442‑3549730

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