Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	3,894,998	Δ6,935 bp	IS150‑mediated	rbsD–[yieO]	rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	3893603–3894996	3901932	6937–8330	24 [22]	[0] 77	[insJ‑5]–[yieO]	[insJ‑5],insK‑5,rbsD,rbsA,rbsC,rbsB,rbsK,rbsR,[yieO]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3894997	NA (NA)	76 (1.200) +GT	62/146	0.0	100%	intergenic (+27/‑161)	insK‑5/rbsD	IS150 putative transposase/predicted cytoplasmic sugar‑binding protein
	?	REL606	3901933 =	0 (0.000)					coding (487/1428 nt)	yieO	predicted multidrug or homocysteine efflux system

TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAC  >  REL606/3894924‑3894999                                                                          |   tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:1313339/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:600687/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:1613416/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:2066399/76‑1 (MQ=34)                                                                          |   TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAC  >  REL606/3894924‑3894999

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