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.

AGGTCGGAATAATATTGCCGCAGGTCTTCATCAATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGT  >  NZ_CP009273/2416828‑2417003                                                                                         |                                                                                        tgaggggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccgc                                                                                        <  1:125332‑M1/2‑1 (MQ=255)  gaggggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccgcg                                                                                       <  1:442250‑M1/3‑1 (MQ=255)  gaggggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccgcg                                                                                       <  1:361536‑M1/3‑1 (MQ=255)   aggggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccgcgc                                                                                      >  2:456774‑M1/87‑90 (MQ=255)     gggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccGCGCgg                                                                                    >  1:279565‑M1/85‑90 (MQ=255)     gggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccGCGCgg                                                                                    >  2:221524‑M1/85‑90 (MQ=255)     gggtgacggaacgccatgacaatatcgacaacttccacgccgcatgatgcggtatttaaatcttttttacgccatccagacaccGCGCgg                                                                                    >  2:357077‑M1/85‑90 (MQ=255)       gAATAATATTGCCGCAGGTCTTCATCAATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGa                                                                                  <  2:275504/90‑1 (MQ=255)       caataatATTGCCGCTGGTCTTCATCAATAAAACTATTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGTAGCGGCGCGGGa                                                                                  <  1:482734/89‑1 (MQ=255)                tGCCGCAGGTCTTCATCAATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAata                                                                         <  2:101554/90‑1 (MQ=255)                 gCCGCAGGTCTTCATCAATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAatat                                                                        >  1:428026/1‑90 (MQ=255)                  ccGCAGGTCTTCATCAATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATc                                                                       <  1:193602/90‑1 (MQ=255)                                caATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcg                                                         <  2:401764/90‑1 (MQ=255)                                 aaTAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgc                                                        <  2:583037/90‑1 (MQ=255)                                 aaTAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgc                                                        <  2:550057/90‑1 (MQ=255)                                  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)                                                      tttAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTaaaaa                                   <  2:236482/90‑1 (MQ=255)                                                      tttAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTaaaaa                                   <  2:606881/90‑1 (MQ=255)                                                        gtatttaaatcttttttacgccatccagacaccGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGCTGGCGTAAAAAAg                                 <  1:242554‑M1/57‑1 (MQ=255)                                                                   aaaTCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATAcc                      >  2:116414/1‑90 (MQ=255)                                                                             gTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgcg                                                       <  2:232507/47‑1 (MQ=255)                                                                             gTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCcgcgcg                                                       >  1:232507/1‑47 (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)                                                                                       cggcgCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGt  >  2:315714/1‑90 (MQ=255)                                                                                         |                                                                                        AGGTCGGAATAATATTGCCGCAGGTCTTCATCAATAAAACTGTTTGGTTCCAGTTTAAGCGTCGTTAAATCACACAGTTTGCGCAGCGGCGCGGGAAGATGAATATCAATAAAATCCCGCGCGGTGTCTGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGT  >  NZ_CP009273/2416828‑2417003

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.