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

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT...................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  REL606/1733649‑1733614
                                                                       
  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATGT                                   <  1:5977080/36‑1 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATGTT                                  >  1:1457406/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                  <  1:589376/36‑1 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                  >  1:1375041/1‑36 (MQ=1)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACATGTTG                                 <  1:7625522/36‑1 (MQ=1)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACATGTTG                                 <  1:2354904/36‑1 (MQ=1)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACATGTTG                                 >  1:6911413/1‑36 (MQ=1)
     TTAACTGTCCACCTTTTTGGGGTCAGTACATGTTGA                                <  1:1244634/36‑1 (MQ=0)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                               >  1:5940827/1‑36 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                               >  1:4623959/1‑36 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                               <  1:4061294/36‑1 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                               >  1:6590643/1‑36 (MQ=1)
       AACTTTCCAACTTTTTGGGGTCAGTACATGTTGAGG                              <  1:4471424/36‑1 (MQ=0)
       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                              >  1:1352476/1‑36 (MQ=1)
         CTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCC                            <  1:6963972/36‑1 (MQ=1)
          TGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCT                           >  1:3144897/1‑36 (MQ=1)
          TGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCT                           >  1:3623799/1‑36 (MQ=1)
          TGTCCAACTTTTTGCGGTCAGTACATGTTGAGGCCT                           >  1:6345976/1‑36 (MQ=0)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                          >  1:4922092/1‑36 (MQ=1)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                          >  1:1370259/1‑36 (MQ=1)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                        <  1:6720214/36‑1 (MQ=1)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                        <  1:2873443/36‑1 (MQ=1)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                       <  1:5003982/36‑1 (MQ=32)
               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                      <  1:6183802/36‑1 (MQ=35)
               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                      <  1:4310337/36‑1 (MQ=35)
                ACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGG                     >  1:980089/1‑36 (MQ=38)
                 CTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGT                    <  1:369749/36‑1 (MQ=38)
                  TTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTT                   <  1:4653535/36‑1 (MQ=38)
                  TTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTT                   >  1:5942946/1‑36 (MQ=38)
                   TTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTT                  <  1:604363/36‑1 (MQ=38)
                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT                 >  1:6904157/1‑36 (MQ=38)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC                <  1:2555365/36‑1 (MQ=38)
                     TTGGGGTCAATACATGTTGAGGCCTTCCTGGTTTTC                >  1:4305331/1‑36 (MQ=0)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG               <  1:5959476/36‑1 (MQ=38)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG               <  1:4811789/36‑1 (MQ=38)
                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA              >  1:1213674/1‑36 (MQ=38)
                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA              >  1:1843478/1‑36 (MQ=38)
                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA              >  1:3131583/1‑36 (MQ=38)
                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT             <  1:6054337/36‑1 (MQ=38)
                         GGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATT            <  1:1055340/36‑1 (MQ=38)
                          GTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTT           <  1:4683240/36‑1 (MQ=38)
                          GTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTT           >  1:5670162/1‑36 (MQ=38)
                           TCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTT          <  1:3909576/36‑1 (MQ=38)
                           TCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTT          <  1:7858885/36‑1 (MQ=38)
                            CAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTG         >  1:6852691/1‑36 (MQ=38)
                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGG        <  1:7771319/36‑1 (MQ=35)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA       >  1:6411285/1‑36 (MQ=32)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA       >  1:3944634/1‑36 (MQ=32)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA       <  1:6826474/36‑1 (MQ=32)
                               TACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAG      >  1:7837541/1‑36 (MQ=1)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA     >  1:3188733/1‑36 (MQ=1)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    >  1:4309542/1‑36 (MQ=1)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    <  1:7331300/36‑1 (MQ=1)
                                  ATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG   <  1:4221973/36‑1 (MQ=1)
                                  ATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG   >  1:5005022/1‑36 (MQ=1)
                                                                       
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT...................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  REL606/1733649‑1733614

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