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 222 total aligned reads displayed.

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666096‑666165
                                  |                                   
cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    <  1:5532712‑M1/36‑2 (MQ=255)
cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    >  1:6105976‑M1/1‑35 (MQ=255)
 gtgtTTAACTGTCCATCTTTTTGGGGTCAGTACAta                                   <  1:3394847‑M1/36‑3 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtt                                   >  1:5718519/1‑36 (MQ=35)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtg                                   >  1:4563594‑M1/1‑34 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtg                                   >  1:283299‑M1/1‑34 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAta                                   <  1:7288832‑M1/36‑3 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAta                                   >  1:6597808‑M1/1‑34 (MQ=255)
  tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:4497454/36‑1 (MQ=38)
  tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:5176112/36‑1 (MQ=38)
  tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:3198029/36‑1 (MQ=38)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAtacg                                 >  1:8256105‑M1/1‑32 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTc                                 <  1:864890/36‑1 (MQ=40)
    tttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                                <  1:6684485/36‑1 (MQ=255)
     ttAACTGTCCAACTTTTTGTGGTCAGTACAtacgcg                               <  1:5368482‑M2/36‑7 (MQ=255)
     ttAACTGTCCAACTTTTTGGGGTCAGTACAtacgcg                               >  1:4016436‑M2/1‑30 (MQ=255)
     ttAACTGTCCAACTTTTTGGGGTCAGTACActttct                               <  1:5623797/36‑4 (MQ=11)
     ttAACTGTCCAACTTTTTGGGGTCAGTACATTTCAg                               >  1:6057244/1‑36 (MQ=255)
     ttAACTGTCCAACTTTTTGGGGTCAGTACATTTCAg                               <  1:629383/36‑1 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACAtacgcgt                              >  1:1830227‑M2/1‑29 (MQ=255)
       aaCTGTCCAACTTTTTGGGGTCAGTACAtgttgagg                             <  1:2330569‑M1/36‑9 (MQ=255)
        aCTGTCCAACTTTTTGGGGTCAGTACAtacgcgtaa                            >  1:4532500‑M2/1‑27 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACAtgttgaggcct                          >  1:3084401‑M1/1‑25 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                          >  1:2674192/1‑36 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                          >  1:8045813/1‑36 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                          <  1:8318640/36‑1 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         <  1:5775843‑M1/36‑13 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         >  1:5799340‑M1/1‑24 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         <  1:1749360‑M1/36‑13 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAtacgcgtaaaac                         >  1:3732123‑M2/1‑24 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                         <  1:2603605/36‑1 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                         >  1:7359554/1‑36 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                         <  1:5106800/36‑1 (MQ=255)
           gTCCAACTTTTTGGGATCAATACATTTCAGGGCttt                         >  1:7781328/1‑36 (MQ=25)
             ccAACTTTTTGGGGTCAGTACAtgttgaggccttcc                       <  1:4138268‑M1/36‑15 (MQ=255)
              gaacTTTTTGGGGTCAGTACAtacgcgtaaaacagc                      <  1:7521257‑M2/35‑16 (MQ=255)
              cAACTTTTTGGGGTCAGTACAtgttgaggccttcct                      <  1:6826955‑M1/36‑16 (MQ=255)
              cAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                      <  1:1595579/36‑1 (MQ=255)
              cAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                      >  1:6716786/1‑36 (MQ=255)
               aaCTTTTTGGGGTCAGTCCATTTCAGGGCTTTTTCa                     <  1:4598308/36‑1 (MQ=255)
                 cTTTTTGGGGTCAGTACAtgttgaggccttcctggt                   <  1:1232796‑M1/36‑19 (MQ=255)
                  tttttGGGGTCAGGACAtgttgaggccttcctggtt                  >  1:1031303‑M1/1‑17 (MQ=255)
                   ttttGGGGTCAGTACAtacgagtaaaacagccagcg                 <  1:5223287‑M2/36‑21 (MQ=255)
                    tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                >  1:4801186/1‑36 (MQ=255)
                    tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                >  1:6966086/1‑36 (MQ=255)
                     ttGGGGTCAGTGCAtacgcgtaaaacagccagcgct               >  1:7406886‑M2/1‑14 (MQ=255)
                     ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc               >  1:2246241/1‑36 (MQ=255)
                      tGGGGTCAGTACAtgtttaggccttcctggttttcg              <  1:1650661‑M1/36‑24 (MQ=255)
                      tGGGGTCAGTACAtgttgaggccttcctggttttcg              <  1:3600293‑M1/36‑24 (MQ=255)
                      tGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCg              >  1:1898452/1‑36 (MQ=255)
                      tGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCg              <  1:3004766/36‑1 (MQ=255)
                       ggggTCAGTACAtgttgaggccttcctggttttcga             >  1:4390915‑M1/1‑12 (MQ=255)
                       ggggTCAGTACAtgttgaggccttcctggttttcga             <  1:6186999‑M1/36‑25 (MQ=255)
                       ggggTCAGTACATTTCAGGGCTTTTTCACGAGCCgg             >  1:1532346/1‑36 (MQ=255)
                       ggggTAAGTACAtacgcgtaaaacagccagcgctgg             >  1:8284507‑M2/1‑12 (MQ=255)
                        gggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgc            >  1:1693242/1‑36 (MQ=255)
                         ggTCAGTACAtacgcgtaaaacagccagcgctggcg           >  1:883393‑M2/1‑10 (MQ=255)
                         ggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgcg           <  1:142857/36‑1 (MQ=255)
                          gTCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc          >  1:33542/1‑36 (MQ=255)
                            cAGTACAtgttgaggccttcctggttttcgattttg        >  1:2504170‑M1/1‑7 (MQ=255)
                            cAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCtt        <  1:8001028/36‑1 (MQ=255)
                              gTACAtgttgaggccttcctggttttcgattttgga      <  1:6002454‑M1/36‑32 (MQ=255)
                               tACAtacgcgtaaaacagccagcgctggcgcgattt     <  1:1210858‑M2/36‑33 (MQ=255)
                               tACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTAt     >  1:5481076/1‑36 (MQ=255)
                               gtaatacgcgtaaaacagccagcgctggcgcgattt     >  1:1917533‑M2/4‑4 (MQ=255)
                                taatacgcgtaaaacagccagcgctggcgcgattta    >  1:7197773‑M2/3‑3 (MQ=255)
                                aCATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:3226990/36‑1 (MQ=255)
                                 cAtgttgaggccttcctggttttcgattttggagat   <  1:2273683‑M1/36‑35 (MQ=255)
                                 cAtgttgaggccttcctggttttcgattttggagat   <  1:6806048‑M1/36‑35 (MQ=255)
                                 cAtccgcgtaaaacagccagcgctggcgcgatttag   >  1:5682537‑M2/1‑2 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:3265132‑M2/1‑2 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   <  1:2160940‑M2/36‑35 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:6450833‑M2/1‑2 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:8376787‑M2/1‑2 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   <  1:1224768‑M2/36‑35 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   <  1:1565834‑M2/36‑35 (MQ=255)
                                 cAtacgcgtaaaacagccagcgctggcgcgatttag   <  1:4978556‑M2/36‑35 (MQ=255)
                                 cAtacgcataaaacagccagcgctggcgcgatttag   >  1:1137536‑M2/1‑2 (MQ=255)
                                 aatacgcgtaaaacagccagcgctggcgcgatttag   >  1:5147111‑M2/2‑2 (MQ=255)
                                  atgttgaggccttcctggttttcgattttggagatg  <  1:3777289‑M1/36‑36 (MQ=255)
                                  atacgcgtaaaacatccagcgctggcgcgatttagc  <  1:7820006‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgttttagc  <  1:1772432‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:427608‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:8415686‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:7761092‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:771670‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  >  1:7613531‑M2/1‑1 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:4871962‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:7481210‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:4960886‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  >  1:6940414‑M2/1‑1 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:6892541‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:3297419‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:6553309‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:6306379‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:6211417‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:5871999‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:2355163‑M2/36‑36 (MQ=255)
                                  atacgcgtaaaacagccagcgctggcgcgatttagc  >  1:2470337‑M2/1‑1 (MQ=255)
                                  atacgcgtaaaacagcaagcgctggcgcgatttagc  <  1:4233763‑M2/36‑36 (MQ=255)
                                  |                                   
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666096‑666165

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG