New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 4161725	20 (0.350)	3 (0.000)	3/122	4.0	0.45%	intergenic (+41/‑260)	rrfB/murB	5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
	?	Exported	= 2325	1315 (0.060)					intergenic (+91/–)	fdh/–	formate dehydrogenase from Pseudomonas sp. 101 with an N‑terminal his‑tag/–
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTC  >  CP009273/4161659‑4161743                                                                   |                   tGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGt                   >  1:2887989/1‑68 (MQ=35)   cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGttt                 >  1:5380192/1‑68 (MQ=35)   cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGttt                 >  1:2337058/1‑68 (MQ=35)   cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGttt                 >  1:2719066/1‑68 (MQ=35)   cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGgtt                 >  1:5615683/1‑68 (MQ=25)   cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTACGttt                 >  1:5270768/1‑68 (MQ=25)     agagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTa               <  1:204319/68‑1 (MQ=35)     agagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTa               >  1:3325716/1‑68 (MQ=35)      gagTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAt              <  1:4207977/68‑1 (MQ=35)        gTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCt            >  1:5161463/1‑68 (MQ=35)        gTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCt            >  1:1853044/1‑68 (MQ=35)           gggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgtt         <  1:3622921/68‑1 (MQ=35)            ggAACTGCCAGGCATCAATTAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgttg        >  1:4267235/1‑68 (MQ=255)            ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgttg        >  1:3040731/1‑68 (MQ=35)             gAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgttgt       >  1:2933398/1‑68 (MQ=35)             gAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgttgt       <  1:3814636/68‑1 (MQ=35)             gAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgttgt       >  1:2434248/1‑68 (MQ=35)               aCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTCCGTTTTATCTGttgttt     >  1:4025525/1‑68 (MQ=25)                 tGCCAGGCATCAAATAAAACGAAAGGCACAGTCGAAAGACTGGGCCTTTCGTTTTATCTGttgtttgt   >  1:233477/1‑68 (MQ=18)                  gCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTc  >  1:5228246/1‑68 (MQ=31)                                                                   |                   TGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTC  >  CP009273/4161659‑4161743

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