New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3630716	NA (NA)	65 (0.820)	50/100	0.3	NA	intergenic (‑171/‑305)	ldrD/yhjV	toxic polypeptide, small/predicted transporter
	?	REL606	= 3652533	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region

GCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGCAC  >  REL606/3630697‑3630766                    |                                                   gCAACATTGCGCCTCCCCGGAACCTTTCCGGCTTGCCGCTGATTCTACGAt                     <  1:5167440/51‑1 (MQ=9) gCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGAt                     <  1:372756/51‑1 (MQ=31) gCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGAt                     <  1:6317008/51‑1 (MQ=31)  cAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGAtt                    >  1:1818805/1‑51 (MQ=33)   aaCATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTg                   <  1:2795228/51‑1 (MQ=35)   aaCATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTg                   >  1:2435712/1‑51 (MQ=35)   aaCATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTg                   >  1:3139687/1‑51 (MQ=35)    aCATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGc                  >  1:6654003/1‑51 (MQ=38)     cATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCt                 >  1:3858145/1‑51 (MQ=40)       ttGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCttt               <  1:7067225/51‑1 (MQ=255)        tGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCtttt              <  1:6190500/51‑1 (MQ=255)        tGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCtttt              <  1:1710810/51‑1 (MQ=255)          cgcCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTa            >  1:1924666/1‑51 (MQ=255)          cgcCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTa            <  1:1245441/51‑1 (MQ=255)           gcCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTAc           >  1:5760219/1‑51 (MQ=255)           gcCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTAc           <  1:2433324/51‑1 (MQ=255)              tCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGt        <  1:5432404/51‑1 (MQ=255)               cACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTc       >  1:7251054/1‑51 (MQ=255)               cACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTc       <  1:1054840/51‑1 (MQ=255)                aCCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCa      <  1:5820030/51‑1 (MQ=255)                aCCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCa      <  1:1469008/51‑1 (MQ=255)                  cGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGc    >  1:5211195/1‑51 (MQ=255)                   ggAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGCa   >  1:1426441/1‑51 (MQ=255)                   ggAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGCa   <  1:4372062/51‑1 (MQ=255)                   ggAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGCa   <  1:401643/51‑1 (MQ=255)                    gAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGCAc  <  1:6971700/51‑1 (MQ=255)                    |                                                   GCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGATTGCTTTTTACTGTCAGCAC  >  REL606/3630697‑3630766

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