New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4580793 =	32 (1.370)	3 (0.130)	3/244	NT	8.5%	coding (2679/3321 nt)	treSB	bifunctional trehalose synthase B/maltokinase
	?	NC_002947	4580833 =	34 (1.520)					coding (2639/3321 nt)	treSB	bifunctional trehalose synthase B/maltokinase

CCGGCCTTCCAGCCGCGCCCCAGCGATGCGCAGGACAGCGAGCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4580840‑4580793 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctggAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTG  >  NC_002947/4580833‑4580934                                                                                                                                                        CCGGCCTTCCAGCCGCGCCCCAGCGATGCGCAGGACAGCGAGCGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTG  >  1:247230/1‑150         CCAGCCGCGCCCCAGCGATGCGCAGGACAGCGAGCGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGA                      >  1:164061/1‑122         CCAGCCGCGCCCCAGCGATGCGCAGGACAGCGAGCGCTGGAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGA                      <  2:164061/122‑1                                                                                                                                                        CCGGCCTTCCAGCCGCGCCCCAGCGATGCGCAGGACAGCGAGCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4580840‑4580793 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gctggAAGGCCGGGTCTTCGGTGGGTGCGGCCAGCAGCAGGTGCATCTCGCCCAGGCGCTGGCCAAGCAAGGCGGCAAAGCCGCTGAGTTCGGCCAATGCGTCAGTG  >  NC_002947/4580833‑4580934

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