Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_012967	3,894,998	Δ6,934 bp	100%		rbsD–[mdtD]	rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [mdtD]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_012967	3894998	3901931	6934	13 [0]	[0] 14	rbsD–[mdtD]	rbsD,rbsA,rbsC,rbsB,rbsK,rbsR,[mdtD]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	= 3894997	0 (0.000)	4 (0.470)	3/44	NT	100%	intergenic (+27/‑161)	ECB_RS25215/rbsD	IS3 family transposase/D‑ribose pyranase
	?	NC_012967	3901932 =	0 (0.000)					coding (488/1428 nt)	mdtD	multidrug transporter subunit MdtD

GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCCTTTCTTCATAAGGGGGCGTTTTTGTTT > NC_012967/3894963‑3895032 | gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAGg < 2:92476‑M1/36‑2 (MQ=255) tttttaaCTGTCCAACTTTTTGGGGTCAGTACAGgt < 2:507901‑M1/34‑3 (MQ=255) tgcttaacTGTCCAACTTCTTGGGGTCAGTACAGgt > 2:192327‑M1/4‑34 (MQ=255) gtTTAACTGTCCAACTTTTTGGGGTCAGTACAGgtt > 2:682036‑M1/1‑33 (MQ=255) tttAACTGTCCAACTTTTTGGGGTCAGTACAGgttg < 1:40355‑M1/36‑5 (MQ=255) aaCTGTCCAACTTTTTGGGGTCAGTACAGgttgccc > 2:187745‑M1/1‑29 (MQ=255) aaCTGTCCAACTTTTTGCGGTCAGTACAGgttgccc > 2:295431‑M1/1‑29 (MQ=255) cTTTTTGGGGTCAGTACAGgttgcccaggtgaccag < 1:203217‑M1/36‑18 (MQ=255) ggggTCAGTACAGgttgcccaggtgaccagcacgcc < 1:768231‑M1/36‑24 (MQ=255) aCAGgttgcccaggtgaccagcacgccgccaagaac < 1:404236‑M1/36‑33 (MQ=255) cAGgttgcccaggtgaccagcacgccgccaagaacg > 1:947199‑M1/1‑3 (MQ=255) ggttgcccaggtgaccagcacgccgccaagaacggg > 1:776216‑M1/1‑1 (MQ=255) ggttgcccaggtgaccagcacgccgccaagaacggg > 2:695312‑M1/1‑1 (MQ=255) | GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAGACGCCTCCTTTCTTCATAAGGGGGCGTTTTTGTTT > NC_012967/3894963‑3895032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑