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: ‑