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

AAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTC  >  REL606/1733550‑1733620
                                   |                                   
aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                                     <  1:1897303/36‑1 (MQ=255)
aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                                     >  1:6919044/1‑36 (MQ=255)
aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                                     >  1:3006100/1‑36 (MQ=255)
aaGCGTTCTGACGTTATCGAAATCCGTGAGCACCTg                                     >  1:2979654/1‑36 (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:4681567‑M1/15‑1 (MQ=255)
        atcgaaatccgtgagcacctgGCACCTGAAAGCGCa                             <  1:5012000‑M1/15‑1 (MQ=255)
        atcgaaatccgtgagcacctgGCACCTGAAAGCGCa                             >  1:6730858‑M1/22‑36 (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:5895799‑M1/19‑1 (MQ=255)
            aaatccgtgagcacctgGCACCTGAAAGCGCAcggc                         <  1:3797703‑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:1249450‑M1/16‑36 (MQ=255)
              atccgtgagcacctgGCACCTGAAAGCGCAcggcgg                       >  1:964955‑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:3797703‑M2/36‑20 (MQ=255)
                   aaaTCCGTGAGCACCTGgcacctgaaagcgcacggc                  <  1:5895799‑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:964955‑M2/1‑15 (MQ=255)
                     aTCCGTGAGCACCTGgcacctgaaagcgcacggcgg                >  1:1249450‑M2/1‑15 (MQ=255)
                     aTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAAc                <  1:5350473/36‑1 (MQ=255)
                     aTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAAc                <  1:5910490/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:5350008‑M1/6‑36 (MQ=255)
                        acctgGCACCTGAAAGCGCACGGCGGCGAAAACATc             <  1:6052798‑M1/31‑1 (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:3014077‑M1/3‑36 (MQ=255)
                           tgGCACCTGAAAGCGCACGGCGGCGAAAACATCcac          >  1:6303443‑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:1841187/36‑1 (MQ=255)
                             gCACCTGAAAGCGCACGGCGGCGAAAACATCCAcat        <  1:5183931/36‑1 (MQ=255)
                              caacTGAAAGCGCACGGCGGCGAAAACATCAAcatc       >  1:4350996/4‑36 (MQ=37)
                               acctggcacctgaaagcgcacggcggcgaaaacatc      >  1:5350008‑M2/1‑5 (MQ=255)
                               acctggcacctgaaagcgcacggcggcgaaaacatc      <  1:6052798‑M2/36‑32 (MQ=255)
                               aCCTGAAAGCGCACGGCGGCGAAAACATCCAcatca      >  1:2416504/1‑36 (MQ=255)
                                 cTGAAAGCGCACGGCGGCGAAAACATCCACatcatc    <  1:5119261/36‑1 (MQ=255)
                                  tggcacctgaaagcgcacggcggcgaaaacatccac   >  1:3014077‑M2/1‑2 (MQ=255)
                                  tggcacctgaaagcgcacggcggcgaaaacatccac   >  1:6303443‑M2/1‑2 (MQ=255)
                                  tGAAAGCGCACGGCGGCGAAAACATCCACATCAtct   <  1:5258420/36‑1 (MQ=255)
                                   ggcacctgaaagcgcacggcggcgaaaacatccaca  <  1:5522938‑M2/36‑36 (MQ=255)
                                   gAAAGCGCACGGCGGCGAAAACATCCACATCAtctc  <  1:7089964/36‑1 (MQ=255)
                                   gAAAGCGCACGGCGGCGAAAACATCCACATCAtctc  <  1:2796532/36‑1 (MQ=255)
                                   gAAAGCGCACGGCGGCGAAAACATCCACATCAtctc  <  1:6234814/36‑1 (MQ=255)
                                   gAAAGCCCACGGCAGCGAAAACATCCACATCAtctc  <  1:6860799/36‑1 (MQ=37)
                                   |                                   
AAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTC  >  REL606/1733550‑1733620

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