New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	2416916 =	21 (0.900)	3 (0.140)	3/154	1.8	15.8%	coding (103/891 nt)	rpnB	recombination‑promoting nuclease RpnB
	?	NZ_CP009273	2416950 =	13 (0.620)					coding (69/891 nt)	rpnB	recombination‑promoting nuclease RpnB
Rejected: Frequency below/above cutoff threshold.

ATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCA  >  NZ_CP009273/2416861‑2417037                                                                                          |                                                                                        aTAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgcg                                                                                         <  2:428026/90‑1 (MQ=255)   cttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGgt                                                                                       <  2:439145‑M1/37‑1 (MQ=255)         gTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTgg                                                                                 >  1:129241/1‑90 (MQ=255)          tttGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGa                                                                                >  1:364733/1‑90 (MQ=255)           ttGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGAt                                                                               <  1:449082/90‑1 (MQ=255)           ccgcatgatgcggtatttaaatcttttttacgccatccagacaccgcgcgggattttattgatattcatcttcccgcgcGGTGTCTGGAt                                                                               <  1:49710‑M2/11‑1 (MQ=255)                     tttAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTaaaaa                                                                     <  2:606881/90‑1 (MQ=255)                     tttAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTaaaaa                                                                     <  2:236482/90‑1 (MQ=255)                       gtatttaaatcttttttacgccatccagacaccGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGCTGGCGTAAAAAAg                                                                   <  1:242554‑M1/57‑1 (MQ=255)                                 cttttttacgccagccagacaccgcgcgggattttattgatattcatcttcccgcgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATAc                                                         <  1:242554‑M2/33‑1 (MQ=255)                                  aaaTCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATAcc                                                        >  2:116414/1‑90 (MQ=255)                                            gTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgcg                                                                                         >  1:232507/1‑47 (MQ=255)                                            gTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgcg                                                                                         <  2:232507/47‑1 (MQ=255)                                             atccagacaccGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATgcgg                                             <  1:49710‑M1/79‑1 (MQ=255)                                              ttGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATgcggc                                            >  2:329703/1‑90 (MQ=255)                                                    aGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGaa                                      >  2:93361/1‑90 (MQ=255)                                                     gcggcgCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAg                                     >  2:226306/1‑90 (MQ=255)                                                     accgcgcgggattttattgatattcatcttcccgcgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAg                                     <  2:439145‑M2/53‑1 (MQ=255)                                                      cggcgCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGt                                    >  2:315714/1‑90 (MQ=255)                                                                                    ccgcgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCAccc      >  2:221524‑M2/7‑90 (MQ=255)                                                                                    ccgcgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCAccc      >  2:357077‑M2/7‑90 (MQ=255)                                                                                    ccgcgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCAccc      >  1:279565‑M2/7‑90 (MQ=255)                                                                                      gcgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCt    >  2:456774‑M2/5‑90 (MQ=255)                                                                                       cgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTc   <  1:442250‑M2/87‑1 (MQ=255)                                                                                       cgcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTc   <  1:361536‑M2/87‑1 (MQ=255)                                                                                        gcGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCa  <  1:125332‑M2/88‑1 (MQ=255)                                                                                          |                                                                                        ATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCA  >  NZ_CP009273/2416861‑2417037

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