New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1541977 =	36 (0.640)	3 (0.050)	3/272	5.6	6.2%	coding (614/3741 nt)	narZ	nitrate reductase Z subunit alpha
	?	NC_000913	= 4406770	55 (0.990)					coding (117/2442 nt)	rnr	RNase R
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCCGCTCCAACTGGCAGGAACTAAACCAGCTGATTGCCGCCGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1542115‑1541977 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gCTCATCACGGCTGGCCGGTTTTTCACGTTTGGTTAAATGTTCGAGGATAAATTCCCGACTAGGGATGGGATTCGCGTATTTTTCAGCTTCGCGTTCCTGGAAAG  <  NC_000913/4406770‑4406667                                                                                                                                                                                                                                                     TCCGCTCCAACTGGCAGGAACTAAACCAGCTGATTGCCGCCGCCAACGTCTGGACCATCAAAACCCACGGCCCCGGTCGCGTTGCCGGGTTCTCGCCGATCCCCGCGATGTCGGTGGTTTCTTACGCCCCCGGAACGCGC                                                                                                         >  2:368357/1‑140                                     CCGCCGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGCTCATCACGGCTGGCCGGTTTTTCACGTTTGGTTAAA                                                                     >  1:329361/1‑140                                                                                                        GGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGCTCATCACGGCTGGCCGGTTTTTCACGTTTGGTTAAATGTTCGAGGATAAATTCCCGACTAGGGATGGGATTCGCGTATTTTTCAGCTTCGCGTTCCTGGAAAG  <  2:329361/140‑1                                                                                                                                                                                                                                                     TCCGCTCCAACTGGCAGGAACTAAACCAGCTGATTGCCGCCGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1542115‑1541977 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gCTCATCACGGCTGGCCGGTTTTTCACGTTTGGTTAAATGTTCGAGGATAAATTCCCGACTAGGGATGGGATTCGCGTATTTTTCAGCTTCGCGTTCCTGGAAAG  <  NC_000913/4406770‑4406667

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