| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| MC JC | NC_000913 | 4,583,432 | Δ3,267 bp | coding (64‑3330/3513 nt) | hsdR ← | type I restriction enzyme EcoKI endonuclease subunit |
| Missing coverage evidence... | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| seq id | start | end | size | ←reads | reads→ | gene | description | |||
| * | * | ÷ | NC_000913 | 4583432 | 4586698 | 3267 | 13 [0] | [0] 13 | hsdR | type I restriction enzyme EcoKI endonuclease subunit |
| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 4583431 | 0 (0.000) | 13 (0.780) | 13/292 | 0.2 | 100% | coding (3331/3513 nt) | hsdR | type I restriction enzyme EcoKI endonuclease subunit |
| ? | NC_000913 | 4586699 = | 0 (0.000) | coding (63/3513 nt) | hsdR | type I restriction enzyme EcoKI endonuclease subunit | |||||
GGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTCCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACT > NC_000913/4583195‑4583628 | ggTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcata > 1:102248‑M1/1‑237 (MQ=255) gtCAGGCCAGCTCGTCCCAGATATAAGCGCTGAATTTGCCCAGCAGGGTATCGAGATGATCTTCAAAGGTTCTGTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGATCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggc < 1:247438‑M1/187‑4 (MQ=255) cAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGCTCCCGGTTTTGAAGACATCGTCGTCGAGCACCAATTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcacttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtg > 1:407590‑M1/1‑167 (MQ=255) cccGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggcctgtccaaatatttaaggcccataacatc > 1:207023‑M1/1‑96 (MQ=255) gccgcGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGGCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaaggcgttgacgcccttcagtaattcaaaattggatttattcatcattgttagtaatccattgctgtgcgggcctgttcaaatatttaaggcccataacatctcatcttagctttctgtacctttccgggcaatgac < 1:279182‑M1/246‑154 (MQ=255) ggAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCCCTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggacttcaaaattggctttattcatcattgttattaatccattgctgtgcgggcctgtccaaatatttaaggcccataacatctcatcttagctttctgtacctttccgggcaac > 1:67773‑M1/1‑82 (MQ=255) gtcgtcgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggcctgtccaaatatgtaaggcccataacatctcatcttagc < 1:143007‑M1/185‑129 (MQ=255) gtcgAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtagttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggcctgcccaaatatttaatgcccataacatctcatcttagctttctgtacctttccgggcactgaccacggtcacagcaactgactcatttctaacgtgttcgtctattt > 1:300586‑M1/1‑54 (MQ=255) tctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggnnttattcnnc > 2:25385‑M1/1‑38 (MQ=255) tctTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggcctgtccaaatatttaaggcccataacatctcatcttagctttctgtac < 1:25385‑M1/175‑138 (MQ=255) cGATCGGNCCAGCTTAAggcgatggcataagtgaagtcgntnacgcccttcaggaatgcaaaattggatttattcatcattgttattaatncagtgctgtgngggnntgtcnananatttaaggcccataacatctcatcttagcnttctgtacctttccgggcaatgannacnnncanngcaac < 2:317222‑M1/185‑169 (MQ=255) cGATCGAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggactgtccaaatatttaaggcccataacatctcctcttagctttctgtacctttccgggcaatgaccacggtcacagcaac > 1:317222‑M1/1‑17 (MQ=255) gAGCCAGCTTAAggcgatggcataagtgaagtcgttgacgcccttcaggaattcaaaattggatttattcatcattgttattaatccattgctgtgcgggcctgtccaaatatttaaggcccataacatctcatcttagctttctgtacctttccgggcaatgaccac > 1:3899‑M1/1‑12 (MQ=255) | GGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAAGACATCGTCGTCGAGCACCACTTTCTCTTTCAGCGCCTGCGCTAAACGATCGAGCCAGCTTAATTGCTCGCTGCTCCAGTCGTTTTCGCCCTTAATGCGCGTCAGCGCGTGATCGACACGTTCCTCAAACGGTTTCAGCGCATCGCCCACCGCAGCGCGGCGAATATGACCAATCAGCCGGGCGGCGATATCTTCATTGCGCGTCTCTTTCCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAACCACT > NC_000913/4583195‑4583628 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
| Reads not counted as support for junction |
| read_name Not counted due to insufficient overlap past the breakpoint. |
| read_name Not counted due to not crossing MOB target site duplication. |