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 =6 (0.150)40 (1.040) 34/146 0.2 90.0% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 3 (0.080)coding (331/1413 nt) pykF pyruvate kinase

TCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733214‑1733363
                                                                          |                                                                           
tctcTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTaa                                                                            <  2:660211/76‑1 (MQ=255)
agctggtcagacctttactttcaccactgataaagctgttatcggcaacagcgaaatggttgcggtaacgtatgaa                                                                            >  1:494540‑M1/75‑76 (MQ=255)
 ctctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAAc                                                                           <  2:165698/76‑1 (MQ=255)
  tctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACg                                                                          <  2:1217850/76‑1 (MQ=255)
  ctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacg                                                                          <  1:51434‑M1/4‑1 (MQ=255)
   ctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGt                                                                         >  1:651261/1‑76 (MQ=255)
    tGNAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCCGTAACNTa                                                                        <  2:1032348/76‑1 (MQ=255)
     gtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtat                                                                       <  2:490970‑M1/7‑1 (MQ=255)
      aaaGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATg                                                                      >  2:867814/1‑76 (MQ=255)
       cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatgtttgcgttaacgcatcaacgtatga                                                                     <  1:12636‑M1/9‑1 (MQ=255)
        aGCTGGTCAGACCTTTACTTTCACCACTGATAAAGCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                                                                    >  1:494540‑M2/1‑74 (MQ=255)
         gacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAg                                                                   >  2:901904‑M1/66‑76 (MQ=255)
         gacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAg                                                                   >  2:890841‑M1/66‑76 (MQ=255)
          cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                  <  1:51434‑M2/76‑5 (MQ=255)
             tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                               <  2:777150‑M1/15‑1 (MQ=255)
             tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                               >  2:1205040‑M1/62‑76 (MQ=255)
             gTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                                                               <  2:490970‑M2/76‑8 (MQ=255)
               cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGTTTGCGTTAACGCATCaacgtatga                                                             <  1:12636‑M2/76‑10 (MQ=255)
                actttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCac                                                            >  2:1214682‑M1/59‑76 (MQ=255)
                actttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCac                                                            <  2:971141‑M1/18‑1 (MQ=255)
                 gACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                                                           >  2:890841‑M2/1‑65 (MQ=255)
                 gACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                                                           >  2:901904‑M2/1‑65 (MQ=255)
                 ctttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCact                                                           <  1:270410‑M1/19‑1 (MQ=255)
                  tttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCacta                                                          >  1:825804‑M1/57‑76 (MQ=255)
                  tttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCacta                                                          >  1:744797‑M1/57‑76 (MQ=255)
                    tcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactact                                                        <  1:540158‑M1/22‑1 (MQ=255)
                    tcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactact                                                        <  2:500687‑M1/22‑1 (MQ=255)
                     tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                       >  2:1205040‑M2/1‑61 (MQ=255)
                     tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                       <  2:777150‑M2/76‑16 (MQ=255)
                      accactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGa                                                      >  1:543656‑M1/53‑76 (MQ=255)
                        cactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAcc                                                    >  1:1069047‑M1/51‑76 (MQ=255)
                        aCTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcac                                                    <  2:971141‑M2/76‑19 (MQ=255)
                        aCTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcac                                                    >  2:1214682‑M2/1‑58 (MQ=255)
                         cTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcact                                                   <  1:270410‑M2/76‑20 (MQ=255)
                         actgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACct                                                   <  1:1199041‑M1/27‑1 (MQ=255)
                          tttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcacta                                                  >  1:825804‑M2/1‑56 (MQ=255)
                          tttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcacta                                                  >  1:744797‑M2/1‑56 (MQ=255)
                            tCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactact                                                <  1:540158‑M2/76‑23 (MQ=255)
                            tCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactact                                                <  2:500687‑M2/76‑23 (MQ=255)
                              taaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctg                                              >  2:152289‑M1/45‑76 (MQ=255)
                              accacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactga                                              >  1:543656‑M2/1‑52 (MQ=255)
                               aaatctgttatcggcaacagagaagtggttgcggtaacgtatgAACGTATGAACGTCTGACGACTGACctgtctgt                                             <  1:58323‑M1/33‑1 (MQ=255)
                                cacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacc                                            >  1:1069047‑M2/1‑50 (MQ=255)
                                 atctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTg                                           <  1:476734‑M1/35‑1 (MQ=255)
                                 acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacct                                           <  1:1199041‑M2/76‑28 (MQ=255)
                                    tgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCa                                        >  2:924774‑M1/39‑76 (MQ=255)
                                      ttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAac                                      <  1:902929‑M1/40‑1 (MQ=255)
                                      tAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctg                                      >  2:152289‑M2/1‑44 (MQ=255)
                                       aaaTCTGTTATCGGCAACAGAGAAGTGGTTGCGGTAACGTATGaacgtatgaacgtctgacgactgacctgtctgt                                     <  1:58323‑M2/76‑34 (MQ=255)
                                         tcggcaacngcgtaatggttgcggtaaggtatgGGCGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACAcc                                   >  2:201097‑M1/34‑76 (MQ=255)
                                         aTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttg                                   <  1:476734‑M2/76‑36 (MQ=255)
                                          cggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCg                                  >  2:141460‑M1/33‑76 (MQ=255)
                                            tGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggca                                >  2:924774‑M2/1‑38 (MQ=255)
                                             caacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTAc                               <  2:913414‑M1/47‑1 (MQ=255)
                                              ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaac                              <  1:902929‑M2/76‑41 (MQ=255)
                                              aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                              <  1:160335‑M1/48‑1 (MQ=255)
                                              aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                              <  1:454063‑M1/48‑1 (MQ=255)
                                              aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                              <  1:588035‑M1/48‑1 (MQ=255)
                                                 tCGGCAACNGCGTAATGGTTGCGGTAAGGTATGggcgtatgaaggtttcactactgacctgtctgttggcaacacc                           >  2:201097‑M2/1‑33 (MQ=255)
                                                  cGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccg                          >  2:141460‑M2/1‑32 (MQ=255)
                                                   cgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTg                         <  2:848068‑M1/53‑1 (MQ=255)
                                                     cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtac                       <  2:913414‑M2/76‑48 (MQ=255)
                                                      aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                      <  1:588035‑M2/76‑49 (MQ=255)
                                                      aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                      <  1:454063‑M2/76‑49 (MQ=255)
                                                      aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                      <  1:160335‑M2/76‑49 (MQ=255)
                                                       atggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGa                     <  2:899192‑M1/57‑1 (MQ=255)
                                                          gttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATgg                  <  2:397259‑M1/60‑1 (MQ=255)
                                                          gttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATgg                  <  1:235162‑M1/60‑1 (MQ=255)
                                                           cGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttg                 <  2:848068‑M2/76‑54 (MQ=255)
                                                             gcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCt               >  1:1016349‑M1/14‑76 (MQ=255)
                                                              cggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTg              >  2:352094‑M1/13‑76 (MQ=255)
                                                               ggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGa             >  1:852970‑M1/12‑76 (MQ=255)
                                                               aTGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacga             <  2:899192‑M2/76‑58 (MQ=255)
                                                                  gTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatgg          <  1:235162‑M2/76‑61 (MQ=255)
                                                                  gTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatgg          <  2:397259‑M2/76‑61 (MQ=255)
                                                                   acgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCgg         >  2:1162392‑M1/8‑76 (MQ=255)
                                                                     gCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtct       >  1:1016349‑M2/1‑13 (MQ=255)
                                                                      tatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTAt      >  2:651260‑M1/5‑76 (MQ=255)
                                                                      cGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctg      >  2:352094‑M2/1‑12 (MQ=255)
                                                                       ggTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctga     >  1:852970‑M2/1‑11 (MQ=255)
                                                                       atgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATg     >  2:616941‑M1/4‑76 (MQ=255)
                                                                       atgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATg     <  1:878200‑M1/73‑1 (MQ=255)
                                                                         gAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGa   >  1:291435‑M1/2‑76 (MQ=255)
                                                                          AACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGaa  >  2:303949‑M1/1‑76 (MQ=255)
                                                                          |                                                                           
TCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733214‑1733363

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