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 =3 (0.070)40 (0.930) 33/146 0.3 92.0% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 4 (0.090)coding (331/1413 nt) pykF pyruvate kinase

TTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733213‑1733363
                                                                           |                                                                           
aagctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatga                                                                             >  1:756785‑M1/76‑76 (MQ=255)
 tctcTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTaa                                                                            <  1:982171/76‑1 (MQ=255)
  actggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaac                                                                           >  1:1179234‑M1/74‑76 (MQ=255)
   ctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacg                                                                          <  2:1366050‑M1/4‑1 (MQ=255)
   ctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacg                                                                          <  2:857115‑M1/4‑1 (MQ=255)
    ctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGt                                                                         >  1:1032367/1‑76 (MQ=255)
      gAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTAt                                                                       <  2:810770/76‑1 (MQ=255)
       tcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatg                                                                      >  2:152188‑M1/69‑76 (MQ=255)
       tcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatg                                                                      <  2:58157‑M1/8‑1 (MQ=255)
        aaGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGa                                                                     >  1:756785‑M2/1‑75 (MQ=255)
         agacctttactttcaccactgataaatctgttatcgtcaacagcgaaatggttgcggnaacgtatgaacgtatgaa                                                                    >  1:1340952‑M1/67‑76 (MQ=255)
          actGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaac                                                                   >  1:1179234‑M2/2‑73 (MQ=255)
           cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                  <  2:1366050‑M2/76‑5 (MQ=255)
           cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                  <  2:857115‑M2/76‑5 (MQ=255)
              tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                               <  1:1035936‑M1/15‑1 (MQ=255)
               ttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTc                                                              >  1:144201‑M1/61‑76 (MQ=255)
               ttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggaaaggtgtgAACGTATGAAGGTTTc                                                              <  1:1321484‑M1/16‑1 (MQ=255)
               tCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                                                              >  2:152188‑M2/1‑68 (MQ=255)
               tCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                                                              <  2:58157‑M2/76‑9 (MQ=255)
                ttacttccaccactgataatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCa                                                             >  2:1372808‑M1/60‑76 (MQ=255)
                tactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCa                                                             >  1:921309‑M1/60‑76 (MQ=255)
                 actttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGGTTCac                                                            >  1:487072‑M1/59‑76 (MQ=255)
                 aGACCTTTACTTTCACCACTGATAAATCTGTTATCGTCAACAGCGAAATGGTTGCGGNAACGTATGaacgtatgaa                                                            >  1:1340952‑M2/1‑66 (MQ=255)
                   tttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCacta                                                          >  1:706571‑M1/57‑76 (MQ=255)
                     tcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactact                                                        >  1:744960‑M1/55‑76 (MQ=255)
                      tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                       <  1:1035936‑M2/76‑16 (MQ=255)
                      caccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTg                                                       >  1:283443‑M1/54‑76 (MQ=255)
                       ttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttc                                                      >  1:144201‑M2/1‑60 (MQ=255)
                       ttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGAAAGGTGTGaacgtatgaaggtttc                                                      <  1:1321484‑M2/76‑17 (MQ=255)
                       ttACTTCCACCACTGATAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttca                                                     >  2:1372808‑M2/1‑59 (MQ=255)
                       accactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGa                                                      <  1:1048110‑M1/24‑1 (MQ=255)
                        tACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttca                                                     >  1:921309‑M2/1‑59 (MQ=255)
                         cactcataaatctgttagcggcaacagcgaaatggttgcggtaacgtatgCACGTATGAAGGTTACACTACTGNcc                                                    <  2:1316788‑M1/26‑1 (MQ=255)
                         aCTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaagggttcac                                                    >  1:487072‑M2/1‑58 (MQ=255)
                           tttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcacta                                                  >  1:706571‑M2/1‑56 (MQ=255)
                            tgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgt                                                 <  2:920018‑M1/29‑1 (MQ=255)
                            tgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGAATGAAGGTTTCACTACTGACctgt                                                 <  1:1090395‑M1/29‑1 (MQ=255)
                             tCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactact                                                >  1:744960‑M2/1‑54 (MQ=255)
                              caccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactg                                               >  1:283443‑M2/1‑53 (MQ=255)
                               taaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctg                                              >  2:331652‑M1/45‑76 (MQ=255)
                               accacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactga                                              <  1:1048110‑M2/76‑25 (MQ=255)
                                aaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctgt                                             >  1:559486‑M1/44‑76 (MQ=255)
                                aaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctgt                                             <  1:389685‑M1/33‑1 (MQ=255)
                                 cacTCATAAATCTGTTAGCGGCAACAGCGAAATGGTTGCGGTAACGTATGcacgtatgaaggttacactactgncc                                            <  2:1316788‑M2/76‑27 (MQ=255)
                                  atctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTg                                           >  2:1001226‑M1/42‑76 (MQ=255)
                                  atctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTg                                           >  1:469218‑M1/42‑76 (MQ=255)
                                    tGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgt                                         <  2:920018‑M2/76‑30 (MQ=255)
                                    tGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgaatgaaggtttcactactgacctgt                                         <  1:1090395‑M2/76‑30 (MQ=255)
                                     tgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCa                                        <  1:1048132‑M1/38‑1 (MQ=255)
                                      gttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCATTACTGACCTGTCTGTTGGCaa                                       >  2:60211‑M1/38‑76 (MQ=255)
                                       tAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctg                                      >  2:331652‑M2/1‑44 (MQ=255)
                                        aaaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgt                                     <  1:389685‑M2/76‑34 (MQ=255)
                                        aaaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgt                                     >  1:559486‑M2/1‑43 (MQ=255)
                                          aTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttg                                   >  1:469218‑M2/1‑41 (MQ=255)
                                          aTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttg                                   >  2:1001226‑M2/1‑41 (MQ=255)
                                             tGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggca                                <  1:1048132‑M2/76‑39 (MQ=255)
                                              gTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcattactgacctgtctgttggcaa                               >  2:60211‑M2/1‑37 (MQ=255)
                                               aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                              >  1:893721‑M1/29‑76 (MQ=255)
                                                  agcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCANTACTGACCTGTCTGTTGGCGACACCGTACTGGt                           >  2:1191761‑M1/26‑76 (MQ=255)
                                                     gaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGa                        >  1:1036348‑M1/23‑76 (MQ=255)
                                                     gaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGa                        <  2:142914‑M1/54‑1 (MQ=255)
                                                       aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                      >  1:893721‑M2/1‑28 (MQ=255)
                                                          aGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcantactgacctgtctgttggcgacaccgtactggt                   >  2:1191761‑M2/1‑25 (MQ=255)
                                                           nttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATgg                  <  1:429758‑M1/60‑1 (MQ=255)
                                                             gAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttga                >  1:1036348‑M2/1‑22 (MQ=255)
                                                             gAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttga                <  2:142914‑M2/76‑55 (MQ=255)
                                                               cggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTg              <  1:847246‑M1/64‑1 (MQ=255)
                                                                ggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGa             >  2:1065688‑M1/12‑76 (MQ=255)
                                                                 ttaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGAt            <  2:814935‑M1/66‑1 (MQ=255)
                                                                   nttgCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatgg          <  1:429758‑M2/75‑61 (MQ=255)
                                                                   aacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCg          >  1:1237714‑M1/9‑76 (MQ=255)
                                                                    acgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCgg         <  2:1221011‑M1/69‑1 (MQ=255)
                                                                       tatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTAt      >  2:1034209‑M1/5‑76 (MQ=255)
                                                                       cGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctg      <  1:847246‑M2/76‑65 (MQ=255)
                                                                        ggTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctga     >  2:1065688‑M2/1‑11 (MQ=255)
                                                                        atgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATg     <  1:416790‑M1/73‑1 (MQ=255)
                                                                         ttaACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgat    <  2:814935‑M2/75‑67 (MQ=255)
                                                                          gAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGa   >  1:215039‑M1/2‑76 (MQ=255)
                                                                          gAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCNGATCGGTATGGa   >  1:13051‑M1/2‑76 (MQ=255)
                                                                           aacgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcg  >  1:1237714‑M2/1‑8 (MQ=255)
                                                                           AACGTACGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTCGTTGACGATGGTCTGATCGGTATGGaa  >  1:1380277‑M1/1‑76 (MQ=255)
                                                                           |                                                                           
TTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733213‑1733363

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