| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | 1207790 = | 6 (0.320) | 13 (0.910) | 11/100 | 0.2 | 77.6% | coding (290/630 nt) | stfP | e14 prophage; uncharacterized protein |
| ? | NC_000913 | 1209619 = | 3 (0.210) | pseudogene (1/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTA > NC_000913/1209552‑1209682 | gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac > 2:156492‑M2/68‑68 (MQ=255) ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg > 2:637565‑M2/67‑68 (MQ=255) tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt > 2:404344‑M2/66‑68 (MQ=255) tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt > 2:106485‑M2/66‑68 (MQ=255) tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg > 2:548642‑M2/55‑68 (MQ=255) gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag > 1:683607/1‑68 (MQ=255) gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa > 1:275566‑M2/44‑68 (MQ=255) gtagccgttgctgcattgagttttgctacttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa > 1:185376‑M2/44‑68 (MQ=255) ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa > 2:96188‑M2/40‑68 (MQ=255) ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc < 2:8117‑M2/35‑1 (MQ=255) attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca < 1:512389‑M2/39‑1 (MQ=255) tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt < 1:478022‑M2/41‑1 (MQ=255) gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt > 2:538544‑M2/25‑68 (MQ=255) tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt < 1:644016/64‑1 (MQ=255) tGCAGAGCCTTCTCACAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt > 2:644018/1‑64 (MQ=255) tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc > 2:219436‑M2/23‑68 (MQ=255) caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa < 1:585078‑M2/60‑1 (MQ=255) aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTTCTTTTACCAGGGGGAt < 1:658010‑M2/61‑1 (MQ=255) ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa > 1:116254‑M2/5‑68 (MQ=255) | GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTA > NC_000913/1209552‑1209682 |
| 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. |