New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 514780 =90 (0.880)3 (0.050) 3/28 NT 8.3% coding (727/1146 nt) glxK glycerate kinase II
?REL606 1871739 = 14 (0.230)coding (75/1911 nt) yoaA conserved protein with nucleoside triphosphate hydrolase domain

TCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCG  >  REL606/1871707‑1871773
                                |                                  
tCTATCGCCTGGGTGACGGCTACCGCCATCTGTCg                                  >  1:3716286/1‑35 (MQ=255)
   aTCGCCTGGGTGACGGCTAACGCCATCTGTTGccg                               >  1:11444830/1‑33 (MQ=25)
      gCCTGGGTGACGGCTACCGCCATCTGTCGCCGTgg                            >  1:16780060/1‑35 (MQ=255)
       ccTGGGTGACGGCTACCGCCATCTGTCGCTGGggg                           >  1:9120880/1‑34 (MQ=255)
       ccTGGGTGACGGCTACCGCCATCTGTCGCTGGggg                           >  1:384614/1‑34 (MQ=255)
         tGGGTGACGGCTACCGCCATCTGTCGCTGTgggtt                         >  1:7403825/1‑32 (MQ=255)
         ggggCGACGTCTCCCGCCATCTGTCGCTGTGGTtc                         <  1:2523051/34‑1 (MQ=21)
          gggTGCCGGCTCCCGCCATCTGTCGCTGTGGTtct                        <  1:20381338/35‑1 (MQ=25)
          gggTGCCGGCTCCCGCCATCTGTCGCTGTGGTtct                        <  1:8805171/35‑1 (MQ=25)
          gggTGACGGCTAACGCCATCTGTCGCTGGGGGtct                        >  1:4712347/1‑35 (MQ=25)
           gtggACGGTTCCAGCCATCTGTCGCTGTGGTtctc                       <  1:16630412/32‑1 (MQ=16)
           ggTGCCGCCTCCCGCCATCTGTCGCTGTGGTtctc                       <  1:1852839/35‑1 (MQ=21)
            gTGCCGCCTCCCGCCATCTGTCGCTGTGGTTCTcg                      <  1:17879228/35‑1 (MQ=21)
            gTGCCGCCTACCGCCATCTGTCGCTGTGGTTCTcg                      <  1:12430589/35‑1 (MQ=25)
              gccgcctcccgccatctgTCGCTGTGGTTCTcgcg                    <  1:7131828‑M2/17‑1 (MQ=255)
              gACGGCTCCAGCCATCTGTCGCTGTGGTTCTcgcg                    <  1:9581401/35‑1 (MQ=25)
               gcGGCTACCGCCATCTGTCGCTGTGGTTCTCGCgg                   <  1:317785/34‑1 (MQ=255)
               ccgcctcccgccatctgTCGCTGTGGTTCTCGCgg                   <  1:4716095‑M2/18‑1 (MQ=255)
                cgcctcccgccatctgTCGCTGTGGTTCTCGCGGc                  <  1:14158055‑M2/19‑1 (MQ=255)
                  gctcccgccatctgTCGCTGTGGTTCTCGCGGCtt                <  1:11250287‑M2/21‑1 (MQ=255)
                  gctcccgccatctgTCGCTGTGGTTCTCGCGGCtt                <  1:13858409‑M2/21‑1 (MQ=255)
                    tcccgccatctgTCGCTGTGGTTCTCGCGGCTTaa              <  1:7991570‑M2/23‑1 (MQ=255)
                    tcccgccatctgTCGCTGTGGTTCTCGCGGCTTaa              <  1:7365380‑M2/23‑1 (MQ=255)
                     cccgccatctgTCGCTGTGGTTCTCGCGGCTTaaa             <  1:5110079‑M2/24‑1 (MQ=255)
                     cccgccatctgTCGCTGTGGTTCTCGCGGCTTaaa             <  1:17208693‑M2/24‑1 (MQ=255)
                     cccgccatctgTCGCTGTGGTTCTCGCGGCTTaaa             <  1:6150779‑M2/24‑1 (MQ=255)
                     cccgccatctgTCGCTGTGGTTCTCGCGGCTTaaa             <  1:5605228‑M2/24‑1 (MQ=255)
                     aCCGCCATCTGTCGCTGTGGTTCTCGCGGCTTtaa             >  1:13693236/1‑35 (MQ=255)
                      ccgccatctgTCGCTGTGGTTCTCGCGGCTTAAAg            <  1:20193929‑M2/25‑1 (MQ=255)
                      ccGCCATCTGTCGCTGTGGGTCTCGGGGGTTTAAg            >  1:2991563/1‑35 (MQ=16)
                      cagccatctgTCGCTGTGGTTCTCGCGGCTTAAAg            <  1:241780‑M2/25‑1 (MQ=255)
                       cgccatctgTCGCTGTGGTTCTCGCGGCTTAAAGc           <  1:783318‑M2/26‑1 (MQ=255)
                       cgccatctgTCGCTGTGGTTCTCGCGGCTTAAAGc           <  1:12906909‑M2/26‑1 (MQ=255)
                       cgccatctgTCGCTGTGGTTCTCGCGGCTTAAAGc           <  1:10018240‑M2/26‑1 (MQ=255)
                       cgccatctgTCGCTGTGGTTCTCGCGGCTTAAAGc           <  1:10612336‑M2/26‑1 (MQ=255)
                       cgccatctgTCGCTGTGGTTCTCGCGGCTTAAAGc           <  1:15071661‑M2/26‑1 (MQ=255)
                        gccatctgTCGCTGTGGTTCTCGCGGCTTAAAGcc          <  1:5901218‑M2/27‑1 (MQ=255)
                        gccatctgTCGCTGTGGTTCTCGCGGCTTAAAGcc          <  1:12566607‑M2/27‑1 (MQ=255)
                         ccatctgTCGCTGTGGTTCTCGCGGCTTAAAGCCt         <  1:9552715‑M2/28‑1 (MQ=255)
                         ccatctgTCGCTGTGGTTCTCGCGGCTTAAAGCCt         <  1:8212171‑M2/28‑1 (MQ=255)
                          catctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTg        <  1:13381616‑M2/29‑1 (MQ=255)
                          catctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTg        <  1:1812372‑M2/29‑1 (MQ=255)
                           atctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTgg       <  1:7364137‑M2/30‑1 (MQ=255)
                           atctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTgg       <  1:2027094‑M2/30‑1 (MQ=255)
                            tctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGt      <  1:6332921‑M2/31‑1 (MQ=255)
                            tctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGt      <  1:17717180‑M2/31‑1 (MQ=255)
                             ctgTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTa     <  1:6660471‑M2/32‑1 (MQ=255)
                              tgTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTAt    <  1:5993727‑M2/33‑1 (MQ=255)
                               gTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATc   <  1:9666799‑M2/34‑1 (MQ=255)
                               gTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATc   <  1:1385074‑M2/34‑1 (MQ=255)
                                tCGCTGTGGTTCTCGCGGCTTAAATCCTGGTATCg  >  1:17460341/1‑35 (MQ=255)
                                tCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCg  <  1:12814542/35‑1 (MQ=255)
                                tCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCg  <  1:6923603/35‑1 (MQ=255)
                                tCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCg  <  1:15789288/35‑1 (MQ=255)
                                tCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCg  >  1:5470759/1‑35 (MQ=255)
                                tCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCg  >  1:18608266/1‑35 (MQ=255)
                                tCGCTGTGGTGCTCGCGGCTTAAAGCCTGGTATCg  <  1:14336181/35‑1 (MQ=255)
                                |                                  
TCTATCGCCTGGGTGACGGCTACCGCCATCTGTCGCTGTGGTTCTCGCGGCTTAAAGCCTGGTATCG  >  REL606/1871707‑1871773

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