New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4582041 =	70 (1.170)	4 (0.070)	4/234	NT	5.8%	coding (1431/3321 nt)	treSB	bifunctional trehalose synthase B/maltokinase
	?	NC_002947	4582087 =	62 (1.060)					coding (1385/3321 nt)	treSB	bifunctional trehalose synthase B/maltokinase

CAAGCAGCAAAAGGCCTTTGGCCGCGGCAGCCTGCGTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4582080‑4582041 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑accGCCAGCATGCGCCGGGTCCAGTTCAGCAGTGAATGCGGGTCGTGGCTTTGCGCCTCGACGTTGACAGTCTGGTAACCGTAGAGCGGGTCCATGATCGGCGGCAGCACCAGGCGCTGCGGGTCGGCGCGGGAAAAACCGCCGTTGC  >  NC_002947/4582087‑4582231                                                                                                                                                                                           CAAGCAGCAAAAGGCCTTTGGCCGCGGCAGCCTGCGTACCGCCAGCATGCGCCGGGTCCAGTTCAGCAGTGAATGCGGGTCGTGGCTTTGCGCCTCGACGTTGACAGTCTGGTAACCGTAGAGCGG                                                             >  1:1446122/1‑126 CAAGCAGCAAAAGGCCTTTGGCCGCGGCAGCCTGCGTACCGCCAGCATGCGCCGGGTCCAGTTCAGCAGTGAATGCGGGTCGTGGCTTTGCGCCTCGACGTTGACAGTCTGGTAACCGTAGAGCGG                                                             <  2:1446122/126‑1          AAAGGCCTTTGGCCGCGGCAGCCTGCGTACCGCCAGCATGCGCCGGGTCCAGTTCAGCAGTGAATGCGGGTCGTGGCTTTGC                                                                                                >  1:1290486/1‑82          AAAGGCCTTTGGCCGCGGCAGCCTGCGTACCGCCAGCATGCGCCGGGTCCAGTTCAGCAGTGAATGCGGGTCGTGGCTTTGC                                                                                                <  2:1290486/82‑1                                    gcnCCGCCAGCATGCGCCGGGTCCAGTNCAGCAGTGAATGNGGGTCGNGNCTTTGCGCCTCGACGTTGACAGTCTGGTAACCGTAGAGCGGGTCCATGATCGGCGGCAGCACCAGGCGCTGCGGGTCGGCGCGGGAAAAACCGCCGTTGC  >  1:1136289/4‑150                                                                                                                                                                                           CAAGCAGCAAAAGGCCTTTGGCCGCGGCAGCCTGCGTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/4582080‑4582041 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑accGCCAGCATGCGCCGGGTCCAGTTCAGCAGTGAATGCGGGTCGTGGCTTTGCGCCTCGACGTTGACAGTCTGGTAACCGTAGAGCGGGTCCATGATCGGCGGCAGCACCAGGCGCTGCGGGTCGGCGCGGGAAAAACCGCCGTTGC  >  NC_002947/4582087‑4582231

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