| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,902,939 | C→T | 41.8% | A34A (GCG→GCA) | recB ← | exonuclease V (RecBCD complex), beta subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,902,939 | 0 | C | T | 41.8% | 33.5 / 61.0 | 86 | A34A (GCG→GCA) | recB | exonuclease V (RecBCD complex), beta subunit |
| Reads supporting (aligned to +/- strand): ref base C (32/18); new base T (26/10); total (58/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01 | |||||||||||
GGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAG > minE/1902884‑1902999 | ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:988879/69‑1 (MQ=255)ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:78656/69‑1 (MQ=255)ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:553775/69‑1 (MQ=255)ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:4962/69‑1 (MQ=255)ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:1375592/69‑1 (MQ=255)ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:2791531/69‑1 (MQ=255)ggcgggGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAgg < 1:2509135/69‑1 (MQ=255) ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc > 1:424370/1‑66 (MQ=255) ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc > 1:343303/1‑66 (MQ=255) ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc > 1:3356915/1‑66 (MQ=255) ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc > 1:3327195/1‑66 (MQ=255) ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc > 1:3003246/1‑66 (MQ=255) ggAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGcc > 1:199557/1‑66 (MQ=255) aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:2405765/1‑69 (MQ=255) aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:2476041/1‑69 (MQ=255) aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:2347071/1‑69 (MQ=255) aCCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:1055290/1‑69 (MQ=255) aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:2511078/1‑69 (MQ=255) aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:1465924/1‑69 (MQ=255) aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:1718677/1‑69 (MQ=255) aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:1783011/1‑69 (MQ=255) aCCGCCTAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:311617/1‑69 (MQ=255) aCCGCCGAGCCCAAGTAACAGGCGCAAATAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga > 1:436901/1‑69 (MQ=255) cGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCagag < 1:2555926/68‑1 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATc > 1:2200833/1‑68 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:2942689/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:901665/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:2756831/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:931329/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:1037088/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:2500947/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:236655/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:164080/1‑69 (MQ=255) cAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCa > 1:2157157/1‑69 (MQ=255) aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTc < 1:379912/62‑1 (MQ=255) aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTc < 1:312606/62‑1 (MQ=255) aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTc < 1:2557048/62‑1 (MQ=255) aGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGATTc < 1:2649884/62‑1 (MQ=255) tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg < 1:2090219/69‑1 (MQ=255) tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg < 1:1952624/69‑1 (MQ=255) tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg < 1:645660/69‑1 (MQ=255) tAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGgcg < 1:2483563/69‑1 (MQ=255) ggCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAc < 1:1527602/69‑1 (MQ=255) ggCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAc < 1:2012254/69‑1 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:3378495/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:3236411/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:1161508/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:3017733/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:1597890/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:2036836/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCAccc > 1:979653/1‑69 (MQ=255) cgcAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCCGCAGAGGCTTCAATCAGGCGCTCAccc > 1:3185470/1‑69 (MQ=255) cgcAAATAGAGAGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAAt > 1:1435927/1‑56 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:1474722/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:1289084/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:366484/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:2228451/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:2268229/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:1706986/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:2866754/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:624881/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCaga < 1:1812502/42‑1 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:906316/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:192487/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:166284/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1209308/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1613812/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1229520/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:627866/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1944118/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:3338341/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1396920/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:3302193/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:3301179/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:3117603/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:3077913/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1975495/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:1035632/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:2402818/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:2259970/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTaa > 1:2172295/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGATTCAATCAGGCGCTCACCCTGTaa > 1:2122219/1‑69 (MQ=255) aTAGAGCGCTGCAATCGTAAAGGTCTTGCCTGTGCCGGCaga < 1:148723/42‑1 (MQ=255) tAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAg > 1:1209134/1‑69 (MQ=255) tAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAg > 1:2043541/1‑69 (MQ=255) tAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAg > 1:222217/1‑69 (MQ=255) | GGCGGGGAAAGGCGGCGGAACCGCCTAGTCCAAGTAACAGGCGCAAATAGAGCGCCGCAATCGTAAAGGTTTTGCCTGTGCCGGCAGAGGCTTCAATCAGGCGCTCACCCTGTAAG > minE/1902884‑1902999 |
| 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: ‑ |