Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 4,581,599 | T→C | 31.4% | 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 | 31.4% | 87.2 / 61.2 | 70 | K348E (AAA→GAA) | mcrC | 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
Reads supporting (aligned to +/- strand): ref base T (13/35); new base C (8/14); total (21/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01 |
ATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTA > W3110S.gb/4581537‑4581669 | atgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGata < 1:1832328/69‑1 (MQ=255) atgatgaGCCGCTCCCGATGTGGTGGCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATAtt < 1:2100593/71‑1 (MQ=255) gaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATc > 1:1616497/1‑70 (MQ=255) gCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGaa < 1:3465209/71‑1 (MQ=255) gCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGaa < 1:3281528/71‑1 (MQ=255) ccGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa < 1:4031978/71‑1 (MQ=255) ccGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc < 1:3725947/70‑1 (MQ=255) ccGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc < 1:3632505/70‑1 (MQ=255) ccGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCg < 1:2431745/71‑1 (MQ=255) tgtgGTGTCGGGAGCGGTATTTTCTATAAAAGTTACCGCTTATTTGAGATATTCATCGaa < 1:2815548/60‑1 (MQ=255) gtgtCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTc > 1:2304342/1‑71 (MQ=255) tgtCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc < 1:1569898/61‑1 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaa < 1:569567/53‑1 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa > 1:288049/1‑61 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa > 1:1009361/1‑61 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTg < 1:2463106/71‑1 (MQ=255) tCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTg < 1:3236476/71‑1 (MQ=255) gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa < 1:511199/70‑1 (MQ=255) gggAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa < 1:1296927/70‑1 (MQ=255) cGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:1833107/63‑1 (MQ=255) cGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:4040175/63‑1 (MQ=255) ggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata < 1:1481002/70‑1 (MQ=255) ggTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta < 1:3866010/68‑1 (MQ=255) gTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGt > 1:95414/1‑49 (MQ=255) gTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATac < 1:1179745/70‑1 (MQ=255) tATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta < 1:3837285/66‑1 (MQ=255) tATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta < 1:710170/66‑1 (MQ=255) tATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta < 1:2138538/66‑1 (MQ=255) aTTTTCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACg > 1:3410596/1‑71 (MQ=255) ttttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:2244757/51‑1 (MQ=255) ttttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt < 1:492389/51‑1 (MQ=255) ttttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTaa > 1:808010/1‑53 (MQ=255) ttttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACgg > 1:3396074/1‑71 (MQ=255) ttttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACgg > 1:171393/1‑71 (MQ=255) ttttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACg < 1:868288/70‑1 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCg > 1:386401/1‑37 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:1934040/1‑70 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc < 1:3064705/70‑1 (MQ=255) ttCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc < 1:3870050/70‑1 (MQ=255) ttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:3605459/1‑70 (MQ=255) ttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:1485526/1‑70 (MQ=255) ttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:3855067/1‑70 (MQ=255) ttCTATAAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc > 1:1717611/1‑70 (MQ=255) tCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGcc < 1:3266342/70‑1 (MQ=255) tCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGcc < 1:2678560/70‑1 (MQ=255) tCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTAGAGTAATTCTTGa > 1:1994368/1‑57 (MQ=255) tataAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta < 1:3861013/59‑1 (MQ=255) tataAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta < 1:2140007/59‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTaa < 1:888422/46‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTaa < 1:2985917/46‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTaa < 1:2296835/46‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTaa < 1:1363466/46‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:3065827/51‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:364809/51‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAt < 1:1657415/68‑1 (MQ=255) taAAACTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAt < 1:971310/68‑1 (MQ=255) aataCTTACGGCTTTTTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTc < 1:1674418/66‑1 (MQ=255) aaaCTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc > 1:1468420/1‑69 (MQ=255) cTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTaa > 1:796864/1‑41 (MQ=255) cTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGAcc > 1:2251016/1‑71 (MQ=255) cTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGAcc > 1:3577681/1‑71 (MQ=255) cTTACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGAcc > 1:992982/1‑71 (MQ=255) ttACCGCTTATTCGAGATATTCATCGAAAATGTCGAGTAATTCtt < 1:4061871/45‑1 (MQ=255) aCCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:3798584/70‑1 (MQ=255) aCCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:3070234/70‑1 (MQ=255) ccGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTaa < 1:1040763/70‑1 (MQ=255) ccGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:132907/69‑1 (MQ=255) gCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg < 1:841419/44‑1 (MQ=255) tATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:464091/64‑1 (MQ=255) tATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTa < 1:2058598/64‑1 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGa < 1:177159/37‑1 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg < 1:1157548/39‑1 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg < 1:1066291/39‑1 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg < 1:4063546/39‑1 (MQ=255) tttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACgg < 1:1203056/71‑1 (MQ=255) ttGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACgg < 1:3298041/70‑1 (MQ=255) tGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTa > 1:3746002/1‑71 (MQ=255) tGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTa > 1:437912/1‑71 (MQ=255) tGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTa > 1:12414/1‑71 (MQ=255) | ATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTA > W3110S.gb/4581537‑4581669 |
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: ‑ |