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.

ACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTC  >  NZ_CP009273/4506619‑4506794                                                                                          |                                                                                       aCTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCgg                                                                                        >  2:202525/1‑90 (MQ=255)    gAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCaa                                                                                     <  2:707250/90‑1 (MQ=255)                             gATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGcc                                                            <  1:822540/90‑1 (MQ=255)                                       aGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTaa                                                  <  1:93059/90‑1 (MQ=255)                                        ggATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAAc                                                 <  2:742139/90‑1 (MQ=255)                                         ttcagcgataaaccgctgggtgaggtgatagccacgctaacccgttaccGCAACGGCGGGATCCCAGCGCAGCACGCCGCTGCGGGAACg                                                <  1:627225‑M1/41‑1 (MQ=255)                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCggg                                                                              >  1:36122/1‑57 (MQ=255)                                            cGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCggg                                                                              <  2:36122/57‑1 (MQ=255)                                             gCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGtt                                            >  2:696034/1‑90 (MQ=255)                                                 cGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCg                                        <  2:115054/90‑1 (MQ=255)                                                  ggTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTg                                                        <  1:744047/73‑1 (MQ=255)                                                  ggTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTg                                                        >  2:744047/1‑73 (MQ=255)                                                      tttttCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCt                                   <  2:830738/90‑1 (MQ=255)                                                      tttttCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCt                                   <  1:619566/90‑1 (MQ=255)                                                        tttCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAt                                 <  2:75938/90‑1 (MQ=255)                                                             gCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATTCCAGCGCAGCACGCCGGTGGGGGAACGGGTTAGCGGGGCTAtcaccc                            >  1:93245/1‑89 (MQ=255)                                                             gCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTAtcacct                            <  1:336758/90‑1 (MQ=255)                                                               aggtgatagccacgctaacccgttaccGCAACGGCGGGATCGCAGCGCAGCACGCCCGTGCGGGAACGGGTTAGCGGGGgtatcacccca                          <  1:698783‑M1/63‑4 (MQ=255)                                                                     gTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGc                    <  2:482723/90‑1 (MQ=255)                                                                     gTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGc                    <  2:363981/90‑1 (MQ=255)                                                                         cGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGtt                >  2:578684/1‑90 (MQ=255)                                                                          acgctaacccgttaccGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGGTAGCGTGGCTATCACCTCACCCagcgggtt               <  2:272513‑M1/74‑4 (MQ=255)                                                                                 cgcAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTg        <  2:138186/90‑1 (MQ=255)                                                                                  gcAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGa       <  2:465115/90‑1 (MQ=255)                                                                                  gcAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGa       <  2:702504/90‑1 (MQ=255)                                                                                   cAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGaa      <  1:456104/90‑1 (MQ=255)                                                                                     gCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGGGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGc    >  2:57277/1‑90 (MQ=255)                                                                                       ccGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTc  <  2:481757/90‑1 (MQ=255)                                                                                          |                                                                                       ACTGAATTTTAACGGGAAGCGTTTGCGCGATAACGTTCAGGATCGCATCGGTATTTTTCAGCGGGAACGTCCCGCTCAGGCGCAGCCCGGCAACGGCGGGATCGCAGCGCAGCACGCCGTTGCGGTAACGGGTTAGCGTGGCTATCACCTCACCCAGCGGTTTATCGCTGAAGCTC  >  NZ_CP009273/4506619‑4506794

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