Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	3,894,997	Δ6,934 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	3893609–3894994	3901930	6937–8322	4 [3]	[0] 23	[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	= 3652533	NA (NA)	23 (1.070)	20/146	0.2	100%	noncoding (1443/1443 nt)	IS150	repeat region
	?	REL606	3901931 =	0 (0.000)					coding (489/1428 nt)	yieO	predicted multidrug or homocysteine efflux system

GGGCATATGTTTGCGCGCGTAAAGAAGGCCTGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGG  >  REL606/3901859‑3901993                                                                         |                                                               gactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGtt                                                             <  2:636323‑M2/4‑1 (MQ=255)    tccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGcc                                                          >  1:106627‑M2/70‑76 (MQ=255)    tccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGcc                                                          >  2:376269‑M2/70‑76 (MQ=255)       aattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAg                                                       >  1:582872‑M2/67‑76 (MQ=255)            aatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGgtgcc                                                  >  2:325246‑M2/62‑74 (MQ=255)              tatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCa                                                >  2:669821‑M2/60‑76 (MQ=255)              tatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCa                                                <  1:683802‑M2/17‑1 (MQ=255)               atcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAg                                               <  1:277201‑M2/18‑1 (MQ=255)               atcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAg                                               >  1:423398‑M2/59‑76 (MQ=255)                      tcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTg                                                    >  2:709744‑M2/52‑64 (MQ=255)                      tcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTg                                                    <  1:709744‑M2/13‑1 (MQ=255)                      tcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCAcgccg                                        <  2:699102‑M2/25‑1 (MQ=255)                          acctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCaa                                    <  2:202792‑M2/29‑1 (MQ=255)                          acctatatgcctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCaa                                    <  2:375011‑M2/29‑1 (MQ=255)                                  gcctcgtgtttaactgtccacctttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACgggg                            <  2:553099‑M2/37‑1 (MQ=255)                                   cctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGc                           <  2:312002‑M2/38‑1 (MQ=255)                                   cctcgtgtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGc                           <  2:60745‑M2/38‑1 (MQ=255)                                         gtttaactgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTaaaa                     >  1:174690‑M2/33‑76 (MQ=255)                                              actgtccaactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAATTTGGc                >  2:702613‑M2/28‑76 (MQ=255)                                                     aactttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCccacca         >  1:461326‑M2/21‑76 (MQ=255)                                                       ctttttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGa       >  2:470682‑M2/19‑76 (MQ=255)                                                           ttggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCg   <  2:148658‑M2/62‑1 (MQ=255)                                                            tggggtcagtacaGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCgg  <  2:452230‑M2/63‑1 (MQ=255)                                                                         |                                                               GGGCATATGTTTGCGCGCGTAAAGAAGGCCTGCAATACCTATGGGGATATTGATTAAAAATATCCAGTGCCAGGTTGCCCAGGTGACCAGCACGCCGCCAAGAACGGGGCCTAAAATTGGCCCCACCAGACCCGG  >  REL606/3901859‑3901993

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