Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 98.1% 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)40 (0.750) 32/62 NT 97.6% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 1 (0.020)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)63 (1.080) 33/68 NT 98.4% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 1 (0.020)coding (685/1413 nt) pykF pyruvate kinase

Only 100 of 139 total aligned reads displayed.

GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588460‑588530
                                   |                                   
gccgtcagaggcttcgaggatttcgttgaagttgtt                                     <  1:2098290‑M1/1‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     >  1:6219761/1‑36 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     <  1:5968842/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     <  1:4543664/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     >  1:311282/1‑36 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     >  1:2398427/1‑36 (MQ=255)
 aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTg                                    <  1:4761066/36‑1 (MQ=255)
 aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTg                                    <  1:344076/36‑1 (MQ=255)
  cgtcagaggcttcgaggatttcgtcgaagttgttgt                                   <  1:7077003‑M1/3‑1 (MQ=255)
    tcagaggcttcgaggctttcgtcgaagttgtTGTAc                                 <  1:2359024‑M1/5‑1 (MQ=255)
    tcagaggcttcgaggatttcgtcgaagttgtTGTAc                                 >  1:908203‑M1/32‑36 (MQ=255)
    cATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTAc                                 >  1:6510522/1‑36 (MQ=255)
     cagaggcttcgaggatttcgtcgaagttgtTGTACt                                >  1:3803752‑M1/31‑36 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                                >  1:5285838/1‑36 (MQ=255)
      tGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTg                               <  1:3593218/36‑1 (MQ=255)
      tGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTg                               >  1:5105095/1‑36 (MQ=255)
      agaggcttcgaggatttcgtcgaagttgtTGTACTg                               >  1:5921742‑M1/30‑36 (MQ=255)
       gCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGc                              <  1:36329/36‑1 (MQ=255)
         ggcttcgaggatttcgtcggagttgtTGTACTGCAc                            <  1:1450951‑M1/10‑1 (MQ=255)
         ggcttcgaggatttcgtcgaagttgtTGTACTGCAc                            >  1:2864339‑M1/27‑36 (MQ=255)
         ggctccgaggatttcgtcgaagttgtTGTACTGCAc                            >  1:6310895‑M1/27‑36 (MQ=255)
         aGAAAACAACGGGTTCAGCTTCAAGTTGTACTGCAc                            <  1:7819860/36‑1 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            >  1:5782160/1‑36 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            >  1:4134076/1‑36 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            >  1:5623221/1‑36 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            <  1:1840949/36‑1 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           >  1:3021808/1‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           >  1:2478183/1‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           <  1:2080948/36‑1 (MQ=255)
           cttcgaggatttcgtcgaagttgtTGTACTTCAccc                          <  1:123758‑M1/12‑1 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTTCACCCa                         <  1:6265446‑M1/13‑1 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                         >  1:2745889‑M1/24‑36 (MQ=255)
             aaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAt                        <  1:4288915/36‑1 (MQ=255)
              cgaggatttcgtcgaagttgtTGTACTGCACCCAtt                       <  1:6507697‑M1/15‑1 (MQ=255)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                       >  1:6439704/1‑36 (MQ=255)
               gaggatttcgtcgaagttgtTGTACTGCACCCAttt                      <  1:5380533‑M1/16‑1 (MQ=255)
               gaggatttcgtcgaagttgtTGTACTGCACCCAttt                      <  1:3173666‑M1/16‑1 (MQ=255)
               aaaCGGGTTCAGCTTCAAGGTGTACTGCACCCAttt                      >  1:6850976/1‑36 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                     <  1:7253340/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                     <  1:286961/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                     <  1:7372386/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                     >  1:5775700/1‑36 (MQ=255)
                 ggatttcgtcgaagttgtTGTACTGCACCCATTttg                    <  1:2564286‑M1/18‑1 (MQ=255)
                  gatttcgtcgaagttgtTGTACTGCACCCATTttgt                   >  1:665546‑M1/18‑36 (MQ=255)
                  cGGGTTCAGCTTCAAGTTGTACTGCACCCATTttgt                   >  1:6827931/1‑36 (MQ=255)
                  cGGGTTCAGCTTCAAGTTGTACTGCACCCATTttgt                   <  1:3959249/36‑1 (MQ=255)
                   gggTTCAGCTTCAAGTTGTACTGGACCCATTttgtt                  >  1:2703791/1‑36 (MQ=255)
                   gggTTCAGCTTCAAGTTGTACTGCACCCATTttgtt                  >  1:391034/1‑36 (MQ=255)
                   atttcgtcgaagttgtTGTACTGCACCCATTttgtt                  >  1:4214091‑M1/17‑36 (MQ=255)
                    tttcgtcgaagttgtTGTACTGCACCCATTttgttg                 <  1:7461775‑M1/21‑1 (MQ=255)
                    ggTTCAGCTTCAAGTTGTACTGCACCCATTttgttg                 <  1:1260923/36‑1 (MQ=255)
                     ttcgtcgaagttgtTGTACTGCACCCATTTTGTTgg                <  1:3833232‑M1/22‑1 (MQ=255)
                     ttcgtcgaagttgtTGTACTGCACCCATTTTGTTgg                <  1:5823721‑M1/22‑1 (MQ=255)
                     gTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTgg                >  1:762681/1‑36 (MQ=255)
                     gTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTgg                >  1:7350749/1‑36 (MQ=255)
                      ttCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGa               >  1:7377449/1‑36 (MQ=255)
                      ttCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGa               >  1:6250609/1‑36 (MQ=255)
                      tcgtcgaagttgtTGTACTGCACCCATTTTGTTGGa               >  1:3592864‑M1/14‑36 (MQ=255)
                       tCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGAc              >  1:2271623/1‑36 (MQ=255)
                       tCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGAc              <  1:1219449/36‑1 (MQ=255)
                        gtcgaagttgtTGTACTGCACCCATTTTGTTGGACg             >  1:6665799‑M1/12‑36 (MQ=255)
                         tcgaagttgtTGTACTGCACCCATTTTGTTGGACGa            >  1:1720457‑M1/11‑36 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           >  1:1290947/1‑36 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           >  1:6642168/1‑36 (MQ=255)
                          cgaagttgtTGTACTGCACCCATTTTGTTGGACGAt           <  1:5830176‑M1/27‑1 (MQ=255)
                           gaagttgtTGTACTGCACCCATTTTGTTGGACGATg          <  1:5314379‑M1/28‑1 (MQ=255)
                           gaagttgtTGTACTGCACCCATTTTGTTGGACGATg          <  1:3330933‑M1/28‑1 (MQ=255)
                            ttCAAGTTGTACTGCACCCATTTTGTTGGACGATGa         <  1:5991877/36‑1 (MQ=255)
                            aagttgtTGTACTGCACCCATTTTGTTGGACGATGa         <  1:4953085‑M1/29‑1 (MQ=255)
                              cAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:2194758/1‑36 (MQ=40)
                              cAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:7466480/1‑36 (MQ=40)
                               ttgtTGTACTGCACCCATTTTGTTTGACGATGAAAt      >  1:1667201‑M1/5‑36 (MQ=255)
                               aaGTTGTACTGCACCCATTTTGTTGGACGATGAAAt      <  1:4862207/36‑1 (MQ=38)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:7774606‑M1/33‑1 (MQ=255)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:1006445‑M1/33‑1 (MQ=255)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:1327762‑M1/4‑36 (MQ=255)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:2086126/1‑36 (MQ=35)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:4353655/36‑1 (MQ=35)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:2903998/36‑1 (MQ=35)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg    >  1:259348‑M1/3‑36 (MQ=255)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg    <  1:2859390‑M1/34‑1 (MQ=255)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg    >  1:3218305‑M1/3‑36 (MQ=255)
                                 gTTGTACTGCACCCATTCTGTTGAACGATGAAATgg    <  1:4025793/36‑1 (MQ=11)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:924248‑M1/35‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:1012524‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:2420141‑M1/35‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:6988553‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:3756294‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:3911571‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:4002176‑M1/35‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:4151592‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:4852909‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:554621‑M1/35‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:7065093‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:7309921‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:4593760‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:1564602‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:1006864‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:5676209‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:6094769‑M1/36‑1 (MQ=255)
                                   |                                   
GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588460‑588530

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