New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 487811	36 (0.620)	3 (0.050)	3/330	NT	8.0%	coding (491/990 nt)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase
	?	AE015451	= 487820	34 (0.610)					coding (482/990 nt)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase

AGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/487583‑487811 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgtggCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGAACGTGTAACCCGTATCCTGCATGCCGACATGCGCGACAAGTTTGGCATTGCCAACCCGCGTATCCTCGTGTGCGGCCTCAACCCGCATGCGGGCGAGG  <  AE015451/487820‑487659                                                                                                                                                                                                                                                                                                                                                                                                         AGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGG                                                                                                                                                                    >  2:851470/1‑229                                                                                                   GCCGCACACGAGGATCCGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGGCCCTGGTGCCCACGCACCTGCCACTGCGTGATGTTGCCGA                                                                                                                            <  1:248162/171‑1                                                                                                                                                                                                            AACATCACGCAGTGGCAGGTGCGTGGCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGA                                                                                                    >  2:997768/1‑90                                                                                                                                                                                                                       GTGGCAGGTGCGTGGCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGAACGTGTAACCCGTATCCTGCATGCCGACATGCGCGACAAGTTTGGCATTGCCAACCCGCGTATCCTCGTGTGCGGCCTCAACCCGCATGCGGGCGAGG  <  1:997771/177‑1                                                                                                                                                                                                                                  GTGGCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGAACGTGTAACCCGTATCCTGCATGCCGACATGCGCGACCAGTTTCGCATTGCCAACCCGCGTATCCTCGTG                              >  2:72967/1‑138                                                                                                                                                                                                                                  GTGGCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGAACG                                                                                                 <  1:855501/71‑1                                                                                                                                                                                                                                  GTGGCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGAACG                                                                                                 >  2:855498/1‑71                                                                                                                                                                                                                                                                                                                                                                                                         AGGTCCATGCCTTCGGTGCGCAGGCGGGCCAGGGTCGGCTCGATGATGTCGATTTCTTCGCGGCCCAGGTGGCCGCCCTCGCCCGCATGCGGGTTGAGGCCGCACACGAGGATACGCGGGTTGGCAATGCCAAACTTGTCGCGCATGTCGGCATGCAGGATACGGGTTACACGTTCGACGCGCTCTGCAGTGATAGCGTCGGCAACATCACGCAGTGGCAGGTGCGTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/487583‑487811 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgtggCCCTGGTGACCACGCACCTGCCACTGCGTGATGTTGCCGACGCTATCACTGCAGAGCGCGTCGAACGTGTAACCCGTATCCTGCATGCCGACATGCGCGACAAGTTTGGCATTGCCAACCCGCGTATCCTCGTGTGCGGCCTCAACCCGCATGCGGGCGAGG  <  AE015451/487820‑487659

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