Predicted mutation
evidence seq id position mutation annotation gene description
JC REL606 1,733,297 (CGTATGAA)1→2 coding (333/1413 nt) pykF → pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1733288 =44 (0.260)167 (1.060) 62/68 0.1 83.4% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 25 (0.160)coding (331/1413 nt) pykF pyruvate kinase

Only 100 of 368 total aligned reads displayed.

TAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGC  >  REL606/1733252‑1733324
                                    |                                    
  aaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAAc                                    <  2:10504316/37‑1 (MQ=255)
  aaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAAc                                    <  1:185817/37‑1 (MQ=255)
   tcggcaacagcgaaatggttgcggtaacgtatgaacg                                   <  2:7592355‑M1/4‑1 (MQ=255)
   tcggcaaaagcgaaatggttgcggtaacgtatgaacg                                   <  1:9005068‑M1/4‑1 (MQ=255)
    tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACgg                                  >  1:4438472/1‑36 (MQ=255)
    tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGt                                  >  1:10010119/1‑37 (MQ=255)
    tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGt                                  >  1:5831826/1‑37 (MQ=255)
    tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGt                                  <  2:5661497/37‑1 (MQ=255)
    cggcaacagcgaaatggttgcggtaacgtatgaacgt                                  >  1:1641789‑M1/33‑37 (MQ=255)
    cggcaacagcgaaatggttgcggtaacgtatgaacgt                                  >  2:555802‑M1/33‑37 (MQ=255)
     ggcaacagcgaaatggttgcggtaacgtatgaacgta                                 >  1:10178986‑M1/32‑37 (MQ=255)
     ggcaacagcgaaatggttgcggtaacgtatgaacgta                                 >  2:10433854‑M1/32‑37 (MQ=255)
     ggcaacagcgaaatggttgcggtaacgtatgaacgta                                 >  2:1165708‑M1/32‑37 (MQ=255)
     ggcaacagcgaaatggttgcggtaacgtatgaacgta                                 <  1:693371‑M1/6‑1 (MQ=255)
     ggcaacagcgaaatggttgaggtaacgtatgaacgta                                 <  1:249950‑M1/6‑1 (MQ=255)
     cTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTa                                 <  1:3660394/37‑1 (MQ=255)
      tGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTAt                                >  1:7408390/1‑37 (MQ=255)
      tGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTAt                                <  1:2587753/37‑1 (MQ=255)
      gcaacagcgaaatggttgcggtaacgtatgaacgtat                                >  2:1458690‑M1/31‑37 (MQ=255)
      gcaacagcgaaatggttgcggtaacgtatgaacgtat                                <  2:5267228‑M1/7‑1 (MQ=255)
      gcaacagcgaaatggttgcggtaacgtatgaacgtat                                >  2:4546330‑M1/31‑37 (MQ=255)
       gTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATg                               >  2:7795733/1‑37 (MQ=255)
       caacagcgaaatggttgcggtaacgtatgaacgtatg                               >  1:5088783‑M1/30‑37 (MQ=255)
       caacagagaaatggttgcggtaacgtatgaacgtatg                               <  2:1934179‑M1/8‑1 (MQ=255)
        ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGa                              <  1:4538174‑M2/37‑2 (MQ=255)
        ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGa                              >  1:7533943‑M2/1‑36 (MQ=255)
        ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGa                              >  2:2824058‑M2/1‑36 (MQ=255)
        aacagcgaaatggttgcggtaacgtatgaacgtatga                              <  2:3657511‑M1/9‑1 (MQ=255)
        aacagcgaaatggttgcggtaacgtatgaacgtatga                              <  1:9039265‑M1/9‑1 (MQ=255)
         tATCGGCAACAGCGAAATGGTTGCGGTAACGTATTaa                             >  1:2937640/1‑37 (MQ=255)
         tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                             <  1:11220662‑M2/37‑3 (MQ=255)
         tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                             <  2:9176085‑M2/37‑3 (MQ=255)
         tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                             <  2:8788273‑M2/37‑3 (MQ=255)
         tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                             >  2:3094624‑M2/1‑35 (MQ=255)
         acagcgaaatggttgcggtaacgtatgaacgtatgaa                             <  1:6601678‑M1/10‑1 (MQ=255)
         acagcgaaatggttgcggtaacgtatgaacgtatgaa                             >  1:3444770‑M1/28‑37 (MQ=255)
          cagcgaaatggttgcggtaacgtatgAACGTATGAAg                            >  1:6323871‑M1/27‑37 (MQ=255)
          cagcgaaatggttgcggtaacgtatgAACGTATGAAg                            >  2:659777‑M1/27‑37 (MQ=255)
          aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaac                            <  2:3273009‑M2/37‑4 (MQ=255)
          aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaac                            >  2:3982839‑M2/1‑34 (MQ=255)
          aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaac                            <  2:6889966‑M2/37‑4 (MQ=255)
          aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaac                            <  2:9756950‑M2/37‑4 (MQ=255)
          aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaac                            <  1:10999226‑M2/37‑4 (MQ=255)
           tCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                           <  2:5816846‑M2/37‑5 (MQ=255)
           tCGGCAACAGCGAAATGGTTGCGGGAACGCACGaacg                           <  1:1031916‑M2/37‑5 (MQ=255)
           agcgaaatggttgcggtaacgtatgAACGTATGAAgg                           >  2:2527260‑M1/26‑37 (MQ=255)
             cgaaatgggtgcggtaacgtatgAAAGTATGGAGGtt                         <  1:9982849‑M1/14‑1 (MQ=255)
              gaaatggttgcggtaacgtatgAACGTATGAAGGttt                        <  1:8915953‑M1/15‑1 (MQ=255)
              gaaatggttgcggtaacgtatgAACGTATGAAGGttt                        >  2:4707994‑M1/23‑37 (MQ=255)
              gaaatggttgcggtaacgtatgAACGTATGAAGGttt                        <  2:3404503‑M1/15‑1 (MQ=255)
              gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                        >  2:7214879‑M2/1‑30 (MQ=255)
              gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                        <  2:203874‑M2/37‑8 (MQ=255)
              gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                        <  2:10099834‑M2/37‑8 (MQ=255)
              gCAACAGCGAAATGGTGGCGGTAACGTATGaacgtat                        >  1:5311400‑M2/1‑30 (MQ=255)
               cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                       <  1:7661643‑M2/37‑9 (MQ=255)
               cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                       <  1:8862834‑M2/37‑9 (MQ=255)
               cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                       <  1:8226252‑M2/37‑9 (MQ=255)
               aaattgttgcggtaacgtatgAACGTATGAAGGGTTc                       <  1:6794646‑M1/16‑1 (MQ=255)
               aaatggttgcggtaacgtatgAACGTATGAAGGTTTc                       <  1:4623682‑M1/16‑1 (MQ=255)
                aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatga                      >  2:413987‑M2/1‑28 (MQ=255)
                 atggttgcggtaacgtatgAACGTATGAAGGTTTCac                     <  2:6948048‑M1/18‑1 (MQ=255)
                  tggttgcggtaacgtatgAACGTATGACGGTTTCact                    <  2:2110293‑M1/19‑1 (MQ=255)
                  cAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                    >  2:1536981‑M2/1‑26 (MQ=255)
                  cAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                    <  1:6836907‑M2/37‑12 (MQ=255)
                  cAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                    <  1:1161049‑M2/37‑12 (MQ=255)
                    gCGAAATGGTTGCGGTAACGTATGaacgtatgaaggt                  <  2:10971723‑M2/37‑14 (MQ=255)
                    gCGAAATGGTTGCGGTAACGTATGaacgtatgaaggt                  <  1:9546815‑M2/37‑14 (MQ=255)
                     ttgcggtaacgtatgAACGTATGAAGGTTTCactact                 <  1:4382645‑M1/22‑1 (MQ=255)
                      tgcggtaacgtatgAACGTATGAAGGTTTCACTACTg                >  2:3574600‑M1/15‑37 (MQ=255)
                      gAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                >  1:8773272‑M2/1‑22 (MQ=255)
                       gcggtaacgtatgAACGTATGAAGGTTTCACTACTGa               >  2:4582642‑M1/14‑37 (MQ=255)
                        aaTGGTTGCGGTAACGTATGaacgtatgaaggtttca              >  1:410851‑M2/1‑20 (MQ=255)
                        aaTGGTTGCGGTAACGTATGaacgtatgaaggtttca              <  2:5520566‑M2/37‑18 (MQ=255)
                         ggtaacgtatgAACGTATGAAGGTTTCACTACTGAcc             >  2:10970601‑M1/12‑37 (MQ=255)
                         ggtaacgtatgAACGTATGAAGGTTTCACTACTGAcc             >  1:8177997‑M1/12‑37 (MQ=255)
                          tGGTTGCGGTAACGTATGaacgtatgaaggtttcact            <  2:8532045‑M2/37‑20 (MQ=255)
                           taacgtatgAACGTATGAAGGTTTCACTACTGACctg           >  2:10057260‑M1/10‑37 (MQ=255)
                           taacgtatgAACGTATGAAGGTTTCACTACTGACctg           <  1:7027845‑M1/28‑1 (MQ=255)
                            gTTGCGGTAACGTATGaacgtatgaaggtttcactac          <  2:9587296‑M2/37‑22 (MQ=255)
                            gTTGCGGTAACGTATGaacgtatgaaggtttcactac          <  2:8093708‑M2/37‑22 (MQ=255)
                            gTTGCGGTAACGTATGaacgtatgaaggtttcactac          >  1:4272196‑M2/1‑16 (MQ=255)
                            aacgtatgAACGTATGAAGGTTTCACTACTGACctgt          <  1:4804349‑M1/29‑1 (MQ=255)
                              cgtatgAACGTATGAAGGTTTCACTACTGACctgtct        >  2:5154531‑M1/7‑37 (MQ=255)
                              cgtatgAACGTATGAAGGTTTCACTACTGACctgtct        >  1:9259183‑M1/7‑37 (MQ=255)
                               gtatgAACGTATGAAGGTTTCACTACTGACctgtctg       <  1:3824335‑M1/32‑1 (MQ=255)
                               gCGGTAACGTATGaacgtatgaaggtttcactactga       <  1:2424569‑M2/37‑25 (MQ=255)
                                tatgAACGTATGAAGGTTTCACTACTGACctgtctgt      >  2:4860502‑M1/5‑37 (MQ=255)
                                tatgAACGTATGAAGGTTTCACTACTGACctgtctgt      <  1:2988229‑M1/33‑1 (MQ=255)
                                tatgAACGTATGAAGGTTTCACTACTGACctgtctgt      <  2:7104240‑M1/33‑1 (MQ=255)
                                cGGTAACGTATGaacgtatgaaggtttcactactgac      <  1:7910199‑M2/37‑26 (MQ=255)
                                cGGTAACGTATGaacgtatgaaggtttcactactgac      >  2:3815429‑M2/1‑12 (MQ=255)
                                cGGTAACGTATGaacgtatgaaggtttcactactgac      >  1:5472403‑M2/1‑12 (MQ=255)
                                cGGTAACGTATGaacgtatgaaggtttcactactgac      <  1:9939857‑M2/37‑26 (MQ=255)
                                 atgAACGTATGAAGGTTTCACTACTGACCTGTCTGtt     <  2:2231288‑M1/34‑1 (MQ=255)
                                  tgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTg    >  2:8354281‑M1/3‑37 (MQ=255)
                                  gTAGCGTATGaacgtatgaaggtttcactactgacct    >  2:6758301‑M2/1‑10 (MQ=255)
                                  gTAACGTATGaacgtatgaaggtttcactactgacct    >  1:9076634‑M2/1‑10 (MQ=255)
                                    aacgtatgaacgtatgaaggtttcactactgacctgt  <  1:619465‑M2/37‑30 (MQ=255)
                                    AACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGc  <  1:2145993‑M1/37‑1 (MQ=255)
                                    AACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGc  <  1:712118‑M1/37‑1 (MQ=255)
                                    |                                    
TAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGC  >  REL606/1733252‑1733324

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