Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 98.1% 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)40 (0.750) 32/62 NT 97.6% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 1 (0.020)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)63 (1.080) 33/68 NT 98.4% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 1 (0.020)coding (685/1413 nt) pykF pyruvate kinase

ATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588531‑588495
...................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683
                                                                          
    CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    <  1:1006445/36‑1 (MQ=1)
    CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    <  1:7774606/36‑1 (MQ=1)
    CATTTCATCGTCCAACAAAATGGGTGCAGTACAACA                                    >  1:1327762/1‑36 (MQ=1)
     ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   <  1:7612709/36‑1 (MQ=1)
     ATTTCATCGTCAAACAAAATGGGTGCAGTACAACAA                                   >  1:1667201/1‑36 (MQ=0)
        TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:978595/36‑1 (MQ=1)
        TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                <  1:4953085/36‑1 (MQ=1)
         CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:5314379/36‑1 (MQ=1)
         CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:3330933/36‑1 (MQ=1)
          ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              <  1:5830176/36‑1 (MQ=1)
           TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             >  1:1720457/1‑36 (MQ=1)
            CGTCCAACAAAATGGGTGCAGTACAACAACTTCGAC                            >  1:6665799/1‑36 (MQ=1)
              TCCAACAAAATGGGTGCAGTACAACAACTTCGACGA                          >  1:3592864/1‑36 (MQ=1)
               CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:5823721/36‑1 (MQ=1)
               CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:3833232/36‑1 (MQ=1)
                CAACAAAATGGGTGCAGTACAACAACTTCGACGAAA                        <  1:7461775/36‑1 (MQ=1)
                 AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       >  1:4214091/1‑36 (MQ=1)
                 AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       >  1:2695741/1‑36 (MQ=1)
                  ACAAAATGGGTGCAGTACAACAACTTCGACGAAATC                      >  1:665546/1‑36 (MQ=1)
                   CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     <  1:2564286/36‑1 (MQ=1)
                    AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    >  1:1282515/1‑36 (MQ=1)
                     AAATGGGTGCAGTACAACAACTTCGACGAAATCCac                   >  1:663926/1‑34 (MQ=1)
                     AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   <  1:3173666/36‑1 (MQ=1)
                     AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   <  1:5380533/36‑1 (MQ=1)
                     AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   >  1:7104671/1‑36 (MQ=1)
                      AATGGGTGCAGTACAACAACTTCGACGAAATCCTCG                  >  1:1677531/1‑36 (MQ=32)
                      AATGGGTGCAGTACAACAACTTCGACGAAATCCTCG                  <  1:6507697/36‑1 (MQ=32)
                        TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:2745889/1‑36 (MQ=38)
                        TGGGTGAAGTACAACAACTTCGACGAAATCCTCGAA                <  1:6265446/36‑1 (MQ=255)
                         GGGTGAAGTACAACAACTTCGACGAAATCCTCGAAG               <  1:123758/36‑1 (MQ=255)
                          GGTGCAGTACAACAACTTCGACGAAATCCTCGAAGC              >  1:308886/1‑36 (MQ=255)
                           GTGCAGTACAACAACTTCGACGAAATCCTCGGAGCC             >  1:6310895/1‑36 (MQ=255)
                           GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             >  1:2864339/1‑36 (MQ=255)
                           GTGCAGTACAACAACTCCGACGAAATCCTCGAAGCC             <  1:1450951/36‑1 (MQ=255)
                              CAGTACAACAACTTCGACGAAATCCTCGAAGCCTCT          >  1:5921742/1‑36 (MQ=255)
                               AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         >  1:3803752/1‑36 (MQ=255)
                                GTACAACAACTTCGACGAAATCCTCGAAGCCTCTGA        >  1:908203/1‑36 (MQ=255)
                                GTACAACAACTTCGACGAAAGCCTCGAAGCCTCTGA        <  1:2359024/36‑1 (MQ=255)
                                 TACAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       <  1:6711578/36‑1 (MQ=255)
                                  ACAACAACTTCGACGAAATCCTCGAAGCCTCTGACG      <  1:7077003/36‑1 (MQ=255)
                                                                          
ATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588531‑588495
...................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG