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

GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGA  >  REL606/588460‑588529
                                   |                                  
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                    <  1:3122256/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                    <  1:2264738/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                    >  1:1679524/1‑36 (MQ=255)
 ccgtcagaggcttcgaggatttcgtcgaagctgttg                                   >  1:6022668‑M1/35‑36 (MQ=255)
 aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTg                                   >  1:1001353/1‑36 (MQ=255)
 aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTg                                   >  1:3368953/1‑36 (MQ=255)
 aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTg                                   >  1:5866364/1‑36 (MQ=255)
  gTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGt                                  >  1:2878117/1‑36 (MQ=255)
  gTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGt                                  >  1:2535278/1‑36 (MQ=255)
  cgtcagaggcttcgaggatttcgtcgaagttgttgt                                  <  1:7620501‑M1/3‑1 (MQ=255)
  cgtcagaggcttcgaggatttcgtcgaagttgttgt                                  >  1:2563‑M1/34‑36 (MQ=255)
  cgtcagaggcttcgaggatttcgtcgaagttgttgt                                  >  1:3456069‑M1/34‑36 (MQ=255)
  cgtcagaggcttcgaggattccgtcgaagttgttgt                                  >  1:4831023‑M1/34‑36 (MQ=255)
   taaTGAAGAAAAAAACGGGTTCAGCTTCAAGTTGTa                                 <  1:1124097/34‑1 (MQ=37)
   gtcagaggcttcgaggatttcgtcgaagttgtTGTa                                 >  1:7042640‑M1/33‑36 (MQ=255)
    cATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTAc                                <  1:6742456/36‑1 (MQ=255)
     cagaggcttcgaggatttcgtcgaagttgtTgtagt                               <  1:1137337‑M1/6‑3 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                               >  1:6549941/1‑36 (MQ=255)
      tGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTg                              <  1:5761315/36‑1 (MQ=255)
      agaggcttcgaggatttcgtcgaagttgtTGTACTg                              >  1:7141649‑M1/30‑36 (MQ=255)
        cAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCa                            >  1:7146160/1‑36 (MQ=255)
        cAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCa                            >  1:3351544/1‑36 (MQ=255)
        aggcttcgaggatttcgtcgaagttgtTGTACTGCa                            >  1:5568504‑M1/28‑36 (MQ=255)
         ggcttcgaggatttcgtcgaagttgtTGTACTGCAc                           >  1:5562216‑M1/27‑36 (MQ=255)
         ggcttcgaggatttcgtcgaagttgtTGTACTGCAc                           >  1:5807516‑M1/27‑36 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                           <  1:4396319/36‑1 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                           <  1:437332/36‑1 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGAAc                           >  1:4583402/1‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                          <  1:5464274/36‑1 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                        <  1:4441841‑M1/13‑1 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                        >  1:4344987‑M1/24‑36 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                        >  1:4273253‑M1/24‑36 (MQ=255)
            aaaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCa                        >  1:4918437/1‑36 (MQ=255)
            aaaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCa                        >  1:3212247/1‑36 (MQ=255)
             tcgaggatttcgtcgaagttgtTGTACTGCACCCAt                       >  1:2528471‑M1/23‑36 (MQ=255)
             aaaaaCGGGTTCAGCTTCAAGTTGCACTGCACACAt                       >  1:761857/1‑36 (MQ=37)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                      <  1:6492300/36‑1 (MQ=255)
               gaggatttcgtcgaagttgtTGTACTGCACCCAttt                     >  1:90215‑M1/21‑36 (MQ=255)
               gaggatttcgtcgaagttgtTGTAATGCACCCAttt                     <  1:519542‑M1/16‑1 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                     <  1:4425533/36‑1 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                     <  1:5184174/36‑1 (MQ=255)
                aggatttcgtcgaagttgtTGTACTGCACCCAtttt                    >  1:3806704‑M1/20‑36 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                    <  1:8072603/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                    <  1:5528051/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                    >  1:1339828/1‑36 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                    <  1:4249935/36‑1 (MQ=255)
                 ggatttcgtcgaagttgtTGTACTGCACCCATTttg                   >  1:1051747‑M1/19‑36 (MQ=255)
                  gatttcgtcgaagttgtTGTACTGCACCcatttttt                  <  1:2128551‑M1/19‑3 (MQ=255)
                   gggTTCAGCTTCAAGTTGTACTGCACCCATTttgtt                 <  1:181351/36‑1 (MQ=255)
                   gggTTCAGCTTAAAGTTGGACTGCACCAATTttgtt                 >  1:4564237/1‑36 (MQ=37)
                   atttcgtcgaagttgtTGTACTGCACCCATTttgtt                 <  1:1830335‑M1/20‑1 (MQ=255)
                    tttcgtcgaagttgtTGTACTGCACCCATTtttttg                <  1:1799255‑M1/21‑1 (MQ=255)
                    ggTTCAGCTTCAAGTTGTACTGCACCCATTttgttg                <  1:4255206/36‑1 (MQ=255)
                     ttcgtcgaagttgtTGTACTGCACCCATTTTGTTgg               >  1:5758703‑M1/15‑36 (MQ=255)
                      tcgtcgaagttgtTGTACTGCACCCATTTTGTTGGa              >  1:7695062‑M1/14‑36 (MQ=255)
                       tCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGAc             >  1:7783588/1‑36 (MQ=255)
                       cgtcgaagttgtTGTACTGCACCCATTTTGTTGGAc             <  1:5062364‑M1/24‑1 (MQ=255)
                       cgtcgaagttgtTGTACTGCACCCATTTTGTTGGAc             >  1:1862846‑M1/13‑36 (MQ=255)
                       cgtcgaagttgtTGTACTGCACCCATTTTGTTGGAc             <  1:2657975‑M1/24‑1 (MQ=255)
                       cgtcgaagttgtTGTACTGCACCCATTTTGTTGGAc             <  1:2836694‑M1/24‑1 (MQ=255)
                        gtcgaagttgtTGTACTGCACCCATTTTGTTGGACg            <  1:4506375‑M1/25‑1 (MQ=255)
                        cAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACg            <  1:3489464/36‑1 (MQ=255)
                         tcgaagttgtTGTACTGCACCCATTTTGTTGGACGa           <  1:8284533‑M1/26‑1 (MQ=255)
                         aGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGa           >  1:3143191/1‑36 (MQ=255)
                         aGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGa           >  1:3540420/1‑36 (MQ=255)
                         aGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGa           >  1:5005775/1‑36 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt          <  1:1248883/36‑1 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt          >  1:3026364/1‑36 (MQ=255)
                          cgaagttgtTGTACTGCACCCATTTTGTTGGACGAt          >  1:2178080‑M1/10‑36 (MQ=255)
                          cgaagttgtTGTACTGCACCCATTTTGTTGGACGAt          >  1:2086658‑M1/10‑36 (MQ=255)
                          cgaagttgtTGTACTGCACCCATTTTGTTGGACGAt          <  1:2745405‑M1/27‑1 (MQ=255)
                           gaagttgtTGTACTGCACCCATTTTGTTGGACGATg         <  1:5301234‑M1/28‑1 (MQ=255)
                           gaagttgtTGTACTGCACCCATTTTGTTGGACGATg         <  1:3555139‑M1/28‑1 (MQ=255)
                           gaagttgtTGTACTGCACCCATTTTGTTGGACGATg         <  1:1689617‑M1/28‑1 (MQ=255)
                           cTTCAAGTTGTACTGCACCCATTTTGTTGGACGATg         >  1:2828371/1‑36 (MQ=255)
                            cagttgtTGTACTGCACCCATTGTGTTGGACGATGa        >  1:1213416‑M1/8‑36 (MQ=255)
                            aagttgtTGTACTGCGCCCATTTTGTTGGACGATGa        >  1:393194‑M1/8‑36 (MQ=255)
                            aagttgtTGTACTGCACCCATTTTGTTGGACGATGa        >  1:4643769‑M1/8‑36 (MQ=255)
                            aagttgtTGTACTGCACCCATTTTGTTGGACGATGa        <  1:7745737‑M1/29‑1 (MQ=255)
                             tCAAGTTGTACTGCACCCATTTTGTTGGACGATGaa       >  1:4591236/1‑36 (MQ=255)
                             tCAAGTTGTACTGCACACATTTTGTTGGACGATGaa       >  1:5763557/1‑36 (MQ=255)
                             agttgtTGTACTGCACCCATTTTGTTGGACGATGaa       >  1:3739072‑M1/7‑36 (MQ=255)
                             agttgtTGTACTGCACCCATTTTGTTGGACGATGaa       <  1:8167947‑M1/30‑1 (MQ=255)
                              gttgtTGTACTGCACCCATTTTGTTGGACGATGaaa      <  1:1509913‑M1/31‑1 (MQ=255)
                               ttgtTGTACTGCACCCATTTTGTTGGACGATGAAAt     <  1:7950660‑M1/32‑1 (MQ=255)
                               ttgtTGTACTGCACCCATTTTGTTGGACGATGAAAt     >  1:6297801‑M1/5‑36 (MQ=255)
                               ttgtTGTACTGCACCCATTTTGTTGGACGATGAAAt     >  1:2428670‑M1/5‑36 (MQ=255)
                               ttgtTGTACTGCACCCATTTTGTTGGACGATGAAAt     >  1:5029481‑M1/5‑36 (MQ=255)
                               aaGTTGTACTGCACCCATTTTGTTGGACGATGAAAt     <  1:3229041/36‑1 (MQ=38)
                               aaGTTGTACTGCACCCATTTTGTTGGACGATGAAAt     >  1:5203687/1‑36 (MQ=38)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg    >  1:1476208‑M1/4‑36 (MQ=255)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg    <  1:7199417‑M1/33‑1 (MQ=255)
                                gctTGTACTGCACCCATTTTGTTGGACGATGAAATg    <  1:6403297‑M1/33‑1 (MQ=255)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg    <  1:8307745/36‑1 (MQ=35)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg   >  1:138058‑M1/3‑36 (MQ=255)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg   <  1:7610980‑M1/34‑1 (MQ=255)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg   <  1:597741‑M1/34‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa  >  1:2642965‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa  <  1:4176100‑M1/35‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa  >  1:2806009‑M1/2‑36 (MQ=255)
                                   |                                  
GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGA  >  REL606/588460‑588529

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG