New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5827863 =	33 (1.500)	3 (0.140)	3/244	NT	8.9%	intergenic (+9/+190)	mtgA/rpoH	monofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
	?	NC_002947	5827987 =	30 (1.430)					intergenic (+133/+66)	mtgA/rpoH	monofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32

GTGAGCGGCTTTTTTGTTTGGCTTTGCCGGCCTCTTCGCGGGTAAGCCCGCTCCTACAGGGGATGCGGTTGCCTGTAGGCGCGGGTTTGCCCGCGAGGAGGCCGGGGATGCGGTTGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5828012‑5827863 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACT  >  NC_002947/5827987‑5828049                                                                                                                                                                                                                       GTGAGCGGCTTTTTTGTTTGGCTTTGCCGGCCTCTTCGCGGGTAAGCCCGCTCCTACCAGGGATTCGGTTGCCTGTAGGCGCGGGTTTGCCCGCGAGGAGGCCGGGGATGCGGTTGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGG                                                                 <  2:17155/150‑1                                          GTAAGCCCGCTCCTACAGGGGATGCGGTTGCCTGTAGGCGCGGGTTTGCCCGCGAGGAGGCCGGGGATGCGGTTGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCG                        >  2:183152/1‑150                                                                                             GCGAGGAGGCCGGGGATGCGGTTGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACT  >  1:553967/1‑121                                                                                             GCGAGGAGGCCGGGGATGCGGTTGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACT  <  2:553967/121‑1                                                                                                                                                                                                                       GTGAGCGGCTTTTTTGTTTGGCTTTGCCGGCCTCTTCGCGGGTAAGCCCGCTCCTACAGGGGATGCGGTTGCCTGTAGGCGCGGGTTTGCCCGCGAGGAGGCCGGGGATGCGGTTGCCCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5828012‑5827863 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACT  >  NC_002947/5827987‑5828049

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