New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3033617 =	35 (1.090)	4 (0.130)	5/244	2.9	11.9%	intergenic (+4/+40)	idi/lysS	isopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive
	?	NC_000913	3033639 =	26 (0.860)					intergenic (+26/+18)	idi/lysS	isopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive
Rejected: Frequency below/above cutoff threshold.

GCCGACAGCAGGTCTGGGAATTGGTATCGACCGTATGGTAATGCTGTTCACCAACAGCCATACCATCCGCGACGTTATTCTGTTCCCGGCGATGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGGCAAATGTCGGGGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3033765‑3033617 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggttGTGAGCATAACGTAATGCTTATTTTACCGGAC  >  NC_000913/3033639‑3033670                                                                                                                                                                                       GCCGACAGCAGGGCTGGGAATTGGTATCGACCGTATGGTAATGCTGGTCCCCAACAGCCATACCATCCGCGACGTTATTCTGTTCCCGGCGAGGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGACATTTGCCGGGGTT                                  <  3:27295/149‑1                        GTATCGACCGTATGGTAATGCTGTTCACCAACAGCCATGCCATCCGCGACGTTATTCTGTTCCCGGCGATGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATT           >  7:71334/1‑149                                CGTATGGTAATGCTGTTCACCAACAGCCATACCATCCGCGACGTTATTCTGTTCCCGGCGATGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGA   <  8:71334/149‑1                                                                     GCGACGTTATTCTGTTCCCGGCGATGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGAC  >  1:144518/1‑113                                                                     GCGACGTTATTCTGTTCCCGGCGATGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGAC  <  2:144518/113‑1                                                                                                                                                                                       GCCGACAGCAGGTCTGGGAATTGGTATCGACCGTATGGTAATGCTGTTCACCAACAGCCATACCATCCGCGACGTTATTCTGTTCCCGGCGATGCGTCCGGTAAAATAAGCATTACGTTATGCTCACAACCCCGGCAAATGTCGGGGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3033765‑3033617 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggggttGTGAGCATAACGTAATGCTTATTTTACCGGAC  >  NC_000913/3033639‑3033670

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