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

TGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGC  >  NC_000913/3547853‑3548119                                                                                                                                      |                                                                                                                                      tgctgctgAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTActt                                                                                                                                  >  2:634185/1‑139 (MQ=255)  gctgctgAATATGGTGATTTGATAACACTGACTCCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCGGGCATCAACCGCACTTTActtc                                                                                                                                 >  1:295027/1‑139 (MQ=255)              ggTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCa                                                                                                                     >  2:794901/1‑139 (MQ=255)              ggTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCa                                                                                                                     >  2:524906/1‑139 (MQ=255)                     tGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTgc                                                                                                              >  1:99861/1‑139 (MQ=255)                          acacTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                         >  1:16195/1‑139 (MQ=255)                           cacTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTc                                                                                                        <  2:16195/139‑1 (MQ=255)                            actgactgCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCt                                                                                                       <  2:32662/139‑1 (MQ=255)                                     cGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                         >  2:697848/1‑128 (MQ=255)                                     cGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGt                                                                                                         <  1:697848/128‑1 (MQ=255)                                          gcggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCaaa                                                                                                >  1:602095/1‑132 (MQ=255)                                          gcggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCaaa                                                                                                <  2:602095/132‑1 (MQ=255)                                            ggcgTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGc                                                                                       <  2:1424449/139‑1 (MQ=255)                                                             cGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcatc                                                                      <  2:99861/139‑1 (MQ=255)                                                                     gATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCaa                                                              >  1:395854/1‑139 (MQ=255)                                                                          gggATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATc                                                         >  2:424433/1‑139 (MQ=255)                                                                                    agtagGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCg                                               <  1:833009/139‑1 (MQ=255)                                                                                          ccGGATAAGACGCGTCACGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcatc                                                                      >  1:958845/1‑110 (MQ=255)                                                                                          ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCcatc                                                                      <  2:958845/110‑1 (MQ=255)                                                                                          ccGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCgg                                         >  2:1279795/1‑139 (MQ=255)                                                                                           cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGc                                                                                       <  1:109788/92‑1 (MQ=255)                                                                                           cGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGc                                                                                       >  2:109788/1‑92 (MQ=255)                                                                                              aTAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt                                              <  1:1100541/130‑1 (MQ=255)                                                                                              aTAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGt                                              >  2:1100541/1‑130 (MQ=255)                                                                                               tAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                  >  2:128346/1‑77 (MQ=255)                                                                                               tAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCa                                                                                                  <  1:128346/77‑1 (MQ=255)                                                                                                   aCGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                   <  2:494918/56‑1 (MQ=255)                                                                                                   aCGCGTCAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGc                                                                                                                   >  1:494918/1‑56 (MQ=255)                                                                                                        tCAAGCATCGCATCCGGCATCTACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTACAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCATTGCGGAAAGCTTGcgccgc                           >  2:581759/1‑139 (MQ=255)                                                                                                         cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                          <  2:363509/139‑1 (MQ=255)                                                                                                         cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                          <  2:395854/139‑1 (MQ=255)                                                                                                         cAAGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGcc                          <  1:1279795/139‑1 (MQ=255)                                                                                                          aaGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                           <  1:879463/89‑1 (MQ=255)                                                                                                          aaGCATCGCATCCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCATCAGCAACTTGTTCACg                                                                           >  2:879463/1‑89 (MQ=255)                                                                                                                  catcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTtat                 <  2:1051733/139‑1 (MQ=255)                                                                                                                  catcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTtat                 >  1:834362/1‑139 (MQ=255)                                                                                                                   atcCGGCATCAACCGCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTtata                <  2:834362/139‑1 (MQ=255)                                                                                                                             aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTgc                                                                                                              <  1:1447982/35‑1 (MQ=37)                                                                                                                             aaCCGCACTTTACTTCTTCTTCGCTGCAGCTCTgc                                                                                                              >  2:1447982/1‑35 (MQ=37)                                                                                                                                 gCACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGc  <  1:634185/139‑1 (MQ=255)                                                                                                                                  cACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                           <  2:739178/65‑1 (MQ=255)                                                                                                                                  cACTTTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACg                                                                           >  1:739178/1‑65 (MQ=255)                                                                                                                                      |                                                                                                                                      TGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAGCATCGCATCCGGCATCAACCGCACTCTACTTCTTCTTCGCTGCAGCTCTGCGCCGTCTGTCCAAATCCTTCAGCAACTTGTTCACGCCATCATCGGCAAACATCGACTCAAGCGTTGCGGAAAGCTTGCGTCGCCAGTTTTTATACTGGTAACTGGTGC  >  NC_000913/3547853‑3548119

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