Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 816,173 A→G 100% E82E (GAA→GAG hemA → glutamyl tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE816,1730AG100.0% 130.3 / NA 36E82E (GAA→GAGhemAglutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/0);  total (36/0)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCG  >  minE/816162‑816232
           |                                                           
gttAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCAATTTAATGCg  >  1:605551/2‑71 (MQ=255)
ctaaaTGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:488931/4‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAacga                       >  1:645477/1‑50 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGAcgc                    >  1:625169/1‑53 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:554495/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:372129/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:372910/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:404085/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:442904/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:4701/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:496147/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:525261/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:34031/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:601866/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:623748/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:636438/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:67919/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:87560/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:15785/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:112010/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:119727/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:126985/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:128848/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:139904/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:141715/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:142131/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:146501/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:107829/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:203604/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:22352/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:241315/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:2612/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:264964/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:296748/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:302676/1‑71 (MQ=255)
cTTAATGAAGAGGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCg  >  1:319923/1‑71 (MQ=255)
           |                                                           
CTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCG  >  minE/816162‑816232

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