| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,408,035 | A→G | 32.7% | V63V (GTA→GTG) | hemC → | hydroxymethylbilane synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,408,035 | 0 | A | G | 32.7% | 46.9 / 32.6 | 64 | V63V (GTA→GTG) | hemC | hydroxymethylbilane synthase |
| Reads supporting (aligned to +/- strand): ref base A (24/19); new base G (14/7); total (38/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGATGT > minE/2407968‑2408100 | tGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTaa < 1:1491838/69‑1 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:1195822/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:1705742/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:3394564/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGaa > 1:2419660/1‑69 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGa > 1:515762/1‑68 (MQ=255) ggcgATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGa > 1:1684217/1‑68 (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:2017387/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2084233/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2414843/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:2986871/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:3131207/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:450020/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCgcgc > 1:651920/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:1537272/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:1440814/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:3174039/1‑69 (MQ=255) tgATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCgcgc > 1:2071058/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:3155138/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:578662/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:791033/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:2291699/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:958998/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:1160376/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:1911672/1‑69 (MQ=255) aTTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctc > 1:1803238/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:1280522/1‑68 (MQ=255) cTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGctcct > 1:1313100/1‑68 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGATAAGGCTTATTTGTAAAAGAGCTTGAAGTCGCGCTCCTCGaaaa > 1:1975214/1‑69 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGaa < 1:1890085/67‑1 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGaa < 1:950039/67‑1 (MQ=255) tACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGaa < 1:2393160/67‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:193596/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:940017/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:1347287/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:369163/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:1488376/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGGAAAGGCTTATTTGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:175819/68‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:2127417/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:1557467/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:781452/69‑1 (MQ=255) aCGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAAt < 1:1294795/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:2926709/69‑1 (MQ=255) tAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCg < 1:2930602/69‑1 (MQ=255) ttATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCg < 1:1317956/55‑1 (MQ=255) ttATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGaaa < 1:2063858/69‑1 (MQ=255) ttATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGaaa < 1:2714988/69‑1 (MQ=255) aTTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCg < 1:727748/53‑1 (MQ=255) tttGTGAAAGAGCTGGAAGTCGCGCTCCTGGaaaa < 1:1356718/35‑1 (MQ=37) tGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGAtgt > 1:2925033/1‑69 (MQ=255) tGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGAtgt > 1:3131714/1‑69 (MQ=255) tGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGAtgt > 1:1733348/1‑69 (MQ=255) tGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGAtg > 1:1664303/1‑68 (MQ=255) | TGGTGACGCGCGGCGATGTGATTCTTGATACGCCGCTGGCGAAAGTAGGCGGAAAAGGCTTATTTGTAAAAGAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGATGT > minE/2407968‑2408100 |
| 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: ‑ |