Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,737,646	G→A	T50T (ACG→ACA)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,737,646	0	G	A	100.0%	167.7 / NA	47	T50T (ACG→ACA)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/22);  total (25/22)

GTTTCCCGCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACGCTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAACCTT  >  NC_000913/3737568‑3737719                                                                               |                                                                          gTTTCCCGCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACt                                                                         <  2:731111/81‑1 (MQ=255)  tttCCCGCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTa                                                                        <  2:731112/81‑1 (MQ=255)    tCCCGCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTaaa                                                                      >  1:731167/1‑81 (MQ=255)        gCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAAtttt                                                                  <  1:731114/81‑1 (MQ=255)        gCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAAtttt                                                                  <  1:731113/81‑1 (MQ=255)          cTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGa                                                                <  1:731115/81‑1 (MQ=255)          cTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGa                                                                >  1:731165/1‑81 (MQ=255)                 gAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAg                                                         >  1:731151/1‑81 (MQ=255)                     aGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCt                                                     >  2:731183/1‑81 (MQ=255)                     aGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCt                                                     >  2:731185/1‑81 (MQ=255)                      ggggACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTg                                                    <  1:731116/81‑1 (MQ=255)                             ggAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCt                                             <  2:731117/81‑1 (MQ=255)                              gAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTg                                            >  2:731194/1‑81 (MQ=255)                               aaCATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGc                                           >  1:731175/1‑81 (MQ=255)                                aCATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCg                                          <  1:731118/81‑1 (MQ=255)                                  aTTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGa                                        >  1:731181/1‑81 (MQ=255)                                   tttGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGAt                                       <  2:731119/81‑1 (MQ=255)                                      gTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGc                                     >  2:731187/1‑80 (MQ=255)                                        cAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGa                                   >  1:731189/1‑80 (MQ=255)                                         aGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGAt                                  >  1:731188/1‑80 (MQ=255)                                          gCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACAAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATaa                                <  2:731120/81‑1 (MQ=255)                                            cACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATaaaa                              >  2:731192/1‑81 (MQ=255)                                            cACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATaaaa                              <  1:731121/81‑1 (MQ=255)                                              cgcgCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACt                            <  2:731122/81‑1 (MQ=255)                                              cgcgCAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACt                            <  2:731123/81‑1 (MQ=255)                                               gcgcAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACt                            >  1:731196/1‑80 (MQ=255)                                                 gcAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCaa                         >  1:731260/1‑81 (MQ=255)                                                 gcAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCaa                         >  2:731195/1‑81 (MQ=255)                                                 gcAGTTGCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCaa                         >  2:731212/1‑81 (MQ=255)                                                       gCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGa                    <  2:731124/80‑1 (MQ=255)                                                       gCCCAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGa                    <  2:731125/80‑1 (MQ=255)                                                        cccAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATg                  <  1:731126/81‑1 (MQ=255)                                                        cccAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATg                  <  1:731127/81‑1 (MQ=255)                                                         ccAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGc                 >  1:731197/1‑81 (MQ=255)                                                         ccAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGc                 <  1:731128/81‑1 (MQ=255)                                                         ccAAAGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGc                 >  2:731198/1‑81 (MQ=255)                                                           aaaGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCtt               >  1:731209/1‑81 (MQ=255)                                                             aGGTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTc             <  2:731129/81‑1 (MQ=255)                                                               gTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAGCTCAATCAGATGCTTTccc           >  2:731190/1‑81 (MQ=255)                                                               gTACGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTccc           <  2:731130/81‑1 (MQ=255)                                                                 aCGCGTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTg         >  2:731207/1‑81 (MQ=255)                                                                  cgcgTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGc        <  2:731131/81‑1 (MQ=255)                                                                   gcgTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGc        >  1:731158/1‑80 (MQ=255)                                                                   gcgTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGc        >  1:731205/1‑80 (MQ=255)                                                                     gTCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAAc     >  2:731225/1‑81 (MQ=255)                                                                      tCCACTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAAcc    <  2:731132/81‑1 (MQ=255)                                                                        cactaacactaaATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAACCtt  <  1:731135/81‑1 (MQ=255)                                                                        cactaacactaaATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAACCtt  <  2:731133/81‑1 (MQ=255)                                                                        cactaacactaaATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAACCt   <  2:731134/80‑1 (MQ=255)                                                                          cTAACACTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAActtg  >  1:731204/1‑76 (MQ=255)                                                                               |                                                                          GTTTCCCGCCTTTAGCGAACAGGGGACAGGAACATTTGTCAGCCACGCGCAGTTGCCCAAAGGTACGCGTCCACTAACGCTAAATTTTGACCAACAGTGCTGGCAGCCTGCGGATGCGATAAAACTCAATCAGATGCTTTCCCTGCAACCTT  >  NC_000913/3737568‑3737719

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