New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 2387701	62 (0.850)	4 (0.060)	4/238	NT	6.1%	coding (418/1023 nt)	pyrD	dihydroorotate dehydrogenase
	?	NC_002947	= 2387735	63 (0.900)					coding (452/1023 nt)	pyrD	dihydroorotate dehydrogenase

CGCAGCCTGGCAACCCCAAGCCACGGTTGTTCCGCCTGCCGGAAGCCACGGCCATCATCAACCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/2387552‑2387701 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑catcgACAGCACGTTCGACCGGGGTATCGAAGTTCTTGCCGATGTTGATGCCCAGCACCCCC  <  NC_002947/2387735‑2387679                                                                                                                                                                                                                 CGCAGCCTGGCAACCCCAAGCCACGGTTGTTCCGCCTGCCGCAAGCCACTCCCATCATCAACCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCG                                                           <  2:743372/150‑1                                                          TCAACCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCGACAGCACGTTCGACCGGGGTATCGAAGTTCTTGCCGATGTTGATGCCCAGCACCCCC  >  1:484185/1‑150                                                          TCAACCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCGACAGCACGTTCGACCGGGGTATCGAAGTTCTTGCCGATGTTGATGCCCAGCACCCCC  <  2:484185/150‑1                                                              CCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCGACAGCACGTTCGACCGGGGTATCGAAGTTCTTGCCGATGTTGA                <  1:1924608/132‑1                                                              CCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCGACAGCACGTTCGACCGGGGTATCGAAGTTCTTGCCGATGTTGA                >  2:1924608/1‑132                                                                                                                                                                                                                 CGCAGCCTGGCAACCCCAAGCCACGGTTGTTCCGCCTGCCGGAAGCCACGGCCATCATCAACCGCATGGGCTTCAACAACCTGGGCGTTGATCACCTGCTGGGGCGCGTACGCGCCTCGCGCTACAAGGGGGTGCTGGGCATCAACATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/2387552‑2387701 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑catcgACAGCACGTTCGACCGGGGTATCGAAGTTCTTGCCGATGTTGATGCCCAGCACCCCC  <  NC_002947/2387735‑2387679

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