| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 1207805 | 2 (0.080) | 12 (0.620) | 8/100 | 0.7 | 87.2% | coding (305/630 nt) | stfP | e14 prophage; uncharacterized protein |
| ? | NC_000913 | = 1209602 | 2 (0.100) | pseudogene (18/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA > NC_000913/1209538‑1209669 | tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcct < 2:200506‑M2/68‑4 (MQ=255) gcctgcctGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct > 2:234010‑M2/1‑62 (MQ=255) ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACGCCAACAGGTAATGCAGGGccttctc < 2:808329‑M2/68‑8 (MQ=255) ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc < 1:749196‑M2/68‑8 (MQ=255) gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc > 1:703164‑M2/1‑39 (MQ=255) aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg > 1:596900‑M2/1‑37 (MQ=255) aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg > 1:740446‑M2/1‑37 (MQ=255) aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg > 1:961471‑M2/1‑37 (MQ=255) aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtcttcgagagccg > 1:890623‑M2/1‑37 (MQ=255) aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc > 1:220365‑M2/1‑32 (MQ=255) gACACCTACAGGTAATGCAGAGccttctcccaaatcaaggttttcgagagccgtttgcacagtgccat > 2:914534‑M2/1‑22 (MQ=255) ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga > 1:927169‑M2/1‑18 (MQ=255) aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca < 2:278217/68‑1 (MQ=255) tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc < 1:286905‑M2/68‑62 (MQ=255) agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccgg < 2:952959‑M2/68‑65 (MQ=255) agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa < 1:308699‑M2/68‑65 (MQ=255) gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg < 1:870156‑M2/68‑68 (MQ=255) gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa > 2:763506/1‑68 (MQ=255) | TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA > NC_000913/1209538‑1209669 |
| 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. |