| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 1207805 | 2 (0.060) | 29 (1.010) | 25/128 | 0.2 | 95.7% | coding (305/630 nt) | ycfK | e14 prophage; protein StfP |
| ? | NC_000913 | = 1209602 | 1 (0.030) | pseudogene (54/537 nt) | stfE | e14 prophage; putative side tail fiber protein fragment | |||||
ACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAT > NC_000913/1209527‑1209679 | aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttc < 1:231508‑M2/81‑6 (MQ=255) aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaac < 2:231509‑M2/81‑14 (MQ=255) agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag > 1:231510‑M2/1‑63 (MQ=255) gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt > 2:231511‑M2/1‑62 (MQ=255) ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt < 1:231512‑M2/81‑21 (MQ=255) cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttt > 1:231513‑M2/1‑60 (MQ=255) agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt < 1:231514‑M2/81‑23 (MQ=255) gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc > 2:231515‑M2/1‑58 (MQ=255) gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga < 2:231516‑M2/79‑26 (MQ=255) ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag < 1:231269‑M2/81‑29 (MQ=255) gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg > 1:231517‑M2/1‑50 (MQ=255) gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg < 1:231268‑M2/81‑32 (MQ=255) cGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt > 1:231518‑M2/1‑49 (MQ=255) gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt > 2:231519‑M2/1‑47 (MQ=255) gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttg < 1:231267‑M2/80‑36 (MQ=255) tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt < 2:231520/81‑1 (MQ=255) ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca < 1:231266‑M2/80‑38 (MQ=255) gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccataccaaggttttcgagagccgtttgcacagtg < 2:231264‑M2/81‑43 (MQ=255) gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg < 1:231265‑M2/81‑43 (MQ=255) aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg > 2:231521‑M2/1‑32 (MQ=255) aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcc < 1:231263‑M2/80‑49 (MQ=255) aaCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt < 2:231260‑M2/81‑54 (MQ=255) aaCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt < 1:231261‑M2/81‑54 (MQ=255) cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga < 2:231258‑M2/80‑56 (MQ=255) cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga < 2:231259‑M2/80‑56 (MQ=255) ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat > 2:231522‑M2/1‑23 (MQ=255) ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat < 1:231257‑M2/81‑59 (MQ=255) acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc > 1:231523‑M2/1‑19 (MQ=255) aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaa > 1:231524‑M2/1‑16 (MQ=255) aaTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattc < 2:231255‑M2/81‑73 (MQ=255) aaTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattc < 2:231254‑M2/81‑73 (MQ=255) agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg < 2:231253‑M2/81‑78 (MQ=255) | ACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAT > NC_000913/1209527‑1209679 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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. |