New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AP009048	= 2033398	44 (0.750)	3 (0.060)	3/92	NT	7.2%	coding (1057/1419 nt)	dcm	DNA cytosine methylase
	?	AP009048	= 2033414	38 (0.720)					coding (1041/1419 nt)	dcm	DNA cytosine methylase

TCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/2033325‑2033398 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttcggTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACAAAGAT  <  AP009048/2033414‑2033345                                                                                                                                                  TCGATTAATATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTT                                                                           <  1:1860777/71‑1                     ATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGTTATGGAATGTTTTATC                                                       <  1:685635/71‑1                              ATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGGTTATGGAATGGTTTATCCGAA                                                   >  1:2218021/1‑66                                                ACGCGTGACGCTTTGCGGATTGTTCGGTTATGGAATGGTTTATCCG                                                     <  1:4388668/46‑1                                                                      TTCGGTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGT                             >  1:2134548/1‑48                                                                        CGGTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACA       >  1:263484/1‑68                                                                        CGGTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTAAGCTGTCTGCGCGTTATTAC        >  1:4090979/1‑67                                                                         GGTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACAAAG    <  1:2645015/70‑1                                                                          GTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACAAAGAT  <  1:2277160/71‑1                                                                          GTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACAAAGAT  <  1:2283295/71‑1                                                                          GTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACAAAGAT  <  1:4875977/71‑1                                                                                                                                                  TCGATTAAAATTTCCGCGCCATCTTTGTAATAACGCGCAGACAGCGTACGCGTGACGCTTTGCGGATTGTTCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/2033325‑2033398 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttcggTTATGGAATGGTTTATCCGAACAATCCGCAAAGCGTCACGCGTACGCTGTCTGCGCGTTATTACAAAGAT  <  AP009048/2033414‑2033345

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.