Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_002947	3,866,735	Δ8 bp	6.0%	coding (180‑187/1020 nt)	PP_3415 ←	LacI family transcriptional regulator

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3866734	66 (1.260)	4 (0.080)	4/240	NT	6.0%	coding (188/1020 nt)	PP_3415	LacI family transcriptional regulator
	?	NC_002947	3866743 =	60 (1.170)					coding (179/1020 nt)	PP_3415	LacI family transcriptional regulator

TTCGGCAAGGTCGTAATGATAGTTGCCGATCAGCACCTCTAGCCCACGAGGGCGCATGACCTCGTGAATGGCCTCCAGGGTGTCGATGAACAGTTGGTTGGACAACGACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3866583‑3866734 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGAC  >  NC_002947/3866743‑3866786                                                                                                                                                                                                      TTCGGCAAGGACGTAAAGATAGTTGCCGATCAGCACCTATAGCCCACGAGGGCGCATGACCTCGAGAATGGCCTCAAGGGTGTCGATGAACAGTTGGTTGGACAACGACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAG                                                >  2:413095/1‑150                                             CACGAGGGCGCATGACCTCGTGAATGGCCTCCAGGGTGTCGATGAACAGTTGGTTGGACAACGACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTG            >  1:577820/1‑142                                             CACGAGGGCGCATGACCTCGTGAATGGCCTCCAGGGTGTCGATGAACAGTTGGTTGGACAACGACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTG            <  2:577820/142‑1                                                                                                           GACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGAC  <  1:71105/90‑1                                                                                                           GACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGAC  >  2:71105/1‑90                                                                                                                                                                                                      TTCGGCAAGGTCGTAATGATAGTTGCCGATCAGCACCTCTAGCCCACGAGGGCGCATGACCTCGTGAATGGCCTCCAGGGTGTCGATGAACAGTTGGTTGGACAACGACGGGATCAACACCACCACCGACTGGCTGCGTGCCGAGGCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3866583‑3866734 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGAC  >  NC_002947/3866743‑3866786

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