New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2037554 =	68 (1.090)	3 (0.050)	3/136	6.6	4.9%	coding (577/1176 nt)	wzy	O antigen polymerase (Wzy)
	?	REL606	2037718 =	54 (0.950)					coding (413/1176 nt)	wzy	O antigen polymerase (Wzy)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAGTTAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTAAAA  >  REL606/2037479‑2037628                                                                            |                                                                           aaGTTAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAAt                                                                            <  1:1109934/76‑1 (MQ=255)  aGTTAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATa                                                                           <  1:1810871/76‑1 (MQ=255)   gTTAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATaa                                                                          <  2:1587461/76‑1 (MQ=255)    ttAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATaaa                                                                         <  1:1313821/76‑1 (MQ=255)     tAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATaaaa                                                                        <  1:672393/76‑1 (MQ=255)     tAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATaaaa                                                                        <  1:1640069/76‑1 (MQ=255)      aaTTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAAc                                                                       >  1:1629192/1‑76 (MQ=255)         tGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACata                                                                    >  1:1332630/1‑76 (MQ=255)          gCGTGCTGATAGCTTCGCCACTTCTGATATAATTCTTAATCTGTTCTACAAATGGCAGAGCCACAATAAAACatat                                                                   >  1:1502549/1‑76 (MQ=255)          gCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACatat                                                                   >  2:842595/1‑76 (MQ=255)          gCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACatat                                                                   <  2:958654/76‑1 (MQ=255)           cGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATAtt                                                                  >  2:1329262/1‑76 (MQ=255)            gTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTa                                                                 <  2:1404990/76‑1 (MQ=255)            gTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTa                                                                 <  2:1078807/76‑1 (MQ=255)            gTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTa                                                                 <  1:867061/76‑1 (MQ=255)            gTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTa                                                                 <  1:1003274/76‑1 (MQ=255)                tGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAAtac                                                             >  2:1323460/1‑76 (MQ=255)                tGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAAtac                                                             <  1:566021/76‑1 (MQ=255)                 gATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAAtact                                                            >  2:1079528/1‑76 (MQ=255)                 gATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAAtact                                                            >  2:167352/1‑76 (MQ=255)                  aTAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAAtacta                                                           <  2:126506/76‑1 (MQ=255)                    aGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTAcc                                                         >  1:1323817/1‑76 (MQ=255)                      cTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTAcccc                                                       >  1:108262/1‑76 (MQ=255)                        tCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAg                                                     >  2:1308087/1‑76 (MQ=255)                        tCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAg                                                     <  2:1935897/76‑1 (MQ=255)                         cGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGc                                                    >  2:660888/1‑76 (MQ=255)                           ccACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCtt                                                  >  1:460091/1‑76 (MQ=255)                           ccACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCtt                                                  >  1:1375072/1‑76 (MQ=255)                            cACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCttt                                                 >  1:1229325/1‑76 (MQ=255)                             aCTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTc                                                >  1:1751961/1‑76 (MQ=255)                             aCTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTc                                                >  1:1832798/1‑76 (MQ=255)                             aCTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTc                                                >  1:1559215/1‑76 (MQ=255)                               ttCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCa                                              <  2:1535353/76‑1 (MQ=255)                                 cTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAAt                                            <  2:1555111/76‑1 (MQ=255)                                   gATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAAttt                                          <  1:438410/76‑1 (MQ=255)                                      aGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAAtttttt                                       >  1:1717935/1‑76 (MQ=255)                                      aGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAAtttttt                                       >  1:1752881/1‑76 (MQ=255)                                       gAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAAttttttt                                      <  1:1959981/76‑1 (MQ=255)                                          ttCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGt                                   <  2:266782/76‑1 (MQ=255)                                           tCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGtt                                  >  1:1731355/1‑76 (MQ=255)                                             ttAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGtttt                                <  2:809839/76‑1 (MQ=255)                                             ttAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGtttt                                <  2:649203/76‑1 (MQ=255)                                              tAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGttttt                               >  1:630617/1‑76 (MQ=255)                                                  cTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAAt                           <  1:776886/76‑1 (MQ=255)                                                   tGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCGATc                          <  1:542624/76‑1 (MQ=255)                                                    gTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATcg                         >  2:1636445/1‑76 (MQ=255)                                                     ttCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATcga                        >  1:1208748/1‑76 (MQ=255)                                                      tCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATcgac                       >  1:676663/1‑76 (MQ=255)                                                        tACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATcgacga                     <  2:611677/76‑1 (MQ=255)                                                           aaaTGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGAcgatcg                  >  2:304495/1‑76 (MQ=255)                                                             aTGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGAcgatcgat                <  2:153656/76‑1 (MQ=255)                                                             aTGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGAcgatcgat                <  1:1276543/76‑1 (MQ=255)                                                             aTGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGAcgatcgat                <  1:824693/76‑1 (MQ=255)                                                             aTGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGAcgatcgat                <  1:824683/76‑1 (MQ=255)                                                              tGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATa               >  1:1144130/1‑76 (MQ=255)                                                              tGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATa               <  1:1278516/76‑1 (MQ=255)                                                               ggCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATaa              <  1:70138/76‑1 (MQ=255)                                                                 cAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATNGACGATCGATAAAt            <  1:1079528/76‑1 (MQ=255)                                                                 cAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGCTAAAt            >  2:1990994/1‑76 (MQ=255)                                                                    agCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATgag         >  1:402634/1‑76 (MQ=255)                                                                      ccAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGAt       >  1:7441/1‑76 (MQ=255)                                                                        aGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTa     <  2:375912/76‑1 (MQ=255)                                                                        aGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTa     <  2:1917932/76‑1 (MQ=255)                                                                        aGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTa     <  2:812727/76‑1 (MQ=255)                                                                        aGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTa     <  1:528121/76‑1 (MQ=255)                                                                        aGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTa     <  2:892885/76‑1 (MQ=255)                                                                          aaTAAACCATATTAATCCTCCCCAAGCTTTCCAA‑TTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTaaaa  <  2:524124/76‑1 (MQ=255)                                                                          aaTAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTaaa   <  2:385121/76‑1 (MQ=255)                                                                            |                                                                           AAGTTAATTGCGTGCTGATAGCTTCGCCACTTCTGATAGAATTCTTAATCTGTTCTACAAATGGCAGAGCCAGAATAAAACATATTAATACTACCCCAGCTTTCCAATTTTTTTTGTTTTTCAATCGACGATCGATAAATGAGATTAAAA  >  REL606/2037479‑2037628

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.