New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	3303564 =	53 (1.030)	3 (0.060)	3/144	5.3	5.9%	coding (798/1368 nt)	dcuD	predicted transporter
	?	REL606	= 4448682	44 (0.890)					coding (776/1422 nt)	treB	fused trehalose(maltose)‑specific PTS enzyme: IIB component/IIC component
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/3303640‑3303564 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcTGGTGGTGGTGCCCGTCTGTTCGCTGATCCTCGCGGTGTTCCTCGCCCATGCGCTGATTGGTCCGTTTGGTCGCA  <  REL606/4448682‑4448608                                                                                                                                                          TGCTGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCAT                                                                              >  2:9419/1‑76 TGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCAT                                                                              >  2:1325122/1‑76 TGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGAAGTTCTGACTGCAT                                                                              >  2:477438/1‑76  GCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCATC                                                                             <  2:1358507/76‑1                                               CACCACCACCAGATGCAGTTCTGACTGCATCTGGTGGTGGTG                                                                  >  2:879499/1‑42                                                  CACCACCAGATGCAGTTCTGACTGCATCTGGTGGTG                                                                     <  2:1281438/36‑1                                                  CACCACCAGATGCAGTTCTGACTGCATCTGGCGGTG                                                                     >  1:1281438/1‑36                                                                             CTGGTGGTGGTGCCCGTCTGTTCGCTGATCCTCGCGGTGTTCCTCGCCCATGCGCTGATTGGTCCGTTTGGTCGCA  <  2:1441969/76‑1                                                                                                                                                          TGCGGAAGAACTCAACAAACATCGTCACAGTCAAACTCAGTAACATCACCACCACCAGATGCAGTTCTGACTGCATC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/3303640‑3303564 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcTGGTGGTGGTGCCCGTCTGTTCGCTGATCCTCGCGGTGTTCCTCGCCCATGCGCTGATTGGTCCGTTTGGTCGCA  <  REL606/4448682‑4448608

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