New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2037842 =	57 (0.970)	3 (0.060)	3/134	6.2	5.5%	coding (289/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	2037944 =	52 (0.990)					coding (187/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATGGCGATCGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAAC  >  REL606/2037871‑2038017                                                                          |                                                                          ttGGCGATCGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATataata                                                                         >  2:862235/2‑76 (MQ=255)  tGGCGATCGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATAc                                                                        <  1:357293/76‑1 (MQ=255)         cGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGa                                                                 <  2:102211/76‑1 (MQ=255)         cGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGa                                                                 >  2:1506369/1‑76 (MQ=255)         cGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGa                                                                 <  2:724336/76‑1 (MQ=255)         cGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGa                                                                 <  2:724332/76‑1 (MQ=255)          ggTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGaa                                                                <  2:970735/76‑1 (MQ=255)           gTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTCTATAATACCGATAGaaa                                                               >  2:97900/1‑76 (MQ=255)            tCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaat                                                              <  1:775608/76‑1 (MQ=255)            tCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaat                                                              <  1:230059/76‑1 (MQ=255)             ctactaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaata                                                             <  2:742589/76‑1 (MQ=255)              tactaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaataa                                                            >  1:753676/1‑76 (MQ=255)              tactaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaataa                                                            <  1:12697/76‑1 (MQ=255)               actaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAaataat                                                           >  1:1876129/1‑76 (MQ=255)                ctaTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGACAGAAataata                                                          <  2:1800785/76‑1 (MQ=255)                 naTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAAtaataa                                                         >  1:1266347/2‑76 (MQ=255)                  aTTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAg                                                        >  1:318274/1‑76 (MQ=255)                      gAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTaa                                                    >  2:1583355/1‑76 (MQ=255)                          cAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAgaga                                                >  2:347559/1‑76 (MQ=255)                           ctcTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGat                                               >  2:1549113/3‑76 (MQ=255)                                aTTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGAtatata                                          <  1:1189443/76‑1 (MQ=255)                                aTTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGAtatata                                          <  1:378158/76‑1 (MQ=255)                                  tttGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATaaa                                        <  2:1488249/76‑1 (MQ=255)                                   ttgttgTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAAt                                       <  2:560226/76‑1 (MQ=255)                                   ttgttgTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAAt                                       <  1:531681/76‑1 (MQ=255)                                    tgttgtAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAAtt                                      >  1:996610/1‑76 (MQ=255)                                       tgtACGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTaga                                   >  2:1658644/1‑76 (MQ=255)                                        gtAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTagag                                  <  1:1549113/76‑1 (MQ=255)                                         tAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAgaga                                 <  1:701072/76‑1 (MQ=255)                                         tAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAgaga                                 <  1:701073/76‑1 (MQ=255)                                         tAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAgaga                                 >  1:47008/1‑76 (MQ=255)                                           aGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGAtt                               <  1:372543/76‑1 (MQ=255)                                             tGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGa                             <  1:1526393/76‑1 (MQ=255)                                               gTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAat                           >  1:141303/1‑76 (MQ=255)                                               gTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAat                           >  1:141308/1‑76 (MQ=255)                                                 tttATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAatta                         >  2:730954/1‑76 (MQ=255)                                                    aTTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAAttatta                      <  1:421344/76‑1 (MQ=255)                                                      tttttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACg                    <  1:347559/76‑1 (MQ=255)                                                      tttttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACg                    <  2:222191/76‑1 (MQ=255)                                                        tttCCAAGCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGAc                  <  1:140559/76‑1 (MQ=255)                                                         ttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACa                 >  1:515042/1‑76 (MQ=255)                                                         ttCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACa                 >  1:1193528/1‑76 (MQ=255)                                                          tCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                >  2:121851/1‑76 (MQ=255)                                                          tCCAATCCATTATATAACACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAg                <  1:935814/76‑1 (MQ=255)                                                             aaTCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGttt             >  1:1170418/1‑76 (MQ=255)                                                             aaTCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGttt             >  1:137168/1‑76 (MQ=255)                                                             aaTCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGttt             >  1:137154/1‑76 (MQ=255)                                                              aTCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTc            <  2:917726/76‑1 (MQ=255)                                                               tCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCt           <  1:670009/76‑1 (MQ=255)                                                                ccATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTg          <  2:1810144/76‑1 (MQ=255)                                                                ccATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTg          <  2:363852/76‑1 (MQ=255)                                                                 cATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGa         >  1:1139335/1‑76 (MQ=255)                                                                   ttATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGAtt       >  2:1686152/1‑76 (MQ=255)                                                                      tataATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTaa    <  1:1323149/76‑1 (MQ=255)                                                                      tataATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTaa    >  2:1524022/1‑76 (MQ=255)                                                                       ataataCCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTaaa   >  1:1001730/1‑76 (MQ=255)                                                                        taataCCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAAc  >  1:593082/1‑76 (MQ=255)                                                                          |                                                                          ATGGCGATCGGTCTACTATTAGAATCAGCTAATTTTGTTGTAAGTGGTTTTATTTTTTCCAATCCATTATATAATACCGATAGAAATAATAAGCTAAGAGATATATAAATTAGAGATTGAATTATTACGACAGTTTCTGATTTAAAC  >  REL606/2037871‑2038017

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.