New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2036832 =	54 (0.920)	7 (0.130)	6/136	5.5	12.4%	coding (127/1113 nt)	wbbC	putative glycosyltransferase WbbC
	?	REL606	2037060 =	50 (0.940)					coding (1071/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAGCATTGT  >  REL606/2036987‑2037131                                                                          |                                                                        tgtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCAGGAAATaaaaa                                                                       <  1:1489994/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                      <  1:317751/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                      <  1:979216/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                      <  1:317733/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                      <  1:1491756/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                      <  1:580291/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAac                                                                      <  1:317754/76‑1 (MQ=255)  gtACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCCGTGAATTGGAATATCCAGCCCGGAAATAAAAac                                                                      <  1:317739/76‑1 (MQ=255)   tACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAaca                                                                     <  1:1359642/76‑1 (MQ=255)   tACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAaca                                                                     <  1:891698/76‑1 (MQ=255)   tACTTTAAT‑AAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAacac                                                                    <  1:453702/76‑1 (MQ=255)    aCTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAacac                                                                    <  1:453685/76‑1 (MQ=255)    aCTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAacac                                                                    <  1:936637/76‑1 (MQ=255)    aCTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAacac                                                                    <  1:453679/76‑1 (MQ=255)     cTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACt                                                                   >  1:382729/1‑76 (MQ=255)       ttAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTaa                                                                 <  1:1857753/76‑1 (MQ=255)       ttAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTaa                                                                 <  2:1163736/76‑1 (MQ=255)        taataaAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAAc                                                                >  2:957307/1‑76 (MQ=255)        taataaAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAAc                                                                <  1:952209/76‑1 (MQ=255)         aataaaataaATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAAcc                                                               <  1:403917/76‑1 (MQ=255)            aaaataaaTACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCaaa                                                            <  1:1745884/76‑1 (MQ=255)              aataaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaac                                                          <  1:1375456/76‑1 (MQ=255)              aataaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaac                                                          <  1:1375452/76‑1 (MQ=255)                aaaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAAAACCCAACCCGGAAATAAAAACACAAACCAAaacaa                                                        <  1:677162/75‑1 (MQ=255)                 aaataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAaacaac                                                       <  1:1862734/76‑1 (MQ=255)                  aataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAAcc                                                      <  2:1201291/76‑1 (MQ=255)                   ataCTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCa                                                     >  2:1502381/1‑76 (MQ=255)                       tAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCAtctc                                                 <  1:1103913/76‑1 (MQ=255)                         aaaGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCat                                               <  1:819680/76‑1 (MQ=255)                         aaaGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCat                                               <  1:819684/76‑1 (MQ=255)                           aGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCatat                                             <  2:555176/76‑1 (MQ=255)                             cGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTc                                           <  1:118548/76‑1 (MQ=255)                              gAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCa                                          <  2:1087781/76‑1 (MQ=255)                                 aaTATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTg                                       <  2:1774663/76‑1 (MQ=255)                                      aaCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCatat                                  >  2:1218614/1‑76 (MQ=255)                                      aaCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCatat                                  >  1:1371833/1‑76 (MQ=255)                                      aaCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCatat                                  >  1:1371839/1‑76 (MQ=255)                                       aCTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCAtata                                 <  2:1315462/76‑1 (MQ=255)                                        ctgctgTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAg                                <  1:912151/76‑1 (MQ=255)                                        ctgctgTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAg                                >  2:516702/1‑76 (MQ=255)                                           ctgTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTc                             >  1:1672889/1‑76 (MQ=255)                                               gAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTc                         >  1:1606864/1‑76 (MQ=255)                                                  ttGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCa                      <  1:445203/76‑1 (MQ=255)                                                  ttGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCa                      <  2:1545954/76‑1 (MQ=255)                                                   tGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCaa                     >  2:1120567/1‑76 (MQ=255)                                                     gAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCaata                   <  2:731248/76‑1 (MQ=255)                                                         atCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAAttt               <  2:187753/76‑1 (MQ=255)                                                          tCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAAtttt              <  2:217915/76‑1 (MQ=255)                                                          tCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAAtttt              <  2:217916/76‑1 (MQ=255)                                                          tCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAAtttt              <  2:217935/76‑1 (MQ=255)                                                                    aaaTAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAGCAtt    <  1:1501633/76‑1 (MQ=255)                                                                      aTAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAGCATTGt  <  1:701875/76‑1 (MQ=255)                                                                          |                                                                        TGTACTTTAATAAAATAAATACTAAAAGCGAAAATATAACTGCTGTGAATTGGAATATCCAACCCGGAAATAAAAACACTAACCAAAACAACCATCTCATATTCATTGCATATAGGTCATTCGCAATAATTTTAGATAGCATTGT  >  REL606/2036987‑2037131

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