New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NZ_CP009273 | = 3512449 | 2 (0.220) | 12 (1.320) | 10/166 | 0.1 | 85.7% | intergenic (‑26/+375) | aroK/hofQ | shikimate kinase AroK/DNA uptake porin HofQ |
? | NZ_CP009273 | 3576788 = | NA (NA) | coding (391/417 nt) | BW25113_RS17825 | YrhA family protein |
TGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAATTATC > NZ_CP009273/3512369‑3512538 | gtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccGACTATTCGTTAATGATa > 1:64059‑M2/73‑90 (MQ=255) ttttCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgcca > 2:221139‑M1/1‑77 (MQ=255) ggCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttac < 2:102928‑M1/90‑20 (MQ=255) tGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg < 2:261821‑M1/90‑49 (MQ=255) gCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggt < 1:221139‑M1/90‑50 (MQ=255) aaaacactatcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGaa < 1:136558‑M2/59‑1 (MQ=255) ctatcaataagttggagtcattaccGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCa > 1:137921‑M2/26‑90 (MQ=255) ttttttCGGTACTACTAAGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttc < 1:221789‑M1/90‑64 (MQ=255) aGACTATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtc < 1:90670‑M1/90‑81 (MQ=255) tATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctc > 1:13422‑M1/1‑6 (MQ=255) tATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg < 1:162620‑M1/54‑49 (MQ=255) tATTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg > 2:162620‑M1/1‑6 (MQ=255) tATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGaaaa > 1:88076/1‑90 (MQ=255) aTTCGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcc > 1:194473‑M1/1‑5 (MQ=255) ttCGggtgatgctgccaacttactgatttagtgtat > 1:149165‑M1/1‑4 (MQ=255) ttCGggtgatgctgccaacttactgatttagtgtat < 2:149165‑M1/36‑33 (MQ=255) gTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAATTAtc < 2:88076/90‑1 (MQ=255) | TGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAATTATC > NZ_CP009273/3512369‑3512538 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ 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. |