Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 96.3% 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.000)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)43 (0.720) 31/68 NT 91.5% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 4 (0.060)coding (685/1413 nt) pykF pyruvate kinase

Only 100 of 126 total aligned reads displayed.

TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684
                                   |                                   
tcatcTCCAAAATCGAAAACCAGGAAGGCCTcacca                                     <  1:1967239/36‑1 (MQ=255)
tcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACA                                     >  1:5380310‑M2/1‑36 (MQ=255)
tcataTACAAAATCGAAAACCAGGAAGGCCTcaaca                                     <  1:2555852/36‑1 (MQ=37)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     >  1:1476208‑M2/34‑36 (MQ=255)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     <  1:7199417‑M2/3‑1 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                    <  1:3777289‑M2/36‑2 (MQ=255)
 atttcatcgtccaacaaaatgggtgcagtacaacaa                                    <  1:7950660‑M2/4‑1 (MQ=255)
 atttcatcgtccaacaaaatgggtgcagtacaacaa                                    >  1:5029481‑M2/33‑36 (MQ=255)
 atttcatcgtccaacaaaatgggtgcagtacaacaa                                    >  1:6297801‑M2/33‑36 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   >  1:4906455‑M2/1‑34 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   <  1:6550055‑M2/36‑3 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   <  1:6806048‑M2/36‑3 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   <  1:2273683‑M2/36‑3 (MQ=255)
   tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                  >  1:8330184‑M2/1‑33 (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:5710148‑M2/8‑1 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:5301234‑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:6990722‑M2/28‑36 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               >  1:2178080‑M2/28‑36 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               >  1:2086658‑M2/28‑36 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               <  1:2745405‑M2/9‑1 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               <  1:7541698‑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:5062364‑M2/12‑1 (MQ=255)
         gtccaacaaaatgggtgcagtacaACAACTTcgacg                            <  1:2657975‑M2/12‑1 (MQ=255)
         aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                            >  1:5802670‑M2/1‑27 (MQ=255)
          aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                           <  1:4255917‑M2/36‑11 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          >  1:5758703‑M2/23‑36 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          >  1:7723644‑M2/23‑36 (MQ=255)
            tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                         <  1:6186999‑M2/36‑13 (MQ=255)
            caacaaaatgggtgcagtacaACAACTTCGACGaaa                         >  1:5166251‑M2/22‑36 (MQ=255)
            caaaaaaatgggtgcagtacaACAACTTCGACGaaa                         <  1:1799255‑M2/15‑1 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        >  1:2998714‑M2/1‑23 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        <  1:3600293‑M2/36‑14 (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)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      >  1:1051747‑M2/19‑36 (MQ=255)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      >  1:3520749‑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)
                 aaatgggtgcattacaACAACTTCGACGAAATCCTc                    <  1:519542‑M2/20‑1 (MQ=255)
                 aaCCAGGAAGGCCTCAACAtgtcctgaccccaaaaa                    >  1:1031303‑M2/1‑19 (MQ=255)
                 aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                    >  1:3159756‑M2/1‑19 (MQ=255)
                  aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                   >  1:3945251‑M2/1‑18 (MQ=255)
                  aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                   <  1:1232796‑M2/36‑19 (MQ=255)
                   ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                  >  1:3298903‑M2/1‑17 (MQ=255)
                   ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                  <  1:1375605‑M2/36‑20 (MQ=255)
                   atgggtgcagtacaACAACTTCGACGAAATCCTCGa                  >  1:2528471‑M2/15‑36 (MQ=255)
                   atgggtgcagtacaACAACTTCGACGAAATCCTCGa                  >  1:4414119‑M2/15‑36 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  1:4441841‑M2/23‑1 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 >  1:4273253‑M2/14‑36 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 >  1:4344987‑M2/14‑36 (MQ=255)
                    tgggtgcagtacaACAACTCCGACGAAATCCTCGaa                 <  1:2392110‑M2/23‑1 (MQ=255)
                      ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:527267‑M2/36‑23 (MQ=255)
                       gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              >  1:5562216‑M2/11‑36 (MQ=255)
                        tgcagtacaACAACTTCGACGAAATCCTCGAAGCct             >  1:5568504‑M2/10‑36 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             <  1:5775843‑M2/36‑25 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             >  1:5799340‑M2/1‑12 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             <  1:1749360‑M2/36‑25 (MQ=255)
                         gcagtacaACAACTTCGACGAAATCCTCGAAGCctc            <  1:3277748‑M2/28‑1 (MQ=255)
                         aGGCCTCAACAtgtactgaccccaaaaagttggaca            >  1:3084401‑M2/1‑11 (MQ=255)
                          cagtacaACAACTTCGACGAAATCCTCGAAGCctct           >  1:7141649‑M2/8‑36 (MQ=255)
                           actacaACAACTTCGACGAAATCCTCGAAGCCTCTg          <  1:1137337‑M2/30‑1 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggacagtt         <  1:2330569‑M2/36‑29 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggacagtt         <  1:1559359‑M2/36‑29 (MQ=255)
                             tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc        >  1:7042640‑M2/5‑36 (MQ=255)
                              acaACAACTTCGACGGAATCCTCGAAGCCTCTGACg       >  1:4831023‑M2/4‑36 (MQ=255)
                              acaACAACTTCGACGAAATCCTCGAAGCCTCTGACg       <  1:7620501‑M2/33‑1 (MQ=255)
                               caACAGCTTCGACGAAATCCTCGAAGCCTCTGACgg      >  1:6022668‑M2/3‑36 (MQ=255)
                                cgtatgtactgaccccaaaaagttggacagttaaac     >  1:8256105‑M2/4‑4 (MQ=255)
                                 gtatgtactgaccccaaaaagttggacagttaaaca    <  1:3301287‑M2/34‑34 (MQ=255)
                                 ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa    >  1:7731320‑M2/1‑36 (MQ=255)
                                 ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa    >  1:1407272‑M2/1‑36 (MQ=255)
                                  tatgtactgaccccaaaaagttggacagttaaacac   >  1:6597808‑M2/2‑2 (MQ=255)
                                  tatgtactgaccccaaaaagttggacagttaaacac   >  1:1783122‑M2/2‑2 (MQ=255)
                                  tatgtactgaccccaaaaagttggacagttaaacac   <  1:7288832‑M2/35‑35 (MQ=255)
                                  tatgtactgaccccaaaaagatggacagttaaacac   <  1:3394847‑M2/35‑35 (MQ=255)
                                  cAtgtactgaccccaaaaagttggacagttaaacac   >  1:283299‑M2/1‑2 (MQ=255)
                                  CAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAt   <  1:5882214‑M2/36‑1 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  <  1:5791152‑M2/36‑36 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  <  1:3037417‑M2/36‑36 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  >  1:6105976‑M2/1‑1 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  <  1:1084075‑M2/36‑36 (MQ=255)
                                   aacTTCGACTAAATCCTCGAAGCCTCTGACGGcatc  <  1:7510979/36‑1 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  >  1:3462464/1‑36 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  <  1:6407300/36‑1 (MQ=255)
                                   |                                   
TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684

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