New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2698091 =73 (0.710)6 (0.070) 3/40 NT 9.3% coding (1992/2145 nt) nrdE ribonucleotide‑diphosphate reductase alpha subunit
?REL606 3668825 = 56 (0.650)coding (873/2031 nt) malS periplasmic alpha‑amylase precursor

Only 100 of 141 total aligned reads displayed.

GGACGAATCTTGATGCCAATATGGGCAACGAAGCCGATCTACGGACGCTGGTTGATAGCGCACATCAGC  >  REL606/3668791‑3668859
                                  |                                  
ggACGAATCTTGATGCCAATATGGGCAACGAAGcc                                    <  1:6921791/35‑1 (MQ=255)
ggACGAATCTTGATGCCAATATGGGCAACGAAGcc                                    >  1:16984949/1‑35 (MQ=255)
ggACGAATCTTGATGCCAATATGGGCAACGAAGcc                                    <  1:20545773/35‑1 (MQ=255)
ggACGAATCTTGATGCCAATATGGGCAACGAAGcc                                    <  1:10037538/35‑1 (MQ=255)
ggACGAATCTTGATGCCAATATGGGCAACGAAGcc                                    >  1:16393811/1‑35 (MQ=255)
 gACGAATCTTGATGCCAATATGGGCAACGAAGCCg                                   >  1:6340327/1‑35 (MQ=255)
 gACGAATCTTGATGCCAATATGGGCAACGAAGCCg                                   >  1:9523896/1‑35 (MQ=255)
 gACGAATCTTGATGCCAATATGGGCAACGAAGCCg                                   >  1:2359645/1‑35 (MQ=255)
  aCGAATCTTGATGCCAATATGGGCAACGAAGCCGa                                  >  1:17953516/1‑35 (MQ=255)
   cGAATCTTTATGCCAATATGGGCAACGAAGCCGAt                                 >  1:1407956/1‑35 (MQ=255)
   cGAATCTTGATGCCAATATGGGCAACGAAGCCGAt                                 >  1:3156607/1‑35 (MQ=255)
    gAATCTTGATGCCAATATGGGCAACGAAGCCGATc                                >  1:11744264/1‑35 (MQ=255)
    gAATCTTGATGCCAATATGGGCAACGAAGCCGATc                                >  1:15191869/1‑35 (MQ=255)
    gAATCTTGATGCCAATATGGGCAACGAAGCCGATc                                >  1:12519541/1‑35 (MQ=255)
    gAATCTTGATGCCAATATGGGCAACGAAGCCGATc                                >  1:3301106/1‑35 (MQ=255)
    gAATCTTGATGCCAATATGGGCAACGAAGCCGATc                                >  1:14659329/1‑35 (MQ=255)
     aaTCTTGATGCCAATATGGGCAACGAAGCCGATCt                               >  1:1668317/1‑35 (MQ=255)
     aaTCTTGATGCCAATATGGGCAACGAAGCCGATCt                               >  1:17398500/1‑35 (MQ=255)
      aTCTTGATGCCAATATGGGCAACTAAGCCGATCTa                              >  1:14687929/1‑35 (MQ=255)
      aTCTTGATGCCAATATGGGCAACGAAGCCGATCTa                              >  1:13516688/1‑35 (MQ=255)
      aTCTTGATGCCAATATGGGCAACGAAGCCGATCTa                              >  1:6445089/1‑35 (MQ=255)
      aTCTTGATGCCAATATGGGCAACGAAGCCGATCTa                              >  1:2782184/1‑35 (MQ=255)
      aTCTTGATGCCAATATGGGCAACGAAGCCGATCTa                              >  1:5006007/1‑35 (MQ=255)
      aTCTTGATGCCAATATGGGCAACGAAGCCGATCTa                              >  1:10092886/1‑35 (MQ=255)
       tCTTGATGCCAATATGGGCAACGAAGCCGATCTAc                             >  1:17739484/1‑35 (MQ=255)
        cTTGATGCCAATATGGGCAACGAAGCCGATCTACg                            >  1:1497711/1‑35 (MQ=255)
        cTTGATGCCAATATGGGCAACGAAGCCGATCTACg                            >  1:10312632/1‑35 (MQ=255)
         ttGATGCCAATATGGGCAACGAAGCCGATCTACgg                           >  1:13239681/1‑35 (MQ=255)
         ttGATGCCAATATGGGCAACGAAGCCGATCTACgg                           >  1:5871324/1‑35 (MQ=255)
          tGATGCCAATATGGGCAACGAAGCCGATCTACGGa                          >  1:1401440/1‑35 (MQ=255)
          tGATGCCAATATGGGCAACGAAGCCGATCTACGGa                          >  1:9253060/1‑35 (MQ=255)
           gATGCCAATATGGGCAACGAAGCCGATCTACTGAc                         >  1:7091320/1‑35 (MQ=255)
             tGCCAATATGGGCAACGAAGCCGATCTACGGACGc                       >  1:19957683/1‑35 (MQ=255)
              gCCAATATGGGCAATGAAGCCGATCTACGGACGCt                      >  1:1237589/1‑35 (MQ=255)
              gCCAATATGGGCAACGAAGCCGATCTACGGACGCt                      >  1:3886005/1‑35 (MQ=255)
              gCCAATATGGGCAACGAAGCCGATCTACGGACGCt                      >  1:1626024/1‑35 (MQ=255)
              gCCAATATGGGCAACGAAGCCGATCTACGGACGCt                      >  1:8126321/1‑35 (MQ=255)
              gCCAATATGGGCAACGAAGCCGATCTACGGACGCt                      >  1:10127855/1‑35 (MQ=255)
               ccAATATGGGCAACGAAGCCGATCTACGGACGCTg                     >  1:20589119/1‑35 (MQ=255)
               caaaTAGGGGAAAAAAAGCCGATCTACGGACGCTg                     <  1:13388490/33‑1 (MQ=16)
                cAATATGGGCAACGAAGCCGATCTACTGACGCTgg                    >  1:8716379/1‑35 (MQ=255)
                cAATATGGGCAACGAAGCCAATCTACTGACGCTgg                    >  1:7965817/1‑35 (MQ=37)
                 aaTATGGGCAACGAAGCCGATCTACGGACGCTGGt                   >  1:3696333/1‑35 (MQ=255)
                   tatGGGCAACGAAGCCGATCTACGGACGCTGGTTg                 >  1:4739450/1‑35 (MQ=255)
                   tatGGGCAACGAAGCCGATCTACGGACGCTGGGTg                 >  1:1344024/1‑35 (MQ=255)
                   taggggCAAAAAAGCCGATCTACGGACGCTGGTTg                 <  1:7383189/32‑1 (MQ=25)
                   taggggCAAAAAAGCCGATCTACGGACGCTGGTTg                 <  1:19053312/32‑1 (MQ=25)
                    gggggcaaaaaagcCGATCTACGGACGCTGGTTGa                <  1:15571958‑M2/21‑1 (MQ=255)
                    atGGGCAACGAAGCCGATCTACGGACGCTGGCTGa                >  1:20185594/1‑35 (MQ=255)
                     tgggaaaaaaagcCGATCTACGGACGCTGGTTGAt               <  1:2190693‑M2/22‑1 (MQ=255)
                     ggggaaaaaaagcCGATCTACGGACGCTGGTTGAt               <  1:2921361‑M2/22‑1 (MQ=255)
                     ggggCAACAAAGCCGATCTACGGACGCTGGTTGAt               <  1:19954489/34‑1 (MQ=255)
                      gggcaaaaaagcCGATCTACGGACGCTGGTTGATa              <  1:9937222‑M2/23‑1 (MQ=255)
                      gggCAACGAAGCCGATCTACGGACGCTGGTTGATa              >  1:1896980/1‑35 (MQ=255)
                      gggCAACGAAGCCGATCTACGGACGCTGGTTGATa              >  1:13147693/1‑35 (MQ=255)
                       ggCAACGAAGCCGATCTACGGACGCTGGTTGATAg             >  1:1001736/1‑35 (MQ=255)
                       ggCAACGAAGCCGATCTACGGACGCTGGTTGATAg             >  1:49546/1‑35 (MQ=255)
                       ggCAACGAAGCCGATCTACGGACGCTGGTTGATAg             >  1:4848132/1‑35 (MQ=255)
                       ggCAACAAAGCCGATCTACGGACGCTGGTTGATAg             <  1:4510543/35‑1 (MQ=255)
                        gaaaaaaagcCGATCTACGGACGCTGGTTGATAgc            <  1:13022735‑M2/25‑1 (MQ=255)
                        gCAACGAAGCCGATCTACGGACGCTGGTTGATAgc            >  1:18104596/1‑35 (MQ=255)
                        caaagaaagcCGATCTACGGACGCTGGTTGATAgc            <  1:16148801‑M2/25‑1 (MQ=255)
                        acaaaaaagcCGATCTACGGACGCTGGTTGATAgc            <  1:18402055‑M2/25‑1 (MQ=255)
                        aaaaaaaagcCGATCTACGGACGCTGGTTGATAgc            <  1:17736070‑M2/25‑1 (MQ=255)
                        aaaaaaaagcCGATCTACGGACGCTGGTTGATAgc            <  1:16328580‑M2/25‑1 (MQ=255)
                         cAACGAAGCCGATCTACGGACGCTGGTTGATAggg           >  1:3884148/1‑33 (MQ=255)
                         cAACGAAGCCGATCTACGGACGCTGGTCGATAgcg           >  1:5316958/1‑35 (MQ=255)
                         cAACAAAGCCGATCTACGGACGCTGGTTGATAgcg           <  1:4404002/35‑1 (MQ=255)
                          aaCGAAGCCGATCTACGGACGCTGGTTGATAgcgc          >  1:4152275/1‑35 (MQ=255)
                          aaCGAAGCCGATCTACGGACGCTGGTTGATAgcgc          >  1:9765939/1‑35 (MQ=255)
                          aaCGAAGCCGATCTACGGACGCTGGTTGATAgcgc          >  1:2568177/1‑35 (MQ=255)
                           acaaagcCGATCTACGGACGCTGGTTGATAGCGca         <  1:7022865‑M2/28‑1 (MQ=255)
                           aCGAAGCCGATCTACGGACGCTGGTTGATAGCGca         >  1:9598018/1‑35 (MQ=255)
                           aCGAAGCAGATCTACGGACGCTGGTTGATAGCGca         <  1:15940671/35‑1 (MQ=255)
                            caaagcCGATCTACGGACGCTGGTTGATAGCGcac        <  1:5676354‑M2/29‑1 (MQ=255)
                            cGAAGCCGATCTACGGACGCTGGTTGATAGCGcac        >  1:8665726/1‑35 (MQ=255)
                            aaaagcCGATCTACGGACGCTGGTTGATAGCGcac        <  1:5344334‑M2/29‑1 (MQ=255)
                            aaaagcCGATCTACGGACGCTGGTTGATAGCGcac        <  1:8654633‑M2/29‑1 (MQ=255)
                            aaaagcCGATCTACGGACGCTGGTTGATAGCGcac        <  1:13597825‑M2/29‑1 (MQ=255)
                             gAAGCCGATCTACGGACGCTGGTTGATAGCGcaca       >  1:11473085/1‑35 (MQ=255)
                             aaagcCGATCTACGGACGCTGGTTGATAGCGcaca       <  1:4659174‑M2/30‑1 (MQ=255)
                             aaagcCGATCTACGGACGCTGGTTGATAGCGcaca       <  1:6989485‑M2/30‑1 (MQ=255)
                             aaagcCGATCTACAGACGCTGGTTGATAGCGcaca       <  1:5730130‑M2/30‑1 (MQ=255)
                              aagcCGATCTACGGACGCTGGTTGATAGCGCACAt      >  1:8590450‑M2/5‑35 (MQ=255)
                               agcCGATCTACGGACGCTGGTTGATAGCGCACATc     >  1:14541082‑M2/4‑35 (MQ=255)
                               agcCGATCTACGGACGCTGGTTGATAGCGCACATc     <  1:6930891‑M2/32‑1 (MQ=255)
                               agcAGATCTACGGACGCTGGTTGATAGCGCACATc     <  1:2210592‑M2/32‑1 (MQ=255)
                                gcCGATCTACGGACGCTGGTTGATAGCGCACCTCa    >  1:15895738‑M2/3‑35 (MQ=255)
                                gcCGATCTACGGACGCTGGTTGATAGCGCACATCa    >  1:1268475‑M2/3‑35 (MQ=255)
                                gcCGATCTACGGACGCTGGTTGATAGCGCACATCa    <  1:11179803‑M2/33‑1 (MQ=255)
                                 cCGATCTACGGACGCTGGTTGATAGCGCACATCAg   <  1:3730529‑M2/34‑1 (MQ=255)
                                 cCGATCTACGGACGCTGGTTGATAGCGCACATCAg   <  1:18424523‑M2/34‑1 (MQ=255)
                                 cCGATCTACGGACGCTGGTTGATAGCGCACATCAg   >  1:3822356‑M2/2‑35 (MQ=255)
                                 cCGATCTACGGACGCTGGTTGATAGCGCACATCAg   <  1:13839499‑M2/34‑1 (MQ=255)
                                 cCGATCTACGGACGCTGGTTGATAGCGCACATCAg   <  1:9120566‑M2/34‑1 (MQ=255)
                                 cCGATCTACGGACGCTGGTTGATAGCGCACATCAg   >  1:16468031‑M2/2‑35 (MQ=255)
                                  cGATCTACGGACGCTGGTTGATAGCGCACATCAgc  >  1:736260/1‑35 (MQ=255)
                                  cGATCTACGGACGCTGGTTGATAGCGCACATCAgc  >  1:10895575/1‑35 (MQ=255)
                                  cGATCTACGGACGCTGGTTGATAGCGCACATCAgc  <  1:18941850/35‑1 (MQ=255)
                                  cGATCTACGGACGCTGGTTGATAGCGCACATCAgc  <  1:2582577/35‑1 (MQ=255)
                                  |                                  
GGACGAATCTTGATGCCAATATGGGCAACGAAGCCGATCTACGGACGCTGGTTGATAGCGCACATCAGC  >  REL606/3668791‑3668859

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