New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	= 1733649	0 (0.000)	7 (0.730)	7/50	NT	100%	coding (685/1413 nt)	pykF	pyruvate kinase PykF
	?	NC_012967	2774435 =	NA (NA)					intergenic (‑67/+26)	ECB_RS13765/ECB_RS25125	type I toxin‑antitoxin system Hok family toxin/IS3 family transposase

GAATTGTTTGGGTTCGAACGTAAGCCTCAGGTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC > NC_012967/2774403‑2774466 | gAATTGTTTGGGTTCGAACGTAAGCCTCAGGTTGTa < 2:210461/36‑1 (MQ=255) aTTGTTTGGGTTCGAACGTAAGCCTCAGGTTGTACt > 2:636788/1‑36 (MQ=255) aTTGTTTGGGTTCGAACGTAAGCCTCAGGTTGTACt > 2:820090/1‑36 (MQ=255) aTTGTTTGGGTTCGAAAGTAAGCCTCAGGTTGTACt > 2:610854/1‑36 (MQ=38) ttGTTTGGGTTCGAACGTAAGCCTCAGGTTGTACtc > 2:372372/1‑35 (MQ=255) tttGGGTTCGAACGTAAGCCTCAGGTTGTACTGAcc < 1:324600/36‑1 (MQ=255) gggTTCGAACGTAAGCCTCAGGTTGTACTGACCCCa < 2:805314/36‑1 (MQ=255) ggTTCGAACGTAAGCCTCAGGTTGTACTGACCCCaa < 2:870892/36‑1 (MQ=255) gTTCGAACGTAAGCCTCAGGTTGTACTGACCCCaaa < 2:478185/36‑1 (MQ=255) gAACGTAAGCCTCAGGTTGTACTGACCCCAAAAAGt > 2:747940/1‑36 (MQ=255) aaCGTAAGCCTCAGGTTGTACTGACCCCAAAAAGtt > 1:915363/1‑36 (MQ=255) cGTAAGCCTCAGGTTGTACTGACCCCAAAAAGTTgg < 1:453574/36‑1 (MQ=255) gTAAGCCTCAGGTTGTACTGACCCCAAAAAGTTGGa > 2:230400/1‑36 (MQ=255) tAAGCCTCAGGTTGTACTGACCCCAAAAAGTTGGAc < 1:907820/36‑1 (MQ=255) gCCTCAGGTTGTACTGACCCCAAAAAGTTGGACAGt < 2:151088/36‑1 (MQ=255) aGGTTGTACTGACCCCAAAAAGTTGGACAGTTAAac > 2:671232/1‑36 (MQ=40) | GAATTGTTTGGGTTCGAACGTAAGCCTCAGGTTGTACTGACCCCAAAAAGTTGGACAGTTAAAC > NC_012967/2774403‑2774466

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 32 ≤ 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.