New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	87921 =	65 (1.000)	4 (0.060)	4/258	NT	5.6%	coding (1397/1518 nt)	betC	choline‑sulfatase
	?	NC_002947	88027 =	71 (1.120)					coding (1291/1518 nt)	betC	choline‑sulfatase

TGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCATCGGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/88035‑87921 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caaATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGGTCGTCTTCGCTGTACACGAACTTGTAGGGCCCGCGGCGAATCATCATCAGCGGGCCGA  >  NC_002947/88027‑88173                                                                                                                                                                                                                                                                        TGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCATCGGCAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCC                                                                                                                   >  1:372625/1‑149       CATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCTGCCGGCTGGTGGCCGAAGCGCTGGCCATCGGCAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAG                                                                                                             >  5:193445/1‑149                                       ATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCTGCCGGCTGGTGGCCGAAGCGCTGGCCATCGGCAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGT                                                                             <  6:193445/149‑1                                               CCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCATCGGCAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGT                                                                     <  2:372625/149‑1                                                                                                                  AAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGGTCGTCTTCGCTGTACACGAACTTGTAGGGCCCGCGGCGAATCATCATCAGCGGGCCGA  >  1:18089/1‑149                                                                                                                                                                                                                                                                        TGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCATCGGCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/88035‑87921 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caaATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGGTCGTCTTCGCTGTACACGAACTTGTAGGGCCCGCGGCGAATCATCATCAGCGGGCCGA  >  NC_002947/88027‑88173

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