New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | REL606 | 2038432 = | 51 (1.030) | 11 (0.250) | 6/134 | 4.6 | 22.1% | coding (517/822 nt) | wbbB | putative glycosyltransferase WbbB |
? | REL606 | 2038941 = | 32 (0.720) | coding (8/822 nt) | wbbB | putative glycosyltransferase WbbB | |||||
Rejected: Coverage evenness skew score above cutoff. |
TTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGA > REL606/2038869‑2039015 | ttATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAAt < 2:1157168/76‑1 (MQ=255) tatataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATc < 1:771402/76‑1 (MQ=255) tatataATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATc < 1:771403/76‑1 (MQ=255) ttttATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATtat > 2:1311233/2‑76 (MQ=255) tttATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATtata < 1:284800/76‑1 (MQ=255) tattaATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATaaa > 1:501429/1‑76 (MQ=255) taattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatca < 1:1417540/76‑1 (MQ=255) taattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatca < 1:1417522/76‑1 (MQ=255) aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat < 1:1040481/76‑1 (MQ=255) aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat < 1:786060/76‑1 (MQ=255) aattaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcat < 1:1040480/76‑1 (MQ=255) attaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcatc < 1:1163524/76‑1 (MQ=255) attaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAatcatc < 1:1163516/76‑1 (MQ=255) ttaTGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATcc < 2:351600/76‑1 (MQ=255) aTTCTGGGTTTAAAGTATAAGAACGGGCCACAATGCAGAATTTTTTTTTACAATCCATTATAAATCATCCCAattg < 1:607041/76‑1 (MQ=255) tCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgat > 1:379637/1‑76 (MQ=255) cTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgatt < 1:471297/76‑1 (MQ=255) cTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgatt < 1:471293/76‑1 (MQ=255) cTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAattgatt < 1:471287/76‑1 (MQ=255) ggCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAAttgattga < 1:1313439/76‑1 (MQ=255) cTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAg > 1:334561/1‑76 (MQ=255) tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg < 1:690703/76‑1 (MQ=255) tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg < 1:690714/76‑1 (MQ=255) tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg < 1:1082147/76‑1 (MQ=255) tAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTg < 1:1133245/76‑1 (MQ=255) aaGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACaaa > 1:869962/1‑76 (MQ=255) gAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAaaga > 2:532975/1‑76 (MQ=255) gggcacaATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCa > 1:1496483/4‑76 (MQ=255) gCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCAt < 1:1269629/76‑1 (MQ=255) tttttttttACAAACCAATAAAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGaaa < 1:1256453/76‑1 (MQ=255) ttttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAag > 1:1518727/1‑76 (MQ=255) ttttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAag > 1:1518715/1‑76 (MQ=255) ttttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAagag > 1:1882/1‑76 (MQ=255) tttttACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGc > 1:104456/1‑76 (MQ=255) tACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCgggg > 1:1375135/1‑76 (MQ=255) aCAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGa > 1:1374818/1‑76 (MQ=255) aCAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGa > 1:1374800/1‑76 (MQ=255) | TTATATAATTTATTAATTATGAATTCTGGCTTTAAAGTATAAGAACAGGCCACAATGCATAATTTTTTTTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGA > REL606/2038869‑2039015 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ 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. |