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)	292 (1.230)	110/146	0.1	99.6%	noncoding (1/1443 nt)	IS150	repeat region
	?	REL606	3015771 =	1 (0.000)					coding (1192/1677 nt)	ECB_02816	KpsD protein
*	?	REL606	= 3015773	1 (0.000)	370 (1.540)	125/148	0.0	99.7%	coding (1194/1677 nt)	ECB_02816	KpsD protein
	?	REL606	= 3894996	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region

TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACG  >  REL606/3894924‑3895000                                                                         |     tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:1441825/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:960730/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:753043/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:6856693/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:6126019/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:43171/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:3293411/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:2691143/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:2636182/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:2472202/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:2119891/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  2:1450692/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:1300455/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:7774136/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:7634006/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:6049823/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:4173713/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:4124921/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:3738925/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:3092100/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:2442653/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:2074622/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:1940633/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:15597/76‑1 (MQ=34) tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGAc   <  1:1315342/76‑1 (MQ=34)  gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgttt  <  2:6548988/76‑4 (MQ=14)                                                                         |     TGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACG  >  REL606/3894924‑3895000

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