New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 4551446NA (NA)9 (0.090) 8/46 NT 7.0% noncoding (768/768 nt) IS1 repeat region
?REL606 4560889 = 119 (1.210)coding (135/1590 nt) hsdM DNA methylase M

ATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/4551414‑4551446
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cCTGTCCGGTTTCTTTACACATTTTCAAAAACAGC  >  REL606/4560889‑4560922
                                                                   
ATAAAACACTATCAATAAGTTGGAGTCATTACCCT                                  <  1:6231550/35‑1
ATAAAACACTATCAATAAGTTGGAGTCATTACCCT                                  <  1:19111899/35‑1
ATAAAACACTATCAATAAGTTGGAGTCATTACCCT                                  >  1:12358203/1‑35
ATAAAACACTATCAATAAGTTGGAGTCATTACCCT                                  >  1:10453141/1‑35
  AAAACACTATCAATAAGTTGGAGTCATTACCCTGT                                >  1:15714295/1‑35
    AACACTATCAATAAGTTGGAGTCATTACCCTGTCC                              <  1:13608055/35‑1
       ACTATCAATAAGTTGGAGTCATTACCCTGTCCGGT                           <  1:13688259/35‑1
       ACTATCAATAAGTTGGAGTCATTACCCTGTCCGGT                           <  1:7566326/35‑1
            CAATAAGTTGGAGTCATTACCCTGTCCGGTTTCTT                      <  1:20232665/35‑1
                 AGTTGGAGTCATTACCCTGTCCGGTTTCTTTACAC                 <  1:12150954/35‑1
                  GTTGGAGTCATTACCCTGTCCGGTTTCTTTACACA                <  1:2175040/35‑1
                    TGGAGTCATTACCCTGTCCGGTTTCTTTACACATT              <  1:6375011/35‑1
                     GGAGTCATTACCCTGTCCGGTTTCTTTACACATTT             <  1:17649380/35‑1
                     GGAGTCATTACCCTGTCCGGTTTCTTTACACATTT             >  1:2276032/1‑35
                                CCTGTCCGGTTTCTTTACACATTTTCAAAAACAGC  <  1:1471295/35‑1
                                CCTGTCCGGTTTCTTTACACATTTTCAAAAACAGC  >  1:18329536/1‑35
                                CCTGTCCGGTTTCTTTACACATTTTCAAAAACAGC  >  1:6861336/1‑35
                                CCTGTCCGGTTTCTTTACACATTTTCAAAAACAGC  <  1:8393940/35‑1
                                                                   
ATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/4551414‑4551446
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cCTGTCCGGTTTCTTTACACATTTTCAAAAACAGC  >  REL606/4560889‑4560922

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