New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 204200334 (1.540)4 (0.190) 3/144 2.7 11.2% coding (818/1428 nt) wzx O antigen translocase (Wzx)
?REL606 = 2042213 31 (1.460)coding (608/1428 nt) wzx O antigen translocase (Wzx)
Rejected: Frequency below/above cutoff threshold.

TTGGAAAAACTACAACGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGCTAGAGATCTGAT  >  REL606/2041929‑2042078
                                                                          |                                                                           
ttGGAAAAACTACAACGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAat                                                                            <  2:444694/76‑1 (MQ=255)
        aCTACAACGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACaaa                                                                    <  1:692253/76‑1 (MQ=255)
           aCAACGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACaaaaaa                                                                 >  2:641825/1‑76 (MQ=255)
             aaCGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGt                                                               >  2:105302/1‑76 (MQ=255)
              aCGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTg                                                              <  2:216991/76‑1 (MQ=255)
                 ccTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTGCAAAAAAGTGGAt                                                           <  1:548678/76‑1 (MQ=255)
                   tGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATaa                                                         <  2:167608/76‑1 (MQ=255)
                       aTTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTa                                                     <  2:692442/76‑1 (MQ=255)
                       aTTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTa                                                     >  1:701181/1‑76 (MQ=255)
                              gtgtgtTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTg                                              >  2:395943/1‑76 (MQ=255)
                               tgtgtTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGa                                             >  1:14804/1‑76 (MQ=255)
                                 tgtTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGaat                                           >  1:161827/1‑76 (MQ=255)
                                  gtTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGaata                                          >  1:244324/1‑76 (MQ=255)
                                     tAAATAAGCATCTTCTTTTTTCCTTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAataata                                       <  1:324723/76‑1 (MQ=255)
                                     tAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAataata                                       <  1:61097/76‑1 (MQ=255)
                                          aaGCATCATCTATTTTCCCTTTGCTAT‑‑GAAAcattcaaaaagctattttcatttagcagtaacttactaattgctg                                >  2:89325‑M1/1‑31 (MQ=255)
                                          aaGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGt                                  >  1:89325/1‑76 (MQ=255)
                                             catcatCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAg                               >  1:439563/1‑76 (MQ=255)
                                               tcatcTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAgtg                             <  2:327236/76‑1 (MQ=255)
                                               tcatcTATTTTCCCNTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAgtg                             >  1:390464/1‑76 (MQ=255)
                                                 atcTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTc                           >  1:115373/1‑76 (MQ=255)
                                                 atcTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTc                           >  1:115367/1‑76 (MQ=255)
                                                     aTTTTCCCTTTGCCNTTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTg                       >  1:391182/1‑76 (MQ=255)
                                                     aTTTTCCCTTTGCCNTTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTg                       >  1:391175/1‑76 (MQ=255)
                                                     aTTTTCCCTTTGCCATTATAAAcattcaaaaagctattttcatttagcagtaacttactaattgctggagtacttg                       >  2:391182‑M1/1‑22 (MQ=255)
                                                     aTTTTCCCTTTGCCATTATAAAcattcaaaaagctattttcatttagcagtaacttactaattgctggagtacttg                       >  2:391175‑M1/1‑22 (MQ=255)
                                                      ttttCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGc                      >  1:156583/1‑76 (MQ=255)
                                                        ttCCCTTTGCCATTATAAAcattcaaaaagctattttcatttagcagtaacttactaattgctggagtacttgaag                    >  2:375743‑M1/1‑19 (MQ=255)
                                                        ttCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAg                    >  1:375743/1‑76 (MQ=255)
                                                          cccTTTGCCATTATAAATATTACAAAAAAGTGGGTAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGt                  >  1:651745/1‑76 (MQ=255)
                                                            cTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTAc                <  2:459581/76‑1 (MQ=255)
                                                             tttGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACg               <  1:335008/76‑1 (MQ=255)
                                                              ttGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGc              >  2:32698/1‑76 (MQ=255)
                                                               tGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGCt             >  1:357960/1‑76 (MQ=255)
                                                                   aTTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGCTagag         <  1:271234/76‑1 (MQ=255)
                                                                     tataAATATTACAAAAAAGCGGATAAGTTACTCTCTGAATAATATTAGCAAGTGTCATTGCAGCTACGCTAgagac       >  2:105193/1‑75 (MQ=255)
                                                                      ataAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGCTAGAGATc      <  1:94840/76‑1 (MQ=255)
                                                                          atatTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGCTAGAGATCTGAt  <  1:50260/76‑1 (MQ=255)
                                                                          acattcaaaaagctattttcatttagcagtaacttactaattgctggagtacttgaagctttttattcaaatatat  >  1:566264‑M1/1‑1 (MQ=255)
                                                                          |                                                                           
TTGGAAAAACTACAACGCCTGCAATTTTAAGTGTGTTTAAATAAGCATCATCTATTTTCCCTTTGCCATTATAAATATTACAAAAAAGTGGATAAGTTACTCTCTGAATAATATTAGTAAGTGTCATTGCAGGTACGCTAGAGATCTGAT  >  REL606/2041929‑2042078

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