New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2671662	79 (0.990)	3 (0.040)	3/488	NT	3.7%	coding (438/813 nt)	hcaB	2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
	?	NC_000913	= 2671671	81 (1.040)					coding (447/813 nt)	hcaB	2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
Rejected: Frequency below/above cutoff threshold.

TTGGCGAACATATTCTGGCGGTGGAAGGTAACGTGACCTGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2671373‑2671662 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgccAGGATACCAGGCGGCATTTGACAGTGTGAAAATCATGCTGCCTTCACTGGCGATTAA  <  NC_000913/2671671‑2671615                                                                                                                                                                                                                                                                                                                                                             TTGGCGAACATATTCTGGCGGTGGAAGGTAACGTGACCTGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTACCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCG                                                             >  1:95821/1‑288  TGGCGAACATATTCTGGCGGTGGAAGGTAACGTGACCTGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCC                                                           <  1:170322/289‑1  TGGCGAACATATTCTGGCGGTGGAAGGTAACGTGACCTGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGC                                                            >  2:53236/1‑288                                       TGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCC                                                              >  1:62146/1‑249                                       TGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCC                                                              <  2:62146/249‑1                                                           CGCGTTCGATCATGTCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCTAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCATTTGACAGTGTGAAAATCATGCTGCCTTCACTGGCGATTAA  <  1:29743/289‑1                                                                                                                                                                                                               GGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCATTTGACAGTGTGAAAATCATGCTGCCTTCACTG          <  1:577275/133‑1                                                                                                                                                                                                               GGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCATTTGACAGTGTGAAAATCATGCTGCCTTCACTG          >  2:577275/1‑133                                                                                                                                                                                                                                       TGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCC                                                           >  1:40761/1‑60                                                                                                                                                                                                                                       TGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCC                                                           <  2:40761/60‑1                                                                                                                                                                                                                                                                                                                                                             TTGGCGAACATATTCTGGCGGTGGAAGGTAACGTGACCTGTTATGCCGATTATCAACGCGCGGTCGATCAGATCCTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2671373‑2671662 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgccAGGATACCAGGCGGCATTTGACAGTGTGAAAATCATGCTGCCTTCACTGGCGATTAA  <  NC_000913/2671671‑2671615

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