Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
MC JC | REL606 | 3,894,997 | Δ6,460 bp | 100% | IS150‑mediated | rbsD–[yieO] | rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO] |
Missing coverage evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | start | end | size | ←reads | reads→ | gene | description | |||
* | * | ÷ | REL606 | 3893577–3894996 | 3901456 | 6461–7880 | 51 [46] | [0] 180 | insJ‑5–[yieO] | insJ‑5, insK‑5, rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO] |
New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | REL606 | = 3652533 | NA (NA) | 95 (0.960) | 38/46 | NT | 100% | noncoding (1443/1443 nt) | IS150 | repeat region |
? | REL606 | 3901457 = | 0 (0.000) | coding (963/1428 nt) | yieO | predicted multidrug or homocysteine efflux system |
TCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATG > REL606/3893541‑3893586 | tCTAACGCC‑GACTGTACTGCACCCATTTTGTTGGa > 1:7205251/1‑35 (MQ=17) cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa < 1:4759200/35‑1 (MQ=255) cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa < 1:9375308/35‑1 (MQ=255) cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa < 1:12577486/35‑1 (MQ=255) cTAACGCCAGACTGTACTGCACCCATTTTGTTGGa < 1:7564924/35‑1 (MQ=255) cTAACGCC‑GACTGTACTGCACCCATTTTGTTGGAc < 1:3646778/35‑1 (MQ=17) cTAACGCC‑GACTGTACTGCACCCATTTTGTTGGAc < 1:5101731/35‑1 (MQ=17) tAACGCCAGACTGTACTGCACCCATTTTGTGGGAc < 1:10046849/35‑1 (MQ=255) aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg < 1:19520888/35‑1 (MQ=255) aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg > 1:15502888/1‑35 (MQ=255) aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg > 1:15434413/1‑35 (MQ=255) aaCGCCAGACTGTACTGCACCCATTTTGTTGGACg > 1:11718126/1‑35 (MQ=255) aaCGCCAGACTGTACGGCACCCATTTTGTTGGACg > 1:11160342/1‑35 (MQ=255) aaCGCC‑GACTGTACTGCACCCATTTTGTTGGACGa > 1:7885979/1‑35 (MQ=14) aaCGCC‑GACTGTACTGCACCCATTTTGTTGGACGa > 1:6006317/1‑35 (MQ=14) aCGCCAGACTGTACTGCACCCATTTTGTTGGACGa < 1:3365432/35‑1 (MQ=255) aCGCCAGACTGTACTGCACCCATTTTGTTGGACGa > 1:678807/1‑35 (MQ=255) aCGCCAGACTGTACTGCACCCATTTTGTTGGACGa > 1:7307590/1‑35 (MQ=255) cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt < 1:1275610/35‑1 (MQ=255) cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt > 1:19050878/1‑35 (MQ=255) cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt < 1:7521486/35‑1 (MQ=255) cGCCAGACTGTACTGCACCCATTTTGTTGGACGAt < 1:15207157/35‑1 (MQ=255) gccgACTGTACTGCACCCATTTTGTTGGACGATGa > 1:19070967/4‑35 (MQ=255) gccgACTGTACTGCACCCATTTTGTTGGACGATGa < 1:11701875/32‑1 (MQ=255) gccgACTGTACTGCACCCATTTTGTTGGACGATGa > 1:3688206/4‑35 (MQ=255) cccgACTGTACTGCACCCATTTTGTTGGACGATGa < 1:11748627/32‑1 (MQ=255) ccAGACTGTACTGCACCCATTTTGTTGGACGATGa > 1:11600026/1‑35 (MQ=255) ccAGACTGTACTGCACCCATTTTGTTGGACGATGa > 1:11225130/1‑35 (MQ=255) ccAGACTGTACTGCACCCATTTTGTTGGACGATGa < 1:13220283/35‑1 (MQ=255) ccgACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:18274383/3‑35 (MQ=25) cAGACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:3007875/1‑35 (MQ=40) cAGACTGTACTGCACCCATTTTGTTGGACGATGaa > 1:18691505/1‑35 (MQ=40) gtaCTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:14084420/3‑35 (MQ=17) cgACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:3903288/2‑35 (MQ=25) cgACTGTACTGCACCCATTTTGTTGGACGATGaaa < 1:3473095/34‑1 (MQ=25) cgACTGTACTGCACCCATTTTGTTGGACGATGaaa < 1:1657623/34‑1 (MQ=25) aaaCTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:8413227/3‑35 (MQ=17) aGACTGTACTGCACCCATTTTGTTGGACTATGaaa > 1:5659756/1‑35 (MQ=255) aGACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:514955/1‑35 (MQ=38) aGACTGTACTGCACCCATTTTGTTGGACGATGaaa > 1:12036791/1‑35 (MQ=38) gACTGTACTGCACCCATTTTGTTGGACGATGAAAt > 1:4768112/1‑35 (MQ=35) gACTGTACTGCACCCATTTTGTTGGACGATGAAAt < 1:14295949/35‑1 (MQ=35) gACTGTACTGCACCCATTTTGTTGGACGATGAAAt > 1:15295427/1‑35 (MQ=35) gACTGTACTGCACACATTTTGTTGGACGATGAAAt < 1:5519752/35‑1 (MQ=255) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg < 1:993263/35‑1 (MQ=32) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg > 1:7352109/1‑35 (MQ=32) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg > 1:7309178/1‑35 (MQ=32) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg > 1:16830497/1‑35 (MQ=32) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg > 1:6324007/1‑35 (MQ=32) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg > 1:2436271/1‑35 (MQ=32) aCTGTACTGCACCCATTTTGTTGGACGATGAAATg > 1:3816741/1‑35 (MQ=32) | TCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATG > REL606/3893541‑3893586 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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. |