New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	3991398 =	23 (1.050)	3 (0.150)	3/156	1.6	12.5%	coding (56/2163 nt)	uvrD	DNA helicase II
	?	NZ_CP009273	3991437 =	21 (1.050)					coding (95/2163 nt)	uvrD	DNA helicase II
Rejected: Frequency below/above cutoff threshold.

GTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC  >  NZ_CP009273/3991349‑3991524                                                                                         |                                                                                        gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                                                        <  2:64593/90‑1 (MQ=255) gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                                                        <  2:26506/90‑1 (MQ=255) ctcatcaaccaggcgatacgatgcaccagtacgcgcgtcttaccactgcccgcgcccgccagcaccagaaggttgctgcgtggcgcggGc                                                                                        <  2:496264‑M2/2‑1 (MQ=255)          cTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                                                               <  1:147786/90‑1 (MQ=255)          cTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                                                               <  1:230697/90‑1 (MQ=255)            gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGGAGTGGGAAGAc                                                                             >  1:541292/1‑90 (MQ=255)            gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAc                                                                             >  2:13635/1‑90 (MQ=255)                   aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                                                        <  2:48974/72‑1 (MQ=255)                   aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                                                        >  1:48974/1‑72 (MQ=255)                      tgcaccagtacgcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCCGTGGGAAGACGCGCGTACTg                                                                   <  2:276450‑M1/62‑1 (MQ=255)                      cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGGACTg                                                                   >  2:559210/1‑90 (MQ=255)                       gcaccagtccgcgcgtcttacccctgcccgcgcccgccagcaccagaaggttgctgcgtggcgcggGCAGTGGTAAGACGCGCGTACTgg                                                                  <  1:146677‑M2/24‑1 (MQ=255)                          ccagtacgcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGGACTGGTGc                                                               <  1:146677‑M1/66‑1 (MQ=255)                           cagtacgcgcgtcttcccacggcccgcgcccgccagcaccagaaggttgctgcgtggcgcggGCAGTGGTAAGACGCGCGTACTGGTGCa                                                              <  2:276450‑M2/28‑1 (MQ=255)                                    gcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGc                                                     >  2:85435/1‑90 (MQ=255)                                    gcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAAACGCGCGTACTGGTGCATCGTATCGc                                                     >  2:200060/1‑90 (MQ=255)                                    gcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGGAAGACGCGCGGACTGGTGCATCGTATCGc                                                     >  1:68151/1‑90 (MQ=255)                                    gcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGGAAGACGCGCGGACTGGTGCATCGTATCGc                                                     >  1:554076/1‑90 (MQ=255)                                      gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCt                                                   <  2:581340/90‑1 (MQ=255)                                        gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGCGCTGGCGGGCGGGGGGAGGGGTAAAACGCGCCGACTGGGGCATCGGATCGCCTgg                                                 >  1:229157/1‑90 (MQ=255)                                                gcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAg                                         <  2:496264‑M1/88‑1 (MQ=255)                                                      GCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGa                                   <  2:113683‑M1/90‑1 (MQ=255)                                                              gcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCt                           <  1:113683/90‑1 (MQ=255)                                                                      cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGGGCGTGGGAAACTGCTCGCCATAc                   >  2:421199/1‑90 (MQ=255)                                                                         gTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGGTGATGAGCGTGGAAAACTGCTCGCCATACTCg                >  1:423266/1‑90 (MQ=255)                                                                           gCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg                                                                  >  2:302019/1‑38 (MQ=255)                                                                           gCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg                                                                  <  1:302019/38‑1 (MQ=255)                                                                            cTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAtt             >  1:33463/1‑90 (MQ=255)                                                                                   gcgcggGCAGTGGGAAAACGCGCGTACTGGTGCATCGGATCGCCCGGTTGATGAGGGTGGAAAACTGCTCGCCATACTCGATTATGGCgg      >  1:515017‑M2/7‑90 (MQ=255)                                                                                       ggGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAc  <  2:146677‑M2/88‑1 (MQ=255)                                                                                         |                                                                                        GTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAC  >  NZ_CP009273/3991349‑3991524

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