New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4846196 =	61 (1.060)	3 (0.050)	3/328	NT	4.9%	coding (946/2475 nt)	PP_4261	cation‑transporting P‑type ATPase
	?	AE015451	4846222 =	57 (1.030)					coding (972/2475 nt)	PP_4261	cation‑transporting P‑type ATPase
Rejected: Frequency below/above cutoff threshold.

CGCGCCGGGCAGCACCTGCACGGTGTCGCCGCACTGTAACTCGCTGAGCAGGATGCGTTCGCTGCGCCCGATGGCGTCCAGGCGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4846282‑4846196 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcagCGAACGCATCCTGCTCAGCGAGTTACAGTGCGGCGACACCGTGCAGGTGCTGCCCGGCGCGGTCATCCCTGCCGACGGGCGCATCGTCGAGGGCCGCTCCAGCGTCGA  >  AE015451/4846222‑4846329                                                                                                                                                                                                     CGCGCCGGGCAGCACCTGCACGGTGTCGCCGCACTGTAACTCGCTGAGCAGGATGCGTTCGCTGCGCCCGATGGCGTCCAGGCGCAGCGAACGCATCCTGCTC                                                                                              <  2:141107/103‑1                                                GCAGGATGCGTTCGCTGCGCCCGATGGCGTCCAGGCGCAGCGAACGCATCCTGCTCAGCGAGTTACAGTGCGGCGACACCGTGCAGGTGCTGCCCGGCGCGGCCATCCCTGCCGACGGGCGCATCGTCGAGGGCCGCTCCAGCGTCGA  >  2:276192/1‑148                                                  AGGATGCGTTCGCTGCGCCCGATGGCGTCCAGGCGCAGCGAACGCATCCTGCTCAGCGAGTTACAGTGCGGCGACACCGTGCAGGTGCTGCCCGGCGCG                                                 >  2:653539/1‑99                                                                                     CAGCGAACGCATCCTGCTCAGCGAGTTACAGTGCGGCGACACCGTGCAGGTGCTGCCCGGCGCGGTCATCCCTGCC                                     <  1:110271/76‑1                                                                                     CAGCGAACGCATCCTGCTCAGCGAGTTACAGTGCGGCGACACCGTGCAGGTGCTGCCCGGCGCGGTCATCCCTGCC                                     >  2:110269/1‑76                                                                                                                                                                                                     CGCGCCGGGCAGCACCTGCACGGTGTCGCCGCACTGTAACTCGCTGAGCAGGATGCGTTCGCTGCGCCCGATGGCGTCCAGGCGCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4846282‑4846196 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgcagCGAACGCATCCTGCTCAGCGAGTTACAGTGCGGCGACACCGTGCAGGTGCTGCCCGGCGCGGTCATCCCTGCCGACGGGCGCATCGTCGAGGGCCGCTCCAGCGTCGA  >  AE015451/4846222‑4846329

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