New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2670137	103 (1.190)	3 (0.040)	3/416	NT	3.9%	coding (1106/1362 nt)	hcaE	3‑phenylpropionate dioxygenase, large (alpha) subunit
	?	NC_000913	= 2670349	50 (0.600)					coding (1318/1362 nt)	hcaE	3‑phenylpropionate dioxygenase, large (alpha) subunit
Rejected: Frequency below/above cutoff threshold.

CCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCGCCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTGAAGTGTGGGCGTTCTGTATTACTGACAAAGCCGCCTCCGATGAAGTTAAAGCCGCTTTTGAAAACAGCGCCACTCGTGCTTTTGGTCCTGCTGGTTTTCTCGAGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2669897‑2670137 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgagcaCTTCCTGCCAGCTTTCGCTACTCAGAAGATCGGCCCAGCGTTGGTACATTCCACGAGCGGCAGTTTCTGAAAAGATATAGTTAGTAATGCCAGGAATGCCGTCGTCGCGGCGCTTTTCCTGACCAAGCCCCATTTCCAGACACAGTTTGCTGTTGCGGGCGCGGTGGCCTTTAAGCAATTTCTGGATTTCACACCAGTTCTCCGAGTCATCCTGCTCGAGAAAACCAGCAGGACCAA  <  NC_000913/2670349‑2670113                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                CCGAAGCAGAACAAAGTTTAGGTGGATGTCCCGGCCTGCGCCGTGCGGGTCATAACATGATTTCCCCCACGCTTTCAGTGCAAAAAGTCCCTGCCACACGCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTGAAGTGTGGGCGTTCTGTATTACTGACAAAGCCGCCTCCGATGAAGTTAAAGCCGCTTTTGAAAACAGCGCCACTCGTGCTTTTGGTCCTGCTGGTTTTCTCGAG                                                                                                                                                                                                                                                 <  2:657059/239‑1  CGAAGCAGACCAACGTTTAGGTGAAGTTCGCGCCCTGCGCCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTGAAGTGTGGGCGTTCTGTATTACTGACAAAGCCGCCTCCGATGAAGTTAAAGCCGCTTTTGAAAACAGCGCCACTCGTGCTTTTGGTCCTGCTGGTTTTCTCGAG                                                                                                                                                                                                                                                 <  1:334791/238‑1                                                                                                                                    GTTGAAGTGTGGGCGTTCTGTATTACTGACAAAGCCGCCTCCGATGAAGTTAAAGCCGCTTTTGAAAACAGCGCCACTCGTGCTTTTGGTCCTGCTGGTTTTCTCGAGCACTTCCTGCCAGCTTTCGCTACTCAGAAGATCGGCCCAGCGTTGGTACATTCCACGAGCGGCAGTTTCTGAAAAGATATAGTTAGTAATGCCAGGAATGCCGCCTTCGCGGCGCTTTTCCTGACCAAcc                                                                                                               >  2:480695/1‑236                                                                                                                                                                              GATGAAGTTAAAGCCGGTTTTGAAAACAGCGCCACTCGTGCTTTTGGTCCTGCTGGTTTTCTCGAGCATTTCCTGCCAGTTTTCGCTGCTCAGAAGATCGGCCCAGCGTTGGTACATTCCACGAGCGGCAGTTTCTGAAAAGATATAGTTAGTAATGCCAGGAATGCCGTCGTCGCGGCGCTTTTCCTGACCAAGCCCCATTTCCAGACACAGTTTGCTGTTGCGGGCGCGGTGGCCTT                                                                    <  1:480695/239‑1                                                                                                                                                                                                                              agggggGGTTTTTTCGGGTACTTCCGGCCAGTGTTGGGTACTCAGAAGATGGGACCAGCGTGGGTACATTCCACGAGCGGCAGTGTCTGAAAAGATATAGTTAGTAATGCCAGGAATGCCGTCGTCGCGGCGCTTTTCCTGACCAAGCCCCATTTCCAGACACAGTTTGCTGTTGCGGGCGCGGGGGCCTTTAAGCAATTTCTGGATTTCACACCAGTTCTCCGAGTCATCCTGCTCG                     <  2:994613/232‑1                                                                                                                                                                                                                                                CACTTCCTGCCAGCTTTCGCTACTCAGAAGATCGGCCCAGCGTTGGTACATTCCACGAGCGGCAGTTTCTGAAAAGATATAGTTAGTAATGCCAGGAATGCCGTCGTCGCGGCGCTTTTCCTGACCAAGCCCCATTTCCAGACACAGTTTGCTGTTGCGGGCGCGGTGGCCTTTAAGCAATTTCTGGATTTCAAACCAGATCTCCGAGTGATCCTTCTCGGGAAAACCTGCAGGACCAA  >  1:874792/1‑239                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                CCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCGCCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTGAAGTGTGGGCGTTCTGTATTACTGACAAAGCCGCCTCCGATGAAGTTAAAGCCGCTTTTGAAAACAGCGCCACTCGTGCTTTTGGTCCTGCTGGTTTTCTCGAGCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2669897‑2670137 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgagcaCTTCCTGCCAGCTTTCGCTACTCAGAAGATCGGCCCAGCGTTGGTACATTCCACGAGCGGCAGTTTCTGAAAAGATATAGTTAGTAATGCCAGGAATGCCGTCGTCGCGGCGCTTTTCCTGACCAAGCCCCATTTCCAGACACAGTTTGCTGTTGCGGGCGCGGTGGCCTTTAAGCAATTTCTGGATTTCACACCAGTTCTCCGAGTCATCCTGCTCGAGAAAACCAGCAGGACCAA  <  NC_000913/2670349‑2670113

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