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

ATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588531‑588495
...................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683
                                                                          
    CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    >  1:4465788/1‑36 (MQ=1)
    CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    <  1:4125444/36‑1 (MQ=1)
     ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   <  1:5861480/36‑1 (MQ=1)
      TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  >  1:7620310/1‑36 (MQ=1)
      TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  <  1:7656018/36‑1 (MQ=1)
      TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  >  1:3847570/1‑36 (MQ=1)
      TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  >  1:3767338/1‑36 (MQ=1)
       TTCATCGTCCAACAAAATGGGTGCAGTACAACAACT                                 <  1:4918491/36‑1 (MQ=1)
       TTCATCGTCCAACAAAATGGGTGCAGTACAACAACT                                 <  1:3901854/36‑1 (MQ=1)
        TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:1420488/1‑36 (MQ=1)
         CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:2807239/36‑1 (MQ=1)
         CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:317594/36‑1 (MQ=1)
          ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              <  1:499118/36‑1 (MQ=1)
          ATCGTCCAACAAAATGGGGGCAGTACAACAACTTCG                              >  1:5236970/1‑36 (MQ=0)
           TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             <  1:2738075/36‑1 (MQ=1)
           TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             >  1:6824615/1‑36 (MQ=1)
           TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             >  1:6760903/1‑36 (MQ=1)
           TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             >  1:6679846/1‑36 (MQ=1)
            CGTCCAACAAAATTGGTGCAGTACAACAACTTCGAC                            <  1:6486538/36‑1 (MQ=0)
             GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           >  1:4853182/1‑36 (MQ=1)
               CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:3950285/36‑1 (MQ=1)
               CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:1788330/36‑1 (MQ=1)
                 AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       >  1:3589778/1‑36 (MQ=1)
                 AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       <  1:2227349/36‑1 (MQ=1)
                 AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       <  1:5401105/36‑1 (MQ=1)
                  ACAAAATGGGTGCAGTACAACAACTTCGACGAAATC                      <  1:447752/36‑1 (MQ=1)
                  ACAAAATGGGTGCAGTACAACAACTTCGACGAAATC                      <  1:6075137/36‑1 (MQ=1)
                   CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     >  1:5386949/1‑36 (MQ=1)
                    AAAATGGGTGCAGTACAACAACTTCGCCGAAATCCT                    <  1:7797036/36‑1 (MQ=0)
                    AAAATGGGTGCAGTACAACAACTTCGACGAAATTCT                    >  1:3501915/1‑36 (MQ=0)
                    AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    >  1:3742230/1‑36 (MQ=1)
                     AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   >  1:4003030/1‑36 (MQ=1)
                     AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   <  1:3849578/36‑1 (MQ=1)
                      AATGGGTGCAGTACAACAACTTCGACGAAATCCTCG                  >  1:5815710/1‑36 (MQ=32)
                      AATGGGTGCAGTACAACAACTTCGACGAAATCCTCG                  >  1:5073817/1‑36 (MQ=32)
                       ATGGGTGCAGTACAACAACTTCGACGAAATCCTCGA                 <  1:2260805/36‑1 (MQ=35)
                       ATGGGTGCAGTACAACAACTTCGACGAAATCCTCGA                 >  1:1986145/1‑36 (MQ=35)
                        TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                <  1:487928/36‑1 (MQ=38)
                        TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                <  1:4201173/36‑1 (MQ=38)
                        TGGGTGCAGTAAAACAACTTCGACGAAATCCTCGAA                >  1:4560388/1‑36 (MQ=255)
                          GGTGCAGTACAATAACTTCGACGAAATCCTCGAAGC              <  1:5576135/36‑1 (MQ=255)
                          GGTGCAGTACAACAACTTCGACGAAATCCTCGAAGC              <  1:3194686/36‑1 (MQ=255)
                            aGCAGTACAACAACTGCGACGAAATCCTCGAAGCCT            <  1:2460156/35‑1 (MQ=255)
                             GCAGTACAACAACTTCGACGAAATCCTCGAAGCCTC           <  1:1383247/36‑1 (MQ=255)
                              CAGTAAAACAACTTCGACGAAATCCTCGAAGCCTCT          <  1:361535/36‑1 (MQ=255)
                               AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         <  1:7306423/36‑1 (MQ=255)
                               AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         >  1:5229963/1‑36 (MQ=255)
                               AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         >  1:7490265/1‑36 (MQ=255)
                               AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         <  1:3965450/36‑1 (MQ=255)
                                GTACAACAACTTCGACGAAATCCTCGAAGCCTCTGA        >  1:4505912/1‑36 (MQ=255)
                                 TACAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       >  1:7458637/1‑36 (MQ=255)
                                 TACAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       >  1:785115/1‑36 (MQ=255)
                                                                          
ATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588531‑588495
...................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG