Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,006,526	A→G	D124G (GAC→GGC)	amyA →	cytoplasmic alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,006,526	0	A	G	100.0%	136.0 / NA	39	D124G (GAC→GGC)	amyA	cytoplasmic alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/21);  total (18/21)

CTGCTGGCCGCCATTGATGCTCTGAAACGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGACCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGTATCGACCATATCGAAAACC  >  NC_000913/2006402‑2006660                                                                                                                             |                                                                                                                                       ctgctgGCCGCCATTGATGCTCTGAAACGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTg                                                                                                                          <  2:27008/139‑1 (MQ=255)   gctgGCCGCCATTGATGCTCTGAAATGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGAc                                                                                                                        <  1:302608/139‑1 (MQ=255)     tgGCCGCCATTGATGCTCTGAAACGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACga                                                                                                                      <  2:164837/139‑1 (MQ=255)              ttGATGCTCTGAAACGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCAt                                                                                                             >  2:863850/1‑139 (MQ=255)                gATGCTCTGAAACGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTg                                                                                                           <  2:37587/139‑1 (MQ=255)                               aaTGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCg                                                                                                                                     >  2:285440/1‑98 (MQ=255)                               aaTGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCg                                                                                                                                     <  1:285440/98‑1 (MQ=255)                               aaTGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGa                                                                                            <  1:475668/139‑1 (MQ=255)                                 tGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGAcg                                                                                          >  2:108601/1‑139 (MQ=255)                                       tGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTac                                                                                    >  1:367555/1‑139 (MQ=255)                                                  tGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCcccgc                                                                         <  2:367555/139‑1 (MQ=255)                                                         ggTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGcc                                                                  <  2:366123/139‑1 (MQ=255)                                                          gTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCg                                                                 >  2:714148/1‑139 (MQ=255)                                                           tAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACgaagaa                                                                                                                  <  2:527376/89‑1 (MQ=255)                                                           tAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACgaagaa                                                                                                                  >  1:527376/1‑89 (MQ=255)                                                              tCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCa                                                             >  1:397107/1‑139 (MQ=255)                                                                   cacaAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACt                                                        <  2:544940/139‑1 (MQ=255)                                                                    acaAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTc                                                       <  1:108601/139‑1 (MQ=255)                                                                     caAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCAtt                                                                                                            <  2:13560/85‑1 (MQ=255)                                                                     caAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCAtt                                                                                                            >  1:13560/1‑85 (MQ=255)                                                                     caAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCg                                                      >  1:759355/1‑139 (MQ=255)                                                                       aaaTGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCa                                                    <  2:759355/139‑1 (MQ=255)                                                                           gggCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCggg                                                               <  1:646653/124‑1 (MQ=255)                                                                           gggCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCggg                                                               >  2:646653/1‑124 (MQ=255)                                                                                             agCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGt                              >  2:118592/1‑139 (MQ=255)                                                                                             agCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGt                              >  2:37452/1‑139 (MQ=255)                                                                                             agCTATTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGt                              >  1:34464/1‑139 (MQ=255)                                                                                                 aTTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTa                          <  1:118592/139‑1 (MQ=255)                                                                                                 aTTCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTa                          <  1:37452/139‑1 (MQ=255)                                                                                                  ttCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAg                         <  1:733767/139‑1 (MQ=255)                                                                                                  ttCGCGTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAg                         <  1:863850/139‑1 (MQ=255)                                                                                                     gcgTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCaa                                  <  1:74780/127‑1 (MQ=255)                                                                                                     gcgTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCaa                                  >  2:74780/1‑127 (MQ=255)                                                                                                      cgTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGt                     <  1:714148/139‑1 (MQ=255)                                                                                                      cgTGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGt                     <  2:397107/139‑1 (MQ=255)                                                                                                        tGCAGCGTGTAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGTAt                   >  2:72232/1‑139 (MQ=255)                                                                                                               tgtAAATGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGTATCGACAat            >  2:207051/1‑139 (MQ=255)                                                                                                                     tGCTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGTATCGACCATATCGaa      >  2:59228/1‑139 (MQ=255)                                                                                                                       cTGATGGCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTAACTGGGATTTCAAATGTTTTAGCGGTATCGACCATATCGaaaa    >  2:614236/1‑139 (MQ=255)                                                                                                                         gatggCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTaca                                                                                   <  2:200383/58‑1 (MQ=255)                                                                                                                         gatggCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTaca                                                                                   >  1:200383/1‑58 (MQ=255)                                                                                                                         gatggCCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGTATCGACCATATCGAAAAcc  >  2:757502/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                       CTGCTGGCCGCCATTGATGCTCTGAAACGTAATGACATTGCGGTGCTGTTGGATGTGGTAGTCAACCACAAAATGGGCGCGGATGAAAAAGAAGCTATTCGCGTGCAGCGTGTAAATGCTGATGACCGTACGCAAATTGACGAAGAAATCATTGAGTGTGAAGGCTGGACGCGTTACACCTTCCCCGCCCGTGCCGGGCAATACTCGCAGTTTATCTGGGATTTCAAATGTTTTAGCGGTATCGACCATATCGAAAACC  >  NC_000913/2006402‑2006660

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