New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3341848	29 (0.610)	4 (0.090)	4/252	NT	12.4%	coding (583/987 nt)	kdsD	D‑arabinose 5‑phosphate isomerase KdsD
	?	NC_000913	= 3341856	29 (0.640)					coding (591/987 nt)	kdsD	D‑arabinose 5‑phosphate isomerase KdsD

AGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3341721‑3341848 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acccagTGCGCCGCCTGGGTGTGAGA  <  NC_000913/3341856‑3341837                                                                                                                                                      AGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGTGCGCCGCCTGGGTG       <  2:41028/143‑1      GTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGTGCGCCGCCTGGGTGTGAGA  <  1:154670/143‑1                                                                          GCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGTGCGCCGCCTGGGTGTG     >  1:724135/1‑72                                                                          GCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGTGCGCCGCCTGGGTGTG     <  2:724135/72‑1                                                                                                                                                      AGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3341721‑3341848 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acccagTGCGCCGCCTGGGTGTGAGA  <  NC_000913/3341856‑3341837

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