New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3652533	NA (NA)	46 (1.250)	34/146	0.1	100%	noncoding (1443/1443 nt)	IS150	repeat region
	?	REL606	3901931 =	0 (0.000)					coding (489/1428 nt)	yieO	predicted multidrug or homocysteine efflux system

TCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCA  >  REL606/3652465‑3652605                                                                     |                                                                         tCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                                                                   >  1:831097/1‑76 (MQ=39)  ccAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGa                                                                  >  1:862715/1‑76 (MQ=40)   cAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                                                                 >  2:290429/1‑76 (MQ=255)    aaTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAtt                                                                >  1:743263/1‑76 (MQ=255)         aaTATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAg                                                           <  2:510369/76‑1 (MQ=255)          atatCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                                                          >  2:236474/1‑76 (MQ=255)                gAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTg                                                    >  1:768959/1‑76 (MQ=255)                    cAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGc                                                >  2:815184/1‑76 (MQ=255)                       aCCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCttt                                             <  1:339018/76‑1 (MQ=255)                         cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc                                           >  1:978972/1‑76 (MQ=255)                          tataTGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtca                                          <  2:1030516/76‑2 (MQ=255)                          natatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtct                                          >  1:218610/2‑76 (MQ=255)                             atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctg                                       >  1:800899/1‑76 (MQ=255)                                ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtct                                    >  1:183128/1‑76 (MQ=255)                                 cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtctg                                   <  1:592728/76‑1 (MQ=255)                                 cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtctg                                   <  1:173402/76‑1 (MQ=255)                                   cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGcc                                 >  2:586347/1‑76 (MQ=255)                                   cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGcc                                 >  2:586342/1‑76 (MQ=255)                                     tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCtt                               <  2:1088632/76‑1 (MQ=255)                                      gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt                              <  2:779422/76‑1 (MQ=255)                                       tttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc                             >  1:519392/1‑76 (MQ=255)                                         tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg                           >  1:1104171/1‑76 (MQ=255)                                          aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt                          <  2:665730/76‑1 (MQ=255)                                          aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt                          >  1:761204/1‑76 (MQ=255)                                                 cAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATAcc                   >  2:217373/1‑76 (MQ=255)                                                 cAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATAcc                   >  2:217368/1‑76 (MQ=255)                                                  aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCt                  <  2:378977/76‑1 (MQ=255)                                                     tttttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTgag               <  2:471722/76‑1 (MQ=255)                                                     tttttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTgag               <  1:626669/76‑1 (MQ=255)                                                        ttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACa            <  2:3000/76‑1 (MQ=255)                                                        tgggggTCAGTACAAATTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACa            <  2:2999/74‑1 (MQ=255)                                                          ggggTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACaaa          <  2:1040553/76‑1 (MQ=255)                                                           gggTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAAc         >  1:392379/1‑76 (MQ=255)                                                              tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGt      <  1:1048447/76‑1 (MQ=255)                                                                 gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGtctc   <  2:482469/76‑1 (MQ=255)                                                                  tACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCa  <  1:279891/76‑1 (MQ=255)                                                                  tACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCa  <  2:717573/76‑1 (MQ=255)                                                                     |                                                                         TCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCA  >  REL606/3652465‑3652605

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