| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| MC JC | REL606 | 3,894,998 | Δ6,935 bp | IS150‑mediated | rbsD–[yieO] | rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO] |
| Missing coverage evidence... | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| seq id | start | end | size | ←reads | reads→ | gene | description | |||
| * | * | ÷ | REL606 | 3893607–3894996 | 3901932 | 6937–8326 | 85 [77] | [1] 217 | [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 | = 3894997 | NA (NA) | 208 (1.070) +GT |
106/146 | 0.1 | 99.5% | intergenic (+27/‑161) | insK‑5/rbsD | IS150 putative transposase/predicted cytoplasmic sugar‑binding protein |
| ? | REL606 | 3901933 = | 1 (0.010) | coding (487/1428 nt) | yieO | predicted multidrug or homocysteine efflux system | |||||
TCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAC > REL606/3894922‑3894999 | tCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg > 1:2061809/1‑76 (MQ=31)tCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg > 1:3201989/1‑76 (MQ=31)tCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg > 1:547039/1‑76 (MQ=31)tCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg < 1:6107878/76‑1 (MQ=31)tCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAg > 2:3347671/1‑76 (MQ=31) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgta < 2:3167630/76‑3 (MQ=14) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:5144309/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:1417130/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:504391/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:4508625/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:425132/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:3584081/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:1948275/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 2:1828162/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:863312/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:662545/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:5660644/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:368797/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:3210390/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc < 1:2998459/76‑1 (MQ=34) | TCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAC > REL606/3894922‑3894999 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ 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. |