New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	87928 =	62 (1.030)	4 (0.070)	4/248	NT	6.9%	coding (1390/1518 nt)	betC	choline‑sulfatase
	?	NC_002947	87973 =	50 (0.880)					coding (1345/1518 nt)	betC	choline‑sulfatase

CGAGCGGGAGAACCTCACCGGCAGCCCGGACCACCAGGTGCTGCTGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/88079‑87928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctggccaGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGG  >  NC_002947/87973‑88115                                                                                                                                                                                                                                                                                                         CGAGCGGGAGAACCTCACCGGCAGCCCGGACCACCAGGTGCTGCTGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGG                                                                                                                                                    >  4:57661/1‑149    GCGGGAGAACCTCACCGGCAGCCCGGACCACCAGGTGCTGCTGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCA                                                                                                                                                 >  3:56859/1‑149                                                                                                       GGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCT                                              >  7:6234/1‑149                                                                                                            TGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGCCAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGC                                         >  3:78688/1‑149                                                                                                              GCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGC                                                                                                              <  7:296459/78‑1                                                                                                              GCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGC                                                                                                              >  8:296459/1‑78                                                                                                                                                  CTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGG   <  5:141529/149‑1                                                                                                                                                  CTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGG   <  7:110360/149‑1                                                                                                                                                   TGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGG  >  5:218963/1‑149                                                                                                                                                   TGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGG  >  5:61397/1‑149                                                                                                                                                                                                                                                                                                         CGAGCGGGAGAACCTCACCGGCAGCCCGGACCACCAGGTGCTGCTGCAGGCATTTGTCGATGAAGCGCAACAGCGCTGGGACATCACCAGCCTGCGCCAGCAGGTACTGGCCAGCCAGCGCCGCCGCCGGCTGGTGGCCGAAGCGCTGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/88079‑87928 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctggccaGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATAGAGTAGGCATGGG  >  NC_002947/87973‑88115

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.