Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 96.3% 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)60 (1.100) 35/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)43 (0.720) 31/68 NT 91.5% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 4 (0.060)coding (685/1413 nt) pykF pyruvate kinase

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT...................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  REL606/1733649‑1733614
                                                                       
 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATG                                    >  1:4563594/1‑36 (MQ=32)
 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATG                                    >  1:283299/1‑36 (MQ=32)
       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                              <  1:1559359/36‑1 (MQ=255)
       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                              <  1:2330569/36‑1 (MQ=255)
        ACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGC                             >  1:6033688/1‑36 (MQ=255)
          TGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCT                           >  1:3084401/1‑36 (MQ=255)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                          >  1:5799340/1‑36 (MQ=255)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                          <  1:5775843/36‑1 (MQ=255)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                          <  1:1749360/36‑1 (MQ=255)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                        <  1:4138268/36‑1 (MQ=255)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                        <  1:527267/36‑1 (MQ=255)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                       <  1:6826955/36‑1 (MQ=255)
               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                      >  1:5191603/1‑36 (MQ=255)
                ACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGG                     <  1:1375605/36‑1 (MQ=255)
                ACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGG                     >  1:3298903/1‑36 (MQ=255)
                 CTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGT                    >  1:5764128/1‑36 (MQ=255)
                 CTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGT                    <  1:1232796/36‑1 (MQ=255)
                 CTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGT                    >  1:3945251/1‑36 (MQ=255)
                  TTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTT                   >  1:3159756/1‑36 (MQ=255)
                  TTTTTGGGGTCAGGACATGTTGAGGCCTTCCTGGTT                   >  1:1031303/1‑36 (MQ=255)
                   TTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTT                  <  1:1916619/36‑1 (MQ=255)
                   TTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTT                  <  1:168482/36‑1 (MQ=255)
                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT                 <  1:6302874/36‑1 (MQ=255)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC                <  1:4855464/36‑1 (MQ=255)
                      TGGGGTCAGTACATGTTTAGGCCTTCCTGGTTTTCG               <  1:1650661/36‑1 (MQ=255)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG               <  1:3600293/36‑1 (MQ=255)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG               >  1:4551023/1‑36 (MQ=255)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG               >  1:2998714/1‑36 (MQ=255)
                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA              >  1:4390915/1‑36 (MQ=255)
                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA              <  1:6186999/36‑1 (MQ=255)
                         GGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATT            <  1:4255917/36‑1 (MQ=255)
                          GTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTT           >  1:5802670/1‑36 (MQ=255)
                            CAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTG         >  1:2504170/1‑36 (MQ=255)
                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGG        <  1:5149931/36‑1 (MQ=255)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA       <  1:6002454/36‑1 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA     >  1:8330184/1‑36 (MQ=255)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    <  1:1747236/36‑1 (MQ=255)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    <  1:2273683/36‑1 (MQ=255)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    >  1:4906455/1‑36 (MQ=255)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    <  1:6550055/36‑1 (MQ=255)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT    <  1:6806048/36‑1 (MQ=255)
                                  ATGTTGAGGCCTTCCTGGTTTTCGATTTTTGAGATG   >  1:3205738/1‑36 (MQ=255)
                                  ATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG   <  1:3777289/36‑1 (MQ=255)
                                                                       
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT...................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATGA  <  REL606/1733649‑1733614

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