Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 94.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 (1.010) 26/62 NT 100% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 0 (0.000)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)56 (1.280) 35/68 NT 90.3% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 6 (0.130)coding (685/1413 nt) pykF pyruvate kinase

Only 100 of 107 total aligned reads displayed.

CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682
                                   |                                   
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:1416733‑M2/1‑1 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:4671371‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:4065497‑M2/1‑1 (MQ=255)
tccatttcatcgtccaacaaaatcggtgcagtacaa                                     <  1:4498989‑M2/1‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:2752458‑M2/2‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    >  1:4679702‑M2/35‑36 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    >  1:3956152‑M2/35‑36 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    <  1:3381411/36‑1 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    <  1:488046/36‑1 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  >  1:3997097‑M2/33‑36 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaacaa                                  >  1:4471506‑M2/33‑36 (MQ=255)
    tttcatcgtccaacaaaatgggtgcagtacaacaac                                 >  1:936180‑M2/32‑36 (MQ=255)
    atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                 <  1:3749872‑M2/36‑3 (MQ=255)
     ttcatcgtccaacaaaatgcgtgcagtacaACAACt                                >  1:4190529‑M2/31‑36 (MQ=255)
     tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                <  1:1789080‑M2/36‑4 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               >  1:4625495‑M2/30‑36 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               >  1:5134107‑M2/30‑36 (MQ=255)
      tcatcgtccaacaaaatgggtgcagtacaACAACtt                               >  1:2531073‑M2/30‑36 (MQ=255)
      ctcCAAAATCGAAAACCAGGAAGGCCTCAACGtgta                               >  1:1229949‑M2/1‑32 (MQ=255)
      caccaaAATCGAAAACCAGGAAGGCCTCAACAtgta                               >  1:3820527‑M2/3‑32 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:1937326‑M2/8‑1 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:5376590‑M2/8‑1 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaACAACTTc                              <  1:734171‑M2/8‑1 (MQ=255)
        ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                             >  1:2339346‑M2/1‑30 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaACAACTTcg                             >  1:289791‑M2/28‑36 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaACAACTTcga                            <  1:25505‑M2/10‑1 (MQ=255)
           gtccaacaaaatgggtgcagtacaACAACTTcgacg                          >  1:5166559‑M2/25‑36 (MQ=255)
           aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                          >  1:2898888‑M2/1‑27 (MQ=255)
            aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                         <  1:2686415‑M2/36‑11 (MQ=255)
            aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                         <  1:3381031‑M2/36‑11 (MQ=255)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        <  1:3649576‑M2/14‑1 (MQ=255)
             ccaacaaaatgggtgcagtacaACAACTTCGACGaa                        <  1:2992076‑M2/14‑1 (MQ=255)
             aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                        >  1:1886920‑M2/1‑25 (MQ=255)
             aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                        <  1:1447377‑M2/36‑12 (MQ=255)
               cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                      >  1:4555365‑M2/1‑23 (MQ=255)
               aacaaaatgggtgcagtacaACAACTTCGACGAAAt                      >  1:5531616‑M2/21‑36 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                     >  1:941291‑M2/1‑22 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                     <  1:5259580‑M2/36‑15 (MQ=255)
                acaaaatgggtgcagtacaACAACTTCGACGAAATc                     >  1:238384‑M2/20‑36 (MQ=255)
                 caaaatgggtgcagtacaACAACTTCGACGAAATcc                    >  1:3213291‑M2/19‑36 (MQ=255)
                 aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                    <  1:4651056‑M2/36‑16 (MQ=255)
                 aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                    <  1:568833‑M2/36‑16 (MQ=255)
                 aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                    <  1:5221215‑M2/36‑16 (MQ=255)
                  aaaatgggtgcagtacaACAACTTTGACGAAATCCt                   >  1:409887‑M2/18‑36 (MQ=255)
                  aaaatgggtgcagtacaACAACTTCGACGAAATCCt                   <  1:2107502‑M2/19‑1 (MQ=255)
                  aaaatgggtgcagtacaACAACTTCGACGAAATCCt                   <  1:4492684‑M2/19‑1 (MQ=255)
                  aaaatgggtgcagtacaACAACTTCGACGAAATCCt                   <  1:3029847‑M2/19‑1 (MQ=255)
                  aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                   <  1:675657‑M2/36‑17 (MQ=255)
                  aaaCCAGGAAGGCATCAACAtgtactgaccccaaaa                   >  1:3441423‑M2/1‑20 (MQ=255)
                  aaaCCAGGAAGGCATCAACAtgtactgaccccaaaa                   <  1:3732215‑M2/36‑17 (MQ=255)
                   aaatgggtgcagtacaACAACTTCGACGAAATCCTc                  >  1:3507338‑M2/17‑36 (MQ=255)
                   aaatgggtgcagtacaACAACTTCGACGAAATCCTc                  >  1:1580148‑M2/17‑36 (MQ=255)
                   aaatgggtgcagtacaACAACTTCGACGAAATCCTc                  >  1:1804980‑M2/17‑36 (MQ=255)
                   aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                  <  1:678913‑M2/36‑18 (MQ=255)
                    aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                 <  1:5000256‑M2/36‑19 (MQ=255)
                     atgggtgcagtacaACAACTTCGACGAAATCCTCGa                <  1:485502‑M2/22‑1 (MQ=255)
                     atgggtgcagtacaACAACTTCGACGAAATCCTCGa                >  1:5920655‑M2/15‑36 (MQ=255)
                     atgggtgcagtacaACAACTTCGACGAAAGCCTCGa                >  1:5448741‑M2/15‑36 (MQ=255)
                      tgggtgcagtacaACAACTTCGACGAAATCCTCGaa               >  1:1181169‑M2/14‑36 (MQ=255)
                      tgggtgcagtacaACAACTTCGACGAAATCCTCGaa               <  1:3201918‑M2/23‑1 (MQ=255)
                       aGGAAGGCCTCAACAtgtactgaccccaaaaagttg              <  1:5912946‑M2/36‑22 (MQ=255)
                       aGGAAGGCCTCAACAtgtactgaccccaaaaagttg              <  1:2270580‑M2/36‑22 (MQ=255)
                       aGGAAGGCCTCAACAtgtactgaccccaaaaagttg              <  1:291479‑M2/36‑22 (MQ=255)
                       aGGAAGGCCTCAACAcgtactgaccccaaaaagttg              <  1:4367657‑M2/36‑22 (MQ=255)
                        ggAAGGCCTCAACAtgtactgaccccaaaaagttgg             <  1:4480476‑M2/36‑23 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc            >  1:1069879‑M2/11‑36 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc            >  1:4214161‑M2/11‑36 (MQ=255)
                         gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc            <  1:3417889‑M2/26‑1 (MQ=255)
                         gAAGGCCTCAACAtgtactgaccccaaaaagttgga            <  1:224995‑M2/36‑24 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           <  1:5642823‑M2/36‑25 (MQ=255)
                          aaGGCCTCAACAtgtactgaccccaaaaagttggac           >  1:1755164‑M2/1‑12 (MQ=255)
                           aGGCCTCAACAtgtactgaccccaaaaagttggaca          <  1:905139‑M2/36‑26 (MQ=255)
                            ggCCTCAACAtgtactgaccccaaaaagttggacag         >  1:1127081‑M2/1‑10 (MQ=255)
                            cagtacaACAACTTCGACGAAATCCTCGAAGCctct         <  1:509864‑M2/29‑1 (MQ=255)
                             gCCTCAACAtgtactgaccccaaaaagttggacagt        <  1:1864177‑M2/36‑28 (MQ=255)
                             gCCTCAACAtgtactgaccccaaaaagttggacagt        <  1:2055804‑M2/36‑28 (MQ=255)
                             agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg        >  1:2195698‑M2/7‑36 (MQ=255)
                              gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa       >  1:4038898‑M2/6‑36 (MQ=255)
                              ccTCAACAtgtactgaccccaaaaagttggacagtt       >  1:3983917‑M2/1‑8 (MQ=255)
                              ccTCAACAtgtactgaccccaaaaagttggacagtt       >  1:2972698‑M2/1‑8 (MQ=255)
                               tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc      <  1:57344‑M2/32‑1 (MQ=255)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      <  1:884641‑M2/36‑30 (MQ=255)
                               cTCAACAtgtactgaccccaaaaagttggacagtta      >  1:414481‑M2/1‑7 (MQ=255)
                                tCAACAtgtactgaccccaaaaagttggacagttaa     >  1:4447495‑M2/1‑6 (MQ=255)
                                tCAACAtgtactgaccccaaaaagttggacagttaa     <  1:4891742‑M2/36‑31 (MQ=255)
                                tCAACAtgtactgaccccaaaaagttggacagttaa     <  1:5311735‑M2/36‑31 (MQ=255)
                                gcaACCtgtactgaccccaaaaagttggacagttaa     <  1:1131‑M2/35‑31 (MQ=255)
                                gcaACCtgtactgaccccaaaaagttggacagttaa     <  1:3929976‑M2/35‑31 (MQ=255)
                                 cAACAtgtactgaccccaaaaagttggacagttaaa    >  1:1897360‑M2/1‑5 (MQ=255)
                                 cAACAtgtactgaccccaaaaagttggacagttaaa    <  1:3132601‑M2/36‑32 (MQ=255)
                                 cAACAtgtactgaccccaaaaagtcggacagttaaa    <  1:5144045‑M2/36‑32 (MQ=255)
                                 cAACAtgtactgacaccaaaaagttggacagttaaa    <  1:1540506‑M2/36‑32 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   <  1:3560027‑M2/36‑33 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   <  1:775052‑M2/36‑33 (MQ=255)
                                  aaCCtgtactgaccccaaaaagttggacagttaaac   <  1:1088683‑M2/36‑33 (MQ=255)
                                  aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc   <  1:2083104‑M2/35‑1 (MQ=255)
                                  aACAACTTCGACGAAACCCTCGAAGCCTCTGACGGc   <  1:2199506‑M2/35‑1 (MQ=255)
                                   acaacTTCGACGAAATCCTCGAAGCCTCTGTCggaa  >  1:4462478/1‑34 (MQ=38)
                                   aCAtgtactgaccccaaaaagttggacagttaaaca  <  1:5700712‑M2/36‑34 (MQ=255)
                                   ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa  >  1:2263186‑M2/1‑36 (MQ=255)
                                   |                                   
CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG