New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 2846300 | 0 (0.000) | 24 (0.940) | 21/90 | 0.1 | 100% | coding (172/1710 nt) | hycE | formate hydrogenlyase subunit HycE |
? | NC_000913 | 2867269 = | 0 (0.000) | coding (283/993 nt) | rpoS | RNA polymerase, sigma S (sigma 38) factor |
GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA > NC_000913/2867222‑2867316 | gcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatc < 1:1338542‑M2/2‑1 (MQ=255) cgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatca > 1:1771169‑M2/47‑49 (MQ=255) gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcat > 1:1441831‑M2/46‑49 (MQ=255) cattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc > 1:294269‑M2/38‑49 (MQ=255) ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc > 1:58076‑M2/36‑47 (MQ=255) cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCga > 1:1787802‑M2/34‑47 (MQ=255) cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCga > 1:2021656‑M2/34‑47 (MQ=255) tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc < 1:253139‑M2/19‑1 (MQ=255) gcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGa > 1:195113‑M2/28‑48 (MQ=255) gttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAt > 1:870919‑M2/26‑49 (MQ=255) cgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc > 1:1548915‑M2/23‑49 (MQ=255) cgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtct < 1:250209‑M2/26‑1 (MQ=255) gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc > 1:2020788‑M2/22‑48 (MQ=255) gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc > 1:392546‑M2/22‑48 (MQ=255) cgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCAc > 1:152286‑M2/20‑49 (MQ=255) gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCAc < 1:785302‑M2/30‑1 (MQ=255) gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg > 1:1457886‑M2/19‑49 (MQ=255) gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg > 1:1378421‑M2/19‑49 (MQ=255) taccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg > 1:34960‑M2/17‑47 (MQ=255) ccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt < 1:773073‑M2/35‑1 (MQ=255) acagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcgc > 1:1636963‑M2/11‑49 (MQ=255) agcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa > 1:56729‑M2/9‑49 (MQ=255) gcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc > 1:1870676‑M2/8‑49 (MQ=255) gcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc < 1:438479‑M2/42‑1 (MQ=255) aTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaa < 1:200518‑M2/48‑1 (MQ=255) tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaa > 1:1554893/1‑48 (MQ=255) | GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA > NC_000913/2867222‑2867316 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ 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. |