| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,408,020 | A→G | 35.5% | G58G (GGA→GGG) | hemC → | hydroxymethylbilane synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,408,020 | 0 | A | G | 35.5% | 35.6 / 33.0 | 59 | G58G (GGA→GGG) | hemC | hydroxymethylbilane synthase |
| Reads supporting (aligned to +/- strand): ref base A (22/16); new base G (15/6); total (37/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.02e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACCGATGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCG > minE/2407963‑2408081 | aCCGATGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTAtt > 1:904420/1‑69 (MQ=255) tGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTaa < 1:1491838/69‑1 (MQ=255) ggTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGt > 1:2104766/1‑66 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:1195822/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:1705742/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:2419660/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:3394564/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGa > 1:1684217/1‑68 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGa > 1:515762/1‑68 (MQ=255) gtgATTCTTGATACGCCGGTGGCGAAAGTAGGCGGa > 1:186892/1‑36 (MQ=255) gtgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTcgcg > 1:2593899/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTgaa > 1:2321172/1‑52 (MQ=39) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTcgcg > 1:1448212/1‑68 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2986871/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2017387/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:3131207/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2414843/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:450020/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2084233/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:651920/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:2071058/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:3174039/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:1537272/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:1440814/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:2291699/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:3155138/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:1160376/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:578662/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:791033/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:958998/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:1803238/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:1911672/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGAGctc > 1:3343124/1‑69 (MQ=255) tCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctcc > 1:1575858/1‑68 (MQ=255) tCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctcc > 1:1049912/1‑68 (MQ=255) cTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctcct > 1:2151/1‑68 (MQ=255) cTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctcct > 1:1313100/1‑68 (MQ=255) cTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctcct > 1:1280522/1‑68 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGATAAGGCTTATTTGTAAAAGAGCTTGAAGTCGCGCTCCTCGaaaa > 1:1975214/1‑69 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGaa < 1:950039/67‑1 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGaa < 1:1890085/67‑1 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGaa < 1:2393160/67‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:940017/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:193596/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:175819/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:369163/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:1488376/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:1347287/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:1294795/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:781452/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:2127417/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:1557467/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:1645925/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:2872492/69‑1 (MQ=255) cgccgcTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:2964727/68‑1 (MQ=255) cTGGCTAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGcc < 1:1350519/69‑1 (MQ=255) cTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGcc < 1:3261425/69‑1 (MQ=255) tAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCg < 1:464260/69‑1 (MQ=255) tAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCg < 1:2930602/69‑1 (MQ=255) tAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCg < 1:2926709/69‑1 (MQ=255) | ACCGATGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCG > minE/2407963‑2408081 |
| 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: ‑ |