Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 100% coding (683‑685/1413 nt) pykF → pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 588495 =NA (NA)60 (1.100) 35/62 NT 100% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 0 (0.030)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)43 (0.720) 31/68 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 0 (0.000)coding (685/1413 nt) pykF pyruvate kinase

Only 100 of 114 total aligned reads displayed.

CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682
                                   |                                   
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:2642965‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:2806009‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:4176100‑M2/1‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:597741‑M2/2‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    >  1:138058‑M2/35‑36 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:7610980‑M2/2‑1 (MQ=255)
  tcatcTCCAAAATCGAAAACCAGGAAGGCCTcacca                                   <  1:1967239/36‑1 (MQ=255)
  tcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACA                                   >  1:5380310‑M2/1‑36 (MQ=255)
  tcataTACAAAATCGAAAACCAGGAAGGCCTcaaca                                   <  1:2555852/36‑1 (MQ=37)
  catttcatcgtccaacaaaatgggtgcagtacaagc                                   <  1:6403297‑M2/3‑3 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   <  1:7199417‑M2/3‑1 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                  <  1:3777289‑M2/36‑2 (MQ=255)
   catcTCAAAAATCGAAAACCAGGAAGGCCTCAACAt                                  >  1:3205738‑M2/1‑35 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  >  1:6297801‑M2/33‑36 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  <  1:7950660‑M2/4‑1 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  >  1:5029481‑M2/33‑36 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  >  1:2428670‑M2/33‑36 (MQ=255)
    tttcatcgtccaacaaaatgggtgcagtacaacaac                                 <  1:1509913‑M2/5‑1 (MQ=255)
    atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                 <  1:1747236‑M2/36‑3 (MQ=255)
    atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                 <  1:6550055‑M2/36‑3 (MQ=255)
    atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                 <  1:6806048‑M2/36‑3 (MQ=255)
    atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                 >  1:4906455‑M2/1‑34 (MQ=255)
    atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                 <  1:2273683‑M2/36‑3 (MQ=255)
     ttcatcgtccaacaaaatgggtgcagtacaACAACt                                >  1:3739072‑M2/31‑36 (MQ=255)
     ttcatcgtccaacaaaatgggtgcagtacaACAACt                                <  1:8167947‑M2/6‑1 (MQ=255)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                >  1:8330184‑M2/1‑33 (MQ=255)
      tcatcgtccaacacaatgggtgcagtacaACAActg                               >  1:1213416‑M2/30‑35 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               >  1:4643769‑M2/30‑36 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               <  1:2633944‑M2/7‑1 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               <  1:7745737‑M2/7‑1 (MQ=255)
      tcatcgtccaacaaaatgggcgcagtacaACAACtt                               >  1:393194‑M2/30‑36 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              <  1:6002454‑M2/36‑6 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:5301234‑M2/8‑1 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:5710148‑M2/8‑1 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:3555139‑M2/8‑1 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:1689617‑M2/8‑1 (MQ=255)
        ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                             <  1:5149931‑M2/36‑7 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             <  1:7541698‑M2/9‑1 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             >  1:6990722‑M2/28‑36 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             >  1:2086658‑M2/28‑36 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             <  1:2745405‑M2/9‑1 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             >  1:1935021‑M2/28‑36 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaACAACTTcga                            <  1:8284533‑M2/10‑1 (MQ=255)
         cAAAATCGAAAACCAGGAAGGCCTCAACAtgtactg                            >  1:2504170‑M2/1‑29 (MQ=255)
          cgtccaacaaaatgggtgcagtacaACAACTTcgac                           <  1:4506375‑M2/11‑1 (MQ=255)
           gtccaacaaaatgggtgcagtacaACAACTTcgacg                          <  1:2836694‑M2/12‑1 (MQ=255)
           gtccaacaaaatgggtgcagtacaACAACTTcgacg                          >  1:1862846‑M2/25‑36 (MQ=255)
           gtccaacaaaatgggtgcagtacaACAACTTcgacg                          <  1:5062364‑M2/12‑1 (MQ=255)
           gtccaacaaaatgggtgcagtacaACAACTTcgacg                          <  1:2657975‑M2/12‑1 (MQ=255)
           aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                          >  1:5802670‑M2/1‑27 (MQ=255)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        >  1:5758703‑M2/23‑36 (MQ=255)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        >  1:7723644‑M2/23‑36 (MQ=255)
              tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                       >  1:4390915‑M2/1‑24 (MQ=255)
              tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                       <  1:6186999‑M2/36‑13 (MQ=255)
              caacaaaatgggtgcagtacaACAACTTCGACGaaa                       >  1:5166251‑M2/22‑36 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      >  1:2998714‑M2/1‑23 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      >  1:4551023‑M2/1‑23 (MQ=255)
               cGAAAACCAGGAAGGCCTAAACAtgtactgacccca                      <  1:1650661‑M2/36‑14 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      <  1:1830335‑M2/16‑1 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                     <  1:4855464‑M2/36‑15 (MQ=255)
                aaaaaatgggtgcagtacaACAACTTCGACGAAATc                     <  1:2128551‑M2/17‑1 (MQ=255)
                 caaaatgggtgcagtacaACAACTTCGACGAAATcc                    >  1:3520749‑M2/19‑36 (MQ=255)
                 caaaatgggtgcagtacaACAACTTCGACGAAATcc                    >  1:1051747‑M2/19‑36 (MQ=255)
                 aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                    <  1:6302874‑M2/36‑16 (MQ=255)
                  aaaatgggtgcagtacaACAACTTCGACGAAATCCt                   >  1:3806704‑M2/18‑36 (MQ=255)
                  aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                   <  1:168482‑M2/36‑17 (MQ=255)
                  aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                   <  1:1916619‑M2/36‑17 (MQ=255)
                   aaCCAGGAAGGCCTCAACAtgtcctgaccccaaaaa                  >  1:1031303‑M2/1‑19 (MQ=255)
                   aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                  >  1:3159756‑M2/1‑19 (MQ=255)
                    aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                 <  1:1232796‑M2/36‑19 (MQ=255)
                    aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                 >  1:3945251‑M2/1‑18 (MQ=255)
                     ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                <  1:1375605‑M2/36‑20 (MQ=255)
                     ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                >  1:3298903‑M2/1‑17 (MQ=255)
                     atgggtgcagtacaACAACTTCGACGAAATCCTCGa                >  1:4414119‑M2/15‑36 (MQ=255)
                     atgggtgcagtacaACAACTTCGACGAAATCCTCGa                >  1:2528471‑M2/15‑36 (MQ=255)
                      tgggtgcagtacaACAACTTCGACGAAATCCTCGaa               >  1:4344987‑M2/14‑36 (MQ=255)
                      tgggtgcagtacaACAACTTCGACGAAATCCTCGaa               >  1:4273253‑M2/14‑36 (MQ=255)
                      tgggtgcagtacaACAACTTCGACGAAATCCTCGaa               <  1:4441841‑M2/23‑1 (MQ=255)
                      tgggtgcagtacaACAACTCCGACGAAATCCTCGaa               <  1:2392110‑M2/23‑1 (MQ=255)
                      cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt               >  1:5191603‑M2/1‑16 (MQ=255)
                        ggAAGGCCTCAACAtgtactgaccccaaaaagttgg             <  1:527267‑M2/36‑23 (MQ=255)
                        ggAAGGCCTCAACAtgtactgaccccaaaaagttgg             <  1:4138268‑M2/36‑23 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc            >  1:5562216‑M2/11‑36 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc            >  1:5807516‑M2/11‑36 (MQ=255)
                          tgcagtacaACAACTTCGACGAAATCCTCGAAGCct           >  1:5568504‑M2/10‑36 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           >  1:5799340‑M2/1‑12 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           <  1:5775843‑M2/36‑25 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           <  1:1749360‑M2/36‑25 (MQ=255)
                           gcagtacaACAACTTCGACGAAATCCTCGAAGCctc          <  1:3277748‑M2/28‑1 (MQ=255)
                            cagtacaACAACTTCGACGAAATCCTCGAAGCctct         >  1:7141649‑M2/8‑36 (MQ=255)
                             actacaACAACTTCGACGAAATCCTCGAAGCCTCTg        <  1:1137337‑M2/30‑1 (MQ=255)
                              ccTCAACAtgtactgaccccaaaaagttggacagtt       <  1:1559359‑M2/36‑29 (MQ=255)
                              ccTCAACAtgtactgaccccaaaaagttggacagtt       <  1:2330569‑M2/36‑29 (MQ=255)
                               tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc      >  1:7042640‑M2/5‑36 (MQ=255)
                                acaACAACTTCGACGAAATCCTCGAAGCCTCTGACg     <  1:7620501‑M2/33‑1 (MQ=255)
                                acaACAACTTCGACGAAATCCTCGAAGCCTCTGACg     >  1:6377352‑M2/4‑36 (MQ=255)
                                acaACAACTTCGACGAAATCCTCGAAGCCTCTGACg     >  1:3456069‑M2/4‑36 (MQ=255)
                                 caACAGCTTCGACGAAATCCTCGAAGCCTCTGACgg    >  1:6022668‑M2/3‑36 (MQ=255)
                                   ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa  >  1:7731320‑M2/1‑36 (MQ=255)
                                   ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa  >  1:1407272‑M2/1‑36 (MQ=255)
                                   |                                   
CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG