New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3600745 =41 (0.400)9 (0.120) 3/34 NT 18.7% intergenic (+12/+39) treF/yhjB cytoplasmic trehalase/predicted DNA‑binding response regulator in two‑component regulatory system
?REL606 3792210 = 49 (0.670)coding (316/414 nt) ECB_03532 hypothetical protein

GAAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCGGCTGG‑‑‑‑‑‑‑‑‑‑  <  REL606/3600785‑3600745
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggctggGTGGGGGAAA  >  REL606/3792210‑3792219
                                                   
GAAGAGGATAATTTTACAGCCAGCAGGAAGTTAGC                  >  1:6291826/1‑35
GAAGAGGATAATTTTACAGCCAGCAGGAAGTTAGC                  >  1:5163961/1‑35
GAAGAGGATAATTTTACAGCCAGCAGGAAGTTAGC                  >  1:15175551/1‑35
 AAGAGGATCATTTTACAGCCAGCAGGAAGTTAGCG                 >  1:17022407/1‑35
 AAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCG                 >  1:7017935/1‑35
 AAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCG                 >  1:19078316/1‑35
 AAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCG                 >  1:15206711/1‑35
 AAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCG                 >  1:15656196/1‑35
 AAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCG                 >  1:17042179/1‑35
  AGAGGATAATTTTACAGCCAGCATGAAGTTAGCGG                >  1:20517360/1‑35
  AGAGGATAATTTTACAGCCAGCAGGAAGTTAGCGG                >  1:4395066/1‑35
  AGAGGATAATTTTACAGCCAGCAGGAAGTTAGCGG                >  1:4117178/1‑35
    AGGATAATTTTACAGCCAGCAGGAAGTTAGCGGCT              >  1:18298140/1‑35
     GGATAATTTTACAGCCAGCAGGAAGTTAGCGGCTG             >  1:11300476/1‑35
     GGATAATTTTACAGCCAGCAGGAAGTTAGCGGCTG             >  1:11438645/1‑35
      GATAATTTTACAGCCAGCAGGAAGTTAGCGGCTGG            >  1:3429137/1‑35
       ATAATTTTACAGCCAGCAGGAAGTTAGCGGCTGGG           >  1:11960457/1‑35
       ATAATTTTACAGCCAGCAGGAAGTTAGCGGCTGGG           >  1:14859057/1‑35
         AATTTTACAGCCAGCAGGAAGTTAGCGGCTGGGgg         >  1:13270661/1‑33
          ATTTTACAGCCAGCAGGAAGTTAGCGGCTGGGTGG        >  1:3805542/1‑35
           TTTTACAGCCAGCAGGAAGTTAGCGGCTGGGTGGG       >  1:17291307/1‑35
           TTTTACAGCCAGCAGGAAGTTAGCGGCTGGGTGGG       >  1:17511086/1‑35
           TTTTACAGCCAGCAGGAAGTTAGCGGCTGGGTGGG       >  1:1828605/1‑35
           TTTTACAGCCAGCAGGAAGTTAGCGGCTGGGGGGG       >  1:3284442/1‑35
            TTTACAGCCAGCAGGAAGTTAGCGGCTGGGTGGGG      >  1:15179166/1‑35
            TTTACAGCCAGCAGGAAGTTAGCGGCTGGCTGGGG      >  1:16482691/1‑35
               ACAGCCAGCAGGAAGTTAGCGGCTGGGTGGGGGAA   >  1:18986811/1‑35
               ACAGCCAGCAGGAAGTTAGCGGCTGGGTGGGGGAA   >  1:7128610/1‑35
               ACAGCCAGCAGGAAGTTAGCGGCTGGGTGGGGAAA   >  1:4584055/1‑35
               ACAGCCAGCAGGAAGTTAGCGGCTGGGTGGGGAAA   >  1:11500917/1‑35
                CAGCCAGCAGGAAGTTAGCGGCTGGGTGGGGGAAA  >  1:18383786/1‑35
                CAGCCAGCAGGAAGTTAGCGGCTGGGTGGGGAAAA  >  1:251892/1‑35
                CAGCCAGCAGGAAGTTAGCGGCTGGCTGGGGGAAA  >  1:19324930/1‑35
                                                   
GAAGAGGATAATTTTACAGCCAGCAGGAAGTTAGCGGCTGG‑‑‑‑‑‑‑‑‑‑  <  REL606/3600785‑3600745
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggctggGTGGGGGAAA  >  REL606/3792210‑3792219

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