New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3582944	43 (0.800)	4 (0.080)	4/248	NT	9.4%	coding (1015/1359 nt)	benA	benzoate 1,2‑dioxygenase subunit alpha
	?	NC_002947	= 3582963	36 (0.700)					coding (1034/1359 nt)	benA	benzoate 1,2‑dioxygenase subunit alpha

TTCGGCCAGGCACGTGCCGACTGGATGATCGAAAACTCGCGCAACCTGTGCCTGTACCCCAACGTGTACCTGATGGACCAGTTCAGCTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3582794‑3582944 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcatcgGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATGGTGATTTCGGTCTTGTTCACCGAAATCGGCCGGGCAACGCGGATCTGCGAGCTGAACTGGTCCATCAGGTACACGTTGGGGTACAGGCACAGGTTGCGCGAGTTTTCGA  <  NC_002947/3582963‑3582822                                                                                                                                                                                                                                                                                                       TTCGGCCAGGCACGTGCCGACTGGATGATCGAAAACTCGCGCAACCTGTGCCTGTACCCCAACGTGTACCTGATGGACCAGTTCAGCTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCAT                                                                                                                                                  >  1:305832/1‑149                                                          CCCAACGTGTACCTGATGGACCAGTTCAGCTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCATCGGCGCTCTCGCCTTTCGGCGCGATGCA                                                                                                                      >  5:158243/1‑120                                                          CCCAACGTGTACCTGATGGACCAGTTCAGCTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCATCGGCGCTCTCGCCTTTCGGCGCGATGCA                                                                                                                      <  6:158243/120‑1                                                                                       CTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCATCGGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATG                                                                                                               <  3:232398/98‑1                                                                                       CTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCATCGGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATG                                                                                                               >  4:232398/1‑98                                                                                                                                                     TCGGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATGGTGATTTCGGTCTTGTTCACCGAAATCGGCCGGGCAACGCGGATCTGCGAGCTGAACTGGTCCATCAGGTACACGTTGGGGTACA                          <  5:156395/121‑1                                                                                                                                                     TCGGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATGGTGATTTCGGTCTTGTTCACCGAAATCGGCCGGGCAACGCGGATCTGCGAGCTGAACTGGTCCATCAGGTACACGTTGGGGTACA                          >  6:156395/1‑121                                                                                                                                                      CGGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATGGTGATTTCGGTCTTGTTCACCGAAATCGGCCGGGCAACGCGGATCTGCGAGCTGAACTGGTCCATCAGGTACACGTTGGGGTACAGGCACAGGTTGCGCGAGTTTTCGA  <  5:43767/144‑1                                                                                                                                                      CGGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATGGTGATTTCGGTCTTGTTCACCGAAATCGGCCGGGCAACGCGGATCTGCGAGCTGAACTGGTCCATCAGGTACACGTTGGGGTACAGGCACAGGTTGCGCGAGTTTTCGA  >  6:43767/1‑144                                                                                                                                                                                                                                                                                                       TTCGGCCAGGCACGTGCCGACTGGATGATCGAAAACTCGCGCAACCTGTGCCTGTACCCCAACGTGTACCTGATGGACCAGTTCAGCTCGCAGATCCGCGTTGCCCGGCCGATTTCGGTGAACAAGACCGAAATCACCATCTACTGCATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/3582794‑3582944 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcatcgGCGCTCTCGCCTTTCGGCGCGATGCAGTAGATGGTGATTTCGGTCTTGTTCACCGAAATCGGCCGGGCAACGCGGATCTGCGAGCTGAACTGGTCCATCAGGTACACGTTGGGGTACAGGCACAGGTTGCGCGAGTTTTCGA  <  NC_002947/3582963‑3582822

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