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

Only 100 of 141 total aligned reads displayed.

GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588460‑588530
                                   |                                   
gccgtcagaggcttcgaggatttcgtcgaagttgtt                                     >  1:1312691‑M1/36‑36 (MQ=255)
gccgtcagaggcttcgaggatttcgtcgaagttgtt                                     >  1:2936837‑M1/36‑36 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     >  1:7811613/1‑36 (MQ=255)
  gTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGt                                   <  1:2532885/36‑1 (MQ=255)
  gTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGt                                   <  1:4494747/36‑1 (MQ=255)
   tCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTa                                  <  1:6013764/36‑1 (MQ=255)
   tCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTa                                  >  1:6641137/1‑36 (MQ=255)
   tCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTa                                  <  1:5602776/36‑1 (MQ=255)
   gtcagaggcttcgaggatttcgtcgaagttgtTGTa                                  >  1:7458637‑M1/33‑36 (MQ=255)
    cATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTAc                                 <  1:5250367/36‑1 (MQ=255)
     cagaggcttcgaggatttcgtcgaagttgtTGTACt                                <  1:3965450‑M1/6‑1 (MQ=255)
     cagaggcttcgaggatttcgtcgaagttgtTGTACt                                <  1:7306423‑M1/6‑1 (MQ=255)
     cagaggcttcgaggatttcgtcgaagttgtTGTACt                                >  1:7490265‑M1/31‑36 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                                >  1:1660645/1‑36 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                                >  1:2870172/1‑36 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                                >  1:2771813/1‑36 (MQ=255)
      tgaaGAAAAAAACGGGTTCAGCTTCAAGTTGTACTg                               <  1:2910264/33‑1 (MQ=255)
      tGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTg                               >  1:3384107/1‑36 (MQ=255)
      agaggcttcgaggatttcgtcgaagttgtTTTACTg                               <  1:361535‑M1/7‑1 (MQ=255)
       gaggcttcgaggatttcgtcgaagttgtTGTACTGc                              <  1:1383247‑M1/8‑1 (MQ=255)
       gCAGAAACAAACGGGTTCAGCTTCAAGTTGTACTGc                              <  1:1858768/36‑1 (MQ=255)
        aggcttcgaggatttcgtcgcagttgtTGTACTgct                             <  1:2460156‑M1/9‑2 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            >  1:2220260/1‑36 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            >  1:7096075/1‑36 (MQ=255)
          gcttcgaggatttcgtcgaagttgtTGTACTGCAcc                           <  1:3194686‑M1/11‑1 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           <  1:7344767/36‑1 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           >  1:4506025/1‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           >  1:4429936/1‑36 (MQ=255)
            ttcgaggatttcgtcgaagttgtTTTACTGCACCCa                         >  1:4560388‑M1/24‑36 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                         <  1:487928‑M1/13‑1 (MQ=255)
             tcgaggatttcgtcgaagttgtTGTACTGCACCCAt                        <  1:2260805‑M1/14‑1 (MQ=255)
             aaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAt                        >  1:6933560/1‑36 (MQ=255)
             aaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAt                        >  1:7102607/1‑36 (MQ=255)
             aaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAt                        <  1:4495287/36‑1 (MQ=255)
             aaaaaCGGGTGCAGCTTCAAGTTGTACTGAACCCAt                        >  1:5118988/1‑36 (MQ=37)
              cgaggatttcgtcgaagttgtTGTACTGCACCCAtt                       >  1:5815710‑M1/22‑36 (MQ=255)
              cgaggatttcgtcgaagttgtTGTACTGCACCCAtt                       >  1:5073817‑M1/22‑36 (MQ=255)
              aaaaTGGGTTCAGCTGCAAGTTGTTCTGCACCCAtt                       <  1:2090972/36‑1 (MQ=21)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                       >  1:6052882/1‑36 (MQ=255)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                       <  1:6303157/36‑1 (MQ=255)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                       >  1:1901672/1‑36 (MQ=255)
               gaggatttcgtcgaagttgtTGTACTGCACCCAttt                      >  1:4003030‑M1/21‑36 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                      >  1:5989123/1‑36 (MQ=255)
                aggatttcggcgaagttgtTGTACTGCACCCAtttt                     <  1:7797036‑M1/17‑1 (MQ=255)
                agaatttcgtcgaagttgtTGTACTGCACCCAtttt                     >  1:3501915‑M1/20‑36 (MQ=255)
                aaCGGGTTCAGTTTCAAGTTGTACTGCACCCAtttt                     <  1:6227584/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                     <  1:4244972/36‑1 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAgttt                     >  1:3903866/1‑36 (MQ=255)
                 ggatttcgtcgaagttgtTGTACTGCACCCATTttg                    >  1:5386949‑M1/19‑36 (MQ=255)
                  gatttcgtcgaagttgtTGTACTGCACCCATTttgt                   <  1:447752‑M1/19‑1 (MQ=255)
                  gatttcgtcgaagttgtTGTACTGCACCCATTttgt                   <  1:6075137‑M1/19‑1 (MQ=255)
                  cGGGTTCAGCTTCAAGTTGTACTGCACCCATTttgt                   <  1:7861470/36‑1 (MQ=255)
                   atttcgtcgaagttgtTGTACTGCACCCATTttgtt                  <  1:2227349‑M1/20‑1 (MQ=255)
                    ggTTCAGCTTCAAGTTGTACTGCACCCATTttgttg                 <  1:5549560/36‑1 (MQ=255)
                    ggTTCAGCTTCAAGTTGTACTGCACCCATTttgttg                 >  1:2640310/1‑36 (MQ=255)
                     ttcgtcgaagttgtTGTACTGCACCCATTTTGTTgg                <  1:1788330‑M1/22‑1 (MQ=255)
                     gTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTgg                >  1:1694395/1‑36 (MQ=255)
                      ttCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGa               >  1:2341806/1‑36 (MQ=255)
                      ttCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGa               <  1:4716742/36‑1 (MQ=255)
                       tCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGAc              <  1:7875615/36‑1 (MQ=255)
                        cAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACg             >  1:2450371/1‑36 (MQ=255)
                        cAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACg             >  1:6020062/1‑36 (MQ=255)
                         tcgaagttgtTGTACTGCACCCATTTTGTTGGACGa            >  1:6824615‑M1/11‑36 (MQ=255)
                         tcgaagttgtTGTACTGCACCCATTTTGTTGGACGa            >  1:6679846‑M1/11‑36 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           >  1:3038513/1‑36 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           <  1:4150235/36‑1 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           <  1:6475542/36‑1 (MQ=255)
                          cgaagttgtTGTACTGCCCCCATTTTGTTGGACGAt           >  1:5236970‑M1/10‑36 (MQ=255)
                           gaagttgtTGTACTGCACCCATTTTGTTGGACGATg          <  1:2807239‑M1/28‑1 (MQ=255)
                           cTTCAAGTTGTACTGCACCCATTTTGTTGGACGATg          >  1:6839798/1‑36 (MQ=255)
                            aagttgtTGTACTGCACCCATTTTGTTGGACGATGa         >  1:1420488‑M1/8‑36 (MQ=255)
                             tCAAGTTGTACTGCACCCATTTTGTTGGACGATGaa        >  1:565048/1‑36 (MQ=255)
                             tCAAGTTGTACTGCACCCATTTTGTTGGACGATGaa        >  1:4118065/1‑36 (MQ=255)
                             agttgtTGTACTGCACCCATTTTGTTGGACGATGaa        <  1:4918491‑M1/30‑1 (MQ=255)
                             agttgtTGTACTGCACCCATTTTGTTGGACGATGaa        <  1:3901854‑M1/30‑1 (MQ=255)
                              gttgtTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:3847570‑M1/6‑36 (MQ=255)
                              gttgtTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:3767338‑M1/6‑36 (MQ=255)
                              gttgtTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:7620310‑M1/6‑36 (MQ=255)
                              cAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa       <  1:3200978/36‑1 (MQ=40)
                               aggtTGTACTGCACCCATTTTGTTGGACGATGAAAt      <  1:6548326‑M1/32‑1 (MQ=255)
                               aaGTTGTACTGCACCCATTTTGTTGGACGATGAAAt      <  1:3258885/36‑1 (MQ=38)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:4125444‑M1/33‑1 (MQ=255)
                                ggtTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:4613693‑M1/33‑1 (MQ=255)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:5727353/1‑36 (MQ=35)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:3826503/1‑36 (MQ=35)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg    <  1:5509034‑M1/34‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:3629451‑M1/35‑1 (MQ=255)
                                  tTGTACTGCACCCACTTTGTTGGACGATGAAATGGa   <  1:5767734‑M1/35‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGGCGATGAAATGGaa  >  1:6204572‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:6497183‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:640427‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:3106463‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:2760052‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:2659308‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:2590197‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:2573330‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:4621597‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:1238700‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:1168100‑M1/36‑1 (MQ=255)
                                   TGTACTGCAACCATTTTGTTGGACGATGAAATGGaa  >  1:6679013‑M1/1‑36 (MQ=255)
                                   |                                   
GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588460‑588530

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