New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5433183 =	55 (1.020)	4 (0.080)	4/254	NT	6.9%	coding (2081/2529 nt)	hrpB	ATP‑dependent helicase HrpB
	?	NC_002947	5433238 =	55 (1.040)					coding (2026/2529 nt)	hrpB	ATP‑dependent helicase HrpB

CTGCAGAAGGATGGCCACAGCGAATGGCCTGACCTGGGGGACGAAGCCCTGCTGGCCAGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5433244‑5433183 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAACCTTTGCGCCGGACCAGGCCAAGCAGCGCCCGGGTCCTGGCTTCATCATCCAGCCCAGGCAGCGGCTC  >  NC_002947/5433238‑5433385                                                                                                                                                                                                                    CTGCAGAAGGATGGCCACAGCGAATGGCCTGACCTGGGGGACGAAGCCCTGCTGGCCAGCCTTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAACCTTTGC                                                               >  8:176187/1‑149      GAAGGATGGCCACAGCGAATGGCCTGACCTGGGGGACGAAGCCCTGCTGGCCAGCCTTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAA                                                                      <  7:90701/137‑1      GAAGGATGGCCACAGCGAATGGCCTGACCTGGGGGACGAAGCCCTGCTGGCCAGCCTTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAA                                                                      >  8:90701/1‑137                                 CCTGGGGGACGAAGCCCTGCTGGCCAGCCTTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAACCTTTGCGCCGGACCAGGCCAAGCAGCGCCCGGGTCCTG                               <  7:176187/149‑1                                                              TTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAACCTTTGCGCCGGACCAGGCCAAGCAGCGCCCGGGTCCTGGCTTCATCATCCAGCCCAGGCAGCGGCTC  >  3:138815/1‑149                                                                                                                                                                                                                    CTGCAGAAGGATGGCCACAGCGAATGGCCTGACCTGGGGGACGAAGCCCTGCTGGCCAGCCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5433244‑5433183 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctTCTGCAGGTCGAGCTGACGCAGCAGCGCCACACGTGCCTGCCATTGCCGCAGTTCCGGCGTCCAGGTCAGCAGGTTCAAACCTTTGCGCCGGACCAGGCCAAGCAGCGCCCGGGTCCTGGCTTCATCATCCAGCCCAGGCAGCGGCTC  >  NC_002947/5433238‑5433385

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.