New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	4562960 =	46 (1.150)	3 (0.080)	3/322	NT	6.4%	coding (901/2316 nt)	vgrG‑III	type 6 secretion system protein
	?	AE015451	4563127 =	44 (1.150)					coding (734/2316 nt)	vgrG‑III	type 6 secretion system protein

GCGCTATTGCGAGCAGGTACGTACTGCGCGGCAAGTGCAGCGCGACTACACCTTCACCAACCCGGCCTACCGGCAGGAGCACAAGGCCCAGACCAGTGGCGTCACGCACCAGTCACTGGACTACGAGCGCTTTGATTACCCGGGGCGCTACAAGCGCGATGCGGTGGGCAAGCCCTTCACTGAAACACGGCTGACGGCGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4563162‑4562960 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgctgcACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGC  >  AE015451/4563127‑4563193                                                                                                                                                                                                                                                                                GCGCTATTGCGAGCAGGTACGTACTGCGCGGCAAGTGCAGCGCGACTACACCTTCACCAACCCGGCCTACCGGCAGGAGCACAAGGCCCAGACCAGTGGCGTCACGCACCAGTCACTGGACTACGAGCGCTTTGATTACCCGGGGCGCTACAAGCGCGATGCGGTGGGCAAGCCCTTCACTGAAACACGGCTGACGGC                                                                          <  2:364750/198‑1                                                                                     gcCCCAGACCAGTGGCGTCACGCACCAGTCACTGGACTACGAGCGCTTTGATTACCCGGGGCGCTACAAGCGCGATGCGGTGGGCAAGCCCTTCCCTGAAACACGGCGGACGGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCAT                  <  1:478707/168‑1                                                                                                                               GCGCTTTGATTACCCGGGGCGCTACAAGCGCGATGCGGTGGGCAAGCCCTTCACTGAAACACGGCTGACGGCG                                                                         >  2:686504/1‑73                                                                                                                                                                               CTTCACTGAAACACGGCTGACGGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGC  <  1:601690/96‑1                                                                                                                                                                               CTTCACTGAAACACGGCTGACGGCGCTGCACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGC  >  2:601689/1‑96                                                                                                                                                                                                                                                                                GCGCTATTGCGAGCAGGTACGTACTGCGCGGCAAGTGCAGCGCGACTACACCTTCACCAACCCGGCCTACCGGCAGGAGCACAAGGCCCAGACCAGTGGCGTCACGCACCAGTCACTGGACTACGAGCGCTTTGATTACCCGGGGCGCTACAAGCGCGATGCGGTGGGCAAGCCCTTCACTGAAACACGGCTGACGGCGCTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/4563162‑4562960 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgctgcACTTGCCGCGCAGTACGTACCTGCTCGCAATAGCGCAAGCGCCGCAGGCATGCGCGGGCCTGGTCGC  >  AE015451/4563127‑4563193

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