New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3820609 =84 (0.820)7 (0.100) 3/34 NT 12.1% coding (425/924 nt) yidL predicted DNA‑binding transcriptional regulator
?REL606 = 4248913 42 (0.580)coding (568/714 nt) aphA acid phosphatase/phosphotransferase, class B, non‑specific

CAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCA  >  REL606/4248880‑4248945
                                 |                                
cAATGGCTGCAGGATAAAAATATCCGAATTTTTTa                                 >  1:10258965‑M2/1‑34 (MQ=255)
 aaTGGCTGCAGGATAAAAATATCCGAATTTTTTat                                >  1:9767262‑M2/1‑33 (MQ=255)
 aaTGGCTGCAGGATAAAAATATCCGAATTTTTTat                                >  1:7770036‑M2/1‑33 (MQ=255)
  aTTGCTGCAGGATAAAAATATCCGAATTTTTTatg                               <  1:16605262‑M2/35‑4 (MQ=255)
  aTGGCTGCAGGATAAAAATATCCGAATTTTTTatg                               >  1:721411‑M2/1‑32 (MQ=255)
  aTGGCTGCAGGATAAAAATATCCGAATTTTTTatg                               >  1:3752440‑M2/1‑32 (MQ=255)
  aTGGCTGCAGGATAAAAATATCCGAATTTTTTatg                               >  1:1594552‑M2/1‑32 (MQ=255)
  aTGGCTGCAGGATAAAAATATCCGAATTTTTTatg                               >  1:13956868‑M2/1‑32 (MQ=255)
   tGGCTGCAGGATAAAAATATCCGAATTTTTTatgg                              <  1:9399895‑M2/35‑5 (MQ=255)
   tGGCTGCAGGATAAAAATATCCGAATTTTTTatgg                              <  1:18581216‑M2/35‑5 (MQ=255)
   tGGCTGCAGGATAAAAATATCCGAATTTTTTatgg                              <  1:10018677‑M2/35‑5 (MQ=255)
    ggCTGCAGGATAAAAATATCCGAATTTTTTatggc                             >  1:1973674‑M2/1‑30 (MQ=255)
    ggCTGCAGGATAAAAATATCCGAATTTTTTatggc                             >  1:15086578‑M2/1‑30 (MQ=255)
     gCTGCAGGATAAAAATATCCGAATTTTTTatggcg                            >  1:3972729‑M2/1‑29 (MQ=255)
     gCTGCAGGATAAAAATATCCGAATTTTTTatggcg                            >  1:14941281‑M2/1‑29 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           >  1:11191661‑M2/1‑28 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           <  1:16062957‑M2/35‑8 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           >  1:17078098‑M2/1‑28 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           >  1:5542878‑M2/1‑28 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           <  1:17392060‑M2/35‑8 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           <  1:7842982‑M2/35‑8 (MQ=255)
      cTGCAGGATAAAAATATCCGAATTTTTTatggcga                           >  1:9602989‑M2/1‑28 (MQ=255)
       tGCAGGATAAAAATATCCGAATTTTTTatggcgac                          <  1:9471682‑M2/35‑9 (MQ=255)
       tGCAGGATAAAAATATCCGAATTTTTTatggcgaa                          <  1:19507092‑M2/35‑9 (MQ=255)
       tGCAGGATAAAAATATCCGAATTTTTTatggcgaa                          <  1:14047611‑M2/35‑9 (MQ=255)
       tGCAGGATAAAAATATCCGAATTTTTTATGGCGAt                          <  1:18333425/35‑1 (MQ=255)
       tGCAGGATAAAAATATCCGAATTTTTTATGGCGAt                          <  1:14281419/35‑1 (MQ=255)
       tGCAGGATAAAAATATCCGAATTTTTTATGGCGAt                          <  1:5739229/35‑1 (MQ=255)
        gCAGGATAAAAATATCCGAATTTTTTatggcgaaa                         <  1:9325125‑M2/35‑10 (MQ=255)
        gCAGGATAAAAATATCCGAATTTTTTATGGCGAtt                         <  1:1968877/35‑1 (MQ=255)
        gCAGGATAAAAATATCCGAATTTTTTATGGCGAtt                         <  1:14364024/35‑1 (MQ=255)
        gCAGGATAAAAATATCCGAATTTTTTATGGCGAtt                         <  1:14037092/35‑1 (MQ=255)
         ccggATAAAAATATCCGAATTTTTTATGGCGATTc                        <  1:7251430/33‑1 (MQ=255)
         cAGGATAAAAATATCCGAATTTTTTatggcgaaat                        <  1:838534‑M2/35‑11 (MQ=255)
         cAGGATAAAAATATCCGAATTTTTTATGGCGATTc                        <  1:6216861/35‑1 (MQ=255)
         cAGGATAAAAATATCCGAATTTTTTATGGCGATTc                        <  1:11683076/35‑1 (MQ=255)
         cAGGATAAAAATATCCGAATTTTTTATGGCGATTc                        <  1:15945157/35‑1 (MQ=255)
         cAGGATAAAAATATCCGAATTTTTTATGGCGATTc                        <  1:4196636/35‑1 (MQ=255)
           ggATAAAAATATCCGAATTTTTTATGGCGATTCTg                      <  1:14785839/35‑1 (MQ=255)
           ggATAAAAATATCCGAATTTTTTATGGCGATTCTg                      <  1:8570232/35‑1 (MQ=255)
            gATAAAAATATCCGAATTTTTTATGGCGATTCTGa                     <  1:9047392/35‑1 (MQ=255)
             aTAAAAATATCCGAATTTTTTatggcgaaaattat                    <  1:5860789‑M2/35‑15 (MQ=255)
             aTAAAAATATCCGAATTTTTTatggcgaaaattaa                    <  1:14252996‑M2/35‑15 (MQ=255)
             aTAAAAATATCCGAATTTTTTATGGCGATTCTGAt                    <  1:3375435/35‑1 (MQ=255)
             aTAAAAATATCCGAATTTTTTATGGCGATTCTGAt                    <  1:19875002/35‑1 (MQ=255)
              tAAAAATATCCGAATTTTTTatggcgaaaattata                   <  1:7664345‑M2/35‑16 (MQ=255)
              tAAAAATATCCGAATTTTTTATGGCGATTCTGATa                   <  1:3281761/35‑1 (MQ=255)
              tAAAAATATCCGAATTTTTTATGGCGATTCTGATa                   <  1:5114537/35‑1 (MQ=255)
               aaaaaTATCCGAATTTTTTATGGCGATTCTGATaa                  <  1:18791123/35‑1 (MQ=255)
               aaaaaTATCCGAATTTTTTATGGCGATTCTGATaa                  <  1:11574278/35‑1 (MQ=255)
               aaaaaTATCCGAATTTTTTATGGCGATTCTGATaa                  <  1:729308/35‑1 (MQ=255)
               aaaaaTATCCGAATTTTTTATGGCGATTCTGATaa                  <  1:3144701/35‑1 (MQ=255)
                aaaaTATCCGAATTTTTTatggcgaaaattataat                 <  1:5475791‑M2/35‑18 (MQ=255)
                aaaaTATCCGAATTTTTTatggcgaaaattataat                 <  1:13713380‑M2/35‑18 (MQ=255)
                aaaaTATCCGAATTTTTTATGGCGATTCTGATAAt                 <  1:14676536/35‑1 (MQ=255)
                aaaaTATCCGAATTTTTTATGGCGATTCTGATAAt                 <  1:18666576/35‑1 (MQ=255)
                aaaaTATCCGAATTTTTTATGGCGATTCTGATAAt                 <  1:11634520/35‑1 (MQ=255)
                aaaaTATCCGAATTTTTTATGGCGATTCTGATAAt                 <  1:17022984/35‑1 (MQ=255)
                 aaaTATCCGAATTTTTTatggcgaaacttataatt                <  1:15657991‑M2/35‑19 (MQ=255)
                 aaaTATCCGAATTTTTTatggcgaaaattataatt                <  1:14121080‑M2/35‑19 (MQ=255)
                 aaaTATCCGAATTTTTTATGGCGATTCTGATAATg                <  1:15133103/35‑1 (MQ=255)
                 aaaTATCCGAATTTTTTATGGCGAAAATGATAAta                >  1:16903545/1‑34 (MQ=18)
                   atatCCGAATTTTTTATGGCGATTCTGATAATGaa              <  1:17135076/35‑2 (MQ=255)
                   atatCCGAATATTTTATGGCGATTCTGATAATGat              <  1:13850058/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:1050695/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:14229298/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:15177955/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:15782236/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:3027161/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:4615499/35‑1 (MQ=255)
                    tatCCGAATTTTTTATGGCGATTCTGATAATGata             <  1:5662095/35‑1 (MQ=255)
                     atCCGAATTTTTTATGGCGATTCTGATAATGatat            <  1:17647134/35‑1 (MQ=255)
                     atCCGAATTTTTTATGGCGATTCTGATAATGatat            <  1:20063168/35‑1 (MQ=255)
                     atCCGAATTTTTTATGGCGATTCTGATAATGatat            <  1:9092214/35‑1 (MQ=255)
                     atCCGAATTTTTTATGGCGATTCTGATAATGatat            <  1:3803967/35‑1 (MQ=255)
                      tCCGAATTTTTTATGGCGATTCTGATAATGATAtt           <  1:9652118/35‑1 (MQ=255)
                      tCCGAATTTTTTATGGCGATTCTGATAATGATAtt           <  1:20218982/35‑1 (MQ=255)
                      tCCGAATTTTTTATGGCGATTCTGATAATGATAtt           <  1:545214/35‑1 (MQ=255)
                      tCCGAATTTTTTATGGCGATTCTGATAATGATAtt           <  1:7298983/35‑1 (MQ=255)
                      tCCGAATTTTTTATGGCGATTCTGATAATGATAtt           <  1:1126075/35‑1 (MQ=255)
                       ccGAATTTTTTATGGCGATTCTGATAATGATATTa          <  1:5776831/35‑1 (MQ=255)
                        cGAATTTTTTATGGCGATTCTGATAATGATATTAc         <  1:190287/35‑1 (MQ=255)
                        cGAATTTTTTATGGCGATTCTGATAATGATATTAc         <  1:5258624/35‑1 (MQ=255)
                        cGAATTTTTTATGGCGATTCTGATAATGATATTAc         <  1:16802918/35‑1 (MQ=255)
                        cGAATTTTTTATGGCGATTCTGATAATGATATTAc         <  1:20202648/35‑1 (MQ=255)
                         gAATTTTTTATGGCGATTCTGATAATGATATTaac        <  1:13121873/35‑3 (MQ=255)
                         gAATTTTTTATGGCGATTCTGATAATGATATTAcc        <  1:681046/35‑1 (MQ=255)
                          aaTTTTTTATGGCGATTCTGATAATGATATTAccg       <  1:19873660/35‑1 (MQ=255)
                          aaTTTTTTATGGCGATTCTGATAATGATATTAccg       <  1:16818515/35‑1 (MQ=255)
                               tttATGGCGATTCTGATAATGATATTAACGCCGCa  <  1:10671232/35‑1 (MQ=255)
                               tttATGGCGATTCTGATAATGATATTAACGCCGCa  <  1:5640532/35‑1 (MQ=255)
                                 |                                
CAATGGCTGCAGGATAAAAATATCCGAATTTTTTATGGCGATTCTGATAATGATATTACCGCCGCA  >  REL606/4248880‑4248945

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