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

TGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGA  >  REL606/2038649‑2038790
                                                                           |                                                                  
tGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTAtcttct                                                                    <  1:451678/76‑1 (MQ=255)
  cATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaa                                                                  <  2:1483437/76‑1 (MQ=255)
  cATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTaa                                                                  >  1:47090/1‑76 (MQ=255)
                   aaTGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATtccttcct                                                 >  1:1157168/1‑76 (MQ=255)
                    aTGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTccttcctt                                                <  1:867359/76‑1 (MQ=255)
                        aTCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCT‑‑AAAATCTATTCCTTCCTTATAGGc                                          <  1:1408292/76‑1 (MQ=255)
                            ttCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAg                                        <  1:1554656/76‑1 (MQ=255)
                              cAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGaa                                      <  1:7952/76‑1 (MQ=255)
                                   tGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc                                 >  1:1331504/1‑76 (MQ=255)
                                   tGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTc                                 >  1:1331513/1‑76 (MQ=255)
                                    gagaATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCa                                <  1:1221251/76‑1 (MQ=255)
                                            gTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAgaga                        <  1:453234/76‑1 (MQ=255)
                                               cTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAAtt                     <  2:1456809/76‑1 (MQ=255)
                                                tCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAAttt                    >  1:909041/1‑76 (MQ=255)
                                                 ccTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTg                   >  2:1329647/1‑76 (MQ=255)
                                                   tGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGtn                 <  1:222904/76‑2 (MQ=255)
                                                   tGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGtn                 <  1:222897/76‑2 (MQ=255)
                                                     aTGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTccc               >  1:528253/1‑76 (MQ=255)
                                                               tGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGata     >  1:1211929/1‑76 (MQ=255)
                                                               tGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGata     >  1:1211909/1‑76 (MQ=255)
                                                               tGAGCTATCTTCTAAAAAAACTATTCCTTCCTTATAGGCAGAAAACCCAAAAAAAAAATTTGTCCCCTTTTCCata     >  1:885465/1‑76 (MQ=255)
                                                                 aGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATg   >  1:1516641/1‑76 (MQ=255)
                                                                  gCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGa  <  1:1512454/76‑1 (MQ=255)
                                                                           |                                                                  
TGCATTATGTTTCAAGAACAATGTATCATTCAAAATGAGAATTGGTACTCCTGATGGTAATAATGAGCTATCTTCTAAAAAATCTATTCCTTCCTTATAGGCAGAAAACTCAAAAAGAGAATTTGTCCCGTTTTCGATATGA  >  REL606/2038649‑2038790

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