Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 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 1733647 =0 (0.040)66 (1.270) 36/66 0.1 100% coding (683/1413 nt) pykF pyruvate kinase
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region
* ? REL606 = 17336490 (0.050)47 (0.850) 29/70 0.5 100% coding (685/1413 nt) pykF pyruvate kinase
?REL606 2774435 = NA (NA)noncoding (1443/1443 nt) IS150 repeat region

Only 100 of 115 total aligned reads displayed.

TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684
                                   |                                   
tcatcTCCAAAATCGAAAACCAGGAAGGCCTcaaca                                     <  2:3676893/36‑1 (MQ=255)
 atttcatcgtccaacaaaatgggtgcagtacaacaa                                    <  2:2074278‑M1/4‑1 (MQ=255)
  tttcatcgtccaacaaaatgggtgcagtacaacaac                                   <  1:2312313‑M1/5‑1 (MQ=255)
  aTCTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   >  2:3050603‑M1/1‑34 (MQ=255)
    tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 >  2:2188645‑M1/30‑36 (MQ=255)
    tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 >  2:2190461‑M1/30‑36 (MQ=255)
    tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 <  1:2082809‑M1/7‑1 (MQ=255)
    ctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgta                                 >  1:2616775‑M1/1‑32 (MQ=255)
    atctcgtccaacaaaatgggtgcagtacaACAACtt                                 >  1:2182853‑M1/30‑36 (MQ=255)
     tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                                >  1:802586‑M1/1‑31 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                >  2:181002‑M1/29‑36 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  2:2884993‑M1/8‑1 (MQ=255)
      ccAAGATCGAAAACCAGGAAGGCCTCAACAtgtact                               <  1:1588985‑M1/36‑7 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               <  1:3648480‑M1/9‑1 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               >  1:669229‑M1/28‑36 (MQ=255)
       tcgtccaacaaaatgggtgcagtacaACAACTTcga                              >  1:703351‑M1/27‑36 (MQ=255)
        cgtccaacaaaatgggtgtagtacaACAACTTcgac                             >  1:428739‑M1/26‑36 (MQ=255)
        cgtccaacaaaatgggtgcagtacaACACCTTcgac                             <  2:2387617‑M1/11‑1 (MQ=255)
        aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                             <  2:2740799‑M1/36‑9 (MQ=255)
         gtccaacaaaatgggtgcagtacaACAACTTcgacg                            >  1:3376699‑M1/25‑36 (MQ=255)
         gtccaacaaaatgggtgcagtacaACAACTTcgacg                            <  1:1863032‑M1/12‑1 (MQ=255)
         gtccaacaaaatgggtgcagtacaACAACTTcgacg                            <  1:2894656‑M1/12‑1 (MQ=255)
          tccaacaaaatgggtgcagtacaACAACTTcgacga                           <  2:1449130‑M1/13‑1 (MQ=255)
          tccaacaaaatgggtgcagtacaACAACTTcgacga                           >  1:237377‑M1/24‑36 (MQ=255)
          aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                           >  2:429035‑M1/1‑26 (MQ=255)
          aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                           >  1:1371748‑M1/1‑26 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          <  2:14550‑M1/14‑1 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          >  1:3168735‑M1/23‑36 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          >  2:2006914‑M1/1‑25 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          <  1:691462‑M1/36‑12 (MQ=255)
            tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                         <  1:2723664‑M1/36‑13 (MQ=255)
            caacaaaatgggtgcagtacaACAACTTCGACGaaa                         >  1:1254894‑M1/22‑36 (MQ=255)
            caacaaaatgggtgcagtacaACAACTTCGACGaaa                         >  1:2016766‑M1/22‑36 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                       >  2:3340523‑M1/1‑22 (MQ=255)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      <  2:1060511‑M1/18‑1 (MQ=255)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      >  2:1417424‑M1/19‑36 (MQ=255)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      >  1:2607963‑M1/19‑36 (MQ=255)
               aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      >  2:1760710‑M1/1‑21 (MQ=255)
               aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      <  1:3726680‑M1/36‑16 (MQ=255)
               aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      <  2:323199‑M1/36‑16 (MQ=255)
                aaagtgggtgcagtacaACAACTTCGACGAAATCCt                     >  2:762145‑M1/18‑36 (MQ=255)
                aaaatgggtgcagtacaACAACTTCGACGAAATCCt                     <  1:817429‑M1/19‑1 (MQ=255)
                 atatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  1:1242169‑M1/17‑36 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  2:589228‑M1/17‑36 (MQ=255)
                  aatgggtgcagtacaACAACTTCGACGAAATCCTCg                   <  2:968700‑M1/21‑1 (MQ=255)
                  aatgggtgcagtacaACAACTTCGACGAAATCCTCg                   <  1:2078444‑M1/21‑1 (MQ=255)
                  aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                   >  2:1039357‑M1/1‑18 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCtcgca                 <  1:2228108‑M1/23‑3 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 >  1:1349842‑M1/14‑36 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  2:2372065‑M1/23‑1 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  2:2215131‑M1/23‑1 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 >  2:544944‑M1/14‑36 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  1:1951812‑M1/23‑1 (MQ=255)
                    cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt                 <  1:3356538‑M1/36‑21 (MQ=255)
                     aGGAGGGCCTCAACAtgtactgaccccaaaaagttg                <  1:589780‑M1/36‑22 (MQ=255)
                     aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                <  2:3221748‑M1/36‑22 (MQ=255)
                     aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                >  1:3398128‑M1/1‑15 (MQ=255)
                      ggtgcagtacaACAACTTCGACGAAATCCTCGAAGc               <  2:781527‑M1/25‑1 (MQ=255)
                      ggtgcagtacaACAACTTCGACGAAATCCTCGAAGc               <  1:3097903‑M1/25‑1 (MQ=255)
                      ggagcagtacaACAACTTCGACGAAATCCTCGAAGc               >  2:2471569‑M1/12‑36 (MQ=255)
                      ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:651316‑M1/36‑23 (MQ=255)
                       gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              <  2:314777‑M1/26‑1 (MQ=255)
                       gAAGGCCTCAACAtgtactgaccccaaaaagttgga              >  1:3644252‑M1/1‑13 (MQ=255)
                        tgcagtacaACAACTTCGACGAAATCCTCGAAGCct             >  1:1972381‑M1/10‑36 (MQ=255)
                        tgcagtacaACAACTTCGACGAAATCCTCGAAGCct             <  1:1373527‑M1/27‑1 (MQ=255)
                          cagtacaACAACTTCGACGAAATCCTCGAAGCctct           >  1:2805093‑M1/8‑36 (MQ=255)
                          cagtacaACAACTTCGACGAAATCCTCGAAGCctct           <  2:284243‑M1/29‑1 (MQ=255)
                          cagtaaaACAACTTCGACGAAATCCTCGAAGCctct           >  1:2151486‑M1/8‑36 (MQ=255)
                            gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa         <  2:2383630‑M1/31‑1 (MQ=255)
                            gaacaACAACTTCGACGAAATCCTCGAAGCCTCTGa         >  2:556875‑M1/6‑36 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggacagtt         >  2:2046875‑M1/1‑8 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggacagtt         <  2:1176510‑M1/36‑29 (MQ=255)
                             taaaACAACTTCGACGAAATCCTCGAAGCCTCTGAc        >  2:154502‑M1/5‑36 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        >  1:2387214‑M1/1‑7 (MQ=255)
                              tCAACAtgtactgaccccaaaaagttggacagttaa       <  2:2754198‑M1/36‑31 (MQ=255)
                              tCAACAtgtactgaccccaaaaagttggacagttaa       >  2:2521362‑M1/1‑6 (MQ=255)
                               caACATCTTCGACGAAATCCTCGAAGCCTCTGACgg      >  1:3230652‑M1/3‑36 (MQ=255)
                               cAACCtgtactgaccccaaaaagttggacagttaaa      >  1:40355‑M1/1‑5 (MQ=255)
                               cAACAtgtactgaccgcaaaaagttggacagttaaa      <  2:2258718‑M1/36‑32 (MQ=255)
                               cAACAtctactgaccccaaaaagttggacagttaaa      <  2:2762656‑M1/36‑32 (MQ=255)
                                aaCCtgtactgaccccaaaaagttggacagttaaac     >  2:2986398‑M1/1‑4 (MQ=255)
                                aaCCtgtactgaccccaaaaagttggacagttaaac     <  2:682036‑M1/36‑33 (MQ=255)
                                aaCAtgtactgaccccaaaaagttggacagttaaac     >  1:2608350‑M1/1‑4 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     <  2:95388‑M1/35‑1 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     >  1:2591329‑M1/2‑36 (MQ=255)
                                 gtatgtactgaccccaaaaagttggacagttaaaca    <  1:635345‑M1/34‑34 (MQ=255)
                                 gtatgtactgaccccaaaaagttggacagttaaaca    >  2:591026‑M1/3‑3 (MQ=255)
                                 gtatgtactgaccccaaaaagttggacagttaaaca    <  2:2592512‑M1/34‑34 (MQ=255)
                                 acaacTTCGACGAAATCCTCGAAGCCTCTGACGGca    >  2:715940/1‑36 (MQ=255)
                                 aCAtgtactgaccccaaaaagttggacagttaaaca    >  1:1265789‑M1/1‑3 (MQ=255)
                                  tatgtactgaccccaaaaagttggacagttaaacac   >  1:2026751‑M1/2‑2 (MQ=255)
                                  tatgtactgaccccaaaaagttggacagttaaacac   >  1:3440988‑M1/2‑2 (MQ=255)
                                  caacTTCGACGAAATCCTCGAAGCCTCTGACGGcat   <  1:324202/36‑1 (MQ=255)
                                  caacTGCGACGAAATCCGCGAAGCCGCTGACGGcat   <  2:1997558/36‑1 (MQ=25)
                                  cAtgtactgaccccaaaaagttggacagttaaacac   >  2:3314998‑M1/1‑2 (MQ=255)
                                   atgtgctgaccccaaaaagttggacagttaaacacg  <  2:1264826‑M1/36‑36 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttatacacg  >  2:2847805‑M1/1‑1 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  >  2:632853‑M1/1‑1 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacaag  >  1:2474577‑M1/1‑1 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  >  1:346493/1‑36 (MQ=255)
                                   |                                   
TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG