New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	3616215 =	43 (0.980)	5 (0.120)	5/328	NT	11.9%	coding (550/1266 nt)	codB	Cytosine transporter
	?	AE015451	3616308 =	33 (0.790)					coding (643/1266 nt)	codB	Cytosine transporter

GTGCCGAAGACCATGGTCACTGCCGCCGAAAAGGTCATGGTTTCATGCGGGACGACTGCAGCCAGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/3616282‑3616215 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaggccACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTGGCAACGGCCTGATGGTGGTGGCCGGGGCGTGGTGCGCGATGGTTTACCAGCAGGCCGATATCGTCGACGTGATGATGCTGCAGGGCCTGTCGTTCGCGGCCGTGGTTATGCTGTGCCTAAACCTGTGGACCATTCAGG  >  AE015451/3616308‑3616515                                                                                                                                                                                                                                                                                      GTGCCGAAGACCATGGTCACTGCCGCCGAAAAGGTCATGGTTTCATGCGGGACGACTGCCGCCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGT                                                                                                                                                                               >  2:38242/1‑103                             AAAAGGTCATGGTTTCGTGCGGGACGACTGCAGCCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGTCAGCG                                                                                                                                                              <  2:747327/92‑1                                  GTCATGGTTTCATGCGGGACGACTGCAGCCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTGGCAACGGCCTGATGGTGGTGGCCGGGGCGTGGTGCGCGATGGTTTACCAGCAGGCC                                                                                    >  2:23068/1‑161                                                      GACTGCAGCCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTGG                                                                                                                                           <  1:593076/86‑1                                                      GACTGCAGCCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTGG                                                                                                                                           >  2:593074/1‑86                                                              CCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTG                                                                                                                                            <  1:32681/77‑1                                                              CCAGGCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTG                                                                                                                                            >  2:32681/1‑77                                                                  GCCACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTGGCAACGGCCTGATGGTGGTGGCCGGGGCGTGGTGCGCGATGGTTTACCAGCAGGCCGATATCGTCGACGTGATGATGCTGCAGGGCCTGTCGTTCGCGGCCGTGGTTATGCTGTGCCTAAACCTGTGGACCATTCAGG  <  2:786632/211‑1                                                                                                                                                                                                                                                                                      GTGCCGAAGACCATGGTCACTGCCGCCGAAAAGGTCATGGTTTCATGCGGGACGACTGCAGCCAGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AE015451/3616282‑3616215 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaggccACCAACTGGACACGGCTGGCACGCAGTGGCCGGGTAGCGGTGATCGCCAGCGGGGTGGCATTCCTCCTTGGCAACGGCCTGATGGTGGTGGCCGGGGCGTGGTGCGCGATGGTTTACCAGCAGGCCGATATCGTCGACGTGATGATGCTGCAGGGCCTGTCGTTCGCGGCCGTGGTTATGCTGTGCCTAAACCTGTGGACCATTCAGG  >  AE015451/3616308‑3616515

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