New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2041682 =41 (0.630)3 (0.050) 3/128 7.8 6.3% coding (1139/1428 nt) wzx O antigen translocase (Wzx)
?REL606 2041834 = 55 (1.000)coding (987/1428 nt) wzx O antigen translocase (Wzx)
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCGTTGATACTTTCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCAAACG  >  REL606/2041608‑2041755
                                                                          |                                                                         
tCGTTGATACTTTCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCAt                                                                          <  2:659551/76‑1 (MQ=255)
 cGTTGATACTTTCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCAtt                                                                         >  1:1230890/1‑76 (MQ=255)
  gTTGATACTTTCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACAAATAACAATTg                                                                        >  1:1809266/1‑76 (MQ=255)
   ttGATACTTTCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGc                                                                       <  1:954915/76‑1 (MQ=255)
    tGATACTTTCTGGCTAAATAATCCATTTAAAAACCATGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCa                                                                      >  1:2010086/1‑76 (MQ=255)
          tttCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCatttata                                                                >  1:1891574/1‑76 (MQ=255)
          tttCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCatttata                                                                >  1:1891563/1‑76 (MQ=255)
                    aataatCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCa                                                      >  1:1474110/1‑76 (MQ=255)
                     ataatCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCAt                                                     >  1:952415/1‑76 (MQ=255)
                       aatCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCg                                                   <  1:166738/76‑1 (MQ=255)
                        atCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGt                                                  >  1:1379980/1‑76 (MQ=255)
                            aTTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATg                                              >  1:1793065/1‑76 (MQ=255)
                               tAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACa                                           <  2:989674/76‑1 (MQ=255)
                                 aaaaCCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAg                                         >  2:1679944/1‑76 (MQ=255)
                                    aCCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGaat                                      <  1:748874/76‑1 (MQ=255)
                                    aCCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGaat                                      <  1:748899/76‑1 (MQ=255)
                                     ccACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGaata                                     >  2:1432164/1‑76 (MQ=255)
                                        cGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAAtt                                  <  2:1933037/76‑1 (MQ=255)
                                         gAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTg                                 >  1:1299294/1‑76 (MQ=255)
                                              atatGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGtcatc                            >  2:1937377/1‑76 (MQ=255)
                                               tatGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGtcatca                           <  2:498178/76‑1 (MQ=255)
                                               tatGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGtcatca                           <  2:498181/76‑1 (MQ=255)
                                               tatGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGtcatca                           <  1:2000940/76‑1 (MQ=255)
                                               tatGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGtcatca                           >  2:896094/1‑76 (MQ=255)
                                                atGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCaa                          >  2:416657/1‑76 (MQ=255)
                                                atGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCaa                          >  1:1909864/1‑76 (MQ=255)
                                                   aGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGt                       <  2:1737793/76‑1 (MQ=255)
                                                       tAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCttt                   <  2:1745283/76‑1 (MQ=255)
                                                          aaCTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCtttttt                >  1:1048080/1‑76 (MQ=255)
                                                          aaCTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCtttttt                >  1:284558/1‑76 (MQ=255)
                                                            cTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAg              <  1:843561/76‑1 (MQ=255)
                                                                 ccAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTTATCAATGTCTTTTTTAGTAttt         >  1:1907495/1‑76 (MQ=255)
                                                                 ccAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTAttt         >  1:1698213/1‑76 (MQ=255)
                                                                 ccAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTAttt         >  1:1907518/1‑76 (MQ=255)
                                                                   aaTCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTcc       >  1:905234/1‑76 (MQ=255)
                                                                   aaTCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTcc       >  1:502144/1‑76 (MQ=255)
                                                                      cACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCaaa    >  1:1357351/1‑76 (MQ=255)
                                                                      cACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCaaa    >  2:437581/1‑76 (MQ=255)
                                                                       aCCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCAAAc   <  1:1532084/76‑1 (MQ=255)
                                                                       aCCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCAAAc   <  1:1532094/76‑1 (MQ=255)
                                                                        ccATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCAAACg  >  1:1362651/1‑76 (MQ=255)
                                                                          |                                                                         
TCGTTGATACTTTCTGGCTAAATAATCCATTTAAAAACCACGAAAGATATGAGTGTAAAACTAAACCAATCACCATTGCATTTATATTTATTTGCATCGTAATGACAAGAATAATTGTCATCAATGTCTTTTTTAGTATTTCCAAACG  >  REL606/2041608‑2041755

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