| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913_3_bme_pgi | 1207790 = | 4 (0.130) | 17 (0.640) | 16/254 | 0.7 | 81.8% | coding (290/630 nt) | stfP | e14 prophage; uncharacterized protein |
| ? | NC_000913_3_bme_pgi | 1209619 = | 4 (0.150) | pseudogene (1/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAC > NC_000913_3_bme_pgi/1209482‑1209753 | ggcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa < 2:231502‑M2/12‑1 (MQ=255) gcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt < 2:496407‑M2/13‑1 (MQ=255) ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaataataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAATTga > 2:525912‑M2/135‑149 (MQ=255) tttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcttaatgtagccgttgttgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa < 2:508203‑M2/16‑1 (MQ=255) aTTCTTCAGCAGAAAAAGCGGCACCTATGCATTTCAGCCAGCCTGTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATTCAGAGGCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc < 2:135537/149‑1 (MQ=255) agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCCAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc > 1:213271/1‑149 (MQ=255) ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt > 1:501289‑M2/103‑149 (MQ=255) tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc < 1:382048‑M2/52‑1 (MQ=255) tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttatcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc < 1:112784‑M2/52‑1 (MQ=255) tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctctcaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg > 1:470632‑M2/95‑149 (MQ=255) aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg < 1:451182‑M2/57‑1 (MQ=255) ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa > 1:74972‑M2/86‑149 (MQ=255) caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa > 2:302997‑M2/85‑149 (MQ=255) cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg < 2:278855/149‑1 (MQ=255) accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa > 1:512764‑M2/77‑96 (MQ=255) accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa < 2:512764‑M2/20‑1 (MQ=255) aaccaaaccagcgcgtccaaatgtagccgttgccgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa < 1:55852‑M2/84‑1 (MQ=255) atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc > 1:494068‑M2/65‑149 (MQ=255) cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCTTTTACCAGGGGGATTTAACATACTTATTGGATATGTACGCGTATCAACAAATGACCAGAACACAGATCt > 2:465807/1‑149 (MQ=255) gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc < 2:391870‑M2/118‑1 (MQ=255) ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc < 1:302997‑M2/121‑1 (MQ=255) ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGATTATTGGCTCTGTACGCGTATGAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTgaat > 2:322708‑M2/15‑148 (MQ=255) | ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAAC > NC_000913_3_bme_pgi/1209482‑1209753 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ 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. |