New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5557381	42 (0.710)	3 (0.050)	3/244	NT	6.1%	coding (249/1935 nt)	PP_4888	methyl‑accepting chemotaxis transducer
	?	NC_002947	= 5557420	53 (0.950)					coding (210/1935 nt)	PP_4888	methyl‑accepting chemotaxis transducer

CGAGAGCACCGGCACGTCGTTGTCCCAGATCGCCCGCGGCGCAACCGAAGCCAGCAGCTCGGCCATGTCGTTGGCGGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTGCGCTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5557228‑5557381 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgctcTGGGCTCGACGGCCTGTCCAGCAGTACCCGCCAGCGCCTGAGCGCAGG  <  NC_002947/5557420‑5557373                                                                                                                                                                                                            CGAGAGCACCGGCACGTCGTTGTCCCAGATCGCCCGCGGCGCAACCGAAGCCAGCAGCTCGGCCATGTCGTTGGCGGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTG                                                       >  1:11932/1‑149      GCACCGGCACGTCGTTGTCCCAGATCGCCCGCGGCGCAACCGAAGCCAGCAGCTCGGCCATGTCGTTGGCGGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTGCGCT                                                   >  2:347390/1‑148                                                    CAGCAGCTCGGCCATGTCGTTGGCGGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTGCGCTCTGGGCTCGACGGCCTGTCCAGCAGTACCCGCCAGCGCCTGAGCGCA    <  5:296110/149‑1                                                                            GGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTGCGCTCTGGGCTCGACGGCCTGTCCAGCAGTACCCGCCAGCGCCTGAGCGCAGG  <  3:370605/127‑1                                                                            GGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTGCGCTCTGGGCTCGACGGCCTGTCCAGCAGTACCCGCCAGCGCCTGAGCGCAGG  >  4:370605/1‑127                                                                                                                                                                                                            CGAGAGCACCGGCACGTCGTTGTCCCAGATCGCCCGCGGCGCAACCGAAGCCAGCAGCTCGGCCATGTCGTTGGCGGAATCCTTGAGGTTCTTTTCCAGGGTGCTGCGTAGTTGCTGCTGCTCGCTTTGCAGGCGCTCCGACAACCCTGCGCTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5557228‑5557381 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgctcTGGGCTCGACGGCCTGTCCAGCAGTACCCGCCAGCGCCTGAGCGCAGG  <  NC_002947/5557420‑5557373

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