Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,000,186	G→A	G179D (GGT→GAT)	ydgD →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,000,186	0	G	A	100.0%	186.3 / NA	57	G179D (GGT→GAT)	ydgD	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/26);  total (31/26)

ACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAA  >  minE/1000125‑1000247                                                              |                                                              aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                      >  1:388673/1‑71 (MQ=255) aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                      >  1:1537500/1‑71 (MQ=255) aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                      >  1:3450826/1‑71 (MQ=255) aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGt                                                      >  1:2531906/1‑71 (MQ=255)               aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:1918303/1‑70 (MQ=255)               aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:2580548/1‑70 (MQ=255)               aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:2132398/1‑70 (MQ=255)               aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:71145/1‑70 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:1031657/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:4003007/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:3671314/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:2883568/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:2881248/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:2801653/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:2799585/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:2682028/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:210338/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:1513627/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:1670531/1‑58 (MQ=255)                  tGAGGGAGATAAAGCCGCGCTTAACGCCGCATTAAAAGCGGCAGATCGTAAAGTGACt                                                  >  1:1982174/1‑58 (MQ=39)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:2181638/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:3983216/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:4148774/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:3854385/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:4255034/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:4307128/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:3478250/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:3332448/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:1906701/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:859247/52‑1 (MQ=255)                      ggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGa                                                    <  1:627512/52‑1 (MQ=255)                       gagaTAAAGCCGCGCTTAGCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:4134085/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:3414378/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:1565832/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:2777437/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:260152/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:3809238/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:769555/71‑1 (MQ=255)                       gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                <  1:1411050/71‑1 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:4377926/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:995662/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:743174/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:1071545/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:3937921/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:3885694/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:3627978/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:2206814/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCaggcaggc                                         >  1:2117233/1‑60 (MQ=255)                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAATGACTCaggcaggc                                         >  1:2284429/1‑60 (MQ=255)                          aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTg                                  >  1:693554/1‑66 (MQ=255)                           tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGat                             <  1:3005123/70‑1 (MQ=255)                                               cATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt        <  1:3485159/71‑1 (MQ=255)                                               cATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt        <  1:2102999/71‑1 (MQ=255)                                               cATTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt        <  1:1989816/71‑1 (MQ=255)                                                aTTAAAAGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt        <  1:2267585/70‑1 (MQ=255)                                                     aaGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt        <  1:4109204/65‑1 (MQ=255)                                                     aaGCGGCAGATCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCaa  <  1:636156/71‑1 (MQ=255)                                                              |                                                              ACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAA  >  minE/1000125‑1000247

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