Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 94.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 (1.010) 26/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)56 (1.280) 35/68 NT 90.3% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 6 (0.130)coding (685/1413 nt) pykF pyruvate kinase

TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588529‑588495
.................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683
                                                                        
   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   >  1:4471506/1‑36 (MQ=1)
   ATTTCATCGTCCAACAAAATGGGTGCAGTACAACAA                                   >  1:3997097/1‑36 (MQ=1)
    TTTCATCGTCCAACAAAATGGGTGCAGTACAACAAC                                  >  1:936180/1‑36 (MQ=1)
     TTCATCGTCCAACAAAATGCGTGCAGTACAACAACT                                 >  1:4190529/1‑36 (MQ=0)
      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:5134107/1‑36 (MQ=1)
      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:4625495/1‑36 (MQ=1)
      TCATCGTCCAACAAAATGGGTGCAGTACAACAACTT                                >  1:2531073/1‑36 (MQ=1)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:1937326/36‑1 (MQ=1)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:734171/36‑1 (MQ=1)
       CATCGTCCAACAAAATGGGTGCAGTACAACAACTTC                               <  1:5376590/36‑1 (MQ=1)
        ATCGTCCAACAAAATGGGTGCAGTACAACAACTTCG                              >  1:289791/1‑36 (MQ=32)
         TCGTCCAACAAAATGGGTGCAGTACAACAACTTCGA                             <  1:25505/36‑1 (MQ=35)
           GTCCAACAAAATGGGTGCAGTACAACAACTTCGACG                           >  1:5166559/1‑36 (MQ=40)
             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:2992076/36‑1 (MQ=255)
             CCAACAAAATGGGTGCAGTACAACAACTTCGACGAA                         <  1:3649576/36‑1 (MQ=255)
               AACAAAATGGGTGCAGTACAACAACTTCGACGAAAT                       >  1:5531616/1‑36 (MQ=255)
                ACAAAATGGGTGCAGTACAACAACTTCGACGAAATC                      >  1:238384/1‑36 (MQ=255)
                 CAAAATGGGTGCAGTACAACAACTTCGACGAAATCC                     >  1:3213291/1‑36 (MQ=255)
                  AAAATGGGTGCAGTACAACAACTTTGACGAAATCCT                    >  1:409887/1‑36 (MQ=255)
                  AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    <  1:3029847/36‑1 (MQ=40)
                  AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    <  1:4492684/36‑1 (MQ=40)
                  AAAATGGGTGCAGTACAACAACTTCGACGAAATCCT                    <  1:2107502/36‑1 (MQ=40)
                   AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   >  1:3507338/1‑36 (MQ=38)
                   AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   >  1:1580148/1‑36 (MQ=38)
                   AAATGGGTGCAGTACAACAACTTCGACGAAATCCTC                   >  1:1804980/1‑36 (MQ=38)
                     ATGGGTGCAGTACAACAACTTCGACGAAATCCTCGA                 >  1:5920655/1‑36 (MQ=32)
                     ATGGGTGCAGTACAACAACTTCGACGAAATCCTCGA                 <  1:485502/36‑1 (MQ=32)
                     ATGGGTGCAGTACAACAACTTCGACGAAAGCCTCGA                 >  1:5448741/1‑36 (MQ=12)
                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                <  1:3201918/36‑1 (MQ=1)
                      TGGGTGCAGTACAACAACTTCGACGAAATCCTCGAA                >  1:1181169/1‑36 (MQ=1)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             >  1:1069879/1‑36 (MQ=1)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             >  1:4214161/1‑36 (MQ=1)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             <  1:3417889/36‑1 (MQ=1)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             <  1:2638503/36‑1 (MQ=1)
                         GTGCAGTACAACAACTTCGACGAAATCCTCGAAGCC             <  1:1892253/36‑1 (MQ=1)
                           GCAGTACAACAACTTCGACGAAATCCTCGAAGCCTC           <  1:4845299/36‑1 (MQ=1)
                            CAGTACAACAACTTCGACGAAATCCTCGAAGCCTCT          <  1:509864/36‑1 (MQ=1)
                             AGTACAACAACTTCGACGAAATCCTCGAAGCCTCTG         >  1:2195698/1‑36 (MQ=1)
                              GTACAACAACTTCGACGAAATCCTCGAAGCCTCTGA        >  1:4038898/1‑36 (MQ=1)
                               TACAACAACTTCGACGAAATCCTCGAAGCCTCTGAC       <  1:57344/36‑1 (MQ=1)
                                                                        
TCCATTTCATCGTCCAACAAAATGGGTGCAGTACAAC...................................  <  REL606/588529‑588495
.................................CAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAT  >  REL606/1733647‑1733683

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG