New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 5922049	40 (0.980)	4 (0.100)	4/316	NT	8.2%	coding (2382/2874 nt)	gcvP‑II	glycine dehydrogenase
	?	AE015451	= 5922165	52 (1.350)					coding (2266/2874 nt)	gcvP‑II	glycine dehydrogenase

TCGACAAACCGATCCAGCTCAGCCTTCGATTCACTCTCGGTCGGCTCCACCATCAGCGTGCCCGGCACTGGGAAGGACATGGTCGGTGCGTGGAAGCCATAGTCCATCAGGCGCTTGGCCACGTCCTCTTCGCTGATGCCTGTCAGCGCTTTCAACGGCCGCAGGTCGAGGATGCATTCATGCGCCACGCGCTGGTTGCGCCCGCGGTAGAGCACCGGGAACGCCGCGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5921819‑5922049 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgcgcccTGGGGCAGCGCAAGCATCTTGCCGATCAGCTGGATGTACATCGCCATGATGGGCCCGCAGCTGGCCGACGCCAGCGAGGTGGCTATCCTTTCGGCCAACTACCTGGCCAGCCAGTTGGGCGCGGCGTTCCCGGTGCTCTACCGCGGGCGCAACC  <  AE015451/5922165‑5922012                                                                                                                                                                                                                                                                                                                                                                                                   gCGACAAACCGATCCAGCTCAGCCTTCGATTCACTCTCGGTCGGCTCCACCATCAGCGTGCCCGGCACTGGGAAGGACATGGTCGGTGCGTGGAAGCCATAGTCCATCAGGCGCTTGGCCACGTCCTCTTCGCTGATGCCTGTCAGCGCTTTCAACGGCCGCAGGTCGAGGATGCATTCATGCGCCACGCGCTGGTTGCGCCCGCGGTAGAGCACCGGGAACGCC                                                                                                                                                                  <  1:415100/224‑1      AAACCGATCCAGCTCAGCCTTCGATTCACTCTCGGTCGGCTCCACCATCAGCGTGCCCGGCACTGGGAAGGACATGGTCGGTGCGTGGAAGCCATAGTCCATCAGGCGCTTGGCCACGTCCTCTTCGCTGATGCCTGTCAGCGCTTTCAACGGCCGCAGGTCGAGGATGCATTCATGCGCCACGCGCTGGTTGCGCCCGCGGTAGAGCACCGGGAACGCCGCGCCCTGGGGCAGC                                                                                                                                                   <  1:426452/235‑1                                                                                                              GGCGCTTGGCCACGTCCTCTTCGCTGATGCCTGTCAGCGCTTTCAACGGCCGCAGGTCGAGGATGCATTCATGCGCCACGCGCTGGTTGCGCCCGCGGTAGAGCACCGGGAACGCCGCGCCCTGGGGCAGCGCAAGCATCTTGCCGATCAGCTGGATGTACATCGCCATG                                                                                                            >  1:518112/1‑170                                                                                                              GGCGCTTGGCCACGTCCTCTTCGCTGATGCCTGTCAGCGCTTTCAACGGCCGCAGGTCGAGGATGCATTCATGCGCCACGCGCTGGTTGCGCCCGCGGTAGAGCACCGGGAACGCCGCGCCCTGGGGCAGCGCAAGCATCTTGCCGATCAGCTGGATGTACATCGCCATG                                                                                                            <  2:518111/170‑1                                                                                                                                                                                                            GCGGTAGAGCACCGGGAACGCCGCGCCCTGGGGCAGCGCAAGCATCTTGCCGATCAGCTGGGTGTACATCGCCATGATGGGCCCGCAGCTGGCCGACGCCAGCGAGGTGGCTATCCTTTCGGCCAACTACCTGGCCAGCCAGTTGGGCGCGGCGTTCCCGGTGCTCTACCGCGGGCGCAACC  >  2:518116/1‑182                                                                                                                                                                                                                                                                                                                                                                                                   TCGACAAACCGATCCAGCTCAGCCTTCGATTCACTCTCGGTCGGCTCCACCATCAGCGTGCCCGGCACTGGGAAGGACATGGTCGGTGCGTGGAAGCCATAGTCCATCAGGCGCTTGGCCACGTCCTCTTCGCTGATGCCTGTCAGCGCTTTCAACGGCCGCAGGTCGAGGATGCATTCATGCGCCACGCGCTGGTTGCGCCCGCGGTAGAGCACCGGGAACGCCGCGCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/5921819‑5922049 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgcgcccTGGGGCAGCGCAAGCATCTTGCCGATCAGCTGGATGTACATCGCCATGATGGGCCCGCAGCTGGCCGACGCCAGCGAGGTGGCTATCCTTTCGGCCAACTACCTGGCCAGCCAGTTGGGCGCGGCGTTCCCGGTGCTCTACCGCGGGCGCAACC  <  AE015451/5922165‑5922012

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