| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_017635 | = 2314113 | 13 (0.320) | 16 (0.580) | 16/68 | 0.4 | 52.5% | coding (52/414 nt) | ECW_RS28245 | hypothetical protein |
| ? | NC_017635 | = 2315924 | 20 (0.720) | pseudogene (17/459 nt) | ECW_RS28250 | phage tail protein | |||||
GGCATTGCCGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCAC > NC_017635/2315876‑2315972 | ggCATTGCCGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACc < 2:413255‑M2/50‑2 (MQ=255) cATTGCCGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACctt > 1:143206‑M2/1‑47 (MQ=255) cATTGCCGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACctt > 2:1577319‑M2/1‑47 (MQ=255) gCCGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctc > 2:1481161‑M2/1‑43 (MQ=255) ccGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctccc > 2:1684963‑M2/1‑42 (MQ=255) ccGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcc > 2:1397658‑M2/1‑42 (MQ=255) cgcgGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctccca > 2:1640062‑M2/1‑41 (MQ=255) cgcgGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctccc < 2:2060066‑M2/49‑9 (MQ=255) gcgGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaa > 2:1808168‑M2/1‑40 (MQ=255) gcgGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctccca < 1:493397‑M2/49‑10 (MQ=255) cgGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaaa < 1:1643018‑M2/50‑12 (MQ=255) ggAGGGCCAGAAGAACGGAACGCCAGTAGCTGGAGCACcttctcccaaac > 1:1787566‑M2/1‑38 (MQ=255) ggAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaaac < 2:1967417‑M2/50‑13 (MQ=255) ggAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaaac > 2:1426321‑M2/1‑38 (MQ=255) gAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaaacc > 1:229002‑M2/1‑37 (MQ=255) ggCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaaaccaac < 1:841211‑M2/50‑17 (MQ=255) gCCAGAAGAACGGAACGCCAATAGCTGGAGCACcttctcccaaaccaacg < 1:774519‑M2/50‑18 (MQ=255) agaagaACGGAACGCCAATAGCTGGAGCACcttctcccaaaccaacgtt < 2:711750‑M2/49‑20 (MQ=255) agaagaACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTAt < 1:283615/50‑1 (MQ=255) agaagaACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTAt > 1:212750/1‑50 (MQ=255) agaagaACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTAt < 2:751001/50‑1 (MQ=255) gaagaaCGGAACGCCAATAGCTGGAGCACcttctcccaaaccaacgttta < 2:1465375‑M2/50‑22 (MQ=255) agaaCGGAACGCCAATAGCTGGAGCACcttctcccaaaccaacgtttatg > 2:1748525‑M2/1‑27 (MQ=255) gaacggaacgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCga < 2:1750583/50‑1 (MQ=255) gaacggaacgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCga < 2:1995906/50‑1 (MQ=255) gaacgCCAATAGCTGGAGCACcttctcccaaaccaacgtttatgaaaatg > 1:754851‑M2/1‑21 (MQ=255) acgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGAcc < 1:1411849/50‑1 (MQ=255) cgCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCa < 1:1421064/50‑1 (MQ=255) ccAATAGCTGGAGCACcttctcccaaaccaacgtttatgaaaatgcagag > 1:1912886‑M2/1‑16 (MQ=255) ccAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGc < 1:1047356/50‑1 (MQ=255) cAATAGCTGGAGCACcttctcccaaaccaacgtttatgaaaatgca < 1:492555‑M2/46‑32 (MQ=255) aaTAGCTGGAGCACcttctcccaaaccaacgtttatgaaaatgcagagat > 2:1049632‑M2/1‑14 (MQ=255) aTAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCcac < 1:734039/50‑1 (MQ=255) aGCTGGAGCACcttctcccaaaccaacgtttatgaaaatgcagagataac < 1:564272‑M2/50‑40 (MQ=255) gCTGGAGCACCTTCTCCAAAACAAAGGTATGCGAGAAGACCAGCCACATc < 2:269553/50‑1 (MQ=18) cTGGAGCACcttctcccaaaccaacgtttatgaaaatgcagagataacgg > 2:488095‑M2/1‑9 (MQ=255) cTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCt < 2:277974/50‑1 (MQ=255) cTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCt < 1:1462674/50‑1 (MQ=255) tGGAGCACcttctcccaaaccaacgtttatgaaaatgcagagataacggg < 1:975435‑M2/50‑43 (MQ=255) ggAGCACcttctcccaaaccaacgtttatgaaaatgcagagataacggg < 2:1515576‑M2/49‑43 (MQ=255) ggAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCttt > 2:962680/1‑50 (MQ=255) ggAGCACCTTCTCCCAAACCAAGGAATGCGAGAAGACCAGCCACATCttt > 2:642542/1‑50 (MQ=255) gAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCtttt > 2:1707333/1‑50 (MQ=255) gAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCtttt > 2:371237/1‑50 (MQ=255) aGCACcttctcccaaaccaacgtttatgaaaatgcagagataacgggcaa < 2:1662760‑M2/50‑46 (MQ=255) gCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTcc < 1:1659267/50‑1 (MQ=255) gCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTcc < 1:101414/50‑1 (MQ=255) ccccTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCa < 2:1826922/48‑1 (MQ=255) ccccTTCCCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCa < 2:1569146/48‑1 (MQ=17) cACcttctcccaaaccaacgtttatgaaaatgcagagataacgggcaact < 1:444720‑M2/50‑48 (MQ=255) aCcttctcccaaaccaacgtttatgaaaatgcagagataacgggcaactg < 2:14547‑M2/50‑49 (MQ=255) | GGCATTGCCGCGGAGGGCCAGAAGAACGGAACGCCAATAGCTGGAGCACCTTCTCCCAAACCAAGGTATGCGAGAAGACCAGCCACATCTTTTCCAC > NC_017635/2315876‑2315972 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ 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. |