New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3670807 =	69 (0.670)	3 (0.030)	3/252	NT	3.8%	coding (1216/1650 nt)	treF	cytoplasmic trehalase
	?	NC_000913	3670847 =	86 (0.870)					coding (1256/1650 nt)	treF	cytoplasmic trehalase
Rejected: Frequency below/above cutoff threshold.

ACAATGGCGGCAGCGGAAAACAGCGCCAGTTGTTCGCGTCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3670850‑3670807 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccaTTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCA  >  NC_000913/3670847‑3670992                                                                                                                                                                                                ACAATGGCGGCCGAGGAAAACAGCGCCAGTTGTTCGCGTCGCCATGGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGT                                           >  4:665261/1‑149    ATGGCGGCAGCGGAAAACAGCGCCAGTTGTTCGCGTCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAA                                                                                                                 <  1:213707/76‑1    ATGGCGGCAGCGGAAAACAGCGCCAGTTGTTCGCGTCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAA                                                                                                                 >  2:213707/1‑76                                         GCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAGACCGGTGAACAGTGGGATAAACCCAACGGCTGGGC   >  3:836851/1‑149                                         GCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGC   >  2:195696/1‑149                                          CCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGGACAGTGGGATAAACCCAACGGCTGGGCA  >  2:606029/1‑149                                            ATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTG                                                     >  3:364212/1‑96                                            ATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTG                                                     <  4:364212/96‑1                                                                                                                                                                                                ACAATGGCGGCAGCGGAAAACAGCGCCAGTTGTTCGCGTCGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3670850‑3670807 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccaTTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACGAAACCGGTGAACAGTGGGATAAACCCAACGGCTGGGCA  >  NC_000913/3670847‑3670992

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