New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 373726148 (2.170)4 (0.190) 3/142 2.7 8.8% coding (566/1254 nt) waaL Lipid A‑core, surface polymer ligase
?REL606 = 3737695 37 (1.770)coding (132/1254 nt) waaL Lipid A‑core, surface polymer ligase
Rejected: Frequency below/above cutoff threshold.

AGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGAAAAAATACAT  >  REL606/3737625‑3737770
                                                                      |                                                                           
agaTGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAAt                                                                        >  1:547388/1‑76 (MQ=255)
 gaTGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAAtt                                                                       <  1:84046/76‑1 (MQ=255)
  aTGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTc                                                                      <  1:685506/76‑1 (MQ=255)
   tGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTcc                                                                     <  1:611622/76‑1 (MQ=255)
    gAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATataaaacaa                                                                    >  2:608613‑M2/1‑67 (MQ=255)
    gAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATataaaacaa                                                                    >  2:608607‑M2/1‑67 (MQ=255)
       aaTTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCt                                                                 <  1:154221/76‑1 (MQ=255)
          tAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCa                                                              >  1:325446/1‑76 (MQ=255)
             cAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATga                                                           >  2:53509/1‑76 (MQ=255)
              aaaGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATgaa                                                          >  1:696397/1‑76 (MQ=255)
               aaGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATgaag                                                         <  2:538830/76‑1 (MQ=255)
                aGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATataaaacaacattatatacaa                                                        >  2:365111‑M2/1‑55 (MQ=255)
                 gTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATgaagaa                                                       >  2:313176/1‑76 (MQ=255)
                  tAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGaaa                                                      >  1:490907/1‑76 (MQ=255)
                   aGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAAc                                                     >  2:375581/1‑76 (MQ=255)
                     gATTATTAGTTATTGTTTTTTTTGGTGTCCTTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACtt                                                   >  1:74135/1‑76 (MQ=255)
                      attattAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATataaaacaacattatatacaattaact                                                  >  2:142068‑M2/1‑49 (MQ=255)
                                  tGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAAcatca                                      >  2:27634/1‑76 (MQ=255)
                                      ttttttGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATccc                                  <  2:155413/76‑1 (MQ=255)
                                      ttttttGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATccc                                  <  1:9035/76‑1 (MQ=255)
                                      ttttttGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATccc                                  >  1:698461/1‑76 (MQ=255)
                                       tttttGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCa                                 >  1:347492/1‑76 (MQ=255)
                                         tttGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATa                               <  1:366454/76‑1 (MQ=255)
                                              tgtCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGAc                          <  1:346118/76‑1 (MQ=255)
                                                tCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACaa                        >  1:576032/1‑76 (MQ=255)
                                                   ttAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGt                     >  1:676816/1‑76 (MQ=255)
                                                   ttAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGt                     >  1:676812/1‑76 (MQ=255)
                                                       aaaaTACAATGCAGATATAATTCCAGCTGCTTGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGtcatc                 <  1:328518/76‑1 (MQ=255)
                                                       aaaaTACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGtcatc                 >  2:516278/1‑76 (MQ=255)
                                                        aaaTACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGtcatca                >  2:605505/1‑76 (MQ=255)
                                                        aaaTACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGtcatca                >  2:577653/1‑76 (MQ=255)
                                                         aaTACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCaa               >  2:520830/1‑76 (MQ=255)
                                                          aTACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCaaa              >  2:133107/1‑76 (MQ=255)
                                                           tACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAg             >  1:35120/1‑76 (MQ=255)
                                                            aCAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGa            <  2:708536/76‑1 (MQ=255)
                                                              aaTGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGaaa          >  1:332348/1‑76 (MQ=255)
                                                              aaTGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGaaa          >  1:332340/1‑76 (MQ=255)
                                                                  cAGATATAATTCCNGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGAAAAAAt      >  1:484503/1‑76 (MQ=255)
                                                                    gATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGAAAAAATAc    <  1:23678/76‑1 (MQ=255)
                                                                      tataATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGAAAAAATACat  >  1:333650/1‑76 (MQ=255)
                                                                      tataATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGAAAAAATACat  >  1:333641/1‑76 (MQ=255)
                                                                      |                                                                           
AGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATATAATTCCAGCTGCATGAAGAAACTTTGTTGCAACATCATCCCATAAAGACAAAGTCATCAAAGAAAAAATACAT  >  REL606/3737625‑3737770

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.