New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5299381 =	39 (0.890)	3 (0.070)	3/326	NT	7.7%	coding (721/2076 nt)	recD	ATP‑dependent DNA helicase and dsDNA/ssDNA exonuclease, ssDNA endonuclease
	?	AE015451	5299407 =	35 (0.840)					coding (695/2076 nt)	recD	ATP‑dependent DNA helicase and dsDNA/ssDNA exonuclease, ssDNA endonuclease

TACGAGCGGCGCATCGACCAGGCCCTGCGCCAGCGCCTGAGTCAGGCCGAAGCACCGCCTGCGGACCTGCCAGAGCGTCTGGCGCAACTGTTCGAGGGCGGCGCGCCAGCTGGCCAGGTGGACTGGCAGAAGCTCGCCTGTGCCTTGGCCACCCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5299618‑5299381 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcaggcgcACCACGGTGGTCGTCTTGCCAGTACCGGGGCCGCCAGTGATGATGCTGAAACCGGCACGGGT  >  AE015451/5299407‑5299468                                                                                                                                                                                                                                                                                                              TACGAGCGGCGCATCGACCAGGCCCTGCGCCAGCGCCTGAGTCAGGCCGAAGCACCGCCTGCGGACCTGCCAGAGCGTCTGGCGCAACTGTTCGAGGGCGGCGCGCCAGCTGGCCAGGTGGACTGGCAGAAGCTCGCCTGTGCCTTGGCCACCCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGC                                                                  <  1:474963/236‑1            ccTCGACCAGGCCCTGCGCCAGCGCCTGAGTCAGGCCGAAGCACCGCCTGCGGCCCTGCCAGAGCGTCTGGCGCAACTGTTCGAGGGCGGCGCGCCAGCTGGCCAGGTGGACTGGCAGAAGCTCGCCTGTGCCTTGGCCACCCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGCGC                                                                <  1:774134/225‑1                                                                      CCAGAGCGGCGGGCGCAACTGTTCGAGGGCGGCGCGCCAGCTGGCCAGGTGGACTGGCAGAAGCTCGCCTGTGCCTTGGCCACCCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGC                                                                  <  2:127989/167‑1                                                                                           TTCGGGGGCGGCGCGCCAGCTGGCCAGGTGGACTGGCAGAAGCTCGCCTGTGCCTTGGCCACCCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGCGCACCACGGTGGTCGTCTT                                               <  2:238121/165‑1                                                                                                                                                         CCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGGTGCTGCAGGCGCACCACGGTGGTCGTCTTGCCAGTACCGGGGCCGCCAGTGATGATGCTGAAACCGGCACGGGT  <  1:371576/148‑1                                                                                                                                                                                 gGGCGGCCCCGGTCCTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGCGCACCACGGTGGTCGTCTTGCCAGTACCGGGGCCGCCAGTG                         <  1:438816/100‑1                                                                                                                                                                                                                                          GGCGCACCACGGTGGTCGTCTTGCCAGTACCGGGGCCGCCAGTGATGATGCc                 >  2:739900/1‑51                                                                                                                                                                                                                                                                                                              TACGAGCGGCGCATCGACCAGGCCCTGCGCCAGCGCCTGAGTCAGGCCGAAGCACCGCCTGCGGACCTGCCAGAGCGTCTGGCGCAACTGTTCGAGGGCGGCGCGCCAGCTGGCCAGGTGGACTGGCAGAAGCTCGCCTGTGCCTTGGCCACCCGTGCCGGTTTCAGCATCATCACTGGCGGCCCCGGTACTGGCAAGACGACCACCGTGGTGCGCCTGCTGGCGTTGCTGCAGGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5299618‑5299381 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcaggcgcACCACGGTGGTCGTCTTGCCAGTACCGGGGCCGCCAGTGATGATGCTGAAACCGGCACGGGT  >  AE015451/5299407‑5299468

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