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

TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588529‑588495
.................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683
                                                                        
  CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    <  1:7199417/36‑1 (MQ=1)
  CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    >  1:1476208/1‑36 (MQ=1)
   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   <  1:7950660/36‑1 (MQ=1)
   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   >  1:2428670/1‑36 (MQ=1)
   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   >  1:6297801/1‑36 (MQ=1)
   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   >  1:5029481/1‑36 (MQ=1)
    TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  <  1:1509913/36‑1 (MQ=1)
     TTCATCGTCCAACAAAATGGGTGCAGTACAACAACT                                 >  1:3739072/1‑36 (MQ=1)
     TTCATCGTCCAACAAAATGGGTGCAGTACAACAACT                                 <  1:8167947/36‑1 (MQ=1)
      TCATCGTCCAACACAATGGGTGCAGTACAACAACTg                                >  1:1213416/1‑35 (MQ=0)
      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:4643769/1‑36 (MQ=1)
      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:7745737/36‑1 (MQ=1)
      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:2633944/36‑1 (MQ=1)
      TCATCGTCCAACAAAATGGGCGCAGTACAACAACTT                                >  1:393194/1‑36 (MQ=0)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:3555139/36‑1 (MQ=1)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:5301234/36‑1 (MQ=1)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:5710148/36‑1 (MQ=1)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:1689617/36‑1 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              >  1:6990722/1‑36 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              <  1:2745405/36‑1 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              <  1:7541698/36‑1 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              >  1:2178080/1‑36 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              >  1:1935021/1‑36 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              >  1:2086658/1‑36 (MQ=1)
         TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             <  1:8284533/36‑1 (MQ=1)
          CGTCCAACAAAATGGGTGCAGTACAACAACTTCGAC                            <  1:4506375/36‑1 (MQ=1)
           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           <  1:5062364/36‑1 (MQ=1)
           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           <  1:2836694/36‑1 (MQ=1)
           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           >  1:1862846/1‑36 (MQ=1)
           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           <  1:2657975/36‑1 (MQ=1)
            TCCAACAAAATGGGTGCAGTACAACAACTTCGACGA                          >  1:7695062/1‑36 (MQ=1)
             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         >  1:5758703/1‑36 (MQ=1)
             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         >  1:7723644/1‑36 (MQ=1)
              CAACAAAATGGGTGCAGTACAACAACTTCGACGAAA                        >  1:5166251/1‑36 (MQ=1)
              CAAAAAAATGGGTGCAGTACAACAACTTCGACGAAA                        <  1:1799255/36‑1 (MQ=0)
               AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       <  1:1830335/36‑1 (MQ=1)
                aaAAAATGGGTGCAGTACAACAACTTCGACGAAATC                      <  1:2128551/34‑1 (MQ=1)
                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     >  1:3520749/1‑36 (MQ=1)
                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     >  1:1051747/1‑36 (MQ=1)
                  AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    >  1:3806704/1‑36 (MQ=1)
                   AAATGGGTGCATTACAACAACTTCGACGAAATCCTC                   <  1:519542/36‑1 (MQ=0)
                   AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   >  1:90215/1‑36 (MQ=1)
                     ATGGGTGCAGTACAACAACTTCGACGAAATCCTCGA                 >  1:2528471/1‑36 (MQ=32)
                     ATGGGTGCAGTACAACAACTTCGACGAAATCCTCGA                 >  1:4414119/1‑36 (MQ=32)
                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:4273253/1‑36 (MQ=1)
                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                <  1:4441841/36‑1 (MQ=1)
                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:4344987/1‑36 (MQ=1)
                      TGGGTGCAGTACAACAACTCCGACGAAATCCTCGAA                <  1:2392110/36‑1 (MQ=0)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             >  1:5562216/1‑36 (MQ=1)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             >  1:5807516/1‑36 (MQ=1)
                          TGCAGTACAACAACTTCGACGAAATCCTCGAAGCCT            >  1:5568504/1‑36 (MQ=1)
                           GCAGTACAACAACTTCGACGAAATCCTCGAAGCCTC           <  1:3277748/36‑1 (MQ=1)
                            CAGTACAACAACTTCGACGAAATCCTCGAAGCCTCT          >  1:7141649/1‑36 (MQ=1)
                             acTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         <  1:1137337/34‑1 (MQ=1)
                               TACAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       >  1:7042640/1‑36 (MQ=1)
                                ACAACAACTTCGACGGAATCCTCGAAGCCTCTGACG      >  1:4831023/1‑36 (MQ=0)
                                ACAACAACTTCGACGAAATCCTCGAAGCCTCTGACG      <  1:7620501/36‑1 (MQ=1)
                                ACAACAACTTCGACGAAATCCTCGAAGCCTCTGACG      >  1:6377352/1‑36 (MQ=1)
                                ACAACAACTTCGACGAAATCCTCGAAGCCTCTGACG      >  1:3456069/1‑36 (MQ=1)
                                ACAACAACTTCGACGAAATCCTCGAAGCCTCTGACG      >  1:2563/1‑36 (MQ=1)
                                                                        
TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588529‑588495
.................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683

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