New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4284707 =	29 (0.660)	3 (0.070)	3/324	NT	9.9%	coding (112/849 nt)	hbd	3‑hydroxybutyryl‑CoA dehydrogenase
	?	AE015451	4284752 =	27 (0.650)					coding (157/849 nt)	hbd	3‑hydroxybutyryl‑CoA dehydrogenase

TGGTGGCAATCAGGCAGTCGGCGGCGACGTTGGCTGCCACCTGCTGCAGGATGCGCTGCTTGAGCTGCAGGTTCTCGGTAGCCGCCTCGATCACCAGCTGTGCAGCGCTGAGCTGGGTGTAGTCGGTGCTGGTACGGATGCGGGCTTTGGCGGCTGCGGCCTTGTCGGCGTCGAGAGTGCCCTTGCTGACCTGGCGCTCGAGGTTCTTGCTCAAGGTGGCCACACCGCGCTCGAGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4284944‑4284707 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctcgagcgcCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCTCAGCGCTGCACAGCTGGTGATCGAGGCGGCTACCGAGAACCTGCAGCTCAAGCAGCGCATCCTGCAGCAGGTGGCAGCCAACGTCGCCGCCGACTGCCTGATTGCCACCAATACTTCGTCGCTGTCGGTGACCCAACTGGCCGCCAGCA  >  AE015451/4284752‑4284983                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        TGGTGGCAATCAGGCAGTCGGCGGCGACGTTGGCTGCCACCTGCTGCAGGATGCGCTGCTTGAGCTGCAGGTTCTCGGTAGCCGCCTCGATCACCAGCTGTGCAGCGCTGAGCTGGGTGTAGTCGGTGCTGGTACGGATGCGGGCTTTGGCGGCTGCGGCCTTGTCGGCGTCGAGAGTGCCCTTGCTGACCTGGCGCTCGAGGTTCTTGCTCAAGGTGGCCACACCGCGCTCGAGC                                                                                                                                                                                                                                            <  1:34601/236‑1                                                                                                                                                                                          CTGACCTGGCGCTCGAGGTTCTTGCTCAAGGTGGCCACACCGCGCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCTCAGCGCTGCACAGCTGGTGATCGAGGCGGCTACCGAGAACCTGCAGCTCAAGCAGCGCATCCT                                                                                        >  2:326059/1‑199                                                                                                                                                                                                 GGCGCTCGAGGTTCTGGCTCAAGGTGGCCACACCGCGCTCGAGCGCCAGGTCAGCAAGG                                                                                                                                                                                                                             <  2:226886/59‑1                                                                                                                                                                                                      TCGAGGTTCTTGCTCAAGGTGGCCACACCGCGCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCTCAGCGCTGCACAGCTGGTGATCGAGGCGGCTACCGAGAACCTGCAGCTCAAGCAGCGCATCCTGCAGCAGGTGGCCGCCAACGTCGCCGCCGACTGCCTGA                                                  >  1:225883/1‑225                                                                                                                                                                                                                                       TCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCTCAGCGCTGCACAGCTGGTGATCGAGGCGGCTACCGAGAACCTGCAGCTCAAGCAGCGCATCCTGCAGCAGGTGGCAGCCAACGTCGCCGCCGACTGCCTGATTGCCACCAATACTTCGTCGCTGTCGGTGACCCAACTGGCCGCC      <  1:665242/236‑1                                                                                                                                                                                                                                          AGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCT                                                                                                                                                       <  2:178081/88‑1                                                                                                                                                                                                                                           GCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCTCAGCGCTGCACAGCTGGTGATCGAGGCGGCTACCGAGAACCTGCAGCTCAAGCAGCGCATCCTGCAGCAGGTGGCAGCCAACGTCGCCGCCGACTGCCTGATTGCCACCAATACTTCGTCGCTGTCGGTGACCCAACTGGCCGCCAGCA  >  1:392016/1‑236                                                                                                                                                                                                                                            CGCCAGGTCAGCAAGTGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCG                                                                                                                                                                               <  2:758302/62‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        TGGTGGCAATCAGGCAGTCGGCGGCGACGTTGGCTGCCACCTGCTGCAGGATGCGCTGCTTGAGCTGCAGGTTCTCGGTAGCCGCCTCGATCACCAGCTGTGCAGCGCTGAGCTGGGTGTAGTCGGTGCTGGTACGGATGCGGGCTTTGGCGGCTGCGGCCTTGTCGGCGTCGAGAGTGCCCTTGCTGACCTGGCGCTCGAGGTTCTTGCTCAAGGTGGCCACACCGCGCTCGAGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4284944‑4284707 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctcgagcgcCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGCCGCCAAAGCCCGCATCCGTACCAGCACCGACTACACCCAGCTCAGCGCTGCACAGCTGGTGATCGAGGCGGCTACCGAGAACCTGCAGCTCAAGCAGCGCATCCTGCAGCAGGTGGCAGCCAACGTCGCCGCCGACTGCCTGATTGCCACCAATACTTCGTCGCTGTCGGTGACCCAACTGGCCGCCAGCA  >  AE015451/4284752‑4284983

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