New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 3209533	20 (0.740)	3 (0.130)	3/150	2.1	14.1%	coding (318/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
	?	NZ_CP009273	= 3209595	19 (0.790)					coding (256/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
Rejected: Frequency below/above cutoff threshold.

CTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3209439‑3209533 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccaccgtcgCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGGTGGGGTCGCCAGCGGCTGGGC  <  NZ_CP009273/3209595‑3209507                                                                                                                                                                                          CTTCCGGCCCCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCAC                                                                                                <  2:31939/90‑1                                               ATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCGCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACT                                                  >  1:573671/1‑90                                                                GCATCGCCCAGCCGCTGGCGACCCCACCGTCGCGTGACTATACGCCGCTGTATGACGAACCACGCCCTGAACTGGCGATTGATTTCTTTA                                 >  1:424548/1‑90                                                                                    ACCCCACCGTCGCGTGGCTATACGCCGCTGGATGGCGGACTACGCCATGGACTGGCGATTGATTTCTTTATTCACGACGGGGGGGTCGCC             >  1:654260/1‑90                                                                                          CCGTCGCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGGTGGGGTCGCCAGCGGC       <  1:118396/90‑1                                                                                               GCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGGTGGGGTCGCCAGCGGCTGGGC  >  2:348234/1‑90                                                                                                                                                                                          CTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCGCCTGCATCGCCCAGCCGCTGGCGACCCCACCGTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/3209439‑3209533 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cccaccgtcgCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGGTGGGGTCGCCAGCGGCTGGGC  <  NZ_CP009273/3209595‑3209507

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