New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3683468	111 (1.130)	5 (0.050) +ACAG	5/240	NT	4.4%	intergenic (‑21/+162)	dctA/yhjK	C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
	?	NC_000913	3683469 =	112 (1.110)					intergenic (‑22/+161)	dctA/yhjK	C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Rejected: Frequency below/above cutoff threshold.

ACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3683323‑3683468 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAA  >  NC_000913/3683469‑3683606                                                                                                                                                   ||||                                                                                                                                           ACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCAC                                                                                                                                              <  6:17308/148‑1       CCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGACAG                                                                                                                                        >  8:483062/1‑148                                                           GGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCAC                                                                                                                                              >  5:298823/1‑90                                                           GGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCAC                                                                                                                                              <  6:298823/90‑1                                                                  ACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTACACAGAGAGGTGTTCATAGGGTGTCCTTTAGTAAATCCACAGACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATACCACCCAGACAACATCACAG                                                                              >  5:203925/1‑147                                                                     AATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCA                                                                                           >  5:316404/1‑131                                                                     AATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCA                                                                                           <  6:316404/131‑1                                                                         GCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTA                                                                       <  7:483062/147‑1                                                                                                                                             AAAACCACAGACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAA  >  2:79089/1‑148                                                                                                                                                   ||||                                                                                                                                           ACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3683323‑3683468 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAA  >  NC_000913/3683469‑3683606

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