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.

TCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAG  >  NZ_CP009273/3991309‑3991485                                                                                          |                                                                                        tctcCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGc                                                                                         <  2:446577/90‑1 (MQ=255)  gtcaccgccataatcgagtatggcgagcagttttccacgctcatcaaccaggcgatacgatgcaccagtacgcgcgtcttaccactgccc                                                                                        <  2:146677‑M1/2‑1 (MQ=255)      ccgccataatcgagtatggcgagcagttttccaccctcatcaaccgggcgatccgatgcaccagtacgcgcgttttcccactgcCCgcgc                                                                                    >  1:515017‑M1/85‑90 (MQ=255)        ccAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCa                                                                                  <  1:85437/90‑1 (MQ=255)         cAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCAc                                                                                 <  2:106155/90‑1 (MQ=255)           aCCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgc                                                                               >  2:156297/1‑90 (MQ=255)                     aCGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTc                                                                     >  2:86246/1‑90 (MQ=255)                         ggcggTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                                 >  1:581340/1‑90 (MQ=255)                         ggcggTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                                 >  1:64593/1‑90 (MQ=255)                                       aCGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgg                                                   <  1:489991/90‑1 (MQ=255)                                         gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                 <  2:64593/90‑1 (MQ=255)                                         gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                 <  2:26506/90‑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                                                 >  1:48974/1‑72 (MQ=255)                                                           aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgggc                                                 <  2:48974/72‑1 (MQ=255)                                                              tgcaccagtacgcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCCGTGGGAAGACGCGCGTACTg                            <  2:276450‑M1/62‑1 (MQ=255)                                                              cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGGACTg                            >  2:559210/1‑90 (MQ=255)                                                                  ccagtacgcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGGACTGGTGc                        <  1:146677‑M1/66‑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)                                                                                          |                                                                                        TCTCCGCCCAACCTATTTTTACGCGGCGGTGCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAG  >  NZ_CP009273/3991309‑3991485

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