New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2037848 =	51 (1.040)	28 (0.620)	13/140	2.4	29.5%	coding (283/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	2038306 =	86 (1.890)					coding (643/822 nt)	wbbB	putative glycosyltransferase WbbB

AAAAATAAAACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2037927‑2037848 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑attatGTAACAATTTAGCTTGATACCAAA  >  REL606/2038306‑2038329                                                                                                          AAAAATAAAACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAA                              >  2:217028/1‑76     ATAAAACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTAT                          >  1:82813/1‑76      TAAAACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATG                         <  2:1147672/76‑1      TAAAACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATG                         <  2:1147691/76‑1         AACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAA                      <  2:500481/76‑1         AACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAA                      <  2:500479/76‑1            CACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAA                   <  2:1075673/76‑1             ACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAAT                  <  2:1336429/76‑1                TACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTA               <  2:1157636/76‑1                TACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTA               <  2:1157648/76‑1                TACAACAAAAGTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTA               <  2:1157628/76‑1                 ACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAG              <  2:89573/76‑1                 ACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAG              <  2:499282/76‑1                 ACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAG              <  2:499278/76‑1                 ACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAG              <  2:1569927/76‑1                  CAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGC             <  2:1439001/76‑1                  CAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGC             <  2:350281/76‑1                   AACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCT            <  2:1231345/76‑1                    ACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTT           <  2:882336/76‑1                      AAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGA         <  2:129161/76‑1                        AATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATA       <  2:1445220/76‑1                        AATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATA       <  2:1445237/76‑1                        AATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATA       <  2:415352/76‑1                        AATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATA       <  2:415369/76‑1                        AATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATCTAGCTTGATA       <  2:415357/76‑1                            AGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATACCAA   <  2:1137076/76‑1                            AGCTGATTCTAAGAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATACCAA   <  2:1137094/76‑1                             GCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATACCAAA  <  2:1491564/76‑1                             GCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATACCAAA  <  2:1301233/76‑1                             GCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTATGTAACAATTTAGCTTGATACCAAA  <  2:1301221/76‑1                                                                                                          AAAAATAAAACCACTTACAACAAAATTAGCTGATTCTAATAGTAGACCGATCGCCATTTTTTTAATTTGTTTAGAATTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2037927‑2037848 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑attatGTAACAATTTAGCTTGATACCAAA  >  REL606/2038306‑2038329

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