New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? minE 2409481 =65 (1.250)4 (0.080) 3/96 NT 7.0% coding (697/741 nt) hemD uroporphyrinogen III synthase
?minE 2409516 = 45 (0.920)coding (732/741 nt) hemD uroporphyrinogen III synthase

TCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/2409550‑2409481
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAAG  >  minE/2409516‑2409583
                                                                                                                                          
TCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGC                                                                      <  1:1399125/70‑1
TCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGC                                                                      >  1:560452/1‑70
                TTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCG                                                                        <  1:1186978/52‑1
                                   GCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCATTACAATAACTCTCATAACAGGAAGCCATA                                       <  1:2039692/66‑1
                                   GCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCATTACAATAACTCTCATAACAGGAAGCCATA                                       <  1:2402334/66‑1
                                     CCGTAAAAGCGCATCGTTGTCAGCGTTATCGGCATTACAATAACTCTCATAACAGGAAGCCATAA                                      >  1:1376349/1‑65
                                                  TCGTTGTCAGCGTTATCGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAAC                   <  1:1853263/71‑1
                                                                   GGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAAG  <  1:2386789/71‑1
                                                                   GGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAAG  <  1:2777890/71‑1
                                                                     CATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAA   <  1:195702/68‑1
                                                                                                                                          
TCATTATGGCTTCCTGTTATGAGAGTTATTGTAATGCCCGTAAAAGCGCATCGTTGTCAGCGTTATCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  minE/2409550‑2409481
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAAG  >  minE/2409516‑2409583

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