New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 270568054 (0.920)5 (0.090) 3/94 NT 9.1% coding (498/837 nt) dam DNA adenine methylase
?minE = 2705688 50 (0.920)coding (506/837 nt) dam DNA adenine methylase

ATTTCCCGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/2705607‑2705680
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCT  <  minE/2705688‑2705620
                                                                                                                                               
ATTTCCCGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTAC                                                                          >  1:457061/1‑71
   TCCCGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGC                                                                        >  1:585921/1‑70
                    TATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCATGCTATCGGCGTAAG                                                       <  1:2263610/70‑1
                                GAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTAC                                                                          <  1:699181/39‑1
                                    AAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCATGCTATCGGCGTAAGA                                                      >  1:3285947/1‑55
                                    AAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCATGCTATCGGCGTAAGA                                                      >  1:3713199/1‑55
                                                  TTTTTCTATTGTGAGTCTTACGCCATGCTATCGGCG                                                           <  1:1616332/36‑1
                                                  TTTTTCTATTGTGAGTCTTACGCCATGCTATCGGCG                                                           <  1:3458714/36‑1
                                                                      CGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACT    <  1:2163709/71‑1
                                                                      CGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACT    <  1:308906/71‑1
                                                                       GCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTC   <  1:1531677/71‑1
                                                                       GCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTC   <  1:2328261/71‑1
                                                                       GCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTC   <  1:3192253/71‑1
                                                                        CCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCT  >  1:3447795/1‑71
                                                                                                                                               
ATTTCCCGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/2705607‑2705680
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCT  <  minE/2705688‑2705620

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