New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038724 =32 (0.620)12 (0.260) 6/134 4.6 22.7% coding (225/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2039214 = 53 (1.150)coding (481/747 nt) wbbA putative glycosyltransferase WbbA
Rejected: Coverage evenness skew score above cutoff.

CATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCT  >  REL606/2038651‑2038797
                                                                         |                                                                         
cATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaa                                                                         >  1:473132/1‑76 (MQ=255)
cATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaa                                                                         >  1:473105/1‑76 (MQ=255)
                        cATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGc                                                 >  2:1308283/1‑76 (MQ=255)
                           tCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGa                                              >  1:741868/1‑76 (MQ=255)
                                aTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAAcn                                         <  1:1194683/76‑2 (MQ=255)
                                aTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACt                                         <  1:911725/76‑1 (MQ=255)
                                aTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACt                                         <  1:1321854/76‑1 (MQ=255)
                                aTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACt                                         <  1:605623/76‑1 (MQ=255)
                                aTGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACt                                         <  1:460836/76‑1 (MQ=255)
                                       ttGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAag                                  >  1:575515/1‑76 (MQ=255)
                                          gTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAgaga                               >  1:695000/1‑76 (MQ=255)
                                           tACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGaa                              <  1:1245322/76‑1 (MQ=255)
                                           tACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGaa                              <  1:614370/76‑1 (MQ=255)
                                                 tGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTc                        >  1:1609344/1‑76 (MQ=255)
                                                   aTGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTccc                      >  1:130672/1‑76 (MQ=255)
                                                    tGGTAATAATGAGCTATCTTCTAAAAAATCTATTTCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCg                     >  1:864530/1‑76 (MQ=255)
                                                      gTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGtt                   >  1:1617122/1‑76 (MQ=255)
                                                      gTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGtt                   >  2:1396939/1‑76 (MQ=255)
                                                      gTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGtt                   >  1:1453723/1‑76 (MQ=255)
                                                       taataaTGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGttt                  >  1:459277/1‑76 (MQ=255)
                                                          taatGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCg               >  1:948351/1‑76 (MQ=255)
                                                          taatGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCg               >  1:875992/1‑76 (MQ=255)
                                                            atGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGat             >  1:430526/1‑76 (MQ=255)
                                                            atGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGat             >  1:696068/1‑76 (MQ=255)
                                                            atGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGat             >  1:1532187/1‑76 (MQ=255)
                                                              gAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGatat           <  1:650711/76‑1 (MQ=255)
                                                              gAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGatat           <  1:1154725/76‑1 (MQ=255)
                                                                   aTCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAAt      >  1:1445068/1‑76 (MQ=255)
                                                                     cttctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATcc    >  1:1261971/1‑76 (MQ=255)
                                                                       tctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCt  >  2:1097538/1‑76 (MQ=255)
                                                                       tctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCt  >  2:1097556/1‑76 (MQ=255)
                                                                       tctAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCt  >  2:988603/1‑76 (MQ=255)
                                                                         |                                                                         
CATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGAAATCCCT  >  REL606/2038651‑2038797

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