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.590)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

CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682
                                   |                                   
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     >  1:6026420‑M2/36‑36 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:4243758‑M2/1‑1 (MQ=255)
tccatttcatcgtccaacaaaatgggtgcagtacaa                                     <  1:4135358‑M2/1‑1 (MQ=255)
catcatCTCCAAAATCGAAAACCAGGAAGGCCTcaa                                     <  1:115204/36‑1 (MQ=255)
catcatCTCCAAAATCGAAAACCAGGAAGGCCTcaa                                     <  1:4584580/36‑1 (MQ=255)
 ccatttcatcgtccaacaaaatgggtgcagtacaac                                    <  1:1640207‑M2/2‑1 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    >  1:4685161/1‑36 (MQ=255)
 atcatcTCCAAAATCGAAAACCAGGAAGGCCTcaac                                    <  1:373027/36‑1 (MQ=255)
  tcatcTCCAAAATCGAAAACCAGGAAGGCCTCAACA                                   >  1:957298‑M2/1‑36 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaacg                                   >  1:4396939‑M2/34‑35 (MQ=255)
  catttcatcgtccaacaaaatgggtgcagtacaaca                                   <  1:1508495‑M2/3‑1 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                  <  1:1879710/36‑1 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                  >  1:5773494/1‑36 (MQ=255)
   catcTCCAAAATCGAAAACCAGGAAGGCCTcaacaa                                  <  1:5696605/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:1957012/36‑1 (MQ=255)
        ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTcg                             <  1:6032196/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:4059761/1‑36 (MQ=255)
                gAAAACCAGGAAGGCCTCAACAACTTCGACGAAATc                     >  1:5837232/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:4258695‑M2/14‑36 (MQ=255)
                      tgggtgcagtacaACAACTTCGACGAAATCCTCGaa               <  1:3901353‑M2/23‑1 (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:5504536‑M2/28‑1 (MQ=255)
                           gcagtacaACAACTTCGACGAAATCCTCGAAGCctc          >  1:3549145‑M2/9‑36 (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:333326/36‑1 (MQ=255)
                             gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTg        >  1:6710837/1‑36 (MQ=255)
                             gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTg        <  1:61006/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)
                                   |                                   
CATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCA  >  REL606/1733612‑1733682

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