Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	3,582,206	Δ1,222 bp	100%	IS1‑mediated	[yhhZ]–yrhA	[yhhZ], yrhA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3582206	3584144–3583428	1223–1939	120 [0]	[47] 48	[yhhZ]–[insB1]	[yhhZ], yrhA, insA, [insB1]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3582205	0 (0.000)	107 (1.250)	72/178	NT	100%	coding (343/1179 nt)	yhhZ	putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
	?	NC_000913	3583428 =	NA (NA)					noncoding (1/768 nt)	IS1	repeat region

ATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAA  >  NC_000913/3583331‑3583441                                                                                                  |              aTTTTATAGAAATAATGATGATTTAATAAACCCTGATCTACAAGAAAGGATAGCGATCGAGGATTATAGCATTTCAATATTTACTTAGGACATTAAAgg              >  1:2280696/1‑99 (MQ=255)            ataatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCaa  <  2:2105633/100‑1 (MQ=255)                                                                                                  |              ATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAA  >  NC_000913/3583331‑3583441

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