New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2039335 =	31 (0.540)	5 (0.100)	4/138	6.8	12.4%	coding (360/747 nt)	wbbA	putative glycosyltransferase WbbA
	?	REL606	2039563 =	42 (0.800)					coding (132/747 nt)	wbbA	putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AACTAGTAGATGAAGATCGTCGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2039411‑2039335 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatttgGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGCCG  >  REL606/2039563‑2039637                                                                                                                                                          AACTAGTAGATGAAGATCGTCGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTT                                                                              <  1:1659261/76‑1                 TCGTCGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTGGTCTTTATATTTAGA                                                              <  2:1716571/76‑1                   GTCGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTGGTCTTTATATTTAGATT                                                            <  2:788978/76‑1                     CGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTGGTCTTTATATTTAGATTCT                                                          <  2:1817078/76‑1                     CGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTGGTCTTTATATTTAGATTCT                                                          <  2:60581/76‑1                       AGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTGGTCTTTATATTTAGATTCTTT                                                        <  2:1384582/76‑1                                                                          TTTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGG     >  1:113820/1‑76                                                                           gTGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGC    >  2:1632434/2‑76                                                                            TGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGCC   >  1:1505783/1‑76                                                                            TGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGCC   >  1:272550/1‑76                                                                            TGGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGCC   >  1:272572/1‑76                                                                             GGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGCCG  >  1:1142367/1‑76                                                                                                                                                          AACTAGTAGATGAAGATCGTCGAGTTGGAAAAATTGATAAAAATGATTTAATATATTTTTTAGAAAATGATTATTTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2039411‑2039335 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatttgGTCTTTATATTTAGATTCTTTTATTGTATTTATAAGATTTATAAAACAATTTTCATGACTAAACCACTCAGGCCG  >  REL606/2039563‑2039637

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