New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2038688 =	16 (0.310)	5 (0.100)	3/138	6.8	25.8%	coding (261/822 nt)	wbbB	putative glycosyltransferase WbbB
	?	REL606	2038966 =	14 (0.290)					coding (729/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

CTTTTTGAGTTTTCTGCCTATAAGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2038765‑2038688 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaattGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCT  >  REL606/2038966‑2039037                                                                                                                                                        CTTTTTGAGTTTTCTGCCTAT‑AGGAAGGGATAGGTTTTTTAGGAGATAGCTCATTATTACCATCAGGAGTACCAAT                                                                           >  1:1515957/1‑76 CTTTTTGAGTTTTCTGCCTAT‑AGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAAT                                                                           >  1:1515949/1‑76              CTGCCTATAAGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATTGATTGAAGGTG                                                               <  2:491835/76‑1              CTGCCTATAAGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATTGATTGAAGGTG                                                               <  2:491841/76‑1                     TAAGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATTGATTGAAGGTGATAACAA                                                        <  2:1139386/76‑1                     TAAGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATTGATTGAAGGTGATAACAA                                                        <  2:1139401/76‑1                         GAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATTGATTGAAGGTGATAACAAAGAA                                                    <  2:1050326/76‑1                                                                           gATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCT  <  1:1577939/75‑1                                                                                                                                                        CTTTTTGAGTTTTCTGCCTATAAGGAAGGAATAGATTTTTTAGAAGATAGCTCATTATTACCATCAGGAGTACCAATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2038765‑2038688 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaattGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGGGGATAGGTGAAAGAAGTCTAATTCT  >  REL606/2038966‑2039037

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