New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2041421 =70 (1.380)4 (0.090) 4/132 5.4 6.4% coding (1400/1428 nt) wzx O antigen translocase (Wzx)
?REL606 2041804 = 56 (1.260)coding (1017/1428 nt) wzx O antigen translocase (Wzx)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCCACTCACTTCCT  >  REL606/2041730‑2041879
                                                                          |                                                                           
naTGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTg                                                                            >  1:104245/2‑76 (MQ=255)
aaTGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTg                                                                            <  1:1371802/76‑1 (MQ=255)
 aTGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGc                                                                           >  2:152207/1‑76 (MQ=255)
 aTGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGc                                                                           <  1:1290292/76‑1 (MQ=255)
  tGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCa                                                                          >  1:1087050/1‑76 (MQ=255)
   gTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCAt                                                                         <  2:441406/76‑1 (MQ=255)
     cTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGt                                                                       <  1:941279/76‑1 (MQ=255)
       tttttAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGtat                                                                     <  1:845310/76‑1 (MQ=255)
       tttttAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGtat                                                                     >  2:383776/1‑76 (MQ=255)
          ttAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAgg                                                                  >  2:1151691/1‑76 (MQ=255)
             gTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGAtn                                                               <  1:553009/76‑2 (MQ=255)
             gTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATa                                                               >  2:605617/1‑76 (MQ=255)
              tATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATaa                                                              >  1:1463555/1‑76 (MQ=255)
               aTTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATaaa                                                             <  1:78161/76‑1 (MQ=255)
                tttCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAAt                                                            <  1:1378551/76‑1 (MQ=255)
                ttgccAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAAt                                                            <  1:1378530/73‑1 (MQ=255)
                   ccAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCAt                                                         <  1:1349070/76‑1 (MQ=255)
                   ccAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGNATAGGATAAATCAt                                                         <  1:891348/76‑1 (MQ=255)
                    cAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATa                                                        <  1:1399347/76‑1 (MQ=255)
                    cAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATa                                                        >  1:19589/1‑76 (MQ=255)
                    cAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATa                                                        <  2:147585/76‑1 (MQ=255)
                     aaaCGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATaa                                                       <  1:1051654/76‑1 (MQ=255)
                      aaCGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATCTTAATATTTATTGCATGTATAGGATAAATCATaaa                                                      <  2:1488495/76‑1 (MQ=255)
                      aaCGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAGTCATaaa                                                      <  1:1130964/76‑1 (MQ=255)
                              aaaaGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATAc                                              <  2:281934/76‑1 (MQ=255)
                               aaaGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACa                                             >  1:104487/1‑76 (MQ=255)
                               aaaGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACa                                             <  1:969436/76‑1 (MQ=255)
                                aaGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACaa                                            <  1:1328930/76‑1 (MQ=255)
                                 aGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACaaa                                           <  1:330639/76‑1 (MQ=255)
                                   aTCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGt                                         <  1:663096/76‑1 (MQ=255)
                                   aTCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGt                                         >  1:490502/1‑76 (MQ=255)
                                    tCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTa                                        <  1:1340747/76‑1 (MQ=255)
                                    tCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTa                                        >  2:61782/1‑76 (MQ=255)
                                    tCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTa                                        <  1:513595/76‑1 (MQ=255)
                                     cACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAg                                       <  1:1461023/76‑1 (MQ=255)
                                     cACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAg                                       >  2:1042986/1‑76 (MQ=255)
                                       cTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGaa                                     >  1:69976/1‑76 (MQ=255)
                                        ttCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCTATACAAAGTAGaaa                                    <  1:873373/76‑1 (MQ=255)
                                         tCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAAc                                   <  1:670642/76‑1 (MQ=255)
                                            cccATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGa                                <  1:711445/76‑1 (MQ=255)
                                            cccATGAACCTGTAGGATATTAAAATTTATTGCATGTATAGGATAAACCAAAAAAGCAAAACAAACAACAAACAGa                                >  1:107797/1‑76 (MQ=255)
                                             ccATGAACCTGTAGTATATTAATATTTATTGCATGTATAGTCTAAATCATAAATGCCATACACAGTATAAACAGac                               >  2:84540/1‑75 (MQ=255)
                                             ccACGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGaa                               >  1:80672/1‑76 (MQ=255)
                                              cATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAat                              >  1:250536/1‑76 (MQ=255)
                                                tGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAatat                            <  2:91634/76‑1 (MQ=255)
                                                tGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAatat                            >  2:469063/1‑76 (MQ=255)
                                                tGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAatat                            <  1:435715/76‑1 (MQ=255)
                                                  aaCCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCa                          >  1:1373870/1‑76 (MQ=255)
                                                  aaCCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCa                          >  1:1373872/1‑76 (MQ=255)
                                                   aCCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAg                         <  2:1118645/76‑1 (MQ=255)
                                                      tGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATg                      >  1:1268546/1‑76 (MQ=255)
                                                         aGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTaa                   >  1:734226/1‑76 (MQ=255)
                                                          ggATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTaaa                  <  2:1564680/76‑1 (MQ=255)
                                                          ggATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTaaa                  <  2:289838/76‑1 (MQ=255)
                                                           gATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAAt                 <  1:9015/76‑1 (MQ=255)
                                                               ttAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTcc             <  2:926656/76‑1 (MQ=255)
                                                                  atatTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCcact          >  2:1240229/1‑76 (MQ=255)
                                                                  atatTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCcact          >  2:1240219/1‑76 (MQ=255)
                                                                  atatTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCcact          >  1:952905/1‑76 (MQ=255)
                                                                    atttattGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCcactca        <  1:1617654/76‑1 (MQ=255)
                                                                    atttattGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCcactca        <  1:770840/76‑1 (MQ=255)
                                                                      ttattGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCcactcata      <  2:986717/76‑3 (MQ=255)
                                                                       tagtGCATGTATAGGATGAAGCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCCACTCACtt     <  1:1492079/73‑1 (MQ=255)
                                                                          tGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCCACTCACTTCCt  <  1:683126/76‑1 (MQ=255)
                                                                          |                                                                           
AATGTCTTTTTTAGTATTTCCAAACGCAAAAAAAGATCACTTCGCCCATGAACCTGTAGGATATTAATATTTATTGCATGTATAGGATAAATCATAAATGCAATACAAAGTAGAAACAGAATATCAGATGTAAATTTCCACTCACTTCCT  >  REL606/2041730‑2041879

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