New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5595982	58 (0.970)	3 (0.050)	3/230	NT	5.2%	coding (148/1881 nt)	thiC	phosphomethylpyrimidine synthase
	?	NC_002947	= 5596006	54 (0.940)					coding (124/1881 nt)	thiC	phosphomethylpyrimidine synthase

TCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCGTTGGACTCGCCGCCGAAATCGGTCGGGGTGTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5595831‑5595982 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgtcCCCATGCGTGAAATCAGCCTGGATGACACCCCGACCGATTTCGGCGGCGAGTCCAA  <  NC_002947/5596006‑5595951                                                                                                                                                                                                                  TCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCGTTGGACTCGCCGCCGAAATCGGTCGGGGTG                                                            >  2:778346/1‑150                                   GCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCGTTGGACTCGCCGCCGAAATCGGTCGGGGTGTCCCCATGCGTGAAATCAGCCTGGATGACACCCC                          <  2:1022420/150‑1                                                                                          GCCAGACGTGTCATACACCAGCACCGGGGCGTTGGACTCGCCGCCGAAATCGGTCGGGGTGTCCCCATGCGTGAAATCAGCCTGGATGACAC                             <  1:454990/92‑1                                                                                          GCCAGACGTGTCATACACCAGCACCGGGGCGTTGGACTCACCGCCGAAATCGGTCGGGGTGTCCCCATGCGTGAAATCAGCCTGGATGACAC                             >  2:454990/1‑92                                                                                                                                                       TCCCCATGCGTGAAATCAGCCTGGATGACACCCCGACCGATTTCGGCGGCGAGTCCAA  <  1:319919/58‑1                                                                                                                                                       TCCCCATGCGTGAAATCAGCCTGGATGACACCCCGACCGATTTCGGCGGCGAGTCCAA  >  2:319919/1‑58                                                                                                                                                                                                                  TCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCGTTGGACTCGCCGCCGAAATCGGTCGGGGTGTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5595831‑5595982 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtgtcCCCATGCGTGAAATCAGCCTGGATGACACCCCGACCGATTTCGGCGGCGAGTCCAA  <  NC_002947/5596006‑5595951

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