New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2036917	56 (1.030)	9 (0.180)	7/140	5.1	18.3%	coding (42/1113 nt)	wbbC	putative glycosyltransferase WbbC
	?	REL606	= 2037192	28 (0.550)					coding (939/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGATAGCATTGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGAAATATATCTGGCAAGCAA  >  REL606/2037120‑2037266                                                                         |                                                                           aGATAGCATTGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATaaaa                                                                         >  1:938975/1‑76 (MQ=255)    tAGCATTGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATc                                                                      >  1:1479408/1‑76 (MQ=255)     aGCATTGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAAtat                                                                     <  2:1077241/76‑3 (MQ=255)       cATTGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATccaa                                                                   <  1:636509/76‑1 (MQ=255)          tGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATccaacca                                                                <  1:684179/76‑1 (MQ=255)          tGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATccaacca                                                                >  1:41207/1‑76 (MQ=255)          tGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATccaacca                                                                >  2:232192/1‑76 (MQ=255)             aaaTAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAAtt                                                             <  1:646780/76‑1 (MQ=255)             aaaTAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAAtt                                                             <  1:548288/76‑1 (MQ=255)                tAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAAc                                                          <  1:901137/76‑1 (MQ=255)                   aaaaCACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGAt                                                       >  2:690700/1‑76 (MQ=255)                          aGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTca                                                <  2:343219/76‑1 (MQ=255)                          aGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTca                                                <  2:401162/76‑1 (MQ=255)                           gCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTcac                                               >  1:239249/1‑76 (MQ=255)                            cAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTcaca                                              >  2:183569/1‑76 (MQ=255)                              aaaCAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACAtt                                            <  2:117836/76‑1 (MQ=255)                                  aacCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCAt                                        >  2:92315/1‑76 (MQ=255)                                        tatGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTa                                  >  2:906377/1‑76 (MQ=255)                                         atgatgGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTAt                                 >  1:1561966/1‑76 (MQ=255)                                          tgatgGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATc                                >  1:493228/1‑76 (MQ=255)                                           gatgGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATcc                               >  1:494262/1‑76 (MQ=255)                                            atgGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCa                              >  2:428117/1‑76 (MQ=255)                                            atgGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCa                              >  2:428132/1‑76 (MQ=255)                                                     aCGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGgag                     >  1:1208255/1‑76 (MQ=255)                                                      cGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGgaga                    >  1:1113771/1‑76 (MQ=255)                                                        tataAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGaaa                  <  1:906377/76‑1 (MQ=255)                                                              tCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGAAATATATc            <  2:194791/76‑1 (MQ=255)                                                                aTTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGAAATATATCTg          <  1:1269250/76‑1 (MQ=255)                                                                     gAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGAAATATATCTGgcaag     >  2:670179/1‑76 (MQ=255)                                                                        tAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGAAATATATCTGgcaagcaa  <  2:494262/76‑1 (MQ=255)                                                                         |                                                                           AGATAGCATTGTAAATAGAAAACACAGCAAAACAACCCATATGATGGGTAGAACGTATAAATCATTTTGAATAAAAATCCAACCAATTAACGATGGATTCACATTCCATGAGCTATCCAGAGATGGAGAAATATATCTGGCAAGCAA  >  REL606/2037120‑2037266

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