Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 4,581,599 | T→C | 10.9% | K348E (AAA→GAA) | mcrC ← | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 4,581,599 | 0 | T | C | 10.9% | 143.1 / 8.6 | 55 | K348E (AAA→GAA) | mcrC | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
Reads supporting (aligned to +/- strand): ref base T (23/26); new base C (2/4); total (25/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTA > W3110S.gb/4581529‑4581669 | gCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt > 1:3756323/1‑71 (MQ=255) tcctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAttt > 1:694125/1‑68 (MQ=255) ctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgaga > 1:936136/1‑70 (MQ=255) ctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgaga > 1:2275287/1‑70 (MQ=255) atgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGatat > 1:2426245/1‑70 (MQ=255) atgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATAtt > 1:2108631/1‑71 (MQ=255) gatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa > 1:4333961/1‑71 (MQ=255) ccGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa < 1:776889/71‑1 (MQ=255) ccGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa < 1:4388709/71‑1 (MQ=255) gCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAAt < 1:4116816/71‑1 (MQ=255) cccGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTga > 1:3438904/1‑52 (MQ=255) gATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATAtt > 1:945778/1‑55 (MQ=255) aTGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:3752367/71‑1 (MQ=255) aTGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:1776614/71‑1 (MQ=255) aTGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:622682/71‑1 (MQ=255) tgtgGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTc < 1:4494239/54‑1 (MQ=255) tggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTaa < 1:2250166/70‑1 (MQ=255) ggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa > 1:2271852/1‑52 (MQ=255) ggtgTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa > 1:778583/1‑52 (MQ=255) gtCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCtt > 1:3540215/1‑71 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGGAAATGTCGAGTAATTCtt < 1:3281450/70‑1 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:2204341/70‑1 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:649048/70‑1 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:1213369/70‑1 (MQ=255) gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTg > 1:494803/1‑69 (MQ=255) gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt < 1:4550045/71‑1 (MQ=255) gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt > 1:471925/1‑71 (MQ=255) ggAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg < 1:1651094/71‑1 (MQ=255) gAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg > 1:2595198/1‑70 (MQ=255) gAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGt < 1:726730/71‑1 (MQ=255) ggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata < 1:1098585/70‑1 (MQ=255) gTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt < 1:4284900/64‑1 (MQ=255) gTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt < 1:3374154/64‑1 (MQ=255) aTTTTCTATAAAACTTACCGCTTATTTGAGATAtt > 1:2860220/1‑35 (MQ=255) aTTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa > 1:1622165/1‑61 (MQ=255) aTTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACg > 1:397517/1‑71 (MQ=255) ttttCTATAAAACTTACCGCTTATTCGAGATATTCATCGaaaa > 1:3097301/1‑43 (MQ=255) tttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATg < 1:146590/44‑1 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:2362039/56‑1 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa > 1:4354395/1‑58 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata > 1:1183585/1‑64 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:1096210/1‑70 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATACATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:4419928/1‑70 (MQ=255) ttCTATAAAACTTACCGCTTATTCGAGATATTCATCGaa < 1:1804939/39‑1 (MQ=255) tCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCg > 1:3023943/1‑45 (MQ=255) cTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAt > 1:2252905/1‑71 (MQ=255) tataAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAtt > 1:285438/1‑71 (MQ=255) ataAAACTTACCGCTTATTTGAGATATTCATCGaa < 1:3321150/35‑1 (MQ=255) ataAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATaca < 1:556463/62‑1 (MQ=255) taAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTc < 1:381531/70‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTaa < 1:163904/46‑1 (MQ=255) aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa > 1:128944/1‑51 (MQ=255) aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTg < 1:2349331/71‑1 (MQ=255) aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTg < 1:4231436/71‑1 (MQ=255) aaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata < 1:2749528/56‑1 (MQ=255) aCCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:628796/70‑1 (MQ=255) cTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAAtt < 1:4313888/70‑1 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc > 1:419654/1‑55 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACgg > 1:160713/1‑71 (MQ=255) tGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTa > 1:4574057/1‑71 (MQ=255) | GCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTA > W3110S.gb/4581529‑4581669 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |