| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| MC JC | NC_000913 | 3,582,206 | Δ1,222 bp | 100% | IS1‑mediated | [yhhZ]–yrhA | [yhhZ], yrhA |
| Missing coverage evidence... | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| seq id | start | end | size | ←reads | reads→ | gene | description | |||
| * | * | ÷ | NC_000913 | 3582206 | 3584147–3583428 | 1223–1942 | 107 [5] | [60] 62 | [yhhZ]–[insB1] | [yhhZ], yrhA, insA, [insB1] |
| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 3582205 | 3 (0.030) | 91 (0.880) | 62/178 | NT | 96.8% | coding (343/1179 nt) | yhhZ | putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity |
| ? | NC_000913 | 3583428 = | NA (NA) | noncoding (1/768 nt) | IS1 | repeat region | |||||
TGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGC > NC_000913/3583364‑3583520 | tGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgatg > 2:787277/1‑100 (MQ=255) tACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtgt > 2:420321/1‑99 (MQ=255) gAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGa < 1:787275/100‑1 (MQ=255) aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAg < 2:2149022/100‑1 (MQ=255) gTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa > 1:811758/1‑100 (MQ=255) gATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGt > 1:1127909/1‑100 (MQ=255) aTCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTg < 2:1127911/100‑1 (MQ=255) tATGACATTAAAGGTGATGCTGCCAACTTACTGATTTACTGTATGATGGTGTTTTTGATGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGt < 1:1587030/100‑1 (MQ=255) cATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGc > 1:1436873/1‑100 (MQ=35) | TGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGC > NC_000913/3583364‑3583520 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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. |