| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,142,959 | A→G | 43.2% | 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 | 43.2% | 15.1 / 55.6 | 44 | N468D (AAC→GAC) | katE | hydroperoxidase HPII(III) |
| Reads supporting (aligned to +/- strand): ref base A (16/9); new base G (14/5); total (30/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.34e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 | |||||||||||
ATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAG > minE/1142899‑1143018 | atgcatCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCcgcg > 1:279439/1‑71 (MQ=255) tgcatCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCcgcgc > 1:1855293/1‑71 (MQ=255) gatCTCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:1735991/70‑1 (MQ=255) catCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCcgcgcg > 1:1959108/1‑70 (MQ=255) cGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTTGCCGCGCGaa > 1:1616771/1‑69 (MQ=255) cGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGGGCGaaa > 1:1443208/1‑70 (MQ=255) cGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:3835760/1‑70 (MQ=255) cGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:996557/1‑70 (MQ=255) cGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:3292148/1‑70 (MQ=255) gCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAaca > 1:3719756/1‑71 (MQ=255) aTGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:2087644/65‑1 (MQ=255) aTGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:3141700/65‑1 (MQ=255) aTGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:2820918/65‑1 (MQ=255) tGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:1898199/1‑66 (MQ=255) tGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:3610754/1‑66 (MQ=255) tGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:3464603/1‑66 (MQ=255) tGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:3258620/1‑66 (MQ=255) gggATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACAccgcc > 1:308875/1‑71 (MQ=255) gggATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACAccgcc > 1:3807246/1‑71 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCTGCGCGaaa > 1:423518/1‑59 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGaaa > 1:11561/1‑59 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGaaa > 1:2979903/1‑59 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGaaa > 1:722689/1‑59 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGaaa > 1:2611639/1‑59 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGaaa > 1:4087210/1‑59 (MQ=255) cGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGc > 1:2858470/1‑71 (MQ=255) tAACCCGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAAc > 1:757722/1‑71 (MQ=255) tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAac < 1:70716/54‑1 (MQ=255) aaCCCGGCGAATTACGATCCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:4006172/1‑52 (MQ=255) aaCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:2018042/1‑52 (MQ=255) aaCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa > 1:2424059/1‑52 (MQ=255) cccGGCGAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAAcgcg > 1:624494/1‑71 (MQ=255) cGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGa < 1:520192/46‑1 (MQ=255) gAATTACGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGaa > 1:3098117/1‑60 (MQ=255) aCGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa < 1:1586085/70‑1 (MQ=255) aCGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa < 1:3661704/70‑1 (MQ=255) aCGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa < 1:430379/70‑1 (MQ=255) aCGAACCGAACTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa < 1:1681435/70‑1 (MQ=255) cgaaccgaacTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa < 1:1119539/69‑1 (MQ=255) gaaccgaacTCGATTAACGATGACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt > 1:1578014/1‑70 (MQ=255) ccgaacTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATa < 1:2064891/71‑1 (MQ=255) cgaacTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGTCCGAAACGCGGCGGTTTTGAATCATa < 1:4084596/70‑1 (MQ=255) cGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAg > 1:1578605/1‑71 (MQ=255) gATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATa < 1:2655630/56‑1 (MQ=255) | ATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAG > minE/1142899‑1143018 |
| 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: ‑ |