New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	386633 =	25 (1.540)	3 (0.190)	3/252	NT	12.6%	intergenic (‑95/‑231)	ltaE/glyA‑I	low specificity l‑threonine aldolase/serine hydroxymethyltransferase
	?	NC_002947	386722 =	17 (1.070)					intergenic (‑184/‑142)	ltaE/glyA‑I	low specificity l‑threonine aldolase/serine hydroxymethyltransferase

GTTTGCGCTGGCCCTTTCGCGGGTGAACCCGCTCCTACAGTGATTGGCGTTTGCAGGACCATTGTGGGGGCGGGCGCGCCCGCGAAGGGGCCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/386727‑386633 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcGCAAACAGCACAACAGCGAATGTCGTAAACGCGCCATTGCAAGGCGTGCGCAGGCATCTGACCCGGCCCCACCAGTCATAACATCGCCACAAAGGGCCACAGTGCCCCACGACAAAAAACGATCGCTGCCGGGAGATACACAA  >  NC_002947/386722‑386864                                                                                                                                                                                                                                                GTTTGCGCTGGCCCTTTCGCGGGTGAACCCGCTCCTACAGTGATTGGCGTTTGCAGGACCATTGTGGGGGCGGGCGCGCCCGCGAAGGGGCCGGCGCAAACAGCACAACAGCGAATGTCGTAAACGCGCCATTGCAAGGCGTGCGCAGGC                                                                                          >  1:438138/1‑150                                         TGATTGGCGTTTGCAGGACCATTGTGGGGGCGGGCGCGCCCGCGAAGGGGCCGGCGCAAACAGCACAACAGCGAATGTCGTAAACGCGCCATTGCAAGGCGTGCGCAGGCATCTGACCCGGCCCCACCAGTCATAA                                                                >  1:75146/1‑136                                         TGATTGGCGTTTGCAGGACCATTGTGGGGGCGGGCGCGCCCGCGAAGGGGCCGGCGCAAACAGCACAACAGCGAATGTCGTAAACGCGCCATTGCAAGGCGTGCGCAGGCATCTGACCCGGCCCCACCAGTCATAA                                                                <  2:75146/136‑1                                                                                         GGCCGGCGCAAACAGCACAACAGCGAATGTCGTAAACGCGCCATTGCAAGGCGTGCGCAGGCATCTGACCCGGCCCCACCAGTCATAACATCGCCACAAAGGGCCACAGTGCCCCACGACAAAAAACGATCGCTGCCGGGAGATACACAA  <  2:438138/150‑1                                                                                                                                                                                                                                                GTTTGCGCTGGCCCTTTCGCGGGTGAACCCGCTCCTACAGTGATTGGCGTTTGCAGGACCATTGTGGGGGCGGGCGCGCCCGCGAAGGGGCCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/386727‑386633 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcGCAAACAGCACAACAGCGAATGTCGTAAACGCGCCATTGCAAGGCGTGCGCAGGCATCTGACCCGGCCCCACCAGTCATAACATCGCCACAAAGGGCCACAGTGCCCCACGACAAAAAACGATCGCTGCCGGGAGATACACAA  >  NC_002947/386722‑386864

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