New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 2283591	NA (NA)	29 (1.560)	16/166	0.2	100%	intergenic (‑67/‑388)	BW25113_RS11435/narP	IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP
	?	NZ_CP009273	= 3512441	0 (0.000)					intergenic (‑18/+383)	aroK/hofQ	shikimate kinase AroK/DNA uptake porin HofQ

GGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCC  >  NZ_CP009273/2283507‑2283672                                                                                     |                                                                                  ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATaa                                                                              <  1:249674/90‑1 (MQ=35)      gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCACATCAGGGACTTGTTCGCACCTTCCATAACGCTg                                                                         <  2:426388/90‑1 (MQ=255)            gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                                           >  1:406572/1‑82 (MQ=38)            gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                                           <  2:406572/82‑1 (MQ=38)                    atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACa                                                           <  1:538053/90‑1 (MQ=255)                        aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGata                                                       >  1:385275/1‑90 (MQ=255)                             gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCg                                                  >  2:489670/1‑90 (MQ=255)                                  aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACg                                             >  2:374618/1‑90 (MQ=255)                                     gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACa                                          >  1:461760/1‑90 (MQ=255)                                       tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACaaa                                        >  1:512638/1‑90 (MQ=255)                                          aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                                     <  1:489670/90‑1 (MQ=255)                                              gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                 <  2:168986/90‑1 (MQ=255)                                              gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                 <  2:66631/90‑1 (MQ=255)                                                      tGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAg                         <  1:374618/90‑1 (MQ=255)                                                        aTCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGAt                       <  2:512638/90‑1 (MQ=255)                                                                 cAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAtt              >  1:44658/1‑90 (MQ=255)                                                                      aCTTGTTCGCACCTTCCATAACGCTGTATCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCt         >  1:511794/1‑90 (MQ=255)                                                                        ttGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCttt       >  1:124830/1‑90 (MQ=255)                                                                          gTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCa     <  1:492835/90‑1 (MQ=255)                                                                            tCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATc   <  2:385275/90‑1 (MQ=255)                                                                             cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATcc  <  2:108530/90‑1 (MQ=255)                                                                                     |                                                                                  GGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCC  >  NZ_CP009273/2283507‑2283672

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.