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

GCCCGGGAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAAG  >  minE/2409448‑2409583
                                                                    |                                                                   
gCCCGGGAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCAt                                                                    <  1:503982/70‑1 (MQ=255)
 cccGGGAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTa                                                                  >  1:412278/1‑71 (MQ=255)
   cGGGAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTAc                                                                 >  1:1390908/1‑70 (MQ=255)
   cGGGAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACa                                                                <  1:793139/71‑1 (MQ=255)
     ggAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACATCGATGCGCTTTTACGGGCATTACAAt                                                              <  1:1255615/71‑1 (MQ=255)
         cTGGGCTGGCAAGACATTAAGGTCGCCTATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAAct                                                          >  1:3784163/1‑71 (MQ=255)
         cTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAAct                                                          <  1:204909/71‑1 (MQ=255)
          tGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAAct                                                          >  1:3660032/1‑70 (MQ=255)
            ggCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAActct                                                        >  1:3490992/1‑70 (MQ=255)
            ggCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACtctc                                                       <  1:3585698/71‑1 (MQ=255)
                 gCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATaa                                                   >  1:2069143/1‑70 (MQ=255)
                 gCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAAc                                                  <  1:105873/71‑1 (MQ=255)
                    aGACATTAAGGTCGCCGATAACGCTGACAACGATGCTCTTTTACGGTCATTACAATAACTCTCATAACAg                                                >  1:308206/1‑70 (MQ=255)
                     gACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAgg                                               >  1:2405010/1‑70 (MQ=255)
                     gACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAgg                                               >  1:910550/1‑70 (MQ=255)
                     gACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGa                                              <  1:3048215/71‑1 (MQ=255)
                     gACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGa                                              <  1:2770174/71‑1 (MQ=255)
                         ttAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACa                                                                <  1:3252255/49‑1 (MQ=255)
                         ttAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGc                                           >  1:3845328/1‑70 (MQ=255)
                          tAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATaa                                                   >  1:2886618/1‑61 (MQ=255)
                          tAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATaa                                                   <  1:690350/61‑1 (MQ=255)
                          tAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGcc                                          >  1:2425113/1‑70 (MQ=255)
                           aaGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAAt                                                              >  1:1846145/1‑49 (MQ=255)
                            aGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCAt                                        <  1:2711418/70‑1 (MQ=255)
                            aGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATa                                       <  1:1286980/71‑1 (MQ=255)
                             ggTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACa                                                 <  1:850051/60‑1 (MQ=255)
                              gTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGa                                              <  1:3280231/62‑1 (MQ=255)
                               tCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATg                                    <  1:1002225/71‑1 (MQ=255)
                                cGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATa                                                             >  1:2312090/1‑45 (MQ=255)
                                 gcccgtaaaagcgcatcgttgtcagcgttatcggcATTACAATAACTCTCATAACAGGAAGCCATa                                       <  1:2402334‑M2/31‑1 (MQ=255)
                                 gcccgtaaaagcgcatcgttgtcagcgttatcggcATTACAATAACTCTCATAACAGGAAGCCATa                                       <  1:2039692‑M2/31‑1 (MQ=255)
                                 GCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGa                                   >  1:560452‑M1/1‑70 (MQ=255)
                                 GCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGa                                   <  1:1399125‑M1/70‑1 (MQ=255)
                                   ccgtaaaagcgcatcgttgtcagcgttatcggcATTACAATAACTCTCATAACAGGAAGCCATaa                                      >  1:1376349‑M2/34‑65 (MQ=255)
                                   CGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATaa                                                   <  1:1186978‑M1/52‑1 (MQ=255)
                                     aTAACGCTGACAACGATGCGCTTTTACGGGCATTAc                                                                 <  1:2581076/36‑1 (MQ=255)
                                      tAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGaa                              >  1:3899336/1‑70 (MQ=255)
                                       aaCGCTGACAACGATGCGCTTTTACGGGCATTACaa                                                               >  1:2494102/1‑36 (MQ=255)
                                       aaCGCTGACAACGATGCGCTTTTACGGGCATTACaa                                                               >  1:1717705/1‑36 (MQ=255)
                                         cGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAg                          <  1:2600649/71‑1 (MQ=255)
                                          gCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAg                                                <  1:612493/48‑1 (MQ=255)
                                          gCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACaa                           >  1:269064/1‑69 (MQ=255)
                                          gCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAg                          >  1:3245737/1‑70 (MQ=255)
                                          gCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGa                         <  1:2003155/71‑1 (MQ=255)
                                           cTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGaa                              >  1:2216079/1‑65 (MQ=255)
                                            tGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATaa                                      >  1:977246/1‑56 (MQ=255)
                                              aCAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACaa                           <  1:995732/65‑1 (MQ=255)
                                               cAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGaaaaaa                    <  1:389425/71‑1 (MQ=255)
                                                tcgttgtcagcgttatcggcATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAAc                   <  1:1853263‑M2/51‑1 (MQ=255)
                                                aaCGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGaa                              >  1:769483/1‑60 (MQ=255)
                                                  cGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCt                 >  1:1145517/1‑71 (MQ=255)
                                                   gATGCGCTTTTACGGGCATTACAATAACTCTCATaa                                                   >  1:2753748/1‑36 (MQ=255)
                                                   gATGCGCTTTTACGGGCATTACAATAACTCTCATAACa                                                 >  1:514314/1‑38 (MQ=255)
                                                    aTGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTcc               >  1:1952951/1‑71 (MQ=255)
                                                       cgcTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTccgcc            <  1:2201595/71‑1 (MQ=255)
                                                          ttttACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTg         >  1:267134/1‑71 (MQ=255)
                                                            ttACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAAcc                  >  1:2663694/1‑60 (MQ=255)
                                                              aCGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGaaa                       <  1:3287759/53‑1 (MQ=255)
                                                               cGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTg     <  1:67487/70‑1 (MQ=255)
                                                                gggCATTACAATAACTCTCATAACAGGAAGCCATAAt                                     >  1:1705291/1‑37 (MQ=255)
                                                                 ggcATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAag  <  1:2386789‑M2/68‑1 (MQ=255)
                                                                 ggcATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAag  <  1:2777890‑M2/68‑1 (MQ=255)
                                                                   cATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGaa   <  1:195702‑M2/67‑1 (MQ=255)
                                                                    |                                                                   
GCCCGGGAACTGGGCTGGCAAGACATTAAGGTCGCCGATAACGCTGACAACGATGCGCTTTTACGGGCATTACAATAACTCTCATAACAGGAAGCCATAATGACGGAACAAGAAAAAACCTCCGCCGTGGTTGAAG  >  minE/2409448‑2409583

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.