New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 3652533	NA (NA)	17 (0.260)	13/200	4.5	100%	noncoding (1443/1443 nt)	IS150	repeat region
	?	REL606	= 3697156	0 (0.000)					coding (280/4134 nt)	rhsA	rhsA element core protein RshA
Rejected: Coverage evenness skew score above cutoff.

CAACTTTTTGGGGTCAGTACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/3652513‑3652533 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTCGGTTTC  <  REL606/3697156‑3697057                                                                                                                           CAACTTTTTGGGGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCG                      >  1:504652/1‑101   ACTTTTTGGGGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGG                    >  1:517449/1‑101   ACTTTTTGGGGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGG                    <  2:1180504/101‑1      TTTTGGGGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAG                 <  2:1017477/101‑1      TTTTGGGGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAG                 <  1:161115/101‑1         TGGGGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGC              <  2:679562/101‑1            GGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAG           <  2:638848/101‑1            GGTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAG           <  1:168788/101‑1             GTCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGG          <  1:834712/101‑1              TCAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGGGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGT         >  1:247422/1‑101               CAGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTC        <  1:1466097/101‑1                AGTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTCG       >  1:856289/1‑101                 GTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTCGGAGCGGGAGAGGACGAACGGCAGCGGGCCGGGCAGGGCAAGGTCGG      >  2:934216/1‑101                 GTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTCGG      >  1:106502/1‑101                 GTACAGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGACGG      >  1:515319/1‑101                     AGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTCGGTTTC  >  2:1581706/1‑101                     AGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTCGGTTTC  <  1:1304378/101‑1                                                                                                                           CAACTTTTTGGGGTCAGTACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/3652513‑3652533 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCAAGGTCGGTTTC  <  REL606/3697156‑3697057

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