New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4399677 =	86 (1.720)	4 (0.080)	4/320	NT	4.6%	coding (213/972 nt)	PP_3878	putative Phage minor capsid protein C
	?	AE015451	4399760 =	86 (1.810)					coding (130/972 nt)	PP_3878	putative Phage minor capsid protein C
Rejected: Frequency below/above cutoff threshold.

GCTTGCCGCATGAAACGACATTTGCGCGCTGCCAGCTTGCTGTTCAACCAGCCGCTACTGACTACACCCGACATGCTGGACCTGGCAGTGCGCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4399898‑4399677 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgccgcCCATGTTCAGATTGACGATGTTCAGGCTCATGGTCTGGTTCGCCCAGCGCACTGCCAGGTCCAGCATGTCGGGTGTAGTCAGTAGCGGCTGGTTGAACAGCAAGCTGGCAGCGCGCAAATGTCGTTTCATGCGGCAAGCATCCTCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGT  >  AE015451/4399760‑4399989                                                                                                                                                                                                                                                                                                                                                                                                                                                                      GCTTGCCGCATGAAACGACATTTGCGCGCTGCCAGCTTGCTGTTCAACCAGCCGCTACTGACTACACCCGACATGCTGGACCTGGCAGTGCGCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACA                                                                                                                                                                                                                                               >  2:670250/1‑215                             CTGCCAGCTTGCTGTTCAACCAGCCGCTACTGACTACACCCGACATGCTGGACCTGGCAGTGCGCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCG                                                                                                                                                                                                                                         <  2:846015/193‑1                                           TTCAACCAGCCGCTACTGACTACACCCGACATGCTGGACCTGGCAGTGCGCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAG                                                                                                                                                                                                                                              >  2:807706/1‑174                                                                               gaaCTGGCAGTGCGCGGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCG                                                                                                                                                                                                                                         <  2:51960/140‑1                                                                                            GCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAG                                                                                                                                                                                                                                              <  1:221625/125‑1                                                                                            GCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAG                                                                                                                                                                                                                                              >  2:221623/1‑125                                                                                                              TGAGACTTACCATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCGCCCATGTTCAGATTGACGATGTTCAGGCTCATGGTCT                                                                                                                                                                                                    <  1:47224/149‑1                                                                                                                                                                         TATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCGCCCATGTTCAGATTGACGATGTTCAGGCTCATGGTCTGGTTCGCCCAGCGCACTGCCAGGTCCAGCATGTCG                                                                                                                                                                 <  1:667477/125‑1                                                                                                                                                                         TATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCGCCCATGTTCAGATTGACGATGTTCAGGCTCATGGTCTGGTTCGCCCAGCGCACTGCCAGGTCCAGCATGTCG                                                                                                                                                                 >  2:667474/1‑125                                                                                                                                                                                          CCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCGCCCATGTTCAGATTGACGAGGTTCAGGCTCATGGTCTGGTTCGCCCAGCGCACTGCCAGGTCCAGCATGTCGGGTGTAGTCAGTAGCGGCTGGTTGAACAGCAAGCTGGCAGCGCGCAAATGTCGTT                                                                                                          <  2:38394/163‑1                                                                                                                                                                                                                         CGCCGCCCATGTTCAGATTGACGATGTTCAGGCTCATGGTCTGGTTCGCCCAGCGCACTGCCAGGTCCAGCATGTCGGGTGTAGTCAGTAGCGGCTGGTTGAACAGCAAGCTGGCAGCGCGCAAATGTCGTTTCATGCGGCAAGCATCCTCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGT  >  1:824482/1‑236                                                                                                                                                                                                                                                                                                                                                                                                                                                                      GCTTGCCGCATGAAACGACATTTGCGCGCTGCCAGCTTGCTGTTCAACCAGCCGCTACTGACTACACCCGACATGCTGGACCTGGCAGTGCGCTGGGCGAACCAGACCATGAGCCTGAACATCGTCAATCTGAACATGGGCGGCGCTGCAGCCAATCCATCGATGTTCTATGACGATGAAGACTACCAGGCTGAGCAGGATCGCCGCGAAGAACAGCGCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4399898‑4399677 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cagcgccgcCCATGTTCAGATTGACGATGTTCAGGCTCATGGTCTGGTTCGCCCAGCGCACTGCCAGGTCCAGCATGTCGGGTGTAGTCAGTAGCGGCTGGTTGAACAGCAAGCTGGCAGCGCGCAAATGTCGTTTCATGCGGCAAGCATCCTCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGT  >  AE015451/4399760‑4399989

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