New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2921558	70 (0.710)	4 (0.040)	4/250	NT	5.3%	coding (543/1353 nt)	gudP	putative D‑glucarate transporter
	?	NC_000913	= 2921578	76 (0.790)					coding (523/1353 nt)	gudP	putative D‑glucarate transporter

GTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2921408‑2921558 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gacgCATGAAGTGGGCTGGTCACACGTCTTCTTCTT  <  NC_000913/2921578‑2921547                                                                                                                                                                                         GTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGA                                    <  3:713819/149‑1   TTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGAC                                   >  2:211124/1‑148                                                         TCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGCATGAAGTGGGCTGGTCACACGTCTTCTTCTT  >  3:400049/1‑127                                                         TCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGCATGAAGTGGGCTGGTCACACGTCTTCTTCTT  <  4:400049/127‑1                                                          CTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGCATGAAGTGGGCTGGTCACACGTCTTCTTCT   >  1:425012/1‑125                                                          CTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACGCATGAAGTGGGCTGGTCACACGTCTTCTTCT   <  2:425012/125‑1                                                                                                                                                                                         GTTTTGCTGATCCATATTGATCAGCGCACCACCCGCGGCGATGTACTCCAGCTCTTTCTTATTTACCCCCGGATGTTGATTTGGCTCGTGGATGACTTTCAACCAGATAAAGCTGATGACAATCCCCAGACCGCCCATAAAGAAGAAGACG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2921408‑2921558 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gacgCATGAAGTGGGCTGGTCACACGTCTTCTTCTT  <  NC_000913/2921578‑2921547

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