Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	120,663	C→T	C181Y (TGC→TAC)	speE ←	spermidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	120,663	0	C	T	100.0%	108.3 / NA	33	C181Y (TGC→TAC)	speE	spermidine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/14);  total (19/14)

TCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAG  >  minE/120603‑120720                                                             |                                                          tcttcCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTg                                                   >  1:3166267/1‑69 (MQ=255)         ctgTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCa                                          <  1:877364/70‑1 (MQ=255)         ctgTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCa                                          <  1:136223/70‑1 (MQ=255)          tgTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCAt                                         >  1:88561/1‑70 (MQ=255)                 aGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGcc                                 >  1:1274115/1‑71 (MQ=255)                  gCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGcc                                 >  1:866822/1‑70 (MQ=255)                  gCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGcc                                 >  1:2490086/1‑70 (MQ=255)                      aCGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCtt                                            <  1:1217665/55‑1 (MQ=255)                      aCGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCtt                                            <  1:1569093/55‑1 (MQ=255)                      aCGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCtt                                            <  1:1877313/55‑1 (MQ=255)                       cGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGt                            >  1:2659885/1‑70 (MQ=255)                       cGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGt                            >  1:1960615/1‑70 (MQ=255)                        gCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTg                           >  1:3120030/1‑70 (MQ=255)                          cGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGaa                         >  1:1615444/1‑70 (MQ=255)                           gTTTTGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGaa                         >  1:2821220/1‑69 (MQ=255)                              ttGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCTGAAg                             >  1:551901/1‑62 (MQ=255)                              ttGTGCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAg                             >  1:451699/1‑62 (MQ=255)                               tgtgCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGc                    <  1:1081826/70‑1 (MQ=255)                               tgtgCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGc                    <  1:730890/70‑1 (MQ=255)                                  gCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGcc                                 <  1:642729/54‑1 (MQ=255)                                  gCGACGAAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGcc                                 <  1:1651471/54‑1 (MQ=255)                                   cgacgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCg               >  1:2278036/1‑71 (MQ=255)                                    gacgaAGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGc              <  1:764461/71‑1 (MQ=255)                                         aGATACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGa         >  1:1215331/1‑71 (MQ=255)                                           aTACCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGAc        <  1:593957/70‑1 (MQ=255)                                             aCCGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCg      <  1:1477559/70‑1 (MQ=255)                                              ccGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGa     <  1:298519/70‑1 (MQ=255)                                               cGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGa     >  1:2217323/1‑69 (MQ=255)                                               cGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGa     >  1:181980/1‑69 (MQ=255)                                               cGCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGa     >  1:135109/1‑69 (MQ=255)                                                gCCAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAg  >  1:18367/1‑71 (MQ=255)                                                 ccAGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGAcc       <  1:1411906/65‑1 (MQ=255)                                                   aGGATTCAGGTAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCg               >  1:857716/1‑55 (MQ=255)                                                             |                                                          TCTTCCTGCTGTAAAAAGCAGACGCCGTTTTGTGCGACGAAGATACCGCCAGGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAG  >  minE/120603‑120720

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