| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 2,960,805 | T→C | 37.5% | K94E (AAA→GAA) | recC ← | exonuclease V (RecBCD complex), gamma chain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,960,805 | 0 | T | C | 37.5% | 21.3 / 32.2 | 32 | K94E (AAA→GAA) | recC | exonuclease V (RecBCD complex), gamma chain |
| Reads supporting (aligned to +/- strand): ref base T (11/9); new base C (6/6); total (17/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01 | |||||||||||
GGGTAAAGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGAACATATCC > W3110S.gb/2960757‑2960875 | gggTAAAGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCtgtttgtt > 1:4045499/1‑70 (MQ=255)gggTAAAGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCtgtttgtt > 1:2218258/1‑70 (MQ=255)gggTAAAGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTa < 1:1992716/71‑1 (MQ=255)gggTAAAGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATAAGTTTCCAGCTCATGCTCTGTTTGTTa < 1:277005/71‑1 (MQ=255)gggTAAAGTCTTC‑CGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGc > 1:3201038/1‑59 (MQ=255) aaGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCa < 1:496003/47‑1 (MQ=255) gTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGgc < 1:1022136/69‑1 (MQ=255) cgcgCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGgc > 1:2049236/1‑64 (MQ=255) cgcgCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGctctct > 1:1145014/1‑71 (MQ=255) cgcgCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGctct < 1:986729/69‑1 (MQ=255) cTCCAGCAATTGCGGCAGCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCtt < 1:1393875/71‑1 (MQ=255) ccAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCa > 1:3356981/1‑39 (MQ=255) ccAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTggg > 1:1312148/1‑70 (MQ=255) ccAGCAATTGAG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCa > 1:2618745/1‑39 (MQ=255) gCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGctc < 1:1095878/41‑1 (MQ=255) gCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGAtt > 1:799334/1‑71 (MQ=255) ttGCG‑‑‑GCAGCAGAGTCATCAGTTCCCAGCTCATGCTCtgtttgtt < 1:1422648/45‑1 (MQ=255) tgag‑‑‑gCAGCAGAGTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCgg > 1:2495442/4‑71 (MQ=255) cG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGt < 1:1724109/70‑1 (MQ=255) cG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGt < 1:752765/70‑1 (MQ=255) cG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGt < 1:637558/70‑1 (MQ=255) cG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGt < 1:1579823/70‑1 (MQ=255) agagTCATCAGTTCCCAGCTCATGCTCTGTTTGTTaaa > 1:3058891/1‑38 (MQ=255) agagTCATCAGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCg > 1:2862931/1‑71 (MQ=255) tcatcaGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGt < 1:2261926/59‑1 (MQ=255) atcaGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTaa > 1:2059935/1‑59 (MQ=255) atcaGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTaa > 1:781052/1‑59 (MQ=255) atcaGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGaa > 1:3912955/1‑71 (MQ=255) atcaGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGaa < 1:2335441/71‑1 (MQ=255) atcaGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGaa < 1:310606/71‑1 (MQ=255) atcaGTTCCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGaa > 1:3607491/1‑71 (MQ=255) caGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCttt > 1:2600335/1‑42 (MQ=255) tttCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTg < 1:872369/40‑1 (MQ=255) tttCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGAACatat > 1:1800095/1‑71 (MQ=255) tCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGAt > 1:720380/1‑43 (MQ=255) tCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGAACATATcc > 1:1013971/1‑71 (MQ=255) | GGGTAAAGTCTTCGCGCTCCAGCAATTGCG‑‑‑GCAGCAGAGTCATCAGTTTCCAGCTCATGCTCTGTTTGTTAAAGGCGCTCTCTTTGGGGATTTCCGGTAACACCCGGACGAACATATCC > W3110S.gb/2960757‑2960875 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |