New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	= 2855747	18 (0.540)	10 (0.310)	7/98	NT	25.8%	intergenic (‑285/‑2)	ygbA/mutS	hypothetical protein/methyl‑directed mismatch repair protein
	?	W3110S.gb	4124505 =	40 (1.220)					coding (1/837 nt)	dam	DNA adenine methylase

CGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  W3110S.gb/2855677‑2855747 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGATATTAAA  >  W3110S.gb/4124505‑4124570                                                                                                                                           CGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                                    >  1:1059870/1‑71      TATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                                    <  1:2904546/66‑1                                              TAATATCAGTAAACCGGACATAACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTA                        <  1:1014154/70‑1                                                 TATCAGTAAACCGGACATAACCCATGAAGAAAAATC                                                       >  1:159597/1‑36                                                       taAACCGGACATAACCCATGAAGAAAAATCGCGCTTTTTT                                             <  1:687058/38‑1                                                        aAACCGGACATAACCCATGAAGAAAAATCGCGCTTTTT                                              >  1:2787283/2‑38                                                         AACCGGACATAACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGA            <  1:1247566/71‑1                                                         AACCGGACATAACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGG                                 >  1:1529677/1‑50                                                           CCGGACATAACCCATGAAGAAAAATCGCGCTTTTTT                                             >  1:1546388/1‑36                                                            CGGACATAACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTA                        <  1:2499122/56‑1                                                               ACATAACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGATAT      >  1:1699267/1‑71                                                                 ATAACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGATATTA    <  1:150206/71‑1                                                                    ACCCATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGATATTAAA  <  1:1818491/70‑1                                                                                                                                           CGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  W3110S.gb/2855677‑2855747 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cATGAAGAAAAATCGCGCTTTTTTGAAGTGGGCAGGGGGCAAGTATCCCCTGCTTGATGATATTAAA  >  W3110S.gb/4124505‑4124570

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