| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,770,511 | A→G | 71.1% | E27E (GAA→GAG) | pheA → | fused chorismate mutase P and prephenate dehydratase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,770,511 | 0 | A | G | 71.1% | 27.4 / 10.0 | 31 | E27E (GAA→GAG) | pheA | fused chorismate mutase P and prephenate dehydratase |
| Reads supporting (aligned to +/- strand): ref base A (2/7); new base G (17/5); total (19/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.18e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
ACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTG > minE/1770444‑1770578 | aCCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAAc < 1:2724945/69‑1 (MQ=255) cccGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCATTACTGGCAGAGCg < 1:3328846/69‑1 (MQ=255) cccGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCATTACTGGCAGAGCg < 1:2350360/69‑1 (MQ=255) cccGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCATTACTGGCAGAGCg < 1:523078/69‑1 (MQ=255) cccGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCATTACTGGCAGAGCg < 1:2883202/69‑1 (MQ=255) cccGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCATTACTGGCAGAGCg < 1:2189013/69‑1 (MQ=255) tACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGcgcg < 1:776118/69‑1 (MQ=255) tACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGcgcg < 1:366870/69‑1 (MQ=255) tACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGcgcg < 1:575352/69‑1 (MQ=255) gCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCg > 1:2417314/1‑68 (MQ=255) cgcTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCa < 1:2595173/69‑1 (MQ=255) aCTGTCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTAc > 1:976234/1‑68 (MQ=255) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAGGcc > 1:1404804/1‑38 (MQ=37) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGcc > 1:1621564/1‑38 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:2716281/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:2292537/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:3076294/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:3112178/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:3184839/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:3376304/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:405227/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:2333596/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:2210481/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:2104448/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:1376916/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:96327/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:122964/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:1134093/1‑39 (MQ=38) cAGAGCGACGCGAACTGGCCGTCGAGGTGGGAAAAGCCa > 1:1061600/1‑39 (MQ=38) aGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGat < 1:1078487/67‑1 (MQ=255) aGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGat < 1:673625/67‑1 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCaaa > 1:782620/1‑37 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGctct > 1:2474679/1‑44 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCg > 1:707473/1‑46 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:2628899/1‑68 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:230358/1‑68 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:1974228/1‑68 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:1374229/1‑68 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:1302038/1‑68 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:1301422/1‑68 (MQ=255) aCGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTg > 1:1290528/1‑68 (MQ=255) | ACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAAAAGCCAAACTGCTCTCGCATCGCCCGGTACGTGATATTG > minE/1770444‑1770578 |
| 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: ‑ |