New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 4353280	48 (1.720)	3 (0.120)	3/158	2.2	6.9%	coding (1580/1698 nt)	dsbD	protein‑disulfide reductase DsbD
	?	NZ_CP009273	= 4353321	36 (1.390)					coding (1539/1698 nt)	dsbD	protein‑disulfide reductase DsbD
Rejected: Frequency below/above cutoff threshold.

CGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGCCTTGTCCGTCAAAAAAGAGAATTGTCGGTAGGCCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/4353188‑4353280 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccaaCGACGCACAAGATGTGGCGCTGTTAAAGCATCTTAATGTCCTTGGCCTACCGACAATTCTCTTTTTTGACGGACAAGGCCAGGAG  <  NZ_CP009273/4353321‑4353237                                                                                                                                                                                    CGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGCCTTGTCCGTCAAAAAAGAGAATTGTCGGTAGGC                                                                                          >  1:464517/1‑90        GGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGCCTTGTCCGTCAAAAAAGAGAATTGTCGGTAGGCCAACGAC                                                                                   <  1:352719/90‑1                                                           TTGTCCGTCAAAAAAGAGAATTGTCGGTAGGCCAACGACGCACAAGATGTGGCGCTGTTAAAGC                                                          >  1:322801/1‑64                                                           TTGTCCGTCAAAAAAGAGAATTGTCGGTAGGCCAACGACGCACAAGATGTGGCGCTGTTAAAGC                                                          <  2:322801/64‑1                                                                                         GCCAACGACGCACAAGATGTGGCGCTGTTAAAGCATCTTAATGTCCTTGGCCTACCGACAATTCTCTTTTTTGACGGACAAGGCCAGGAG  <  1:346615/90‑1                                                                                         GCCAACGACGCACAAGATGTGGCGCTGTTAAAGCATCTTAATGTCCTTGGCCTACCGACAATTCTCTTTTTTGACGGACAAGGCCAGGAG  <  1:360562/90‑1                                                                                                                                                                                    CGCTGAAGGTTTCAGCATCCATAAAGCCCGTGACGCGTGCTTGTGGATGCTCCTGGCCTTGTCCGTCAAAAAAGAGAATTGTCGGTAGGCCAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/4353188‑4353280 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccaaCGACGCACAAGATGTGGCGCTGTTAAAGCATCTTAATGTCCTTGGCCTACCGACAATTCTCTTTTTTGACGGACAAGGCCAGGAG  <  NZ_CP009273/4353321‑4353237

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