New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP000731	17008 =	NA (NA)	27 (0.080)	14/148	NT	100%	coding (54/675 nt)	USA300HOU_pUSA300HOUMR0021	IS431 mec transposase
	?	CP000731	17011 =	0 (0.000)					coding (51/675 nt)	USA300HOU_pUSA300HOUMR0021	IS431 mec transposase

GTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAG  >  CP000731/16918‑17018                                                                                           |           atctacCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCa            >  1:1552086‑M2/7‑91 (MQ=255) atctacCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCa            >  1:3485693‑M2/7‑91 (MQ=255) atctacCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCa            >  1:1078814‑M2/7‑91 (MQ=255) atctacCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCa            >  1:2867134‑M2/7‑91 (MQ=255)  tctacCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCaa           >  1:1953023‑M2/6‑91 (MQ=255)   ctacCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAAc          >  1:3402883‑M2/5‑91 (MQ=255)       CGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATATTAGCCAACGGCt      <  1:3189572‑M2/91‑1 (MQ=255)       CGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCt      <  1:2956095‑M2/91‑1 (MQ=255)       CGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCt      <  1:3497575‑M2/91‑1 (MQ=255)       CGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCt      <  1:2955792‑M2/91‑1 (MQ=255)       CGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCt      <  1:1889281‑M2/91‑1 (MQ=255)        GTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTa     <  1:2258702‑M2/91‑1 (MQ=255)        GTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTa     <  1:3624631‑M2/91‑1 (MQ=255)         TTGAATGATTAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    <  1:3963527‑M2/91‑1 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    <  1:3282760‑M2/91‑1 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    >  1:3984313‑M2/1‑91 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    >  1:3982714‑M2/1‑91 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    >  1:3984504‑M2/1‑91 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    >  1:3984517‑M2/1‑91 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    <  1:4063702‑M2/91‑1 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTAc    <  1:1960116‑M2/91‑1 (MQ=255)         TTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCAAACGGCTAc    >  1:3983994‑M2/1‑91 (MQ=255)          TGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACa   >  1:362202‑M2/1‑91 (MQ=255)          TGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACa   >  1:3163888‑M2/1‑91 (MQ=255)          TGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACa   >  1:1988958‑M2/1‑91 (MQ=255)          TGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACa   >  1:492457‑M2/1‑91 (MQ=255)           GAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTCGCCAACGGCTACAg  >  1:3741871‑M2/1‑91 (MQ=255)           GAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAg  >  1:3242372‑M2/1‑91 (MQ=255)           GAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAg  >  1:3034055‑M2/1‑91 (MQ=255)           GAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAg  <  1:2641970‑M2/91‑1 (MQ=255)           GAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAg  <  1:4047810‑M2/91‑1 (MQ=255)                                                                                           |           GTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAG  >  CP000731/16918‑17018

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