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

CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682
                                   |                                   
tccatttcatcgtccaacaaagtgggtgcagtacaa                                     <  1:5767734‑M2/1‑1 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:6024952‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:3629451‑M2/1‑1 (MQ=255)
catcatCTCCAAAATCGAAAACCAGGAAGGCCTcaa                                     <  1:6745298/36‑1 (MQ=255)
catcatCTCCAAAATCGAAAACCAGGAAGGCCTcaa                                     >  1:5717968/1‑36 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    >  1:2911654‑M2/35‑36 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    >  1:3992589/1‑36 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    >  1:3181682/1‑36 (MQ=255)
 atcatcTCAAAAATCGAAAACCAGGAAGGCCTcaac                                    <  1:918605/36‑1 (MQ=255)
  tcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACA                                   >  1:1332448‑M2/1‑36 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaacc                                   <  1:4613693‑M2/3‑2 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   <  1:4125444‑M2/3‑1 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   >  1:4465788‑M2/34‑36 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                  >  1:5005022‑M2/1‑35 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacct                                  <  1:6548326‑M2/4‑3 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  <  1:5861480‑M2/4‑1 (MQ=255)
    tttcatcgtccaacaaaatgggtgcagtacaacaac                                 >  1:3767338‑M2/32‑36 (MQ=255)
    tttcatcgtccaacaaaatgggtgcagtacaacaac                                 >  1:3847570‑M2/32‑36 (MQ=255)
    tttcatcgtccaacaaaatgggtgcagtacaacaac                                 >  1:7620310‑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)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                >  1:3188733‑M2/1‑33 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               >  1:1420488‑M2/30‑36 (MQ=255)
      ctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgta                               >  1:7837541‑M2/1‑32 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              >  1:6411285‑M2/1‑31 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              <  1:6826474‑M2/36‑6 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              >  1:3944634‑M2/1‑31 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:317594‑M2/8‑1 (MQ=255)
        ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                             <  1:7771319‑M2/36‑7 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             <  1:499118‑M2/9‑1 (MQ=255)
        atcgtccaacaaaatgggggcagtacaACAACTTcg                             >  1:5236970‑M2/28‑36 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaACAACTTcga                            <  1:2738075‑M2/10‑1 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaACAACTTcga                            >  1:6679846‑M2/27‑36 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaACAACTTcga                            >  1:6824615‑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)
           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)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        <  1:3950285‑M2/14‑1 (MQ=255)
             aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                        <  1:6054337‑M2/36‑12 (MQ=255)
              tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                       >  1:1213674‑M2/1‑24 (MQ=255)
              tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                       >  1:1843478‑M2/1‑24 (MQ=255)
              tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                       >  1:3131583‑M2/1‑24 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      <  1:4811789‑M2/36‑14 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      <  1:2227349‑M2/16‑1 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      >  1:3589778‑M2/21‑36 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      <  1:5401105‑M2/16‑1 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtattgaccccaa                     >  1:4305331‑M2/1‑22 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                     <  1:2555365‑M2/36‑15 (MQ=255)
                acaaaatgggtgcagtacaACAACTTCGACGAAATc                     <  1:447752‑M2/17‑1 (MQ=255)
                 caaaatgggtgcagtacaACAACTTCGACGAAATcc                    >  1:5386949‑M2/19‑36 (MQ=255)
                 aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                    >  1:6904157‑M2/1‑21 (MQ=255)
                  aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                   <  1:604363‑M2/36‑17 (MQ=255)
                   aaatgggtgcagtacaACAACTTCGACGAAATCCTc                  <  1:3849578‑M2/20‑1 (MQ=255)
                   aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                  <  1:4653535‑M2/36‑18 (MQ=255)
                   aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                  >  1:5942946‑M2/1‑19 (MQ=255)
                    aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                 <  1:369749‑M2/36‑19 (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:6183802‑M2/36‑21 (MQ=255)
                      cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt               <  1:4310337‑M2/36‑21 (MQ=255)
                       aGGAAGGCCTCAACAtgtactgaccccaaaaagttg              <  1:5003982‑M2/36‑22 (MQ=255)
                        ggtgcagtacaATAACTTCGACGAAATCCTCGAAGc             <  1:5576135‑M2/25‑1 (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)
                            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:5229963‑M2/7‑36 (MQ=255)
                             agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg        <  1:7306423‑M2/30‑1 (MQ=255)
                              gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa       >  1:4505912‑M2/6‑36 (MQ=255)
                               tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc      >  1:7458637‑M2/5‑36 (MQ=255)
                               tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc      >  1:785115‑M2/5‑36 (MQ=255)
                               cccaacacCTTCGACGAAAGCCTCGAAGCCTCTGAc      <  1:4021859/34‑1 (MQ=37)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      <  1:4061294‑M2/36‑30 (MQ=255)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      >  1:5940827‑M2/1‑7 (MQ=255)
                                tCAACAtgtactgaccccaaaaaggtggacagttaa     <  1:1244634‑M2/36‑31 (MQ=255)
                                 cAACAtgtactgaccccaaaaagttggacagttaaa    <  1:7625522‑M2/36‑32 (MQ=255)
                                 cAACAtgtactgaccccaaaaagttggacagttaaa    <  1:2354904‑M2/36‑32 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   <  1:589376‑M2/36‑33 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   >  1:1375041‑M2/1‑4 (MQ=255)
                                  aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc   >  1:2936837‑M2/2‑36 (MQ=255)
                                  aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc   >  1:1312691‑M2/2‑36 (MQ=255)
                                  aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc   >  1:6638028‑M2/2‑36 (MQ=255)
                                   aCAtgtactgaccccaaaaagttggacagttaaaca  <  1:5977080‑M2/36‑34 (MQ=255)
                                   |                                   
CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682

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