| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NZ_CP009273 | 3516287 = | 18 (0.870) | 3 (0.160) | 3/154 | 1.6 | 17.0% | coding (58/2553 nt) | mrcA | peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA |
| ? | NZ_CP009273 | 3516322 = | 13 (0.700) | coding (93/2553 nt) | mrcA | peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA | |||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
TTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCAATA > NZ_CP009273/3516236‑3516411 | ttCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACaa < 1:312921/90‑1 (MQ=255) acatcttttaatgtcgccacatccggcagttgtGGCTCGATTTATGGCCTATACCGCTACACCGGACCCCCACTTCCGGGTGTGGGGACa < 1:202317‑M1/57‑1 (MQ=255) ccTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACAtt > 1:399526/1‑90 (MQ=255) tgtccccacacccggaagtgggggtccggtgtagcggtataggccataaatcgagccACAACTGCCGGATGTGGCGACATTAAAAGATGt < 1:202317‑M2/33‑1 (MQ=255) cTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGt > 2:261601/1‑90 (MQ=255) gCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGcc < 1:395112/90‑1 (MQ=255) aTTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTg > 1:193621/1‑90 (MQ=255) ctgctgGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCaa > 2:341763/1‑90 (MQ=255) ctgctgGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCaa > 2:561186/1‑90 (MQ=255) GGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAg < 2:124702‑M1/90‑1 (MQ=255) ccTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGa < 2:241190/90‑1 (MQ=255) tataCCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATg > 1:513006/1‑90 (MQ=255) tataCCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATg > 1:416111/1‑90 (MQ=255) ataggccataaatcgagccACAACTGCCGGATGTGGCGa < 1:44948‑M2/20‑1 (MQ=255) ataggccataaatcgagccACAACTGCCGGATGTGGCGa > 2:44948‑M2/20‑39 (MQ=255) gCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGc > 1:276392/1‑90 (MQ=255) cATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGAt < 2:289814/90‑1 (MQ=255) cATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGAt < 2:329974/90‑1 (MQ=255) gagccACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGAt < 2:99299‑M2/81‑1 (MQ=255) gagccACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGAt > 1:99299‑M2/6‑86 (MQ=255) ccACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCaa < 2:286705‑M2/88‑1 (MQ=255) acaACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCAATa < 1:203625/90‑1 (MQ=255) | TTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGGCTCGATTTATGGCCTATACCGCTACATCGAGCCACAACTGCCGGATGTGGCGACATTAAAAGATGTTCGCCTGCAAATTCCGATGCAGATTTACAGCGCCGATGGCGAGCTGATTGCTCAATA > NZ_CP009273/3516236‑3516411 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ 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. |