| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,087,357 | A→G | 47.2% | G286G (GGT→GGC) | sufS ← | selenocysteine lyase, PLP‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,087,357 | 0 | A | G | 47.2% | 1.5 / 60.1 | 70 | G286G (GGT→GGC) | sufS | selenocysteine lyase, PLP‑dependent |
| Reads supporting (aligned to +/- strand): ref base A (18/19); new base G (21/12); total (39/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.36e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01 | |||||||||||
TTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGT > minE/1087292‑1087422 | ttCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGAcc > 1:1971210/1‑68 (MQ=255)ttCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCa > 1:1247514/1‑69 (MQ=255)ttCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCa > 1:1991718/1‑69 (MQ=255) cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACTTACTCCAGCGCCGCGCCAAGACCAat < 1:2699208/69‑1 (MQ=255) cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat < 1:2063945/69‑1 (MQ=255) cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat < 1:2183758/69‑1 (MQ=255) cTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat < 1:1771701/69‑1 (MQ=255) cATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAATACCAATGATGc > 1:2288325/1‑69 (MQ=255) cTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc > 1:1260866/1‑69 (MQ=255) cTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGccccc > 1:1702736/1‑69 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:1845643/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:456891/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:2167880/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:2207943/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:2560684/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:187646/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:2581662/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:1736633/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:2885617/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:3171479/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:1300905/69‑1 (MQ=255) cGGCAATGTTGTTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCgg < 1:967305/69‑1 (MQ=255) cTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGt < 1:3201161/67‑1 (MQ=255) aTGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTg < 1:2471674/69‑1 (MQ=255) aTGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTg < 1:2641708/69‑1 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTGTTgg > 1:3261112/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:2022915/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:2127022/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:1981867/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:2841767/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:2365854/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:1905795/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:708081/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:1849093/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:546436/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:3327917/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTgg > 1:80228/1‑69 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTg > 1:1989057/1‑68 (MQ=255) tgttgtTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTg > 1:899167/1‑68 (MQ=255) gTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat < 1:2207896/51‑1 (MQ=255) ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat > 1:3374284/1‑50 (MQ=255) ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAat > 1:165908/1‑50 (MQ=255) ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTggg > 1:293727/1‑68 (MQ=255) ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTggg > 1:1599803/1‑68 (MQ=255) ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTgg < 1:2766339/67‑1 (MQ=255) ttattaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTgg < 1:198686/67‑1 (MQ=255) attaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgtgt < 1:1911358/69‑1 (MQ=255) attaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgtgt < 1:1255715/69‑1 (MQ=255) attaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGgtg < 1:292250/68‑1 (MQ=255) ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa > 1:363924/1‑69 (MQ=255) ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa > 1:3320990/1‑69 (MQ=255) ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa > 1:980712/1‑69 (MQ=255) ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa > 1:1780484/1‑69 (MQ=255) ttaAGCACCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTa > 1:442433/1‑69 (MQ=255) aGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCa > 1:412356/1‑69 (MQ=255) aGCGCCGAAACATACTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCa > 1:2101926/1‑69 (MQ=255) aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg > 1:102034/1‑67 (MQ=255) aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg > 1:3345539/1‑67 (MQ=255) aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg > 1:2612804/1‑67 (MQ=255) aCTCCAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg > 1:759155/1‑67 (MQ=255) aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTg < 1:1621249/69‑1 (MQ=255) aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTg < 1:3237951/69‑1 (MQ=255) aCTCAAGCGCCGCGCCAAGGCCAATAATACCTCCGGTATTGGGCGTACCGGCTTCAAACCGCCACggcg > 1:1223514/1‑67 (MQ=255) cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGTTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc > 1:2070141/1‑69 (MQ=255) cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTg < 1:2623709/68‑1 (MQ=255) cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc > 1:826535/1‑69 (MQ=255) cTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATT‑GGTGTACCGGCTTCAAACCGCCATGGTGCt > 1:188798/1‑69 (MQ=255) aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGt < 1:854177/68‑1 (MQ=255) aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGt < 1:522720/68‑1 (MQ=255) aGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGt < 1:2344288/68‑1 (MQ=255) | TTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGT > minE/1087292‑1087422 |
| 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: ‑ |