Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,549,192	Δ39 bp	coding (385‑423/903 nt)	fimH →	type 1 fimbriae D‑mannose specific adhesin

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4549192	4549230	39	7 [0]	[0] 7	fimH	type 1 fimbriae D‑mannose specific adhesin

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4549191	0 (0.000)	7 (0.950)	5/312	0.2	100%	coding (384/903 nt)	fimH	type 1 fimbriae D‑mannose specific adhesin
	?	NC_000913	4549231 =	0 (0.000)					coding (424/903 nt)	fimH	type 1 fimbriae D‑mannose specific adhesin

ACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTG  >  NC_000913/4549007‑4549399                                                                                                                                                                                         |                                                                                                                                                                                                                 aCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttc                                                                                                                                                >  1:6375‑M1/1‑185 (MQ=255)                                                                   tATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttc                                                                               <  1:114272‑M1/250‑132 (MQ=255)                                                                   tATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttc                                                                               <  1:25905‑M1/250‑132 (MQ=255)                                                                   tATCTAATTTTTCCGGGACCGTAAAATATAGGGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttc                                                                               <  1:118803‑M1/250‑132 (MQ=255)                                                                                                                cATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtt                                                                                 <  1:13111‑M1/203‑130 (MQ=255)                                                                                                                                      gccgcGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttatt           >  1:156545‑M1/1‑52 (MQ=255)                                                                                                                                                gttTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTaaagctggctcattaattgccgtgcttattttgcgacagaccaacaactataacagcgatgatttccagtttgtgtggaatatttacgccaataatgatgtggtggtgcctactggcggctgcgatgtttctgctcgtgatgtcaccgttactctgccggactaccctggttcagtgccaattcctcttaccgtttattgtgcgaaaa  >  1:66855‑M1/1‑42 (MQ=255)                                                                                                                                                                                         |                                                                                                                                                                                                                 ACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGGCGCTTTATTTGACGCCTGTGAGCAGTGCGGGCGGGGTGGCGATTAAAGCTGGCTCATTAATTGCCGTGCTTATTTTGCGACAGACCAACAACTATAACAGCGATGATTTCCAGTTTGTGTGGAATATTTACGCCAATAATGATGTGGTGGTGCCTACTGGCGGCTGCGATGTTTCTGCTCGTGATGTCACCGTTACTCTGCCGGACTACCCTG  >  NC_000913/4549007‑4549399

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