Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	REL606	3,015,771	IS150 (–) +3 bp	coding (1192‑1194/1677 nt)	ECB_02816 →	KpsD protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2775877	NA (NA)	407 (1.320)	131/146	0.0	99.5%	noncoding (1/1443 nt)	IS150	repeat region
	?	REL606	3015771 =	2 (0.010)					coding (1192/1677 nt)	ECB_02816	KpsD protein
*	?	REL606	= 3015773	2 (0.010)	466 (1.490)	129/148	0.0	99.6%	coding (1194/1677 nt)	ECB_02816	KpsD protein
	?	REL606	= 3894996	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region

TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAC  >  REL606/3894924‑3894999                                                                         |    tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:1848810/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:5331728/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:5778187/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:9215443/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  1:9279978/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:1273937/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:1570890/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:2572593/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:2654632/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:3680702/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:4821219/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:5059568/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:5936316/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:9231086/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc  <  2:9973634/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGNCCAACTTTTTGGGGTCAGTACAGAc  <  1:3168221/76‑1 (MQ=21)                                                                         |    TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAC  >  REL606/3894924‑3894999

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