New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE = 49427655 (1.060)3 (0.070) 3/90 NT 5.4% coding (1017/1041 nt) galM galactose‑1‑epimerase
?minE = 494284 57 (1.250)coding (1009/1041 nt) galM galactose‑1‑epimerase

GTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/494205‑494276
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tattccAGCCTGACGGAATATCAGTTTATTGCTGAGTAATGTTACGCCCTCGTCATGAGGGCTTTATCTCA  <  minE/494284‑494220
                                                                                                                                         
GTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTC                                                                    >  1:2255091/1‑71
  GATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCC                                                                   <  1:324536/70‑1
            ATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATT                                                                     <  1:8924/58‑1
               TGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCAGCCTGACGGAATA                                                     <  1:557754/71‑1
                 AGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCAGCCTGACGGAATAT                                                    <  1:174508/70‑1
                    TAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCAGCCTGACGGAATATCA                                                  <  1:2151037/69‑1
                             CATGACGAGGGCGTAACATTACTCAGCAATAAACTGATAT                                                                      <  1:2173580/40‑1
                              ATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTC                                                                    >  1:2303077/1‑41
                               TGACGAGGGCGTAACATTACTCAGCAATAAACTGAT                                                                        >  1:2918758/1‑36
                                                                  TATTCCAGCCTGACGGAATATCAGTTTATTGCTGAGTAATGTTACGCCCTCGTCATGAGGGCTTTATCTCA  >  1:2570465/1‑71
                                                                      CCAGCCTGACGGAATATCAGTTTATTGCTGAGTAATGTTACGCCCTCGTCATGAGGGCTTTAT      >  1:161167/1‑63
                                                                      CCAGCCTGACGGAATATCAGTTTATTGCTGAGTAATGTTACGCCCTCGTCATGAGGGCTTTAT      >  1:1875632/1‑63
                                                                                                                                         
GTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  minE/494205‑494276
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tattccAGCCTGACGGAATATCAGTTTATTGCTGAGTAATGTTACGCCCTCGTCATGAGGGCTTTATCTCA  <  minE/494284‑494220

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