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

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT..................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG  <  REL606/1733649‑1733615
                                                                      
 GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATG                                   <  1:5175378/36‑1 (MQ=1)
   GTTTAACTGTCCACCTTTTTGGGGTCAGTACATGTT                                 <  1:1195056/36‑1 (MQ=0)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATGTT                                 >  1:912950/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATGTT                                 >  1:5722040/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATGTT                                 >  1:3099488/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACATGTT                                 >  1:2237493/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                 <  1:4863894/36‑1 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                 >  1:4811070/1‑36 (MQ=1)
   GTTTAACTGTCCAACTTTTTGGGGTCAGTACAGGTT                                 <  1:2534997/36‑1 (MQ=1)
   GTTTAACTGTCCAACCTTTTGGGGTCAGTACATGTT                                 <  1:7415064/36‑1 (MQ=0)
     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                               >  1:4499705/1‑36 (MQ=1)
     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                               <  1:1397773/36‑1 (MQ=1)
     TTAACTGTCCAACTTTTTGGGGTCAGTACATGTTGA                               >  1:6819818/1‑36 (MQ=1)
     TTAACTGTCCAACTTTTTGGGGTCAGTACAGGTTGc                               <  1:2563011/36‑2 (MQ=35)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                              <  1:1024346/36‑1 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                              <  1:5503595/36‑1 (MQ=1)
      TAACTGTCCAACTTTTTGGGGTCAGTACATGTTGAG                              <  1:1498665/36‑1 (MQ=1)
       AACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGG                             >  1:5863688/1‑36 (MQ=1)
        ACTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGC                            >  1:4369694/1‑36 (MQ=1)
         CTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCC                           <  1:618903/36‑1 (MQ=1)
         CTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCC                           >  1:5595487/1‑36 (MQ=1)
         CTGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCC                           <  1:7355592/36‑1 (MQ=1)
          TGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCT                          >  1:6160010/1‑36 (MQ=1)
          TGTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCT                          <  1:5791549/36‑1 (MQ=1)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                         <  1:363499/36‑1 (MQ=1)
           GTCCAACTTTTTGGGGTCAGTACATGTTGAGGCCTT                         <  1:7013053/36‑1 (MQ=1)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                       <  1:5638919/36‑1 (MQ=1)
             CCAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCC                       <  1:1767742/36‑1 (MQ=1)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                      <  1:676429/36‑1 (MQ=1)
              CAACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCT                      <  1:1740802/36‑1 (MQ=1)
               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                     <  1:2208025/36‑1 (MQ=1)
               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                     <  1:124213/36‑1 (MQ=1)
               AACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTG                     >  1:3793884/1‑36 (MQ=1)
                ACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGG                    >  1:3465141/1‑36 (MQ=1)
                ACTTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGG                    >  1:5067228/1‑36 (MQ=1)
                  TTTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTT                  <  1:377793/36‑1 (MQ=1)
                   TTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTT                 >  1:3335562/1‑36 (MQ=1)
                   TTTTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTT                 >  1:5081604/1‑36 (MQ=1)
                     TTGGGGTCAGTGCATGTTGAGGCCTTCCTGGTTTTC               <  1:1149210/36‑1 (MQ=0)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC               >  1:6836706/1‑36 (MQ=1)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC               >  1:4310454/1‑36 (MQ=1)
                     TTGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTC               <  1:5025707/36‑1 (MQ=1)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG              <  1:1480102/36‑1 (MQ=32)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG              <  1:1713749/36‑1 (MQ=32)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG              <  1:7462026/36‑1 (MQ=32)
                      TGGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCG              <  1:7154245/36‑1 (MQ=32)
                       GGGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGA             >  1:4552755/1‑36 (MQ=35)
                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT            >  1:598156/1‑36 (MQ=38)
                        GGGTCAGTACATGTTGAGGCCTTCCTGGTTTTCGAT            >  1:848609/1‑36 (MQ=38)
                          GTCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTT          <  1:886220/36‑1 (MQ=255)
                           TCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTT         <  1:6198348/36‑1 (MQ=255)
                           TCAGTACATGTTGAGGCCTTCCTGGTTTTCGATTTT         <  1:214213/36‑1 (MQ=255)
                             AGTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGG       <  1:7081885/36‑1 (MQ=255)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA      >  1:4231114/1‑36 (MQ=255)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA      <  1:7767779/36‑1 (MQ=255)
                              GTACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGA      <  1:2111451/36‑1 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA    >  1:7193711/1‑36 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA    >  1:1442553/1‑36 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA    <  1:2265101/36‑1 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA    >  1:4958776/1‑36 (MQ=255)
                                ACATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGA    >  1:3962631/1‑36 (MQ=255)
                                  ATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG  >  1:6862688/1‑36 (MQ=255)
                                  ATGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG  >  1:2248579/1‑36 (MQ=255)
                                                                      
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAT..................................  >  REL606/666096‑666130
...................................TGTTGAGGCCTTCCTGGTTTTCGATTTTGGAGATG  <  REL606/1733649‑1733615

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG