Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 18.4% 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)9 (0.180) 9/62 NT 20.9% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 34 (0.680)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)7 (0.130) 7/68 NT 15.9% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 37 (0.640)coding (685/1413 nt) pykF pyruvate kinase

GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666097‑666165
                                 |                                   
gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtt                                   <  1:6203540/36‑1 (MQ=35)
 tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  >  1:7484753/1‑36 (MQ=38)
 tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:1483480/36‑1 (MQ=38)
 tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:1538565/36‑1 (MQ=38)
 tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:7292919/36‑1 (MQ=38)
  gtTTAACTGTCCAACTTTTTGGGGTCAGTACAtttt                                 >  1:2947199/1‑35 (MQ=38)
    ttAACTGTCCAACTTTTTGGGGTCAGTACATTTCAg                               <  1:187007/36‑1 (MQ=255)
     tAACTGTCCAACTTTTTGGGGTCAGTACATTTCAgt                              >  1:4243919/1‑35 (MQ=255)
      aaCTGTCCAACTTTTTGGGGTCAGTACAtgttgagg                             <  1:5148641‑M1/36‑9 (MQ=255)
      aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                             <  1:741369/36‑1 (MQ=255)
      aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                             <  1:6963487/36‑1 (MQ=255)
      aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                             <  1:63450/36‑1 (MQ=255)
       aCTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                            >  1:4681259/1‑36 (MQ=255)
        cTGTCCAACTTTTTGGGGTCAGTACAGTTCAGGGCt                           >  1:6825395/1‑36 (MQ=255)
         tGTCCAACTTTTTGGGGTCAGTACAtgttgaggcct                          >  1:3843735‑M1/1‑25 (MQ=255)
         tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                          <  1:7440747/36‑1 (MQ=255)
          gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         <  1:1608468‑M1/36‑13 (MQ=255)
          gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                         <  1:6384346/36‑1 (MQ=255)
           tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        <  1:4677710/36‑1 (MQ=255)
           tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        <  1:3798524/36‑1 (MQ=255)
             cAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                      <  1:5421185/36‑1 (MQ=255)
              aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCa                     >  1:4701408/1‑36 (MQ=255)
               aCTTTTTGGGGTCAGTACAtgttgaggccttcctgg                    >  1:3884384‑M1/1‑19 (MQ=255)
                  ttttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                 <  1:5377556/36‑1 (MQ=255)
                   tttGGGGTCAGTACAtgttgaggccttcctggtttt                >  1:2674145‑M1/1‑15 (MQ=255)
                   tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                <  1:1061315/36‑1 (MQ=255)
                   tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                <  1:6067438/36‑1 (MQ=255)
                    ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc               <  1:5353920/36‑1 (MQ=255)
                        ggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgcg           <  1:1098718/36‑1 (MQ=255)
                           cAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCtt        >  1:4582459/1‑36 (MQ=255)
                           cAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCtt        >  1:4084208/1‑36 (MQ=255)
                            aGTACAtgttgaggccttcctggttttcgattttgg       >  1:7398551‑M1/1‑6 (MQ=255)
                            aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt       >  1:6154236/1‑36 (MQ=255)
                            aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt       <  1:1758739/36‑1 (MQ=255)
                             gTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTa      >  1:4028150/1‑36 (MQ=255)
                               tcATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:22800/35‑1 (MQ=255)
                               aCATTTTAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:1646392/36‑1 (MQ=255)
                               aCATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:649951/36‑1 (MQ=255)
                                 atgttgaggccttcctggttttcgattttggagatg  >  1:7228221‑M1/1‑1 (MQ=255)
                                 |                                   
GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666097‑666165

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