Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 782,109 T→C Q167R (CAG→CGG)  hlyE ← hemolysin E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE782,1090TC100.0% 121.6 / NA 34Q167R (CAG→CGG) hlyEhemolysin E
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (1/33);  total (1/33)

ACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACT  >  minE/782055‑782148
                                                      |                                       
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:134197/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:855061/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:706162/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:3203627/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:3162799/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:305061/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:3049445/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2878878/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2793347/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2751445/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2661640/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2630271/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2610430/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2608326/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2551816/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2381529/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1156249/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1534309/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1564377/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1608109/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1625095/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1647366/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1671332/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1920146/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1921523/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:1942587/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2147074/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2295328/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2299392/71‑1 (MQ=255)
acacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTCTGACCGGAAATAGCTGCtttt                         <  1:3140874/71‑1 (MQ=255)
 cacCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:580782/70‑1 (MQ=255)
          cGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2686987/61‑1 (MQ=255)
              aCCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCtttt                         <  1:2949256/57‑1 (MQ=255)
                       tGCTTCCTTCCTGATTTTATCTACCTGTGACCGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACt  >  1:886534/1‑71 (MQ=255)
                                                      |                                       
ACACCGGCTGCGGCACCGGCATATGCTTCCTTCCTGATTTTATCTACCTGTGACTGGAAATAGCTGCTTTTTTCTGAAAAATCATTGGTTAACT  >  minE/782055‑782148

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: