New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	= 822945	6 (0.220)	4 (0.160)	4/194	2.0	35.7%	coding (421/885 nt)	h16_B0728	predicted hydrolase/acyltransferase (alpha/beta‑hydrolase superfamily)
	?	AM260480	= 822981	9 (0.360)					coding (457/885 nt)	h16_B0728	predicted hydrolase/acyltransferase (alpha/beta‑hydrolase superfamily)

CAGCGGCTGGACCGGCTGCTGCAGGCACTGGGCGTGCAGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGAGGCGCAGGCAAGGGAAGTGGCGCGGCAGCGGCTGGAGGCGCTGGAATCCCTGGGCGTTGATGGCCTGGCGCAACGCGGCCCG  >  AM260480/822828‑823076                                                                                                                      |                                                                                                                                    cAGCGGCTGGACCGGCTGCTGCAGGCACTGGGCGTGCAGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCcgccgtagccttgcgc                                                                                                                     >  1:152996‑M1/1‑118 (MQ=255)                        ggCACTGGGCGTGCAGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCcgccgtagccttgcgccacgctgagcaacagcaggcgcg                                                                                              <  2:152996‑M1/134‑40 (MQ=255)                                  gTGCAGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGa                                                                                    <  1:807303/134‑1 (MQ=255)                                  gTGCAGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGa                                                                                    <  2:602671/134‑1 (MQ=255)                                      aGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCcgccgtagccttgcgccacgctgagcaacagcaggc                                                                                                 >  1:238906‑M1/1‑81 (MQ=255)                                      aGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCcgccgtagccttgcgccacgctgagcaacagcaggc                                                                                                 <  2:238906‑M1/117‑37 (MQ=255)                                       gCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGAGgcgc                                                                               <  1:22815/134‑1 (MQ=255)                                       gCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGAGgcg                                                                                >  2:22815/1‑133 (MQ=255)                                                                                               ccAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGAGGCGCAGGCAAGGGAAGTGGCGCGGCAGCGGCTGGAGGCGCTGGAATCCCTGGGCGTTGAt                       <  2:563575/134‑1 (MQ=255)                                                                                                        cctcccTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGAGGCGCAGGCAAGGGAAGTGGCGCGGCAGCGGCTGGAGGCGCTGGAATCCCTGGGCGTTGATGGCCTGgcg              >  1:377722/6‑134 (MQ=255)                                                                                                                    cgcgcCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGgcgcgggcagctggcgcgcgttggcaacataggccgccgccatcagcgcgcccagcgaatgcccgacgatcaacgccggctgcacgcccagtgcc  <  2:152996‑M2/134‑96 (MQ=255)                                                                                                                      |                                                                                                                                    CAGCGGCTGGACCGGCTGCTGCAGGCACTGGGCGTGCAGCCGGCGTTGATCGTCGGGCATTCGCTGGGCGCGCTGATGGCGGCGGCCTATGTTGCCAACGCGCGCCAGCTGCCCGCGCGCCTGCTGTTGCTCAGCGTGGCGCAAGGCTACGGCGGCGCGGGCAGGGAGGCGCAGGCAAGGGAAGTGGCGCGGCAGCGGCTGGAGGCGCTGGAATCCCTGGGCGTTGATGGCCTGGCGCAACGCGGCCCG  >  AM260480/822828‑823076

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