New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1551725 =	33 (1.430)	3 (0.140)	3/234	NT	9.4%	coding (1139/1455 nt)	pykA	pyruvate kinase II
	?	NC_002947	1551778 =	27 (1.250)					coding (1086/1455 nt)	pykA	pyruvate kinase II

CAGGCCATGGCACGCATCTGCCTGGGTGCCGAAAAGCACCCGACCAGCCAGAAGTCCAGCCATCGCCTGCACACCACCTTCCAGCGCTGCGACGAAAGCATCGCCTTGGCGGCCATGTACACTGCCAACCACTTCCCGGGCGTAAAGGCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/1551876‑1551725 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaggcgatGCTTTCGTCGCAGCGCTGGAAGGTGGTGTGCAGGCGATGGCTGGACTTCTGGCTGGTCGGGTGCTTTTCGGCACCCAGGCAGATGCGTGCCATGGC  >  NC_002947/1551778‑1551873                                                                                                                                                                                                                                                          CAGGCCATGGCACGCATCTGCCTGGGTGCCGAAAAGCACCCGACCAGCCAGAAGTCCAGCCATCGCCTGCACACCACCTTCCAGCGCTGCGACGAAAGCATCGCCTTGGCGGCCATGTACACTGCCAACCACTTCCCGGGCGTAAAGGCG                                                                                                    >  2:139080/1‑150                                                                                          CGACGAAAGCATCGCCTTGGCGGCCATGTACACTGCCAACCACTTCCCGGGCGTAAAGGCGATGCTTTCGTCGCAGCGCTGGAAGGTGGTGTGCAGGCGATGGCT                                                        <  1:189430/105‑1                                                                                          CGACGAAAGCATCGCCTTGGCGGCCATGTACACTGCCAACCACTTCCCGGGCGTAAAGGCGATGCTTTCGTCGCAGCGCTGGAAGGTGGTGTGCAGGCGATGGCT                                                        >  2:189430/1‑105                                                                                                   CATCGCCTTGGCGGCCATGTACACTGCCAACCACTTCCCGGGCGTAAAGGCGATGCTTTCGTCGCAGCGCTGGAAGGTGGTGTGCAGGCGATGGCTGGACTTCTGGCTGGTCGGGTGCTTTTCGGCACCCAGGCAGATGCGTGCCATGGC  >  2:402293/1‑150                                                                                                                                                                                                                                                          CAGGCCATGGCACGCATCTGCCTGGGTGCCGAAAAGCACCCGACCAGCCAGAAGTCCAGCCATCGCCTGCACACCACCTTCCAGCGCTGCGACGAAAGCATCGCCTTGGCGGCCATGTACACTGCCAACCACTTCCCGGGCGTAAAGGCGAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/1551876‑1551725 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aaggcgatGCTTTCGTCGCAGCGCTGGAAGGTGGTGTGCAGGCGATGGCTGGACTTCTGGCTGGTCGGGTGCTTTTCGGCACCCAGGCAGATGCGTGCCATGGC  >  NC_002947/1551778‑1551873

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