New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3991398 =	23 (1.050)	3 (0.150)	3/156	1.6	12.5%	coding (56/2163 nt)	uvrD	DNA helicase II
	?	NZ_CP009273	3991437 =	21 (1.050)					coding (95/2163 nt)	uvrD	DNA helicase II
Rejected: Frequency below/above cutoff threshold.

TCCACGCTCATCAACCAGGCGATACGATGCACCAGTACGCGCGTCTTACCACTGCCCGCGCCCGCCAGCACCAGAAGGTTGCTGCGTGGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3991491‑3991398 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcggGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC  >  NZ_CP009273/3991437‑3991524                                                                                                                                                                                        TCCACGCTCATCAACCAGGCGATACGATGCACCAGTACGCGCGTCTTACCACTGCCCGCGCCCGCCAGCACCAGAAGGTTGCTGCGTGGC                                                                                              <  2:113683/90‑1       CTCATCAACCAGGCGATACGATGCACCAGTACGCGCGTCTTACCACTGCCCGCGCCCGCCAGCACCAGAAGGTTGCTGCGTGGCGCGGGC                                                                                        <  2:496264/90‑1                             GCACCAGTCCGCGCGTCTTACCCCTGCCCGCGCCCGCCAGCACCAGAAGGTTGCTGCGTGGCGCGGGCAGTGGTAAGACGCGCGTACTGG                                                                  <  1:146677/90‑1                                 CAGTACGCGCGTCTTCCCACGGCCCGCGCCCGCCAGCACCAGAAGGTTGCTGCGTGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCA                                                              <  2:276450/90‑1                                                                                         GCGCGGGCAGTGGGAAAACGCGCGTACTGGTGCATCGGATCGCCCGGTTGATGAGGGTGGAAAACTGCTCGCCATACTCGATTATGGCGG      >  1:515017/1‑90                                                                                             GGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC  <  2:146677/90‑1                                                                                                                                                                                        TCCACGCTCATCAACCAGGCGATACGATGCACCAGTACGCGCGTCTTACCACTGCCCGCGCCCGCCAGCACCAGAAGGTTGCTGCGTGGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3991491‑3991398 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcggGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC  >  NZ_CP009273/3991437‑3991524

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