New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	4506708 =	25 (0.860)	3 (0.120)	3/152	2.0	10.3%	coding (827/954 nt)	fecR	fec operon regulator FecR
	?	NZ_CP009273	4506742 =	30 (1.150)					coding (793/954 nt)	fecR	fec operon regulator FecR
Rejected: Frequency below/above cutoff threshold.

CATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/4506800‑4506708 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgttgcGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAA  >  NZ_CP009273/4506742‑4506790                                                                                                                                                CATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGT                                                      >  1:835868/1‑90   TCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTG                                                    >  2:498148/1‑90   TCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTG                                                    >  2:507175/1‑90                    aacCCGCTGGGTGAGGTGATAGCCACGCTACCCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGT                                   <  2:272513/87‑1                               tggGGTGATACCCCCGCTAACCCGTTCCCGCACGGGCGTGCTGCGCTGCGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCT                        <  1:698783/87‑1                                                 AACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCC                                                        >  1:726991/1‑40                                                 AACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCC                                                        <  2:726991/40‑1                                                     CGTTCCCGCAGCGGCGTGCTGCGCTGGGATCCCGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAA  <  1:627225/90‑1                                                                                                                                                CATCCTGAGCTTCAGCGATAAACCGCTGGGTGAGGTGATAGCCACGCTAACCCGTTACCGCAACGGCGTGCTGCGCTGCGATCCCGCCGTTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/4506800‑4506708 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgttgcGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAA  >  NZ_CP009273/4506742‑4506790

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