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

TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684
                                   |                                   
tcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACA                                     >  1:957298‑M2/1‑36 (MQ=255)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     <  1:1508495‑M2/3‑1 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                    >  1:5773494/1‑36 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                    <  1:5696605/36‑1 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                    <  1:1879710/36‑1 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                    >  1:7228221‑M2/1‑35 (MQ=255)
  atcTCCAAAATCGAAAACCAGGAAGGCCTCaacaac                                   >  1:1690344/1‑36 (MQ=255)
    ctcCAAAATCGAAAACCAGGAAGGCCTCAACAACtt                                 <  1:4372045/36‑1 (MQ=255)
      ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                               >  1:7398551‑M2/1‑30 (MQ=255)
      ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTcg                               <  1:6032196/36‑1 (MQ=255)
      ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTcg                               <  1:1957012/36‑1 (MQ=255)
       cAAAATCGAAAACCAGGAAGGCCTCAACAACTTcga                              >  1:1591585/1‑36 (MQ=255)
       cAAAATCGAAAACCAGGAAGGCCTCAACAACTTcga                              <  1:2392029/36‑1 (MQ=255)
       cAAAATCGAAAACCAGGAAGGCCTCAACAACTTcga                              >  1:1804404/1‑36 (MQ=255)
        aaaaTCGAAAACCAGGAAGGCCTCAGCAACTTcgac                             >  1:2800362/1‑36 (MQ=255)
          aaTCGAAAACCAGGAAGGCCTCAACAACTTcgacga                           >  1:5499198/1‑36 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAACTTCGACGaa                          >  1:6063038/1‑36 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAACTTCGACGaa                          >  1:550331/1‑36 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAACTGCGACGaa                          >  1:7069235/1‑36 (MQ=255)
           aTCGAAAACCAGCAAGGCCTCAACAACTTCGACGaa                          >  1:2358238/1‑36 (MQ=255)
            tCGAAAACCAGGAAGGCCTCAACAACTTCGACGaaa                         >  1:5125060/1‑36 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAAt                        >  1:6268759/1‑36 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAACTTCGACGAAATc                       >  1:5837232/1‑36 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAACTTCGACGAAATc                       >  1:4059761/1‑36 (MQ=255)
              acaaaatgggtgcagtacaACAACTTCGACGAAATc                       <  1:7311324‑M2/17‑1 (MQ=255)
               aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      >  1:2674145‑M2/1‑21 (MQ=255)
               aaaaCCAGGAAGGCCTCAACAACTTCGACGAAATcc                      >  1:2290896/1‑36 (MQ=255)
                aaaatgggtgcagtacaACAACTTAGACGAAATCCt                     >  1:4768149‑M2/18‑36 (MQ=255)
                aaaCCAGGAAGGCCTCAACAACTTCGACGAAATCCt                     >  1:4922321/1‑36 (MQ=255)
                aaaCCAGGAAGGCCTCAACAACTTCGACGAAATCCt                     >  1:3277200/1‑36 (MQ=255)
                 aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTc                    <  1:5193656/36‑1 (MQ=255)
                 aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTc                    >  1:5481224/1‑36 (MQ=255)
                  aCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCg                   >  1:5168515/1‑36 (MQ=255)
                   ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                  >  1:3884384‑M2/1‑17 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  1:3901353‑M2/23‑1 (MQ=255)
                    tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 >  1:4258695‑M2/14‑36 (MQ=255)
                     aGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAg                >  1:2079999/1‑36 (MQ=255)
                      ggtgcagtacaACAACTTCGACGAAATCCTCGAAGc               >  1:4314115‑M2/12‑36 (MQ=255)
                       gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGcc              >  1:3563649/1‑36 (MQ=255)
                       gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGcc              >  1:2068454/1‑36 (MQ=255)
                        aaGGCCTTAACAACTTCGACGAAATCCTCGAAGCct             >  1:1844831/1‑36 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             <  1:1608468‑M2/36‑25 (MQ=255)
                        aaGGCCTCAACAACTTCGACGAAATCCTCGAAGCct             <  1:4862591/36‑1 (MQ=255)
                         gcagtacaACAACTTCGACGAAATCCTCGAAGCctc            >  1:3549145‑M2/9‑36 (MQ=255)
                         gcagtacaACAACTTCGACGAAATCCTCGAAGCctc            <  1:5504536‑M2/28‑1 (MQ=255)
                         aGGCCTCAACAtgtactgaccccaaaaagttggaca            >  1:3843735‑M2/1‑11 (MQ=255)
                          ggCCTCAACAACTTCGACGAAATCCTCGAAGCctct           <  1:1432745/36‑1 (MQ=255)
                          ggCCTCAACAACTACGACGAAATCCTCGAAGCctct           <  1:937434/36‑1 (MQ=255)
                           gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTg          <  1:61006/36‑1 (MQ=255)
                           gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTg          >  1:6710837/1‑36 (MQ=255)
                           gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTg          <  1:333326/36‑1 (MQ=255)
                            gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa         >  1:5748832‑M2/6‑36 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggacagtt         <  1:5148641‑M2/36‑29 (MQ=255)
                              tCAACAACTTCGACGAAATCCTCGAAGCCTCTGACg       <  1:3172826/36‑1 (MQ=255)
                               caACAACTTCGACGAAATCCTCGAAGCCTCTGACgg      >  1:1746479‑M2/3‑36 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     <  1:5196762‑M2/35‑1 (MQ=255)
                                aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     >  1:1685282‑M2/2‑36 (MQ=255)
                                 ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa    >  1:5738751‑M2/1‑36 (MQ=255)
                                 ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa    >  1:4046027‑M2/1‑36 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  >  1:5104695/1‑36 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  <  1:6787779/36‑1 (MQ=255)
                                   aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  >  1:153710/1‑36 (MQ=255)
                                   |                                   
TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684

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