New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2042256 =69 (1.060)5 (0.080) 4/140 7.7 6.7% coding (565/1428 nt) wzx O antigen translocase (Wzx)
?REL606 2042511 = 76 (1.260)coding (310/1428 nt) wzx O antigen translocase (Wzx)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAACCCAAAGC  >  REL606/2042183‑2042331
                                                                         |                                                                           
ttACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAatat                                                                           >  1:1329493/1‑76 (MQ=255)
  aCTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTc                                                                         >  1:1120165/1‑76 (MQ=255)
  aCTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTc                                                                         >  1:1842954/1‑76 (MQ=255)
    tGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAg                                                                       <  2:1604072/76‑1 (MQ=255)
      cTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATAGCATGGATTATAAATATTCAGta                                                                     >  2:1671414/1‑76 (MQ=255)
      cTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGta                                                                     >  2:1671428/1‑76 (MQ=255)
       tAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtac                                                                    >  1:609574/1‑76 (MQ=255)
        aaaTGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtaca                                                                   <  1:2006263/76‑1 (MQ=255)
          aTGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaat                                                                 >  1:151813/1‑76 (MQ=255)
           tGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaata                                                                <  1:2080277/76‑1 (MQ=255)
             aaaaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaataca                                                              <  1:1057531/76‑1 (MQ=255)
             aaaaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaataca                                                              <  1:1057547/76‑1 (MQ=255)
             aaaaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaataca                                                              >  2:102928/1‑76 (MQ=255)
             aaaaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaataca                                                              <  1:1544513/76‑1 (MQ=255)
             aaaaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaataca                                                              <  1:1544512/76‑1 (MQ=255)
              aaaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGtacaatacaa                                                             <  1:583161/76‑1 (MQ=255)
               aaTAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTacaatacaat                                                            <  1:558895/76‑1 (MQ=255)
                 tAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAAttt                                                          <  1:1568694/76‑1 (MQ=255)
                 tAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAAttt                                                          <  1:1568705/76‑1 (MQ=255)
                  aGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATCATAAATATTCAGTACAATACAATTTa                                                         >  1:1187775/1‑76 (MQ=255)
                   gCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTAt                                                        >  1:1057001/1‑76 (MQ=255)
                   gCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTAt                                                        >  1:1056998/1‑76 (MQ=255)
                     tttttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGa                                                      >  1:241773/1‑76 (MQ=255)
                       tttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGaca                                                    >  1:1225139/1‑76 (MQ=255)
                        ttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGacac                                                   >  2:949627/1‑76 (MQ=255)
                        ttGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGacac                                                   >  1:753537/1‑76 (MQ=255)
                         tGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACt                                                  >  1:686885/1‑76 (MQ=255)
                         tGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACt                                                  >  1:739251/1‑76 (MQ=255)
                            aTGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACtata                                               <  1:1681820/76‑1 (MQ=255)
                             tGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAg                                              >  1:1362839/1‑76 (MQ=255)
                              gTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAgn                                             <  1:768260/76‑2 (MQ=255)
                               tttCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGa                                            >  1:535541/1‑76 (MQ=255)
                                  cATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATaa                                         <  2:1804434/76‑1 (MQ=255)
                                  cATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATaa                                         <  2:1804432/76‑1 (MQ=255)
                                  cATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATaa                                         <  1:1498395/76‑1 (MQ=255)
                                    tAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACg                                       >  1:790156/1‑76 (MQ=255)
                                       cAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGAtt                                    <  1:2020541/76‑1 (MQ=255)
                                        aaaaTTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTg                                   <  2:1446750/76‑1 (MQ=255)
                                        aaaaTTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTg                                   >  2:762531/1‑76 (MQ=255)
                                         aaaTTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTgg                                  >  1:592572/1‑76 (MQ=255)
                                          aaTTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGa                                 >  1:1221867/1‑76 (MQ=255)
                                          aaTTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGa                                 <  1:1093860/76‑1 (MQ=255)
                                             ttCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGaaaa                              >  2:2061747/1‑76 (MQ=255)
                                             ttCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGaaaa                              >  1:1125657/1‑76 (MQ=255)
                                              tctcTTTTGGATACCATGGATTATAAATATTAAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAAc                             >  2:1405097/1‑76 (MQ=255)
                                              tctcTTTTGGATACCATGGATTATAAATATTAAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAAc                             >  2:1405078/1‑76 (MQ=255)
                                                  ttttGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACTCTATAGGATAACGATTGGAAAACCAAt                         <  2:306241/76‑1 (MQ=255)
                                                    ttGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGt                       >  1:419443/1‑76 (MQ=255)
                                                      ggATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTcc                     >  1:1636197/1‑76 (MQ=255)
                                                      ggATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTcc                     >  2:587248/1‑76 (MQ=255)
                                                      ggATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTcc                     <  1:1099182/76‑1 (MQ=255)
                                                        aTACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCAt                   >  2:1397782/1‑76 (MQ=255)
                                                         gaCCCTGGATTATAAATATTCAGTACAATACAATGTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATa                  <  2:430690/75‑1 (MQ=255)
                                                          aCCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATAc                 >  1:539414/1‑76 (MQ=255)
                                                           ccATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACt                >  1:7478/1‑76 (MQ=255)
                                                             aTGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTcc              <  2:151813/76‑1 (MQ=255)
                                                             aTGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTcc              <  1:80349/76‑1 (MQ=255)
                                                             aTGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTcc              <  1:80347/76‑1 (MQ=255)
                                                               ggATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCAt            >  2:1975590/1‑76 (MQ=255)
                                                                gATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATa           <  2:1130980/76‑1 (MQ=255)
                                                                   tataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAcc        <  1:1619683/76‑1 (MQ=255)
                                                                    ataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAccc       >  2:1799635/1‑76 (MQ=255)
                                                                    ataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAccc       >  1:1418856/1‑76 (MQ=255)
                                                                    ataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAccc       >  1:2071673/1‑76 (MQ=255)
                                                                    ataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAccc       >  1:1561872/1‑76 (MQ=255)
                                                                    ataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAccc       <  1:1734504/76‑1 (MQ=255)
                                                                    ataAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAAccc       >  2:825894/1‑76 (MQ=255)
                                                                        atatTCAGTACAATACCATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAACCCAAag   >  2:1091965/1‑76 (MQ=255)
                                                                        atatTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAACCCAAag   <  1:396687/76‑1 (MQ=255)
                                                                         tatTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAACCCAAagc  >  2:1074988/1‑76 (MQ=255)
                                                                         |                                                                           
TTACTGCTAAATGAAAATAGCTTTTTGAATGTTTCATAGCAAAATTTCTCTTTTGGATACCATGGATTATAAATATTCAGTACAATACAATTTATGACACTATAGGATAACGATTGGAAAACCAATGTCCATACTCCATAACCCAAAGC  >  REL606/2042183‑2042331

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