New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 51828	56 (0.970)	5 (0.100)	5/310	NT	8.9%	coding (162/1272 nt)	czcC‑I	cobalt‑zinc‑cadmium resistance protein
	?	AE015451	= 51812	51 (0.980)					coding (178/1272 nt)	czcC‑I	cobalt‑zinc‑cadmium resistance protein

CGACATCTGCACGCAGGCCGTTAACGCGACGGTCTAACTCAAGCTGTGCGGCGGTCTGCCCGTATGTGGCGACGTCAACACGAGCCCCCCGCTTACCGCCCAGCTCAAGGGTCTGGCTGAGAGAGACAGTCGTCTGGCTGGTGTTACGACGGGTGTCTTCCACGTCATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGCGGCCAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/51578‑51828 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctggccgcggccaggCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGACGTGGAAGACACCCGTCGTAACACCAGCCAGACGACTGTCTCTCTCAGCCAGACCCTTGAGCTGGGCGGTAAGCGGGGGGCTCGTGTTGACGTCGCCACATACGGGCAGACCGCCGCACAGCTTGAGTTAGACCGTCGCGTTAACGGCCTGC  <  AE015451/51812‑51590                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            CGACATCTGCACGCAGGCCGTTAACGCGACGGTCTAACTCAAGCTGTGCGGCGGTCTGCCCGTATGTGGCGACGTCAACACGAGCCCCCCGCTTACCGCCCAGCTCAAGGGTCTGGCTGAGAGAGACAGTCGTCTGGCTGGTGTTACGACGGGTGTCTTCCACGTCATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGC                                                                                                                                                                                                                                                >  1:849128/1‑236                                GTCTAACTCAAGCTGTGCGGCGGTCTGCCCGTATGTGGCGACGTCAACACGAGCCCCCCGCTTACAGCCCAGCTCAAGGGTCTGGCTGAGAGAGACAGTCGTCTGGCTGGTGTTACGACGGGTGTCTTCCACGTCATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATGCCAATTTCTTGGCTGGCC                                                                                                                                                                                                                                          >  2:556796/1‑211                                                                         CGTCAACACGAGCCCCCCGCTTACCGCCCAGCTCAAGGGTCTGGCTGAGAGAGACAGTCGTCTGGCTGGTGTTACGACGGGTGTCTTCCACGTCATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGCGGCCAGGCAAGAAATTGGGATCGC                                                                                                                                                                                                                <  2:472276/196‑1                                                                                                                                                     cCGGGTGGCTTCCCCGCCATACGAGATTGTGGGGTTGGGGAGAGGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCC                                                                                                                                                                        <  2:287503/159‑1                                                                                                                                                                 CACGTCATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCCCAATCTCGTATGACGTGGAAGACACCCGTCGTAACACCAGC                                                                                                                              >  1:354235/1‑190                                                                                                                                                                      CATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCCGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGACGTGGAA                                                                                                                                                   >  2:909284/1‑164                                                                                                                                                                                            GATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGC                                                                                                                                                                                                                                        <  1:530148/57‑1                                                                                                                                                                                                                                   AATTTCTTGCCTGGCCGCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGACGTGGAAGACACCCGTCGTAACACCAGCCAGACGACTGTCTCTCTCAGCCAGACCCT                                                                                                 >  1:1038168/1‑153                                                                                                                                                                                                                                              TGGCCGCGGCCAGGCAAGAAATTGGGCTCGCTGGTGGTGCTCGCAAGCAGGCCGGTCTTATCCCCAACCCCACAATCTCGTATGACGTGGAAGACACCCGTCGTAACACCAGCCAGACGACTGTCTCTCTCAGCCAGACCCTTGAGCTGGGCGGTAAGCGGGGGGCTCGTGTTGA                                                                <  2:35741/175‑1                                                                                                                                                                                                                                              TGGCCGCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGCCGTCTTCTGCACCAGTAGAAACACCAGCCAGACGACTGACTCTCTAAGCCAGACCCTTGAGCTGGGCGG                                                                                     >  1:808268/1‑154                                                                                                                                                                                                                                                   GCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGACGTGGAAGACACCCGTCGTAACACCAGCCAGACGA                                                                                                                       >  2:388341/1‑115                                                                                                                                                                                                                                                   GCGGCCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATC                                                                                                                                                                                 <  2:530145/57‑1                                                                                                                                                                                                                                                       CCAGGCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGCCGTCTTCTGCTTGAAAAAAAAAACCAGCCAGACGACTGTCTCTCTCAGCCAGACCCTTGAGCTGGGCGGTAAGCGGGGGGCTCGTGTTGACGTCGCCACATACGGGCAGACCGCCGCACAGCTTGAGTTAGACCGTCGCGTTCACGGCCTGC  >  1:782539/1‑228                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            CGACATCTGCACGCAGGCCGTTAACGCGACGGTCTAACTCAAGCTGTGCGGCGGTCTGCCCGTATGTGGCGACGTCAACACGAGCCCCCCGCTTACCGCCCAGCTCAAGGGTCTGGCTGAGAGAGACAGTCGTCTGGCTGGTGTTACGACGGGTGTCTTCCACGTCATACGAGATTGTGGGGTTGGGGATAAGCCCGGCCTGCTTGCGAGCACCATCAGCGATCCCAATTTCTTGCCTGGCCGCGGCCAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/51578‑51828 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cctggccgcggccaggCAAGAAATTGGGATCGCTGATGGTGCTCGCAAGCAGGCCGGGCTTATCCCCAACCCCACAATCTCGTATGACGTGGAAGACACCCGTCGTAACACCAGCCAGACGACTGTCTCTCTCAGCCAGACCCTTGAGCTGGGCGGTAAGCGGGGGGCTCGTGTTGACGTCGCCACATACGGGCAGACCGCCGCACAGCTTGAGTTAGACCGTCGCGTTAACGGCCTGC  <  AE015451/51812‑51590

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