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

CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682
                                   |                                   
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:924248‑M2/1‑1 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:1012524‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:6988553‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:625286‑M2/1‑1 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:4151592‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:3911571‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:3756294‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:2420141‑M2/1‑1 (MQ=255)
catcatCTCCAAAATCGAAAACCAGGAAGGCCTcaa                                     >  1:3290565/1‑36 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:7435840‑M2/2‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:4060134‑M2/2‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    >  1:3218305‑M2/35‑36 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:2859390‑M2/2‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    >  1:259348‑M2/35‑36 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    <  1:6466526/36‑1 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   <  1:7774606‑M2/3‑1 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   <  1:1006445‑M2/3‑1 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   >  1:1327762‑M2/34‑36 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                  >  1:6862688‑M2/1‑35 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                  >  1:2248579‑M2/1‑35 (MQ=255)
   catcTCCAAAATAGAAAAACAGGAAAGCCTcaacat                                  >  1:7854007/1‑35 (MQ=21)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  <  1:7612709‑M2/4‑1 (MQ=255)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                >  1:1442553‑M2/1‑33 (MQ=255)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                >  1:4958776‑M2/1‑33 (MQ=255)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                <  1:2265101‑M2/36‑4 (MQ=255)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                >  1:3962631‑M2/1‑33 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               <  1:4953085‑M2/7‑1 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              <  1:7767779‑M2/36‑6 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              <  1:2111451‑M2/36‑6 (MQ=255)
       tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                              >  1:4231114‑M2/1‑31 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:5314379‑M2/8‑1 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:3330933‑M2/8‑1 (MQ=255)
        ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                             <  1:7081885‑M2/36‑7 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             <  1:5830176‑M2/9‑1 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaACAACTTcga                            >  1:1720457‑M2/27‑36 (MQ=255)
          cgtccaacaaaatgggtgcagtacaACAACTTcgac                           >  1:6665799‑M2/26‑36 (MQ=255)
          aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                           <  1:214213‑M2/36‑9 (MQ=255)
          aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                           <  1:6198348‑M2/36‑9 (MQ=255)
          aaaaTCGAAAACCAGGAAGGACTCAACAACTTcgac                           >  1:6366001/1‑36 (MQ=255)
           aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                          <  1:886220‑M2/36‑10 (MQ=255)
            tccaacaaaatgggtgcagtacaACAACTTcgacga                         >  1:3592864‑M2/24‑36 (MQ=255)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        <  1:5823721‑M2/14‑1 (MQ=255)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        <  1:3833232‑M2/14‑1 (MQ=255)
              tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                       >  1:4552755‑M2/1‑24 (MQ=255)
              caacaaaatgggtgcagtacaACAACTTCGACGaaa                       <  1:7461775‑M2/15‑1 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      <  1:1480102‑M2/36‑14 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      <  1:1713749‑M2/36‑14 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      <  1:7154245‑M2/36‑14 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      >  1:2695741‑M2/21‑36 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      >  1:4214091‑M2/21‑36 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                     >  1:6836706‑M2/1‑22 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                     <  1:5025707‑M2/36‑15 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgcactgaccccaa                     <  1:1149210‑M2/36‑15 (MQ=255)
                acaaaatgggtgcagtacaACAACTTCGACGAAATc                     >  1:665546‑M2/20‑36 (MQ=255)
                 caaaatgggtgcagtacaACAACTTCGACGAAATcc                    <  1:2564286‑M2/18‑1 (MQ=255)
                  aaaatgggtgcagtacaACAACTTCGACGAAATCCt                   >  1:1282515‑M2/18‑36 (MQ=255)
                  aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                   >  1:3335562‑M2/1‑20 (MQ=255)
                  aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                   >  1:5081604‑M2/1‑20 (MQ=255)
                   aaatgggtgcagtacaACAACTTCGACGAAATCCTc                  >  1:7104671‑M2/17‑36 (MQ=255)
                   aaatgggtgcagtacaACAACTTCGACGAAATCCTc                  <  1:5380533‑M2/20‑1 (MQ=255)
                   aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                  <  1:377793‑M2/36‑18 (MQ=255)
                    aatgggtgcagtacaACAACTTCGACGAAATCCTCg                 >  1:1677531‑M2/16‑36 (MQ=255)
                    aatgggtgcagtacaACAACTTCGACGAAATCCTCg                 <  1:6507697‑M2/21‑1 (MQ=255)
                      cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt               <  1:2208025‑M2/36‑21 (MQ=255)
                      cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt               <  1:124213‑M2/36‑21 (MQ=255)
                       gggtgaagtacaACAACTTCGACGAAATCCTCGAAg              <  1:123758‑M2/24‑1 (MQ=255)
                       aGGAAGGCCTCAACAtgtactgaccccaaaaagttg              <  1:676429‑M2/36‑22 (MQ=255)
                       aGGAAGGCCTCAACAtgtactgaccccaaaaagttg              <  1:1740802‑M2/36‑22 (MQ=255)
                        ggtgcagtacaACAACTTCGACGAAATCCTCGAAGc             >  1:308886‑M2/12‑36 (MQ=255)
                        ggAAGGCCTCAACAtgtactgaccccaaaaagttgg             <  1:5638919‑M2/36‑23 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGGAGcc            >  1:6310895‑M2/11‑36 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc            >  1:2864339‑M2/11‑36 (MQ=255)
                         gtgcagtacaACAACTCCGACGAAATCCTCGAAGcc            <  1:1450951‑M2/26‑1 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           <  1:7013053‑M2/36‑25 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           <  1:363499‑M2/36‑25 (MQ=255)
                           aGGCCTCAACAtgtactgaccccaaaaagttggaca          >  1:6160010‑M2/1‑11 (MQ=255)
                            ggCCTCAACAtgtactgaccccaaaaagttggacag         <  1:618903‑M2/36‑27 (MQ=255)
                            ggCCTCAACAtgtactgaccccaaaaagttggacag         <  1:7355592‑M2/36‑27 (MQ=255)
                            ggCCTCAACAtgtactgaccccaaaaagttggacag         >  1:5595487‑M2/1‑10 (MQ=255)
                            cagtacaACAACTTCGACGAAATCCTCGAAGCctct         >  1:5921742‑M2/8‑36 (MQ=255)
                             gCCTCAACAtgtactgaccccaaaaagttggacagt        >  1:4369694‑M2/1‑9 (MQ=255)
                             agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg        >  1:3803752‑M2/7‑36 (MQ=255)
                              gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa       >  1:908203‑M2/6‑36 (MQ=255)
                              gtacaACAACTTCGACGAAAGCCTCGAAGCCTCTGa       <  1:2359024‑M2/31‑1 (MQ=255)
                              ccTCAACAtgtactgaccccaaaaagttggacagtt       >  1:5863688‑M2/1‑8 (MQ=255)
                               tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc      <  1:6711578‑M2/32‑1 (MQ=255)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      <  1:1024346‑M2/36‑30 (MQ=255)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      <  1:1498665‑M2/36‑30 (MQ=255)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      <  1:5503595‑M2/36‑30 (MQ=255)
                                tCAACAtgtactgaccccaaaaagttggacagttaa     >  1:4499705‑M2/1‑6 (MQ=255)
                                gcaACCtgtactgaccccaaaaagttggacagttaa     <  1:2563011‑M2/35‑31 (MQ=255)
                                acaACAACTTCGACGAAATCCTCGAAGCCTCTGACg     <  1:7077003‑M2/33‑1 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   >  1:4811070‑M2/1‑4 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   <  1:4863894‑M2/36‑33 (MQ=255)
                                  aaCAtgtactgaccccaaaaggttggacagttaaac   <  1:7415064‑M2/36‑33 (MQ=255)
                                  aaCAtgtactgaccccaaaaagttggacagttaaac   >  1:3099488‑M2/1‑4 (MQ=255)
                                  aaCAtgtactgaccccaaaaagttggacagttaaac   >  1:5722040‑M2/1‑4 (MQ=255)
                                  aaCAtgtactgaccccaaaaagttggacagttaaac   >  1:912950‑M2/1‑4 (MQ=255)
                                  aACAACTTCAACGAAATCCTCGAAGCCTCTGACGGc   <  1:2098290‑M2/35‑1 (MQ=255)
                                   ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa  >  1:5783622‑M2/1‑36 (MQ=255)
                                   |                                   
CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682

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