| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,142,959 | A→G | 19.0% | N468D (AAC→GAC) | katE → | hydroperoxidase HPII(III) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,142,959 | 0 | A | G | 19.0% | 41.3 / 6.1 | 21 | N468D (AAC→GAC) | katE | hydroperoxidase HPII(III) |
| Reads supporting (aligned to +/- strand): ref base A (4/13); new base G (2/2); total (6/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.44e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01 | |||||||||||
CGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAG > minE/1142890‑1143015 | cGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATaa < 1:1385050/71‑1 (MQ=255) cGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGc < 1:1396483/71‑1 (MQ=255) tCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa < 1:2172401/71‑1 (MQ=255) gCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAaca > 1:1020409/1‑71 (MQ=255) gATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCgg < 1:1333980/71‑1 (MQ=255) acTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:811026/53‑1 (MQ=255) acTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:254150/53‑1 (MQ=255) acTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:2117118/53‑1 (MQ=255) aaCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAAcg > 1:1439503/1‑71 (MQ=255) aaCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAAcg > 1:2167609/1‑71 (MQ=255) aaCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAAcg > 1:1751810/1‑71 (MQ=255) cGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCgggg < 1:981412/59‑1 (MQ=255) gCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGcggcgg < 1:1665592/71‑1 (MQ=255) cGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGcggcgg < 1:1668126/70‑1 (MQ=255) aTTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGa < 1:1223744/57‑1 (MQ=255) ttACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGtttt > 1:2067788/1‑70 (MQ=255) tACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:1123450/38‑1 (MQ=255) aCGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa < 1:443235/70‑1 (MQ=255) cgaaccgaacTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaaa > 1:2061887/1‑55 (MQ=255) cgaaccgaacTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt < 1:2170823/71‑1 (MQ=255) tCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGcggcgg < 1:1206739/53‑1 (MQ=255) aTAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAg < 1:1132301/59‑1 (MQ=255) | CGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAG > minE/1142890‑1143015 |
| 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: ‑ |