New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 4355162	70 (0.980)	3 (0.040)	3/334	NT	4.2%	coding (362/681 nt)	modB	Molybdenum transport system permease protein ModB
	?	AE015451	= 4355227	69 (1.010)					coding (427/681 nt)	modB	Molybdenum transport system permease protein ModB
Rejected: Frequency below/above cutoff threshold.

GCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4355112‑4355162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccagcgGCAGCACCACATGAACGAAGGTATCCCATGGGCTGGCTCGCAGAGTAGCGGCGACTTCCAGGGGGCGCTGGCCGATCGCGCCAAAGGCATTCTGCAGGGGTTGCACCACGAACGGCATGGAGTACACCGTCGAGCCAATCACCAGGCCTGTGAAGCTGAACACCACACTGCCCAGGCCTAGCGCCTGGGTGGCCTGACCAAGCCAGCCGTG  <  AE015451/4355227‑4355017                                                                                                                                                                                                                                                                        GCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGGCAGCACCACATGAACGAAGGTATCCCATGGGCTGGCTCGCAGAGTAGCGGCGACTTCCAGGGGGCGCTGGCCGATCGCGCCAAAGGCATTCTGCAGGGGTTGCACCACG                                                                                                       >  1:36528/1‑161    GTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGGCAGCACCACATGAACGAAGGTATCCCATGGGCTGGCTCGCAGAGTAGCGGCGACTTCCAGGGGGCGCTGGCCGATCGCGCCAAAGGCATTCTGCAGGGGCTGCACCACGAACG                                                                                                   >  1:740924/1‑162        GTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCGGCAGCACCACATGAACGAAGGTATCCCATGGGCTGGCTCGCAGAGGAGC                                                                                                                                                                    >  2:136490/1‑93                                                 GCGGCAGCACCACATGAACGAAGGTATCCCATGGGCTGGCTCGCAGAGTAGCGGCGACTTCCAGGGGGCGCTGGCCGATCGCGCCAAAGGCATTCTGCAGGGGTTGCACCACGAACGGCATGGAGTACACCGTCGAGCCAATCACCAGGCCTGTGAAGCTGAACACCACACTGCCCAGGCCTAGCGCCTGGGTGGCCTGACCAAGCCAGCCGTG  >  2:706020/1‑214                                                                                                                                                                                                                                                                        GCCGTTCGTGGTGCAACCCCTGCAGAATGCCTTTGGCGCGATCGGCCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AE015451/4355112‑4355162 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccagcgGCAGCACCACATGAACGAAGGTATCCCATGGGCTGGCTCGCAGAGTAGCGGCGACTTCCAGGGGGCGCTGGCCGATCGCGCCAAAGGCATTCTGCAGGGGTTGCACCACGAACGGCATGGAGTACACCGTCGAGCCAATCACCAGGCCTGTGAAGCTGAACACCACACTGCCCAGGCCTAGCGCCTGGGTGGCCTGACCAAGCCAGCCGTG  <  AE015451/4355227‑4355017

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