New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	3585335 =	61 (0.990)	3 (0.050)	3/326	NT	5.4%	coding (348/762 nt)	benD	1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase
	?	AE015451	3585374 =	47 (0.800)					coding (387/762 nt)	benD	1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase

GCGTTGACGCGGATGCCCCGGCCGGCGGTTTCGAAGGCCAGGCAGGCAGTGAGCGCGTTGACGCCACCCTTGGCTGCGCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTGGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/3585511‑3585335 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgtggcagTGGCGCCATCGTCAATGTCTCGTCCATCGCTACCCGCGGCGTCAACCGTGTGCCTTACGGCGCAGCCAAGGGTGGCGTCAACGCGCTCACTGCCTGCCTGGCCTTCGAAACCGCCGGCCGGGGCATCCGCGTCAACGCCACTGCGCCCGGCGGCACCGAAGCACCGCCCCGGCGCGTA  >  AE015451/3585374‑3585551                                                                                                                                                                                                                                                                                                                                                                     GCGTTGACGGGGAGGCCCCGGCCGGCGGTTTCGAAGGCCAGGCAGGCAGTGAGCGCGTTGACGCCACCCTTGGCTGCGCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTGGCA                                                                                                                                                                                     <  2:194463/176‑1    TTGACGCGGATGCCCCGGCCGGCGGTTTCGAAGGCCAGGCAGGCAGTGAGCGCGTTGACGCCACCCTTGGCTGCGCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGG                                                                                                                                                                                            <  1:261119/166‑1                         GCGTTTTCGAAAGCGAGGCAGGCAGTGAGCGCGTTGACGCCTCCCTTGGCTGCGCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTG                                                                                                                                                                                        <  2:290481/149‑1                              tgCGAAGGCGAGGCAGGCAGTGAGCGCGTTGACGCCACCCTTGGCTGCGCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTGGCAGTG                                                                                                                                                                                  <  2:174595/148‑1                                                                              GCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTGGCAGTGGCGCCATCGTCAATGTCTCGTCCATCGCTACCCGCGGCGTCAACCGTGTGCC                                                                                                                              <  1:1100291/154‑1                                                                                       CACAGGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGAGGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCG                                                                                                                                                                                          <  2:13093/85‑1                                                                                         CACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTGGCAGTGGCGCCATCGTCAATGTCTCGTCCATCGCTACCCGCGGCGTCAACCGTGTGCCTTACGGCGCA                                                                                                                    <  2:469528/153‑1                                                                                                      GCGGGTGGCGATGGACGAGACATTGACGAGGGCGCCACTGCCACGCTCGAGCATGTCCGGCAGCGCGGCGTGGCAGGGGCGCCATCGTCAATGTCTCGTCCATCGCTACC                                                                                                                                                  <  2:254250/110‑1                                                                                                                                                                 aggcgcgGGCGTGGCAGTGGCGCCATCGTCAATGTCTCGTCCATCGCTACCCGCGGCGTCAACCGTGTGCCTTACGGCGCAGCCAAGGGTGGCGTCAACGCGCTCACTGCCTGCCTGGCCTTCGAAACCGCCGGCCGGGGCATCCGCGTCAACGCCACTGCGCCCGGCGGCACCGAAGCACCGCCCCGGCGCGTA  <  2:722880/188‑1                                                                                                                                                                           GTGGCAGTGGCGCCATCGTCAATGTCTCGTCCATCGCTACCCGCGGCGTCAACCGTGTGCCTTACGGCGCAGCCAAGGGTGGCGTCAACGCGCTCACTGCCTGCCTGGCCTTCGAAACCGCCGGC                                                              >  2:151917/1‑125                                                                                                                                                                                                                                                                                                                                                                     GCGTTGACGCGGATGCCCCGGCCGGCGGTTTCGAAGGCCAGGCAGGCAGTGAGCGCGTTGACGCCACCCTTGGCTGCGCCGTAAGGCACACGGTTGACGCCGCGGGTAGCGATGGACGAGACATTGACGATGGCGCCACTGCCACGCTCGAGCATGTACGGCAGCGCGGCGTGGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/3585511‑3585335 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgtggcagTGGCGCCATCGTCAATGTCTCGTCCATCGCTACCCGCGGCGTCAACCGTGTGCCTTACGGCGCAGCCAAGGGTGGCGTCAACGCGCTCACTGCCTGCCTGGCCTTCGAAACCGCCGGCCGGGGCATCCGCGTCAACGCCACTGCGCCCGGCGGCACCGAAGCACCGCCCCGGCGCGTA  >  AE015451/3585374‑3585551

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