| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 1207805 | 22 (1.100) | 5 (0.330) | 5/100 | 1.0 | 26.6% | coding (305/630 nt) | stfP | e14 prophage; uncharacterized protein |
| ? | NC_000913 | = 1209602 | 11 (0.730) | pseudogene (18/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
TTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT > NC_000913/1209536‑1209666 | tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc > 2:476215‑M2/1‑67 (MQ=255) ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc > 2:388601‑M2/1‑61 (MQ=255) gccTGTTGGCGGAGTGGCTGAAGGCCACGGACACGGGACACAAACAGGTAATGCAGAGccttctccca < 2:172845‑M2/68‑11 (MQ=255) cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa > 1:441207‑M2/1‑56 (MQ=255) ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc > 1:521347/1‑68 (MQ=255) ggCGGAGTGGCTGAAGGCAACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc > 1:179971/1‑68 (MQ=255) gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg > 2:190400/1‑68 (MQ=255) aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg > 2:143890‑M2/1‑37 (MQ=255) aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa > 1:54818/1‑68 (MQ=255) ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa < 1:370783/68‑1 (MQ=255) gACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttatcgagcgccgtttgcacagagccat < 1:204811‑M2/68‑47 (MQ=255) acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg > 2:385925‑M2/1‑19 (MQ=255) aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTAATGAAAATGAAGAAATAACAAGCAAATGGcatca > 1:260485/1‑68 (MQ=255) gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt > 2:513568/1‑68 (MQ=255) gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt < 1:190400/68‑1 (MQ=255) gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt < 1:188021/68‑1 (MQ=255) tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc < 2:222574‑M2/68‑62 (MQ=255) tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc > 1:84339‑M2/1‑7 (MQ=255) tGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg > 2:318447/1‑68 (MQ=255) agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATAATTCCTGCtt > 2:375897/1‑68 (MQ=255) | TTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTT > NC_000913/1209536‑1209666 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 37 ≤ 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. |