New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2037392	34 (0.510)	5 (0.080)	5/134	7.6	14.8%	coding (739/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	= 2037686	27 (0.450)					coding (445/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGT  >  REL606/2037613‑2037759                                                                          |                                                                          atAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATg                                                                          <  2:724579/75‑1 (MQ=255)  tAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGa                                                                         >  1:11981/1‑75 (MQ=255)  tAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGa                                                                         >  2:1595755/1‑75 (MQ=255)  tAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGa                                                                         <  1:236510/75‑1 (MQ=255)  tAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGa                                                                         >  1:252422/1‑75 (MQ=255)         gaTTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATaataat                                                                  <  1:400358/75‑1 (MQ=255)            tAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATggg                                                               >  1:4833/1‑75 (MQ=255)             aaaaTCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATg                                                                 <  2:1257079/72‑1 (MQ=255)             aaaaTCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATg                                                                 >  1:1257079/1‑72 (MQ=255)             aaaaTCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGc                                                              <  1:353360/75‑1 (MQ=255)               aaTCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCtt                                                            <  1:129791/75‑1 (MQ=255)                aTCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTg                                                           <  2:1361221/75‑1 (MQ=255)                   aaaaTGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGtac                                                        <  1:1065989/75‑1 (MQ=255)                        gAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGtacttact                                                   >  2:433874/1‑75 (MQ=255)                            aaaTAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGtt                                               >  1:1028237/1‑75 (MQ=255)                            aaaTAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGtt                                               <  2:1054183/75‑1 (MQ=255)                            aaaTAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACACGACAATAATGGGCTTGTACTTACTTGtt                                               >  1:913096/1‑75 (MQ=255)                                    tatCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACAtg                                       >  2:73255/1‑75 (MQ=255)                                    tatCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACAtg                                       >  2:808179/1‑75 (MQ=255)                                         aaCCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTa                                  <  2:286205/75‑1 (MQ=255)                                            cACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGaa                               <  1:819770/75‑1 (MQ=255)                                              cGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAat                             >  2:702987/1‑75 (MQ=255)                                              cGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAat                             >  1:1073674/1‑75 (MQ=255)                                                cTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAatat                           <  1:1099048/75‑1 (MQ=255)                                                 tGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATAtt                          >  1:1740451/1‑75 (MQ=255)                                                          gaaaTAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGt                 <  2:1654057/75‑1 (MQ=255)                                                             ataatCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACa              >  1:1753176/1‑75 (MQ=255)                                                              taatCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACAt             >  2:1318326/1‑75 (MQ=255)                                                              taatCCACAACATGATAATAATGGGCTTGCACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACAt             >  2:1318347/1‑75 (MQ=255)                                                                atCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCa           >  1:2134296/1‑75 (MQ=255)                                                                        caTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAg   <  2:101419/75‑1 (MQ=255)                                                                         aTGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGt  >  1:347285/1‑75 (MQ=255)                                                                          |                                                                          ATAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACATTCAGT  >  REL606/2037613‑2037759

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