New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2036917	50 (0.790)	3 (0.050)	3/140	8.0	6.8%	coding (42/1113 nt)	wbbC	putative glycosyltransferase WbbC
	?	REL606	= 2037192	36 (0.610)					coding (939/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAGTGTG  >  REL606/2036848‑2036988                                                                      |                                                                        cccATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCt                                                                   >  2:183517/1‑76 (MQ=255)  ccATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTg                                                                  >  1:723652/1‑76 (MQ=255)    aTCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaa                                                                >  1:1248268/1‑76 (MQ=255)    aTCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaa                                                                <  1:366154/76‑1 (MQ=255)     tCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaaa                                                               <  1:507946/76‑1 (MQ=255)       gAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGaaaaa                                                             >  1:855351/1‑76 (MQ=255)         atatTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGa                                                           >  2:1064881/1‑76 (MQ=255)          tatTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGAt                                                          <  1:357915/76‑1 (MQ=255)          tatTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGAt                                                          >  1:535131/1‑76 (MQ=255)           atTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATa                                                         >  1:468582/1‑76 (MQ=255)               cATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTaa                                                     <  2:1844936/76‑1 (MQ=255)               cATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTaa                                                     <  2:1844946/76‑1 (MQ=255)                aTTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTaaa                                                    <  2:468582/76‑1 (MQ=255)                  tAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGATAAAGATAGTAAAtt                                                  <  1:1522764/76‑1 (MQ=255)                  tAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAAtt                                                  <  1:1522791/76‑1 (MQ=255)                    atatTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGa                                                >  2:163545/1‑76 (MQ=255)                     tatTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATCCTGAAAAAGATAGTAAATTGaa                                               >  2:1636474/1‑76 (MQ=255)                        tttttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAAc                                            >  1:1547606/1‑76 (MQ=255)                        tttttGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAAc                                            >  1:45107/1‑76 (MQ=255)                              cAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTAtttt                                      <  1:1234101/76‑1 (MQ=255)                              cAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTAtttt                                      >  2:539407/1‑76 (MQ=255)                              cAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTAtttt                                      >  2:892172/1‑76 (MQ=255)                               aaTTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTc                                     <  1:1256420/76‑1 (MQ=255)                                  tAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTc                                  >  2:345681/1‑76 (MQ=255)                                   aaaTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCa                                 >  1:1173127/1‑76 (MQ=255)                                    aaTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCAt                                >  2:1137225/1‑76 (MQ=255)                                     aTCTCTAAAGTATACTCCAACCCCCTCATAATTTAATTCTGCAAAAGATAGAAAATTGAAAAATCTTTTCGTCATa                               >  1:718964/1‑76 (MQ=255)                                     aTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATa                               <  1:1724227/76‑1 (MQ=255)                                     aTCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATa                               <  1:1724183/76‑1 (MQ=255)                                       ctctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATaaa                             >  1:750376/1‑76 (MQ=255)                                        tctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAg                            >  2:908004/1‑76 (MQ=255)                                         ctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAgg                           >  1:311089/1‑76 (MQ=255)                                         ctAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTACTTTCGTCATAAAgg                           >  2:1437760/1‑76 (MQ=255)                                           naaGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGtt                         >  1:367211/2‑76 (MQ=255)                                               tataCTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTAc                     >  2:1095937/1‑76 (MQ=255)                                                ataCTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTAct                    >  2:723654/1‑76 (MQ=255)                                                 taCTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTActc                   >  1:743040/1‑76 (MQ=255)                                                 taCTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTActc                   >  2:237439/1‑76 (MQ=255)                                                  aCTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTActct                  >  1:140207/1‑76 (MQ=255)                                                   cTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCtt                 <  2:1498740/76‑1 (MQ=255)                                                   cTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCtt                 >  1:899904/1‑76 (MQ=255)                                                     ccAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCtttt               <  2:1199494/76‑1 (MQ=255)                                                     ccAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCtttt               >  1:1633721/1‑76 (MQ=255)                                                        aCCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGt            >  1:1792195/1‑76 (MQ=255)                                                         ccccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTc           >  1:1252207/1‑76 (MQ=255)                                                         ccccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTc           <  1:389924/76‑1 (MQ=255)                                                           ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt         >  2:703686/1‑76 (MQ=255)                                                           ccATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCt         >  2:703691/1‑76 (MQ=255)                                                            cATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTa        >  2:1177862/1‑76 (MQ=255)                                                                  atTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAgtgtg  <  1:1384324/76‑1 (MQ=255)                                                                      |                                                                        CCCATCGAATATTTCATTAATATTTTTTGCAATTAAATCTCTAAAGTATACTCCAACCCCATAATATTTTAATTCTGAAAAAGATAGTAAATTGAAAACTATTTTCGTCATAAAGGTTTTACTCTTTTAGTCCTAAGTGTG  >  REL606/2036848‑2036988

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