| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913_3_pae_tpiA | = 4583053 | 28 (0.970) | 5 (0.180) | 5/274 | 3.2 | 16.4% | intergenic (‑89/+112) | hsdM/hsdR | DNA methyltransferase M/endonuclease R Type I restriction enzyme |
| ? | NC_000913_3_pae_tpiA | = 4583140 | 24 (0.870) | intergenic (‑176/+25) | hsdM/hsdR | DNA methyltransferase M/endonuclease R Type I restriction enzyme | |||||
| Rejected: Coverage evenness skew score above cutoff. | |||||||||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
AAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAA > NC_000913_3_pae_tpiA/4583012‑4583284 | aagaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGtata < 2:359404/149‑1 (MQ=255) gaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACACACGTACGTGACCGGCGTGAACGAGCGCAT‑CCAAAGCAGAGGCAGCCTGAAGGATGAAGTGTATACg > 1:71953/1‑149 (MQ=37) aaTTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCa > 1:440633/1‑149 (MQ=255) agtGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAg < 2:461379/147‑1 (MQ=255) tGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGc < 1:354538/149‑1 (MQ=255) aaTATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGACCGAGCGCAGCCCCACG‑GTAGCCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTcc < 1:497789/149‑1 (MQ=38) tataCACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCa < 1:44408/149‑1 (MQ=255) cacaAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGata < 2:165125/149‑1 (MQ=255) tCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGc > 2:421032/1‑149 (MQ=255) ttCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGa > 1:310466/1‑149 (MQ=255) aTCAAGGCGGCAAGTGAGTGAGTCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcag > 1:194492/1‑149 (MQ=255) aaGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCcagcag > 2:52829/1‑149 (MQ=255) cggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTa < 1:49803/149‑1 (MQ=255) aaaGGGAGTGAATCCCCGGGAGCGTACAGAAGTCCGTGACCGGAGTGACAGAGCGCAG‑CGAAGGCAGAGTAATCCTGAAGGATGAAATTGATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCg < 2:373267/148‑1 (MQ=255) tCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCa > 1:453093/1‑149 (MQ=255) gTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCtt < 2:193018/149‑1 (MQ=255) gAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGCTATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGAAACATACACTCCCcgccgctcc > 2:72842/1‑146 (MQ=255) aCGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACgg < 2:418368/149‑1 (MQ=255) cgcAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGaa < 1:421032/149‑1 (MQ=255) aG‑CGAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGTCGAATTTTCCTAGCAGGTTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGtt < 2:56740/149‑1 (MQ=255) g‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCAAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCCCCTTCCCCCCGCGACCGGGGAAGtgg > 1:133681/1‑147 (MQ=255) aaaggCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGcc < 2:78578/145‑1 (MQ=255) ggCAGTCTGAAGGATGAAGTGGATACGTGTCAGTCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCATCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGaa < 2:310466/149‑1 (MQ=255) ggCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGaa < 1:339538/149‑1 (MQ=255) | AAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAG‑CCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGAA > NC_000913_3_pae_tpiA/4583012‑4583284 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ 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. |