New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	W3110S.gb	= 1499413	31 (1.360)	3 (0.150)	3/92	NT	9.5%	coding (814/1626 nt)	mdoD	glucan biosynthesis protein, periplasmic
	?	W3110S.gb	= 1499426	29 (1.410)					coding (827/1626 nt)	mdoD	glucan biosynthesis protein, periplasmic

AGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  W3110S.gb/1499370‑1499413 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttcaTTAGTACCGCAGCTGAAC  <  W3110S.gb/1499426‑1499409                                                                AGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCATTAGTACCGCAGCTGAAC  >  1:1507302/1‑62     ATTAAACAGCTGGGCATTGCGCCGATGACCAGTATG                        <  1:1799159/36‑1           CAGCTGGGCATTGCGCCGATGACCAGTATGTTCATTAGTACCG           >  1:160412/1‑43                GGGCATTGCGCCGATGACCAGTATGTTCATTAGTACCGC          <  1:432036/39‑1                                                                AGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  W3110S.gb/1499370‑1499413 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttcaTTAGTACCGCAGCTGAAC  <  W3110S.gb/1499426‑1499409

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