New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2038432 =23 (1.040)9 (0.460) 4/134 2.2 24.2% coding (517/822 nt) wbbB putative glycosyltransferase WbbB
?REL606 2038941 = 36 (1.820)coding (8/822 nt) wbbB putative glycosyltransferase WbbB

CTTTTGTATGCTAATAAATGAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2038508‑2038432
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tttttttaCAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGG  >  REL606/2038941‑2039012
                                                                                                                                                     
CTTTTGTATGCTGATTAATGAAAAATCTGCATAAATAATTGTGAGTTGCTACAATGAGCTTCCAAAATATTTTTTT                                                                           <  2:692765/76‑1
    TGTATGCTAATAAATGAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAA                                                                       <  2:345220/76‑1
    TGTATGCTAATAAATGAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAA                                                                       <  2:345221/76‑1
                   GAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGATCTTCCAAAATATTTTTTTACAATCCATTATAAATCAT                                                        <  2:382259/76‑1
                   GAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAATCCATTATAAATCAT                                                        <  2:507917/76‑1
                   GAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAATCCATTATAAATCAT                                                        <  2:507930/76‑1
                   GAAAAATCTGCAAAAAGAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAATCCATTATAAATCAT                                                        <  2:507923/76‑1
                         TCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAATCCATTATAAATCATCCCAAT                                                  <  2:421857/76‑1
                         TCTGCAAAAATAATGGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAATCCATTATAAATCATCCCAAT                                                  <  2:421853/76‑1
                             CAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTACAATCCATTATAAATCATCCCAATTGAT                                              <  2:213200/76‑1
                                                                         TTTACAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGG  >  1:278100/1‑76
                                                                                                                                                     
CTTTTGTATGCTAATAAATGAAAAATCTGCAAAAATAATTGTTAGTTGCTACAATGAGCTTCCAAAATATTTTTTTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2038508‑2038432
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tttttttaCAATCCATTATAAATCATCCCAATTGATTGAAGGTGATAACAAAGAAGCCATACAATGAGTTGAAAGAGCGG  >  REL606/2038941‑2039012

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