New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2037854	32 (0.670)	6 (0.140)	6/132	4.3	18.7%	coding (277/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	= 2037988	24 (0.570)					coding (143/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACTGTTGTTTATAA  >  REL606/2037916‑2038057                                                                         |                                                                      ggTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGaa                                                                    >  1:649118/1‑76 (MQ=255)  gTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAat                                                                   <  1:701837/76‑1 (MQ=255)   ttttATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAatt                                                                  <  1:725417/76‑1 (MQ=255)    tttATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAatta                                                                 >  2:487248/1‑76 (MQ=255)      tATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAattatt                                                               >  1:628391/1‑76 (MQ=255)          ttttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGa                                                           >  1:17668/1‑76 (MQ=255)           tttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAACAGATATTTAAATTAGACATTTCATTATTAcacc                                                          >  1:385246/1‑73 (MQ=255)            ttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACa                                                         >  2:979632/1‑76 (MQ=255)            ttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACa                                                         <  2:446045/76‑1 (MQ=255)             tCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                                                        <  1:1144628/76‑1 (MQ=255)             tCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                                                        <  1:710305/76‑1 (MQ=255)             tCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                                                        <  1:563830/76‑1 (MQ=255)             tCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                                                        <  1:1144618/76‑1 (MQ=255)             tCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                                                        <  1:563849/76‑1 (MQ=255)                         tataATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTaa                                            >  2:636393/1‑76 (MQ=255)                           taataCCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAAc                                          <  1:1223955/76‑1 (MQ=255)                           taataCCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAAc                                          <  1:1223964/76‑1 (MQ=255)                           taataCCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAAc                                          <  1:1223961/76‑1 (MQ=255)                                   aTAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTaa                                  >  1:573915/1‑76 (MQ=255)                                          taataaGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCca                           >  2:1117186/1‑76 (MQ=255)                                          taataaGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCca                           >  2:1117177/1‑76 (MQ=255)                                                gCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCcaatcaat                     <  1:1347627/76‑1 (MQ=255)                                                          atataAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACtgt           >  2:1061455/1‑76 (MQ=255)                                                            ataAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCACTAATTCACtgttg         >  2:1169203/1‑76 (MQ=255)                                                             taAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACtgttgt        >  2:728479/1‑76 (MQ=255)                                                              aaaTTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACTgttgtt       >  2:285169/1‑76 (MQ=255)                                                                 ttAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACTGTTGTTtat    <  1:364144/76‑1 (MQ=255)                                                                 ttAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACTGTTGTTtat    >  1:519046/1‑76 (MQ=255)                                                                   agagATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACTGTTGTTTATaa  <  2:996639/76‑1 (MQ=255)                                                                         |                                                                      GGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAACCCCATTAAATCGCCAATCAATAATTCACTGTTGTTTATAA  >  REL606/2037916‑2038057

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