New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 380111 | 0 (0.000) | 22 (0.830) +CGAGAT |
18/86 | 0.2 | 64.2% | coding (732/774 nt) | yaiO | outer membrane protein YaiO |
? | NC_000913 | 4070444 = | 28 (0.920) | intergenic (‑43/+71) | yihR/yihS | putative aldose 1‑epimerase YihR/sulfoquinovose isomerase |
GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGC > NC_000913/380062‑380156 | gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa < 1:1172584/50‑1 (MQ=255) gAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa > 1:2338467/1‑50 (MQ=255) tCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatg > 1:2364513‑M1/1‑43 (MQ=255) cAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatg > 1:1697561‑M1/1‑42 (MQ=255) aGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatgga > 1:159338‑M1/1‑41 (MQ=255) aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacg > 1:752070‑M1/1‑39 (MQ=255) aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacg < 1:1481071‑M1/50‑12 (MQ=255) cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccc < 1:2201215‑M1/50‑18 (MQ=255) aGACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccc > 1:2078592‑M1/1‑27 (MQ=255) gACGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccc < 1:501216‑M1/50‑25 (MQ=255) cGGGCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccct > 1:1476150‑M1/1‑24 (MQ=255) gggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgccccctt > 1:1264525‑M1/1‑23 (MQ=255) ggCTGCAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttg < 1:2642827‑M1/50‑29 (MQ=255) cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggag < 1:1143684‑M1/50‑34 (MQ=255) cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttgggg < 1:1094846‑M1/48‑32 (MQ=255) cAAGTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggg > 1:2187642‑M1/1‑17 (MQ=255) gTTGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagg > 1:1577121‑M1/1‑14 (MQ=255) tGCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt > 1:3813‑M1/1‑12 (MQ=255) gCAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggt > 1:11692‑M1/1‑11 (MQ=255) cAGACCGTTAcgagatggacggtcccctcgcccccttggggagagggtta < 1:175924‑M1/50‑41 (MQ=255) gACCGTTAcgagatggacggtcccctcgcccccttggggagagggttagg < 1:1379725‑M1/50‑43 (MQ=255) aCCGTTAcgagatggacggtcccctcgcccccttggggagagggttaggg < 1:2089348‑M1/50‑44 (MQ=255) cGTTAcgagatggacggtcccctcgcccccttggggagagggttagggtg < 1:339808‑M1/50‑46 (MQ=255) gTTAcgagatggacggtcccctcgcccccttggggagagggttagggt > 1:978776‑M1/1‑4 (MQ=255) | GAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGC > NC_000913/380062‑380156 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ 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. |