New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 168441711 (0.540)4 (0.220)
+CAGTGACAAC
3/146 1.5 21.6% coding (309/1509 nt) ydgA YdgA family protein
?NZ_CP009273 = 3404280 22 (1.080)coding (642/966 nt) dusB tRNA dihydrouridine synthase DusB
Rejected: Coverage evenness skew score above cutoff.

GTCCCTTCCCGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/1684404‑1684417
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTCGAGCACAGCT  <  NZ_CP009273/3404280‑3404268
              ||||||||||             
GTCCCTTCCCGCTTCAGTGACAACGTCGAGCACAGCT  <  1:234203/37‑1
GTCCCTTCCCGCTTCAGTGACAACGTCGAGCACAGCT  >  2:234203/1‑37
GTCCCTTCCCGCTTCAGTGACAACGTCGAGCACAGC   >  1:293053/1‑36
GTCCCTTCCCGCTTCAGTGACAACGTCGAGCACAGC   <  2:293053/36‑1
              ||||||||||             
GTCCCTTCCCGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/1684404‑1684417
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTCGAGCACAGCT  <  NZ_CP009273/3404280‑3404268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ 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.