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

GCCGTCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTG.................................  <  REL606/1733681‑1733647
...................................TGTACTGCACCCATTTTGTTGGACGATGAAATGGA  <  REL606/2775877‑2775843
                                                                      
   GTCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTTTA                                 >  2:154502/1‑36 (MQ=0)
   GTCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTt                                 >  2:2680159/1‑35 (MQ=1)
    TCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTtc                                >  2:556875/1‑34 (MQ=1)
    TCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTAC                                <  1:3396807/36‑1 (MQ=1)
    TCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTAC                                <  2:2383630/36‑1 (MQ=1)
     CAGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTACT                               >  2:1859020/1‑36 (MQ=1)
     CAGAGGCTGCGAGGATTTCGTCGAAGTTGTTGTACT                               <  1:381863/36‑1 (MQ=0)
      AGAGGCTTCGAGGATTTCGTCGAAGTTGTTTTACTG                              >  1:2151486/1‑36 (MQ=0)
      AGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTTCTt                              >  1:692170/1‑35 (MQ=0)
      AGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTACTG                              <  2:284243/36‑1 (MQ=1)
      AGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTACTG                              >  1:2805093/1‑36 (MQ=1)
      AGAGGCTTCGAGGATTTCGTCGAAGTTGTTGTACTG                              <  2:1344842/36‑1 (MQ=1)
        AGGCTTCGAGGATTTCGTCGAAGTTGTTGTACTGCA                            <  1:1373527/36‑1 (MQ=1)
        AGGCTTCGAGGATTTCGTCGAAGTTGTTGTACTGCA                            >  1:1972381/1‑36 (MQ=1)
         GGCTTCGAGGATTTCGTCGAAGTTGTTGTACTGCAC                           <  2:314777/36‑1 (MQ=1)
          GCTTCGAGGATTTCGTCGAAGTTGTTGTACTGCTCC                          >  2:2471569/1‑36 (MQ=0)
          GCTTCGAGGATTTCGTCGAAGTTGTTGTACTGCACC                          <  1:3097903/36‑1 (MQ=1)
          GCTTCGAGGATTTCGTCGAAGTTGTTGTACTGCACC                          <  2:781527/36‑1 (MQ=1)
            tgCGAGGATTTCGTCGAAGTTGTTGTACTGCACCCA                        <  1:2228108/34‑1 (MQ=1)
            TTCGAGGATTTCGTCGAAGTTGTTGTACTGCACCCA                        <  2:2215131/36‑1 (MQ=1)
            TTCGAGGATTTCGTCGAAGTTGTTGTACTGCACCCA                        <  2:2372065/36‑1 (MQ=1)
            TTCGAGGATTTCGTCGAAGTTGTTGTACTGCACCCA                        >  1:1349842/1‑36 (MQ=1)
            TTCGAGGATTTCGTCGAAGTTGTTGTACTGCACCCA                        <  1:1951812/36‑1 (MQ=1)
            TTCGAGGATTTCGTCGAAGTTGTTGTACTGCACCCA                        >  2:544944/1‑36 (MQ=1)
              CGAGGATTTCGTCGAAGTTGTTGTACTGCACCCATT                      <  2:968700/36‑1 (MQ=35)
              CGAGGATTTCGTCGAAGTTGTTGTACTGCACCCATT                      <  1:2078444/36‑1 (MQ=35)
               GAGGATTTCGTCGAAGTTGTTGTACTGCACCCATat                     >  1:1242169/1‑34 (MQ=17)
               GAGGATTTCGTCGAAGTTGTTGTACTGCACCCATTT                     >  2:589228/1‑36 (MQ=38)
                AGGATTTCGTCGAAGTTGTTGTACTGCACCCATTTT                    <  1:817429/36‑1 (MQ=40)
                AGGATTTCGTCGAAGTTGTTGTACTGCACCCACTTT                    >  2:762145/1‑36 (MQ=12)
                 GGATTTCGTCGAAGTTGTTGTACTGCACCCATTTTG                   <  2:2941915/36‑1 (MQ=255)
                 GGATTTCGTCGAAGTTGTTGTACTGCACCCATTTTG                   >  2:1417424/1‑36 (MQ=255)
                 GGATTTCGTCGAAGTTGTTGTACTGCACCCATTTTG                   >  1:2607963/1‑36 (MQ=255)
                 GGATTTCGTCGAAGTTGTTGTACTGCACCCATTTTG                   <  2:1060511/36‑1 (MQ=255)
                    TTTCGTCGAAGTTGTTGTACTGCACCCATTTTGTTG                >  1:2016766/1‑36 (MQ=255)
                    TTTCGTCGAAGTTGTTGTACTGCACCCATTTTGTTG                >  1:1254894/1‑36 (MQ=255)
                     TTCGTCGAAGTTGTTGTACTGCACCCATTTTGTTGG               >  1:3168735/1‑36 (MQ=255)
                     TTCGTCGAAGTTGTTGTACTGCACCCATTTTGTTGG               <  2:14550/36‑1 (MQ=255)
                     TTCGTCGAAGTTGTTGTACTGCACCCATTTTGTTGG               <  2:649258/36‑1 (MQ=255)
                      TCGTCGAAGTTGTTGTACTGCACCCATTTTGTTGGA              <  2:1449130/36‑1 (MQ=255)
                      TCGTCGAAGTTGTTGTACTGCACCCATTTTGTTGGA              >  1:237377/1‑36 (MQ=255)
                       CGTCGAAGTTGTTGTACTGCACCCATTTTGTTGGAC             >  1:3376699/1‑36 (MQ=255)
                       CGTCGAAGTTGTTGTACTGCACCCATTTTGTTGGAC             <  1:1863032/36‑1 (MQ=255)
                       CGTCGAAGTTGTTGTACTGCACCCATTTTGTTGGAC             <  1:2894656/36‑1 (MQ=255)
                        GTCGAAGTTGTTGTACTACACCCATTTTGTTGGACG            >  1:428739/1‑36 (MQ=255)
                        GTCGAAGGTGTTGTACTGCACCCATTTTGTTGGACG            <  2:2387617/36‑1 (MQ=255)
                         TCGAAGTTGTTGTACTGCACCCATTTTGTTGGACGA           >  1:703351/1‑36 (MQ=255)
                          CGAAGTTGTTGTACTGCACCCATTTTGTTGGACGAT          <  1:3648480/36‑1 (MQ=255)
                          CGAAGTTGTTGTACTGCACCCATTTTGTTGGACGAT          >  1:669229/1‑36 (MQ=255)
                           GAAGTTGTTGTACTGCACCCATTTTGTTGGACGATG         >  2:181002/1‑36 (MQ=255)
                           GAAGTTGTTGTACTGCACCCATTTTGTTGGACGATG         <  2:2884993/36‑1 (MQ=255)
                            AAGTTGTTGTACTGCACCCATTTTGTTGGACGAgat        >  1:2182853/1‑33 (MQ=255)
                            AAGTTGTTGTACTGCACCCATTTTGTTGGACGATGA        <  1:1859173/36‑1 (MQ=255)
                            AAGTTGTTGTACTGCACCCATTTTGTTGGACGATGA        <  1:2082809/36‑1 (MQ=255)
                            AAGTTGTTGTACTGCACCCATTTTGTTGGACGATGA        >  2:2190461/1‑36 (MQ=255)
                            AAGTTGTTGTACTGCACCCATTTTGTTGGACGATGA        >  2:2188645/1‑36 (MQ=255)
                              GTTGTTGTACTGCACCCATTTTGTTGGACGATGAAA      <  1:2312313/36‑1 (MQ=255)
                               TTGTTGTACTGCACCCATTTTGTTGGACGATGAAAT     <  2:2074278/36‑1 (MQ=40)
                                 GTTGTACTGCACCCATTTTGTTGGACGATGAAATGG   <  1:1965067/36‑1 (MQ=35)
                                 GTTGTACTGCACCCATTTTGTTGGACGATGAAATGG   >  2:2239931/1‑36 (MQ=35)
                                 GTTGTACTGCACCCATTTTGTTGGACGATGAAATGG   >  1:518720/1‑36 (MQ=35)
                                  TTGTACTGCACCCATTTTGTTGGACGATGAAATGGA  >  2:2981188/1‑36 (MQ=32)
                                  TTGTACTGCACCCATTTTGTTGGACGATGAAATGGA  <  1:1602596/36‑1 (MQ=32)
                                  TTGTACTGCACCCATTTTGTTGGACGATGAAATGGA  <  1:3096993/36‑1 (MQ=32)
                                  TTGTACTGCACCCATTTTGTTGGACGATGAAATGGA  <  1:3233510/36‑1 (MQ=32)
                                  TTGTACTGCACCCATTTTGTGGGACGATGAAATGGA  >  1:1319931/1‑36 (MQ=12)
                                                                      
GCCGTCAGAGGCTTCGAGGATTTCGTCGAAGTTGTTG.................................  <  REL606/1733681‑1733647
...................................TGTACTGCACCCATTTTGTTGGACGATGAAATGGA  <  REL606/2775877‑2775843

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG