New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | W3110S.gb | = 71265 | 32 (0.860) | 16 (0.600) +AGGGTCCCTAAGTC |
12/72 | NT | 56.9% | coding (879/879 nt) | araC | DNA‑binding transcriptional dual regulator |
? | W3110S.gb | 3464879 = | 2 (0.050) | intergenic (‑255/+46) | murB/rrfB | UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding/5S ribosomal RNA |
TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGC > W3110S.gb/3464818‑3464933 | ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAAc > 1:819810‑M2/62‑67 (MQ=255) ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGa < 1:621464‑M2/8‑1 (MQ=255) ggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGaa < 1:1832753‑M2/9‑1 (MQ=255) ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAgg > 1:1668443‑M2/58‑69 (MQ=255) gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAGGcc < 1:3661064‑M2/14‑1 (MQ=255) gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAGGcc < 1:2388612‑M2/14‑1 (MQ=255) gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAGGcc < 1:3111071‑M2/14‑1 (MQ=255) gaaaaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTc > 1:3311496‑M2/51‑69 (MQ=255) aaagtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTCtt > 1:3208720‑M2/48‑68 (MQ=255) agtgaatgatgtagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAgg > 1:897176‑M2/46‑57 (MQ=255) tagccgtcaagttgtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt > 1:3272500‑M2/35‑69 (MQ=255) aagttgtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt > 1:3557426‑M2/27‑61 (MQ=255) aagttgtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCttt > 1:3027670‑M2/27‑61 (MQ=255) gttgtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAttt < 1:2476669‑M2/45‑1 (MQ=255) gCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAttt > 1:1439526/1‑67 (MQ=35) gtcataaagggtccctaagtcTAAAACGAAAGGccc > 1:98218‑M2/22‑36 (MQ=255) gtcataaagggtccctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTgac > 1:2027634‑M2/22‑68 (MQ=255) aaagggtccctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATttgaccg < 1:1212705‑M2/50‑4 (MQ=255) cGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGc < 1:206146/69‑1 (MQ=35) cctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTg > 1:2391757‑M2/9‑36 (MQ=255) cctaagtcTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAttt < 1:25057‑M2/45‑1 (MQ=255) aGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgg < 1:232480/57‑1 (MQ=35) | TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGC > W3110S.gb/3464818‑3464933 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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. |