New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	5619030 =	39 (0.890)	3 (0.070)	3/328	NT	8.0%	coding (279/606 nt)	PP_4937	toluene tolerance protein
	?	AE015451	5619049 =	32 (0.760)					coding (298/606 nt)	PP_4937	toluene tolerance protein

CACATTGCGTCGAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGCGGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5619270‑5619030 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggcgaACGCTGCTGCAGATGTCGCGCGATGACGG  >  AE015451/5619049‑5619077                                                                                                                                                                                                                                                                                CACATTGCGTCGAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGC                                     <  1:191565/235‑1       GCGTCGAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGC                                     <  2:630808/229‑1         GTCGAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGGAGCGGGCGA                               >  2:688334/1‑233         GTCGAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGAGACATCTGCAGCAGCGTTCGCCCGGGTAGCGGGC                                 >  1:797944/1‑231            GAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGC                                     <  2:687463/224‑1                   GCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGC                                     <  1:489998/217‑1                                 GCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGCGGGCGAACGCTG                         <  2:540671/215‑1                                                                                        GGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGC                                     <  2:609911/148‑1                                                                                                                                  ACCGCGGCGGGGCAGCTGACGGATCAGTGCCAGCAGGCGTGGCAGGTAAGTGTCCCAGCGGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGC                                     <  2:544482/106‑1                                                                                                                                       TGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGCGGGCGAACGCTGCTGCAGATGTCGCGCGATGACGG  >  1:655061/1‑136                                                                                                                                       TGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGCGGGCGAACGCTGCTGCAGATGTCGCGCGATGACGG  <  2:655059/136‑1                                                                                                                                                                                                                                                                                CACATTGCGTCGAATCATGCGGGCCGACAGTGGCGGGCGCAAAAAGCGCATGTCGGCCACGTCGATCAAGCCCAGACGCTGATCGGGGGTGACCACCACGTTGCCCAGGTGCAGCGAGCGGAAATAAACACCGCTGCGGTGCAGCTGACGGATCAGTTCGATCAGGCGTGGCAGGTAAGTGTCCCAGCTGAAGCCGTCATCGCGCGACATCTGCAGCAGCGTTCGCCCGGGTAGCGGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/5619270‑5619030 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgggcgaACGCTGCTGCAGATGTCGCGCGATGACGG  >  AE015451/5619049‑5619077

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