New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5397407	161 (0.980)	8 (0.050)	5/264	NT	5.0%	intergenic (+15/‑86)	hsdR/hsdM	type I restriction modification system endonuclease/type I restriction modification system methylase
	?	NC_002947	= 5397437	143 (0.880)					intergenic (+45/‑56)	hsdR/hsdM	type I restriction modification system endonuclease/type I restriction modification system methylase

AAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5397268‑5397407 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttAGGTGGACAA‑GGCTCCGCCGTTGCCCACCCAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCAGTTCGATGGGTGCAC  <  NC_002947/5397437‑5397301                                                                                                                                                                                                                                                                                        AAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTTAGGTGGACAA                                                                                                                                  >  2:3013654/1‑150                                                                                                                             ATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGTCCACCTAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCA                 <  1:1443570/138‑1                                                                                                                             ATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGTCCACCTAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCA                 <  1:1443665/138‑1                                                                                                                             ATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGTCCACCTAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCA                 >  2:1443570/1‑138                                                                                                                             ATCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGTCCACCTAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCA                 >  2:1443665/1‑138                                                                                                                              TCCGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGCCCACCCAACGGGTAGAGCGGATTAACTGGCGTACATCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCAGTTCGATGGGTGC    >  1:4426836/1‑150                                                                                                                                CGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGTCCACCTAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCAGTTCGATGGGTGCAC  >  1:3013654/1‑150                                                                                                                                CGCTCTACCCGTTGGGTGGGCAACGGCGGAGCC‑TTGTCCACCTAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCAGTTCGATGGGTGCAC  >  1:5185539/1‑150                                                                                                                                                                                                                                                                                        AAATCCTCCAGCTCGATCCCTTTAACCAGCTCGGTGCACCCATCGAACTGGTTAAAGCCTTTGGCGGCAAAGCTGGCTACAACAAAGCTATCCATGAACTGGAAGACGCGCTGTACGCCAGTTAATCCGCTCTACCCGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5397268‑5397407 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttAGGTGGACAA‑GGCTCCGCCGTTGCCCACCCAACGGGTAGAGCGGATTAACTGGCGTACAGCGCGTCTTCCAGTTCATGGATAGCTTTGTTGTAGCCAGCTTTGCCGCCAAAGGCTTTAACCAGTTCGATGGGTGCAC  <  NC_002947/5397437‑5397301

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