New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	727836 =	NA (NA)	3 (0.080)	3/154	2.8	17.2%	coding (2798/4194 nt)	rhsC	rhs element protein RhsC
	?	NZ_CP009273	727859 =	16 (0.400)					coding (2821/4194 nt)	rhsC	rhs element protein RhsC
Rejected: Frequency below/above cutoff threshold.

GACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGA  >  NZ_CP009273/727772‑727947                                                                                        |                                                                                         gACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGt                                                                                        >  1:1120889/1‑90 (MQ=35)  atcatccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTc                                                                                       <  1:994495‑M2/4‑1 (MQ=255)   tcatccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCa                                                                                      >  1:228547‑M2/86‑90 (MQ=255)     atccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACg                                                                                    <  2:693547‑M2/7‑1 (MQ=255)     atccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACg                                                                                    <  1:1060897‑M2/7‑1 (MQ=255)     atccgtgcggataaccccttccgggatgaggtcggttttccctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACg                                                                                    <  1:148062‑M2/7‑1 (MQ=255)      tccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACgg                                                                                   <  1:295804‑M2/8‑1 (MQ=255)      tccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACgg                                                                                   <  2:653424‑M2/8‑1 (MQ=255)      tccgtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACgg                                                                                   >  2:596602‑M2/83‑90 (MQ=255)         gtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACgg                                                                                   >  2:855405‑M2/80‑87 (MQ=255)         gtgcggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACgg                                                                                   <  1:855405‑M2/8‑1 (MQ=255)            gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGcg                                                                             >  1:312545/1‑89 (MQ=18)            gCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCt                                                                             >  2:180914/1‑90 (MQ=33)            cggataaccccttccgggatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACGGCAGGCt                                                                             <  1:955049‑M2/14‑1 (MQ=255)              aCCCGGACAGCGCCCTGAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGa                                                                           <  1:877306/90‑1 (MQ=18)                       gcgcCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGaaaa                                                                  >  2:286463/1‑90 (MQ=33)                             gatgaggtcggttttctctgtcagcctgccgtgacggtcataccggtaaagatagtgcgCGTCACGGCAGGCTGACAGAGAAAACCGAcc                                                            <  1:393493‑M2/31‑1 (MQ=255)                                       gttttccgtgtccgcctgccgggacggtcatcccggtaaagatagggcgCGTCACGGCAGGCTGACAGAGAAACCCGACCTCATCCCGGa                                                  <  1:494519‑M2/41‑1 (MQ=255)                                                                      tcgcTTTACCGGTATGACCGTCCCGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacc                   <  1:917300/87‑1 (MQ=255)                                                                            tACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCgg             >  2:19266/1‑90 (MQ=33)                                                                             aCCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGt            >  2:487293/1‑90 (MQ=33)                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  1:972244/90‑1 (MQ=33)                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  2:4365/90‑1 (MQ=33)                                                                              ccGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTa           <  1:487293/90‑1 (MQ=33)                                                                                 gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:585683/1‑90 (MQ=35)                                                                                 gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  2:801804/1‑90 (MQ=33)                                                                                 gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa        >  1:234944/1‑90 (MQ=33)                                                                                 gTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCCCATCCCGGAAGGGGTTCTCCGCACGGGTGAGGAGCGCACCCACCGGTACCa        >  1:900186/1‑90 (MQ=11)                                                                                       ccGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGa  <  2:951395/90‑1 (MQ=33)                                                                                        |                                                                                         GACCCGGAGCTGCACCCGGACAGCGCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGA  >  NZ_CP009273/727772‑727947

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