Predicted mutation
evidence seq id position mutation freq annotation gene description
JC REL606 1,733,579 7 bp x 2 49.5% duplication pykF → pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1733579 =33 (0.570)27 (0.470) 22/70 NT 49.5% coding (615/1413 nt) pykF pyruvate kinase
?REL606 = 1733585 22 (0.380)coding (621/1413 nt) pykF pyruvate kinase

ATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACAT  >  REL606/1733544‑1733614
                                   |                                   
agcgttctgacgttatcgaactccgtgagcacctgg                                     <  1:4607761‑M1/1‑1 (MQ=255)
aTTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAg                                     <  1:6823146/36‑1 (MQ=255)
 ttCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGc                                    >  1:4508687/1‑36 (MQ=255)
 ttCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGc                                    >  1:6029111/1‑36 (MQ=255)
 ttCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGc                                    <  1:5026253/36‑1 (MQ=255)
 ttCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGc                                    >  1:5291911/1‑36 (MQ=255)
 ttCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGc                                    <  1:3002324/36‑1 (MQ=255)
 gcgttctgacgttatcgaaatccgtgagcacctggc                                    <  1:3185511‑M1/2‑1 (MQ=255)
  tCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCa                                   >  1:6576155/1‑36 (MQ=255)
   cGTAAGCGTTCTTACGTTATCGAAATCCGTGAGCAc                                  <  1:5605066/36‑1 (MQ=255)
   cGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCAc                                  <  1:4150443/36‑1 (MQ=255)
    ttctgacgttatcgaaatccgtgagcacctggcacc                                 >  1:5750833‑M1/32‑36 (MQ=255)
    gTAAGCGTTCTGACGTTATCGAAATCCGTGAGCAcc                                 <  1:5912441/36‑1 (MQ=255)
    gTAAGCGGTCTGACGTTATCGAAATCCGTGAGCAcc                                 >  1:4228988/1‑36 (MQ=255)
     tctgacgttatcgaaatccgtgagcacctggcacct                                <  1:3607571‑M1/6‑1 (MQ=255)
     tAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCt                                <  1:7399515/36‑1 (MQ=255)
     tAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCt                                >  1:5070919/1‑36 (MQ=255)
     tAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCt                                <  1:4691808/36‑1 (MQ=255)
     tAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCt                                <  1:3731234/36‑1 (MQ=255)
      aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                               >  1:2979654/1‑36 (MQ=255)
      aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                               >  1:3006100/1‑36 (MQ=255)
      aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                               >  1:6919044/1‑36 (MQ=255)
      aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                               <  1:1897303/36‑1 (MQ=255)
       aGCGTTCTGACGTTATCGAACTCCGTGAGCACCTGg                              <  1:4607761‑M2/36‑2 (MQ=255)
       aGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGa                              <  1:4395552/36‑1 (MQ=255)
        gacgttatcgaaatccgtgagcacctgGCACCTGaa                             <  1:933320‑M1/9‑1 (MQ=255)
        gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGgc                             <  1:3185511‑M2/36‑3 (MQ=255)
         cGTTCTGACGTTATCGAAATCCGTGAGCACCTGaaa                            <  1:590624/36‑1 (MQ=255)
          gTTCTGACGTTATCTCAATCCTTGAGCACCTTAAAg                           >  1:2870728/1‑36 (MQ=16)
          gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAg                           <  1:402241/36‑1 (MQ=255)
          cgttatcgaaatccgtgagcacctgGCACCTgaaat                           >  1:4585445‑M1/26‑35 (MQ=255)
           ttCTGACGTTATCGAAATCCGTGAGCACCTGgcacc                          >  1:5750833‑M2/1‑31 (MQ=255)
            tCTGACGTTATCGAAATCCGTGAGCACCTGgcacct                         <  1:3607571‑M2/36‑7 (MQ=255)
             tatcgaaatccgtgagcacctgGCACCTGAAAgcgc                        >  1:7467176‑M1/23‑36 (MQ=255)
              atcgaaatccgtgagcacctgGCACCTGAAAGCGCa                       <  1:5012000‑M1/15‑1 (MQ=255)
              atcgaaatccgtgagcacctgGCACCTGAAAGCGCa                       >  1:6730858‑M1/22‑36 (MQ=255)
              atcgaaatccgtgagcacctgGCACCTGAAAGCGCa                       <  1:4681567‑M1/15‑1 (MQ=255)
               tcgaaatccgtgagcacctgGCACCTGAAAGCGCAc                      >  1:5842308‑M1/21‑36 (MQ=255)
               gACGTTATCGAAATCCGTGAGCACCTGgcacctgaa                      <  1:933320‑M2/36‑10 (MQ=255)
                cgaaatccgtgagcacctgGCACCTGAAAGCGCAcg                     <  1:6030752‑M1/17‑1 (MQ=255)
                 cGTTATCGAAATCCGTGAGCACCTGgcacctgaaat                    >  1:4585445‑M2/1‑25 (MQ=255)
                  gTTATCGAAATCCGTGAGCACCTGAAAGCGCAcggc                   <  1:3633588/36‑1 (MQ=255)
                  aaatccgtgagcacctgGCACCTGAAAGCGCAcggc                   <  1:3797703‑M1/19‑1 (MQ=255)
                  aaatccgtgagcacctgGCACCTGAAAGCGCAcggc                   <  1:5895799‑M1/19‑1 (MQ=255)
                   ttATCGAAATCCGTGAGCACCTGAAAGCGCAcggcg                  >  1:5934589/1‑36 (MQ=255)
                   ttATCGAAATCCGTGAGCACCTGAAAGCGCAcggcg                  >  1:1887505/1‑36 (MQ=255)
                   aatccgtgagcacctgGCACCTGAAAGCGCAcggcg                  >  1:3544126‑M1/17‑36 (MQ=255)
                    tATCGAAATCCGTGAGCACCTGgcacctgaaagcgc                 >  1:7467176‑M2/1‑22 (MQ=255)
                    tATCGAAATCCGTGAGCACCTGAAAGCGCAcggcgg                 >  1:1086165/1‑36 (MQ=255)
                    atccgtgagcacctgGCACCTGAAAGCGCAcggcgg                 >  1:964955‑M1/16‑36 (MQ=255)
                    atccgtgagcacctgGCACCTGAAAGCGCAcggcgg                 >  1:1249450‑M1/16‑36 (MQ=255)
                     tccgtgagcacctgGCACCTGAAAGCGCACggcggc                <  1:2829465‑M1/22‑1 (MQ=255)
                     aTCGAAATCCGTGAGCACCTGgcacctgaaagcgca                >  1:6730858‑M2/1‑21 (MQ=255)
                     aTCGAAATCCGTGAGCACCTGgcacctgaaagcgca                <  1:5012000‑M2/36‑16 (MQ=255)
                     aTCGAAATCCGTGAGCACCTGgcacctgaaagcgca                <  1:4681567‑M2/36‑16 (MQ=255)
                      tCGAAATCCGTGAGCACCTGgcacctgaaagcgcac               >  1:5842308‑M2/1‑20 (MQ=255)
                      tCGAAATCCGTGAGCACCTGAAAGCGCACGgcggcg               <  1:3104272/36‑1 (MQ=255)
                      ccgtgagcacctgGCACCTGAAAGCGCACGgcggcg               <  1:196453‑M1/23‑1 (MQ=255)
                       cgtgagcacctgGCACCTGAAAGCGCACGGCGGCGa              >  1:6648471‑M1/13‑36 (MQ=255)
                       cGAAATCCGTGAGCACCTGgcacctgaaagcgcacg              <  1:6030752‑M2/36‑18 (MQ=255)
                        gtgagcacctgGCACCTGAAAGCGCACGGCGGCGaa             <  1:3710981‑M1/25‑1 (MQ=255)
                         aaaTCCGTGAGCACCTGgcacctgaaagcgcacggc            <  1:5895799‑M2/36‑20 (MQ=255)
                         aaaTCCGTGAGCACCTGgcacctgaaagcgcacggc            <  1:3797703‑M2/36‑20 (MQ=255)
                          gagcacctgGCACCTGAAAGCGCACGGCGGCGaaaa           <  1:3566314‑M1/27‑1 (MQ=255)
                          aaTCCGTGAGCACCTGgcacctgaaagcgcacggcg           >  1:3544126‑M2/1‑16 (MQ=255)
                          aaTCCGTGAGCACCTGAAAGCGCACGGCGGCGaaaa           >  1:5904119/1‑36 (MQ=255)
                           aTCCGTGAGCACCTGgcacctgaaagcgcacggcgg          >  1:1249450‑M2/1‑15 (MQ=255)
                           aTCCGTGAGCACCTGgcacctgaaagcgcacggcgg          >  1:964955‑M2/1‑15 (MQ=255)
                           aTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAAc          <  1:5910490/36‑1 (MQ=255)
                           aTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAAc          <  1:5350473/36‑1 (MQ=255)
                            tCCGTGAGCACCTGgcacctgaaagcgcacggcggc         <  1:2829465‑M2/36‑23 (MQ=255)
                             ccGTGAGCACCTGgcacctgaaagcgcacggcggcg        <  1:196453‑M2/36‑24 (MQ=255)
                              cGTGAGCACCTGgcacctgaaagcgcacggcggcga       >  1:6648471‑M2/1‑12 (MQ=255)
                              acctgGCACCTGAAAGCGCACGGCGGCGAAAACATc       <  1:6052798‑M1/31‑1 (MQ=255)
                              acctgGCACCTGAAAGCGCACGGCGGCGAAAACATc       >  1:5350008‑M1/6‑36 (MQ=255)
                               gTGAGCACCTGgcacctgaaagcgcacggcggcgaa      <  1:3710981‑M2/36‑26 (MQ=255)
                               gTGAGCACCTGAAAGCGCACGGCGGCGAAAACATcc      >  1:6399434/1‑36 (MQ=255)
                                tGAGCACCTGAAAGCGCACGGCGGCGAAAACATCca     <  1:533445/36‑1 (MQ=255)
                                 tgGCACCTGAAAGCGCACGGCGGCGAAAACATCcac    >  1:6303443‑M1/3‑36 (MQ=255)
                                 tgGCACCTGAAAGCGCACGGCGGCGAAAACATCcac    >  1:3014077‑M1/3‑36 (MQ=255)
                                 gAGCACCTGgcacctgaaagcgcacggcggcgaaaa    <  1:3566314‑M2/36‑28 (MQ=255)
                                 gAGCACCTGAAAGCGCACGGCGGCGAAAACATCcac    >  1:5591891/1‑36 (MQ=255)
                                  gGCACCTGAAAGCGCACGGCGGCGAAAACATCcaca   <  1:5522938‑M1/35‑1 (MQ=255)
                                   gCACCTGAAAGCGCACGGCGGCGAAAACATCCAcat  <  1:6219446/36‑1 (MQ=255)
                                   gCACCTGAAAGCGCACGGCGGCGAAAACATCCAcat  <  1:5183931/36‑1 (MQ=255)
                                   gCACCTGAAAGCGCACGGCGGCGAAAACATCCAcat  <  1:1841187/36‑1 (MQ=255)
                                   |                                   
ATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACAT  >  REL606/1733544‑1733614

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG