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

Only 100 of 117 total aligned reads displayed.

TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684
                                   |                                   
tcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACA                                     >  1:1332448‑M2/1‑36 (MQ=255)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     <  1:4125444‑M2/3‑1 (MQ=255)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     >  1:4465788‑M2/34‑36 (MQ=255)
 atttcatcgtccaacaaaatgggtgcagtacaacaa                                    <  1:5861480‑M2/4‑1 (MQ=255)
  tttcatcgtccaacaaaatgggtgcagtacaacaac                                   >  1:3767338‑M2/32‑36 (MQ=255)
  tttcatcgtccaacaaaatgggtgcagtacaacaac                                   >  1:7620310‑M2/32‑36 (MQ=255)
  tttcatcgtccaacaaaatgggtgcagtacaacaac                                   >  1:3847570‑M2/32‑36 (MQ=255)
  tttcatcgtccaacaaaatgggtgcagtacaacaac                                   <  1:7656018‑M2/5‑1 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   >  1:4309542‑M2/1‑34 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   <  1:7331300‑M2/36‑3 (MQ=255)
   ttcatcgtccaacaaaatgggtgcagtacaACAACt                                  <  1:3901854‑M2/6‑1 (MQ=255)
   ttcatcgtccaacaaaatgggtgcagtacaACAACt                                  <  1:4918491‑M2/6‑1 (MQ=255)
    tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 >  1:1420488‑M2/30‑36 (MQ=255)
    ctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgta                                 >  1:7837541‑M2/1‑32 (MQ=255)
     tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                                <  1:6826474‑M2/36‑6 (MQ=255)
     tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                                >  1:3944634‑M2/1‑31 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:2807239‑M2/8‑1 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:317594‑M2/8‑1 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               <  1:499118‑M2/9‑1 (MQ=255)
      atcgtccaacaaaatgggggcagtacaACAACTTcg                               >  1:5236970‑M2/28‑36 (MQ=255)
       tcgtccaacaaaatgggtgcagtacaACAACTTcga                              >  1:6824615‑M2/27‑36 (MQ=255)
       tcgtccaacaaaatgggtgcagtacaACAACTTcga                              >  1:6679846‑M2/27‑36 (MQ=255)
       cAAAATCGAAAACCAGGAAGGCCTCAACAtgtactg                              >  1:6852691‑M2/1‑29 (MQ=255)
        cgtccaacaaaattggtgcagtacaACAACTTcgac                             <  1:6486538‑M2/11‑1 (MQ=255)
        aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                             <  1:7858885‑M2/36‑9 (MQ=255)
        aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                             <  1:3909576‑M2/36‑9 (MQ=255)
         gtccaacaaaatgggtgcagtacaACAACTTcgacg                            >  1:4853182‑M2/25‑36 (MQ=255)
         aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                            <  1:4683240‑M2/36‑10 (MQ=255)
         aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                            >  1:5670162‑M2/1‑27 (MQ=255)
          aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                           <  1:1055340‑M2/36‑11 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          <  1:1788330‑M2/14‑1 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          <  1:6054337‑M2/36‑12 (MQ=255)
            tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                         >  1:1213674‑M2/1‑24 (MQ=255)
            tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                         >  1:3131583‑M2/1‑24 (MQ=255)
             aacaaaatgggtgcagtacaACAACTTCGACGAAAt                        >  1:3589778‑M2/21‑36 (MQ=255)
             aacaaaatgggtgcagtacaACAACTTCGACGAAAt                        <  1:5401105‑M2/16‑1 (MQ=255)
             aacaaaatgggtgcagtacaACAACTTCGACGAAAt                        <  1:2227349‑M2/16‑1 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgtattgaccccaa                       >  1:4305331‑M2/1‑22 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                       <  1:2555365‑M2/36‑15 (MQ=255)
              acaaaatgggtgcagtacaACAACTTCGACGAAATc                       <  1:6075137‑M2/17‑1 (MQ=255)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      >  1:5386949‑M2/19‑36 (MQ=255)
               aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      >  1:6904157‑M2/1‑21 (MQ=255)
                aaaatgggtgcagtacaACAACTTCGCCGAAATCCt                     <  1:7797036‑M2/19‑1 (MQ=255)
                aaaatgggtgcagtacaACAACTTCGACGAAATTCt                     >  1:3501915‑M2/18‑36 (MQ=255)
                aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                     <  1:604363‑M2/36‑17 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  1:4003030‑M2/17‑36 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    <  1:3849578‑M2/20‑1 (MQ=255)
                 aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                    >  1:5942946‑M2/1‑19 (MQ=255)
                 aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                    <  1:4653535‑M2/36‑18 (MQ=255)
                  aatgggtgcagtacaACAACTTCGACGAAATCCTCg                   >  1:5815710‑M2/16‑36 (MQ=255)
                  aatgggtgcagtacaACAACTTCGACGAAATCCTCg                   >  1:5073817‑M2/16‑36 (MQ=255)
                  aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                   <  1:369749‑M2/36‑19 (MQ=255)
                   atgggtgcagtacaACAACTTCGACGAAATCCTCGa                  <  1:2260805‑M2/22‑1 (MQ=255)
                   atgggtgcagtacaACAACTTCGACGAAATCCTCGa                  >  1:1986145‑M2/15‑36 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  1:4201173‑M2/23‑1 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  1:487928‑M2/23‑1 (MQ=255)
                    tgggtgcagtaaaACAACTTCGACGAAATCCTCGaa                 >  1:4560388‑M2/14‑36 (MQ=255)
                    cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt                 <  1:4310337‑M2/36‑21 (MQ=255)
                    cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt                 <  1:6183802‑M2/36‑21 (MQ=255)
                     aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                <  1:5003982‑M2/36‑22 (MQ=255)
                      ggtgcagtacaACAACTTCGACGAAATCCTCGAAGc               <  1:3194686‑M2/25‑1 (MQ=255)
                      ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:6720214‑M2/36‑23 (MQ=255)
                      ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:2873443‑M2/36‑23 (MQ=255)
                        agcagtacaACAACTGCGACGAAATCCTCGAAGCct             <  1:2460156‑M2/27‑1 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             >  1:4922092‑M2/1‑12 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             >  1:1370259‑M2/1‑12 (MQ=255)
                         gcagtacaACAACTTCGACGAAATCCTCGAAGCctc            <  1:1383247‑M2/28‑1 (MQ=255)
                         aGGCCTCAACAtgtactgaccgcaaaaagttggaca            >  1:6345976‑M2/1‑11 (MQ=255)
                         aGGCCTCAACAtgtactgaccccaaaaagttggaca            >  1:3623799‑M2/1‑11 (MQ=255)
                         aGGCCTCAACAtgtactgaccccaaaaagttggaca            >  1:3144897‑M2/1‑11 (MQ=255)
                          ggCCTCAACAtgtactgaccccaaaaagttggacag           <  1:6963972‑M2/36‑27 (MQ=255)
                          cagtaaaACAACTTCGACGAAATCCTCGAAGCctct           <  1:361535‑M2/29‑1 (MQ=255)
                           agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg          <  1:3965450‑M2/30‑1 (MQ=255)
                           agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg          >  1:7490265‑M2/7‑36 (MQ=255)
                           agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg          >  1:5229963‑M2/7‑36 (MQ=255)
                            gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa         >  1:4505912‑M2/6‑36 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggaaagtt         <  1:4471424‑M2/36‑29 (MQ=255)
                             tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc        >  1:785115‑M2/5‑36 (MQ=255)
                             tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc        >  1:7458637‑M2/5‑36 (MQ=255)
                             cccaacacCTTCGACGAAAGCCTCGAAGCCTCTGAc        <  1:4021859/34‑1 (MQ=37)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        >  1:6590643‑M2/1‑7 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        >  1:5940827‑M2/1‑7 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        >  1:4623959‑M2/1‑7 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        <  1:4061294‑M2/36‑30 (MQ=255)
                              tCAACAtgtactgaccccaaaaaggtggacagttaa       <  1:1244634‑M2/36‑31 (MQ=255)
                               cAACAtgtactgaccccaaaaagttggacagttaaa      >  1:6911413‑M2/1‑5 (MQ=255)
                               cAACAtgtactgaccccaaaaagttggacagttaaa      <  1:7625522‑M2/36‑32 (MQ=255)
                               cAACAtgtactgaccccaaaaagttggacagttaaa      <  1:2354904‑M2/36‑32 (MQ=255)
                                aaCCtgtactgaccccaaaaagttggacagttaaac     >  1:1375041‑M2/1‑4 (MQ=255)
                                aaCCtgtactgaccccaaaaagttggacagttaaac     <  1:589376‑M2/36‑33 (MQ=255)
                                aaCAtgtactgaccccaaaaagttggacagttaaac     >  1:1457406‑M2/1‑4 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     >  1:1312691‑M2/2‑36 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     >  1:2936837‑M2/2‑36 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     >  1:6638028‑M2/2‑36 (MQ=255)
                                 aCAtgtactgaccccaaaaagttggacagttaaaca    <  1:5977080‑M2/36‑34 (MQ=255)
                                  CAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAt   >  1:5126259‑M2/1‑36 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  >  1:1263016‑M2/1‑1 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  <  1:6720182/36‑1 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  <  1:841172/36‑1 (MQ=255)
                                   aacTCCGACGAAATCCTCGAAGCCTCTGACGGcatc  <  1:4428633/36‑1 (MQ=255)
                                   |                                   
TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG