| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,583,739 | A→G | 12.5% | F94S (TTC→TCC) | mcrB ← | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,583,739 | 0 | A | G | 12.5% | 81.3 / 5.1 | 32 | F94S (TTC→TCC) | mcrB | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB |
| Reads supporting (aligned to +/- strand): ref base A (8/20); new base G (2/2); total (10/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.72e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01 | |||||||||||
TTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAA > W3110S.gb/4583674‑4583804 | ttttAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGc > 1:94003/1‑71 (MQ=255) tttAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGc < 1:653541/70‑1 (MQ=255) tAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGc < 1:136944/68‑1 (MQ=255) ggATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCAt < 1:1759221/69‑1 (MQ=255) taCACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCa < 1:997820/70‑1 (MQ=255) acacCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGGACTGCCATTGGGCAtg > 1:1482138/1‑71 (MQ=255) cccGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTgg < 1:1357627/71‑1 (MQ=255) cGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGtt > 1:1714010/1‑71 (MQ=255) cGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGtt > 1:651631/1‑71 (MQ=255) aGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCAt < 1:951812/71‑1 (MQ=255) aGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCAt > 1:1766673/1‑71 (MQ=255) gTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCAtt < 1:815148/71‑1 (MQ=255) gCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAgaga > 1:473582/1‑43 (MQ=255) gCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGaa > 1:145671/1‑44 (MQ=255) ataaTACTCTGCGAGAGTGTTAGGTATGTCTGAAGAGAACTGCCa < 1:1492411/43‑1 (MQ=255) aaaaTACTCTTCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCgtg < 1:1478742/67‑1 (MQ=255) aaaaTACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCgtg < 1:1645767/67‑1 (MQ=255) aaaaTACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCgtg < 1:958099/67‑1 (MQ=255) aTACTCTGCGATTGTTTTAGGTATGTCTGAAgaga < 1:805300/35‑1 (MQ=255) ctctGCGATTGTTTTAGGTATGTCTGAAGAGGACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTtata < 1:1294736/70‑1 (MQ=255) tGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCa < 1:509966/52‑1 (MQ=255) cGATTGTTTTAGGTATGTCTGAAGAGGACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATa < 1:646902/70‑1 (MQ=255) cGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGt > 1:1198890/1‑47 (MQ=255) cGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATa < 1:333249/70‑1 (MQ=255) aTTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGc < 1:1653181/38‑1 (MQ=255) tGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTc < 1:1269831/45‑1 (MQ=255) tGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTc < 1:336880/45‑1 (MQ=255) tGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTc < 1:1271812/45‑1 (MQ=255) gTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCAtg > 1:785552/1‑38 (MQ=255) gTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCAtg > 1:80712/1‑38 (MQ=255) gTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCa < 1:1638188/70‑1 (MQ=255) ttttACGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCa < 1:703713/44‑1 (MQ=255) tGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACtcatca > 1:13855/1‑70 (MQ=255) aaGAGGACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCaaaa > 1:1654334/1‑71 (MQ=255) aaGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCaa < 1:661027/69‑1 (MQ=255) agaACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATAc > 1:179012/1‑42 (MQ=255) | TTTTAGGATATACACCCGAAGTTGCCTGAAAATACTCTGCGATTGTTTTAGGTATGTCTGAAGAGAACTGCCATTGGGCATGTGGTTCATTCGTGTCGCTTATACCATAAGCCAAAACCAACTCATCAAAA > W3110S.gb/4583674‑4583804 |
| 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: ‑ |