New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2036832 =	50 (0.850)	8 (0.150)	3/136	6.3	15.5%	coding (127/1113 nt)	wbbC	putative glycosyltransferase WbbC
	?	REL606	2037060 =	42 (0.790)					coding (1071/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATTATAGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCC  >  REL606/2036758‑2036906                                                                           |                                                                           aTTATAGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACtatt                                                                           >  1:1449648/1‑76 (MQ=255)    ataGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTattatt                                                                        >  1:598143/1‑76 (MQ=255)     taGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTAttt                                                                       >  2:1274798/1‑76 (MQ=255)      aGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTc                                                                      <  1:858223/76‑1 (MQ=255)       ggTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTca                                                                     <  2:155175/76‑1 (MQ=255)        gTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTcac                                                                    >  1:914885/1‑76 (MQ=255)                 ttGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATcc                                                           >  1:1419251/1‑76 (MQ=255)                   gTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCa                                                         <  1:1309101/76‑1 (MQ=255)                   gTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCa                                                         >  1:1715101/1‑76 (MQ=255)                      ttCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCg                                                      <  1:1701051/76‑1 (MQ=255)                       tCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTGCACAATCAATCCCATCGa                                                     >  2:1814248/1‑76 (MQ=255)                        cTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGaa                                                    >  1:341241/1‑76 (MQ=255)                         tGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAat                                                   >  2:225696/1‑76 (MQ=255)                         tGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAat                                                   <  1:848338/76‑1 (MQ=255)                         tGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAat                                                   <  1:1641/76‑1 (MQ=255)                         tGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAat                                                   >  2:1114346/1‑76 (MQ=255)                          ggAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAata                                                  >  2:56739/1‑76 (MQ=255)                            aaaaTTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATAtt                                                <  2:247840/76‑1 (MQ=255)                             aaaTTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATAttt                                               <  1:357846/76‑1 (MQ=255)                             aaaTTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATAttt                                               >  2:1082143/1‑76 (MQ=255)                                   aaCAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTaa                                         <  2:1168738/76‑1 (MQ=255)                                   aaCAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTaa                                         >  2:369279/1‑76 (MQ=255)                                         tCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATAtttt                                   <  2:592970/76‑1 (MQ=255)                                           tttAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATAtttttt                                 >  2:763217/1‑76 (MQ=255)                                           ttnAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATAtttttt                                 >  2:1324994/1‑76 (MQ=255)                                            ttAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTg                                >  1:1762215/1‑76 (MQ=255)                                              aaCAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCa                              <  1:507818/76‑1 (MQ=255)                                                cAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAAt                            >  2:1778352/1‑76 (MQ=255)                                                cAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAAt                            >  2:344232/1‑76 (MQ=255)                                                 aGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAAtt                           <  2:1313693/76‑1 (MQ=255)                                                    gATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTaaa                        <  1:1138374/76‑1 (MQ=255)                                                      tatTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAAtc                      <  2:152026/76‑1 (MQ=255)                                                         tAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATctct                   <  1:218290/76‑1 (MQ=255)                                                         tAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATctct                   >  2:1868376/1‑76 (MQ=255)                                                          aGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTa                  <  2:194611/76‑1 (MQ=255)                                                          aGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTa                  <  2:1217286/76‑1 (MQ=255)                                                                gTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGtat            <  1:520284/76‑1 (MQ=255)                                                                   aacaacTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACt         <  2:1269375/76‑1 (MQ=255)                                                                    acaacTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTc        >  1:428365/1‑76 (MQ=255)                                                                    acaacTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTc        >  1:1188434/1‑76 (MQ=255)                                                                     caacTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTcc       <  1:1260015/76‑1 (MQ=255)                                                                      aacTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCa      <  1:1226572/76‑1 (MQ=255)                                                                      aacTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCa      <  1:51792/76‑1 (MQ=255)                                                                       acTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCaa     >  2:988214/1‑76 (MQ=255)                                                                       acTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCaa     >  1:594332/1‑76 (MQ=255)                                                                       acTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCaa     <  1:1170900/76‑1 (MQ=255)                                                                        cTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAAc    >  2:1285896/1‑76 (MQ=255)                                                                         tattatTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAAcc   >  1:568177/1‑76 (MQ=255)                                                                          attattTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAAccc  <  2:175060/76‑1 (MQ=255)                                                                          attattTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAAccc  >  1:148058/1‑76 (MQ=255)                                                                           |                                                                           ATTATAGGTAATATTATTGTTTTCTGGAAAATTGAACAATTCTTTAACAGGGATATTAGCCTGGTGAACAACTATTATTTCACAATCAATCCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCC  >  REL606/2036758‑2036906

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