Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,142,959 A→G 100% N468D (AAC→GAC)  katE → hydroperoxidase HPII(III)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,142,9590AG76.7% 50.5 / 13.2 30N468D (AAC→GAC) katEhydroperoxidase HPII(III)
Reads supporting (aligned to +/- strand):  ref base A (7/0);  new base G (23/0);  total (30/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAA  >  minE/1142925‑1142991
                                  |                                
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGTAACACCGCCGGGGCCGaaa  >  1:1455452/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1806999/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:633132/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:586523/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:657285/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:47859/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:473359/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:345896/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:2217429/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1688795/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1647876/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1591050/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:678667/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1547638/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1442451/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:716810/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:732708/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1180730/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1067191/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaa   >  1:1815373/1‑66 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaa   >  1:1927239/1‑66 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaa   >  1:280035/1‑66 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaa   >  1:497073/1‑66 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGtaa  >  1:853553/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:489033/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1590129/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:1375875/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:734955/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGaaa  >  1:108069/1‑67 (MQ=255)
cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGa    >  1:1340196/1‑65 (MQ=255)
                                  |                                
CCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAA  >  minE/1142925‑1142991

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: