New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038724 =51 (0.820)12 (0.220) 6/134 5.8 20.3% coding (225/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2039214 = 49 (0.890)coding (481/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

TGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTT  >  REL606/2038649‑2038798
                                                                           |                                                                          
tGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTAtcttct                                                                            <  1:587003/76‑1 (MQ=255)
 gCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTa                                                                           <  2:524855/76‑1 (MQ=255)
   aTTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaa                                                                         <  1:367347/76‑1 (MQ=255)
     tATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaaaa                                                                       <  1:8821/76‑1 (MQ=255)
     tATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaaaaa                                                                       <  1:1358470/76‑1 (MQ=255)
       tGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAAt                                                                     <  1:1243365/76‑1 (MQ=255)
       tGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAAt                                                                     <  2:835837/76‑1 (MQ=255)
         tttCATGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCt                                                                   >  2:944176/1‑76 (MQ=255)
         tttCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCt                                                                   >  1:163755/1‑76 (MQ=255)
         tctCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCt                                                                   >  1:163747/3‑76 (MQ=255)
           tCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTAt                                                                 >  2:1857549/1‑76 (MQ=255)
           tCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTAt                                                                 >  1:611107/1‑76 (MQ=255)
            cAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTAtt                                                                >  2:291356/1‑76 (MQ=255)
             aaGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTAttc                                                               >  2:488413/1‑76 (MQ=255)
                  cAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTAttccttcc                                                          >  1:1539715/1‑76 (MQ=255)
                  cAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTAttccttcc                                                          >  1:929569/1‑76 (MQ=255)
                   aaTGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATtccttcct                                                         >  2:1210961/1‑76 (MQ=255)
                      gTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTtat                                                      <  1:1291449/76‑1 (MQ=255)
                      gTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTtat                                                      <  2:793927/76‑1 (MQ=255)
                      gTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTtat                                                      >  1:953106/1‑76 (MQ=255)
                      gTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTtat                                                      <  2:1153482/76‑1 (MQ=255)
                           aTTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCa                                                 <  2:544473/76‑1 (MQ=255)
                           aTTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCa                                                 <  2:544461/76‑1 (MQ=255)
                            ttCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAg                                                <  1:614048/76‑1 (MQ=255)
                               aaaaTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGaan                                             <  1:488413/76‑2 (MQ=255)
                                aaaTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGaaaa                                            <  1:1629485/76‑1 (MQ=255)
                                aaaTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGaaaa                                            >  2:833318/1‑76 (MQ=255)
                                 aaTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAAc                                           >  2:307972/1‑76 (MQ=255)
                                   ttagaATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc                                         >  2:490255/3‑76 (MQ=255)
                                        aTTGATACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCaaaaa                                    >  1:1739927/1‑76 (MQ=255)
                                          tGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAaga                                  <  1:1042217/76‑1 (MQ=255)
                                          tGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAaga                                  <  1:1042205/76‑1 (MQ=255)
                                            gTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAgaga                                <  1:87428/76‑1 (MQ=255)
                                                     aTGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTccc                       >  2:1750118/1‑76 (MQ=255)
                                                      tagTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCg                      >  1:1801528/3‑76 (MQ=255)
                                                      tGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCg                      >  1:1560467/1‑76 (MQ=255)
                                                       ggTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGt                     <  1:733530/76‑1 (MQ=255)
                                                       ggTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAACAGAATTTGTCCCGt                     <  1:733513/76‑1 (MQ=255)
                                                         taataaTGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGttt                   <  1:358992/76‑1 (MQ=255)
                                                         taataaTGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGttt                   >  1:308756/1‑76 (MQ=255)
                                                           ataatGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTc                 >  2:1545655/1‑76 (MQ=255)
                                                             aatGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGa               >  1:302258/1‑76 (MQ=255)
                                                               tGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGata             >  1:534728/1‑76 (MQ=255)
                                                                gAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGatat            <  1:1888348/76‑1 (MQ=255)
                                                                 aGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGatata           <  2:656561/76‑2 (MQ=255)
                                                                  gCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGa          >  2:1961048/1‑76 (MQ=255)
                                                                  gCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGa          >  1:200845/1‑76 (MQ=255)
                                                                   cTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGaa         >  2:1180214/1‑76 (MQ=255)
                                                                   cTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGaa         <  1:315977/76‑1 (MQ=255)
                                                                    tATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGaaa        >  2:506282/1‑76 (MQ=255)
                                                                          ctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCtt  <  2:946160/76‑1 (MQ=255)
                                                                           |                                                                          
TGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCTT  >  REL606/2038649‑2038798

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