New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5739225 =	53 (0.890)	4 (0.070)	4/234	NT	7.4%	coding (99/1365 nt)	hutF	formiminoglutamate deiminase
	?	NC_002947	5739261 =	49 (0.840)					coding (135/1365 nt)	hutF	formiminoglutamate deiminase

GGTCTCGCCAGGTCCAGAAACTGTCGTTGGGGTTGCCGGCGACTTCCGCCAGGCCTGCCATGGCGCGCTGGAAGGCATGCGAGTGCAGGTTGGGCATGCCGGGCAGCAGCGGGCCGGCCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5739376‑5739225 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACA  >  NC_002947/5739261‑5739355                                                                                                                                                                                                                                                         GGTCTCGCCAGGTCCAGAAACTGTCGTTGGGGTTGCCGGCGACTTCCGCCAGGCCTGCCATGGCGCGCTGGAAGGCATGCGAGTGCAGGTTGGGCATGCCGGGCAGCAGCGGGCCGGCCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAG                                                                                                   <  2:1117447/150‑1                                                                                                                    GGCCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACA  <  1:143996/132‑1                                                                                                                    GGCCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAGGCCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACA  >  2:143996/1‑132                                                                                                                      CCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAGGCCGGCCCGCTGCTGCCCGGCATGC                                                                          <  1:1225990/58‑1                                                                                                                      CCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAGGCCGGCCCGCTGCTGCCCGGCATGC                                                                          >  2:1225990/1‑58                                                                                                                                                                                                                                                         GGTCTCGCCAGGTCCAGAAACTGTCGTTGGGGTTGCCGGCGACTTCCGCCAGGCCTGCCATGGCGCGCTGGAAGGCATGCGAGTGCAGGTTGGGCATGCCGGGCAGCAGCGGGCCGGCCAGGCGCTCGGCGTCTTCGGCCGAAGCGCCAGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5739376‑5739225 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcCGGCCCGCTGCTGCCCGGCATGCCCAACCTGCACTCGCATGCCTTCCAGCGCGCCATGGCAGGCCTGGCGGAAGTCGCCGGCAACCCCAACGACA  >  NC_002947/5739261‑5739355

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