New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | REL606 | 2038688 = | 16 (0.310) | 5 (0.100) | 3/138 | 6.8 | 25.8% | coding (261/822 nt) | wbbB | putative glycosyltransferase WbbB |
? | REL606 | 2038966 = | 14 (0.290) | coding (729/747 nt) | wbbA | putative glycosyltransferase WbbA | |||||
Rejected: Coverage evenness skew score above cutoff. |
TGAATAATAACTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAA > REL606/2038613‑2038762 | tgAATAATAACTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATgaga > 1:1338177/1‑76 (MQ=255) taataaCTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTg < 1:233220/76‑1 (MQ=255) ataaCTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGt < 2:674237/76‑1 (MQ=255) taaCTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTa > 1:741939/1‑76 (MQ=255) aaaGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAg < 2:269364/76‑1 (MQ=255) aaGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGc < 2:1339760/76‑1 (MQ=255) aGAATTTTGTATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCt > 1:927876/1‑76 (MQ=255) cATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaa < 1:1632647/76‑1 (MQ=255) tATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaaaa > 1:942555/1‑76 (MQ=255) tATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaaaa > 1:942575/1‑76 (MQ=255) tGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAAt > 1:1437320/1‑76 (MQ=255) cAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCtgca < 2:130770/76‑4 (MQ=255) tGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc > 1:592188/1‑76 (MQ=255) tGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc > 1:592191/1‑76 (MQ=255) agaATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCaa < 2:560152/76‑1 (MQ=255) gaATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCaaa < 1:570664/76‑1 (MQ=255) | TGAATAATAACTCACAATTTTCCTTAAAAAGAATTTTGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAA > REL606/2038613‑2038762 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ 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. |