New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3737622 =48 (0.980)3 (0.060) 3/144 5.7 6.1% coding (205/1254 nt) waaL Lipid A‑core, surface polymer ligase
?REL606 3737741 = 47 (1.000)coding (86/1254 nt) waaL Lipid A‑core, surface polymer ligase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATGCATTGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATA  >  REL606/3737547‑3737696
                                                                           |                                                                          
atGCATTGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATaa                                                                            <  1:1378271/76‑1 (MQ=255)
  gCATTGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATg                                                                          <  2:885156/76‑1 (MQ=255)
     ttGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATgaga                                                                       >  2:127081/1‑76 (MQ=255)
     ttGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATgaga                                                                       >  1:1298615/1‑76 (MQ=255)
      tGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGAt                                                                      >  2:1049818/1‑76 (MQ=255)
      tGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGAt                                                                      >  1:1278960/1‑76 (MQ=255)
       ggTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATg                                                                     <  1:1031142/76‑1 (MQ=255)
       ggTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATg                                                                     <  1:1031159/76‑1 (MQ=255)
       ggTATAAACAGAGCCCGAAACTTTAGAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATg                                                                     <  1:1031125/76‑1 (MQ=255)
          ataAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGaaa                                                                  >  2:1171475/1‑76 (MQ=255)
              aCAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTa                                                              <  1:544952/76‑1 (MQ=255)
              aCAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTa                                                              >  2:1294022/1‑76 (MQ=255)
              aCAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTa                                                              <  2:1076186/76‑1 (MQ=255)
                 gaacccGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACa                                                           >  2:1015741/4‑76 (MQ=255)
                   gCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACaaa                                                         <  2:1527974/76‑1 (MQ=255)
                     ccGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGt                                                       <  1:1333698/76‑1 (MQ=255)
                     ccGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGt                                                       <  2:468917/76‑1 (MQ=255)
                      cGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTa                                                      >  2:323396/1‑76 (MQ=255)
                        aaaCTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAgg                                                    >  2:64685/1‑76 (MQ=255)
                        aaaCTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAgg                                                    >  2:108415/1‑76 (MQ=255)
                        aaaCTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAgg                                                    >  2:108412/1‑76 (MQ=255)
                             ttATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGattat                                               >  1:1490345/1‑76 (MQ=255)
                             ttATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGattat                                               >  2:352277/1‑76 (MQ=255)
                             ttATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGattat                                               >  1:1490359/1‑76 (MQ=255)
                             ttATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGattat                                               <  2:1261775/76‑1 (MQ=255)
                             ttATAATGCGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGattat                                               >  1:254265/1‑76 (MQ=255)
                               ataATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGAttatta                                             <  2:1336435/76‑1 (MQ=255)
                                  aTGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGtt                                          >  1:202174/1‑76 (MQ=255)
                                   tgtgAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTa                                         >  2:781342/1‑76 (MQ=255)
                                    gtgAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGtagg                                        <  2:1115956/76‑4 (MQ=255)
                                    gtgAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTAt                                        >  1:1371010/1‑76 (MQ=255)
                                     tgAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTAtt                                       >  2:1316444/1‑76 (MQ=255)
                                         tACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGttt                                   >  2:362109/1‑76 (MQ=255)
                                            cAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGtttttt                                <  1:809603/76‑1 (MQ=255)
                                            cAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGtttttt                                <  1:809591/76‑1 (MQ=255)
                                              gatgatATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGtttttttt                              <  2:1386655/76‑1 (MQ=255)
                                                  atatTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGgtg                          >  1:1139847/1‑76 (MQ=255)
                                                      ttATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCAt                      <  2:1371009/76‑1 (MQ=255)
                                                           aTTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGaa                 >  1:231687/1‑76 (MQ=255)
                                                           aTTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGaa                 >  2:1298831/1‑76 (MQ=255)
                                                              cccAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAAt              <  2:1427279/76‑1 (MQ=255)
                                                              cccAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAAt              <  2:181057/76‑1 (MQ=255)
                                                                cAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATAc            <  1:587556/76‑1 (MQ=255)
                                                                  aaaGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACaa          <  1:389476/76‑1 (MQ=255)
                                                                  aaaGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACaa          >  2:1369197/1‑76 (MQ=255)
                                                                        aTAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGa    >  1:82900/1‑76 (MQ=255)
                                                                        aTAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGa    >  1:193381/1‑76 (MQ=255)
                                                                          aaTGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGata  <  2:1461085/76‑1 (MQ=255)
                                                                           |                                                                          
ATGCATTGGTATAAACAGAGCCCGAAACTTTATAATGTGAATACCAGATGATATTTATGATTCCCAAAAGACATAATGAGATGAAAATTAACAAAGTAGGATTATTAGTTATTGTTTTTTTTGGTGTCATTAGAAAATACAATGCAGATA  >  REL606/3737547‑3737696

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