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

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT.................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT  <  REL606/1733649‑1733616
                                                                     
 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATG                                  >  1:1102539/1‑36 (MQ=32)
 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATG                                  <  1:4864568/36‑1 (MQ=32)
  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATGT                                 <  1:5700712/36‑1 (MQ=35)
  TGTTTAACTGTCCAACTTTTTGGGGTCAGTACATGT                                 <  1:5642141/36‑1 (MQ=35)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                <  1:1088683/36‑1 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                <  1:775052/36‑1 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                <  1:3560027/36‑1 (MQ=1)
    TTTAACTGTCCGACTTTTTGGGGTCAGTACATGTTG                               <  1:5144045/36‑1 (MQ=255)
    TTTAACTGTCCAACTTTTTGGTGTCAGTACATGTTG                               <  1:1540506/36‑1 (MQ=255)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACATGTTG                               >  1:1897360/1‑36 (MQ=38)
    TTTAACTGTCCAACTTTTTGGGGTCAGTACATGTTG                               <  1:3132601/36‑1 (MQ=38)
     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                              <  1:5311735/36‑1 (MQ=40)
     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                              <  1:4891742/36‑1 (MQ=40)
     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                              >  1:4447495/1‑36 (MQ=40)
     TTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGc                              <  1:3929976/36‑2 (MQ=35)
     TTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGc                              <  1:1131/36‑2 (MQ=35)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                             <  1:884641/36‑1 (MQ=255)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                             >  1:414481/1‑36 (MQ=255)
       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                            >  1:2972698/1‑36 (MQ=255)
       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                            >  1:3983917/1‑36 (MQ=255)
        ACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGC                           <  1:1864177/36‑1 (MQ=255)
        ACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGC                           <  1:2055804/36‑1 (MQ=255)
         CTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCC                          >  1:1127081/1‑36 (MQ=255)
          TGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCT                         <  1:905139/36‑1 (MQ=255)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                        <  1:5642823/36‑1 (MQ=255)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                        >  1:1755164/1‑36 (MQ=255)
            TCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTC                       <  1:224995/36‑1 (MQ=255)
            TCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTC                       >  1:4145372/1‑36 (MQ=255)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                      <  1:3436889/36‑1 (MQ=255)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                      <  1:4480476/36‑1 (MQ=255)
              CAACTTTTTGGGGTCAGTACGTGTTGAGGCCTTCCT                     <  1:4367657/36‑1 (MQ=255)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                     <  1:291479/36‑1 (MQ=255)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                     <  1:2270580/36‑1 (MQ=255)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                     <  1:5912946/36‑1 (MQ=255)
                 CTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGT                  <  1:5000256/36‑1 (MQ=255)
                  TTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTT                 <  1:678913/36‑1 (MQ=255)
                   TTTTGGGGTCAGTACATGTTGATGCCTTCCTGGTTT                >  1:3441423/1‑36 (MQ=255)
                   TTTTGGGGTCAGTACATGTTGATGCCTTCCTGGTTT                <  1:3732215/36‑1 (MQ=255)
                   TTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTT                <  1:675657/36‑1 (MQ=255)
                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT               <  1:4651056/36‑1 (MQ=255)
                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT               <  1:5221215/36‑1 (MQ=255)
                    TTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTT               <  1:568833/36‑1 (MQ=255)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC              >  1:941291/1‑36 (MQ=40)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC              <  1:5259580/36‑1 (MQ=40)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC              >  1:983997/1‑36 (MQ=40)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG             >  1:4555365/1‑36 (MQ=38)
                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT           <  1:1447377/36‑1 (MQ=38)
                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT           >  1:1886920/1‑36 (MQ=38)
                         GGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATT          <  1:3381031/36‑1 (MQ=38)
                         GGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATT          <  1:2686415/36‑1 (MQ=38)
                          GTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTT         >  1:2898888/1‑36 (MQ=38)
                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGG      >  1:2339346/1‑36 (MQ=38)
                               TACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGtg    >  1:3820527/1‑34 (MQ=25)
                               TACACGTTGAGGCCTTCCTGGTTTTCGATTTTGGAG    >  1:1229949/1‑36 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA   <  1:1789080/36‑1 (MQ=32)
                                 CATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT  <  1:3749872/36‑1 (MQ=1)
                                                                     
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT.................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGAT  <  REL606/1733649‑1733616

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG