| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 1207805 | 15 (0.770) | 8 (0.550) | 8/100 | 0.4 | 42.9% | coding (305/630 nt) | stfP | e14 prophage; uncharacterized protein |
| ? | NC_000913 | = 1209602 | 10 (0.680) | pseudogene (18/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
GCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA > NC_000913/1209541‑1209669 | gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct > 2:12973‑M2/1‑62 (MQ=255) ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc > 1:419362‑M2/1‑61 (MQ=255) cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa > 2:151619‑M2/1‑56 (MQ=255) ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc > 2:29871/1‑68 (MQ=255) gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg > 2:458158/1‑68 (MQ=255) aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa > 1:104669/1‑68 (MQ=255) tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGGGCCTTCTCCCAAGCCAACGTTTATg > 2:612575/1‑68 (MQ=255) aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg > 2:43144/1‑68 (MQ=255) aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg > 2:362556/1‑68 (MQ=255) aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg > 1:666940‑M2/1‑37 (MQ=255) gCCACGGAACCGGGACACAAACAGGTACTGCAGAGcattctcccaaaccaaggttttcgagagccgtt < 1:294045‑M2/68‑34 (MQ=255) aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt > 2:424795‑M2/1‑27 (MQ=255) cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc > 2:173326‑M2/1‑25 (MQ=255) aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt < 1:156621‑M2/68‑53 (MQ=255) aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat > 1:66639/1‑68 (MQ=255) gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt < 2:662910/68‑1 (MQ=255) aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt > 2:282401/1‑68 (MQ=255) tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc < 1:78134‑M2/68‑62 (MQ=255) gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc < 1:29871/68‑1 (MQ=255) cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca < 2:409828‑M2/68‑64 (MQ=255) gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccggtttgatatcgccaaacg < 2:315333‑M2/68‑68 (MQ=255) | GCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA > NC_000913/1209541‑1209669 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ 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. |