New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2037392	28 (0.590)	9 (0.210)	7/136	4.0	32.5%	coding (739/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	= 2037686	12 (0.280)					coding (445/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.

ATAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACA  >  REL606/2037613‑2037753                                                                          |                                                                    atAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGa                                                                   <  1:1482144/76‑1 (MQ=255)         gaTTAAAATCAAAATGAAAAAATAAGATGTCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGataatatta                                                           <  2:504323/76‑4 (MQ=255)                aTCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGt                                                    >  1:414225/1‑76 (MQ=255)                    aaaTGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGtactt                                                <  2:908021/76‑1 (MQ=255)                          aaaaaTAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGt                                          <  1:1489542/76‑1 (MQ=255)                          aaaaaTAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGt                                          <  1:1489545/76‑1 (MQ=255)                             aaTAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGtttt                                       <  2:360401/76‑1 (MQ=255)                                aaGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAAc                                    >  1:1421362/1‑76 (MQ=255)                                    tatCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACAtgt                                >  1:1158220/1‑76 (MQ=255)                                              cGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAata                      <  1:1283480/76‑1 (MQ=255)                                              cGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAata                      <  1:1336818/76‑1 (MQ=255)                                              cGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAata                      <  1:540757/76‑1 (MQ=255)                                                    agaaaagaaaTAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTaaaa                <  1:1521266/76‑1 (MQ=255)                                                                  ccACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACa  >  2:319658/1‑76 (MQ=255)                                                                          |                                                                    ATAAATGAGATTAAAATCAAAATGAAAAAATAAGATATCCAACCACGCTGAAGAAAAGAAATAATCCACAACATGATAATAATGGGCTTGTACTTACTTGTTTTAACATGTGTTAGAAATATTAAAACAAGTACATCAACA  >  REL606/2037613‑2037753

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