New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	3651091 =	NA (NA)	67 (1.040)	58/200	0.1	NA	noncoding (1/1443 nt)	IS150	repeat region
	?	REL606	3697154 =	NA (NA)					coding (278/4134 nt)	rhsA	rhsA element core protein RshA

GGTGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAA  >  REL606/3697054‑3697166                                                                                                     |             ggTGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCgg              <  1:458120/101‑1 (MQ=40)  gTGAACCCGCCCTTGCCCTGCCCGGCCCGCTGCCGTTCACCCTCCCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGc             <  2:859311/101‑1 (MQ=14)   tGAACCCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc            <  2:354767/101‑1 (MQ=255)   tGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc            >  1:656294/1‑101 (MQ=40)   tGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcc            >  2:720422/1‑101 (MQ=40)    gAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcct           <  2:1581706/101‑2 (MQ=40)    gAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcct           >  1:1304378/1‑100 (MQ=255)        ccGTCCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtac       <  1:515319/101‑6 (MQ=255)        ccGACCTTGCCCTGCCCGGCCCGCTGCCGTTCGTCCTCTCCCGCTCCGACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtac       <  2:934216/101‑6 (MQ=14)        ccGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtac       <  1:106502/101‑6 (MQ=255)         cGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtact      <  1:856289/101‑7 (MQ=255)          gACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactg     >  1:1466097/1‑94 (MQ=255)           aCCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCCCCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactga    <  1:247422/101‑9 (MQ=14)            ccTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactgac   >  1:834712/1‑92 (MQ=255)             cTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactgacc  >  2:638848/1‑91 (MQ=255)             cTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGcctgtactgacc  >  1:168788/1‑91 (MQ=255)                                                                                                     |             GGTGAAACCGACCTTGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAA  >  REL606/3697054‑3697166

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