New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3209524 =	19 (0.650)	3 (0.120)	3/150	2.0	12.8%	coding (327/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
	?	NZ_CP009273	3209606 =	24 (0.930)					coding (245/765 nt)	nfeF	NADPH‑dependent ferric chelate reductase
Rejected: Frequency below/above cutoff threshold.

GGAAGGACCACGCCCACCGTCGCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGGTGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3209617‑3209524 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgacggtgggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAG  >  NZ_CP009273/3209606‑3209685                                                                                                                                                                                GGAAGGACCACGCCCACCGTCGCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGG                                                                                      <  1:441954/90‑1                                   CCGCTGTATGACGAACTACGCCATGACCTGGCGATTGTTTTTTTTTTTCACGACGGTGGGGGTGGTCCTTCCGGCCAGACGATGCCCTCT                                                    <  1:554671/90‑1                                                  CTACGCCATGAACTGGCGATTGATTTCTTTATTCCCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGC                                     <  2:433599/90‑1                                                   tgCGCCATGAACTGGCGATTGATTTCTTTATTCACGGCGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCG                                    <  2:219899/88‑1                                                                                     CGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAG  >  2:264315/1‑90                                                                                     CGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAG  <  2:307154/90‑1                                                                                                                                                                                GGAAGGACCACGCCCACCGTCGCGTGACTATACGCCGCTGTATGACGAACTACGCCATGAACTGGCGATTGATTTCTTTATTCACGACGGTGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3209617‑3209524 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgacggtgggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAG  >  NZ_CP009273/3209606‑3209685

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