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

GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588460‑588530
                                   |                                   
gccgtcagaggcttcgaggatttcgtcgaagttgtt                                     <  1:5196762‑M1/1‑1 (MQ=255)
gccgtcagaggcttcgaggatttcgtcgaagttgtt                                     >  1:1685282‑M1/36‑36 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     >  1:4030469/1‑36 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     <  1:1551076/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     <  1:5170205/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     <  1:2783102/36‑1 (MQ=255)
gAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGtt                                     <  1:5261368/36‑1 (MQ=255)
 ccgtcagaggcttcgaggatttcgtcgaagttgttg                                    >  1:1746479‑M1/35‑36 (MQ=255)
 aGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTg                                    <  1:3391046/36‑1 (MQ=255)
  gTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGt                                   <  1:7164192/36‑1 (MQ=255)
   tCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTa                                  >  1:4007076/1‑36 (MQ=255)
   tCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTa                                  >  1:2093832/1‑36 (MQ=255)
    tcagaggcttcgaggatttcgtcgaagttgtTGTAc                                 >  1:5748832‑M1/32‑36 (MQ=255)
    cATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTAc                                 <  1:4218120/36‑1 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                                <  1:772035/36‑1 (MQ=255)
     aTGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACt                                >  1:5575831/1‑36 (MQ=255)
      tGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTg                               >  1:6983784/1‑36 (MQ=255)
       gaggcttcgaggatttcgtcgaagttgtTGTACTGc                              <  1:5504536‑M1/8‑1 (MQ=255)
       gaggcttcgaggatttcgtcgaagttgtTGTACTGc                              >  1:3549145‑M1/29‑36 (MQ=255)
       gCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGc                              >  1:1936227/1‑36 (MQ=255)
       gCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGc                              >  1:5584980/1‑36 (MQ=255)
       gCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGc                              <  1:4468444/36‑1 (MQ=255)
       gCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGc                              >  1:355781/1‑36 (MQ=255)
       gCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGc                              >  1:4623193/1‑36 (MQ=255)
        cAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCa                             >  1:871454/1‑36 (MQ=255)
        cAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCa                             >  1:2829038/1‑36 (MQ=255)
         aGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAc                            >  1:5255882/1‑36 (MQ=255)
          gcttcgaggatttcgtcgaagttgtTGTACTGCAcc                           >  1:4314115‑M1/26‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTCCTGCAcc                           >  1:3323318/1‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           >  1:6074036/1‑36 (MQ=255)
          gAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCAcc                           >  1:284263/1‑36 (MQ=255)
           aaaaaacCGGGTTCAGCTTCAAGTTGTACTGCAccc                          <  1:6474076/36‑1 (MQ=255)
           aaaaaaaCAGGTTCAGCTTCAAGTTGTACTGCAccc                          >  1:4688463/1‑36 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                         <  1:3901353‑M1/13‑1 (MQ=255)
            ttcgaggatttcgtcgaagttgtTGTACTGCACCCa                         >  1:4258695‑M1/24‑36 (MQ=255)
            aaaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCa                         >  1:2583595/1‑36 (MQ=255)
            aaaaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCa                         >  1:1975267/1‑36 (MQ=255)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                       <  1:972343/36‑1 (MQ=255)
              aaaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtt                       <  1:2851944/36‑1 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                      <  1:6602991/36‑1 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                      <  1:5485063/36‑1 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                      >  1:344041/1‑36 (MQ=255)
               aaaCGGGTTCAGCTTCAAGTTGTACTGCACCCAttt                      >  1:7003247/1‑36 (MQ=255)
                aggatttcgtctaagttgtTGTACTGCACCCAtttt                     >  1:4768149‑M1/20‑36 (MQ=255)
                aaCGGGTTCAGCTTCAAGTTGTACTGCACCCAtttt                     >  1:5903371/1‑36 (MQ=255)
                 aCGGGTTCAGCTTCAAGTTGTACTGCACCCATTttg                    >  1:3414802/1‑36 (MQ=255)
                  gatttcgtcgaagttgtTGTACTGCACCCATTttgt                   <  1:7311324‑M1/19‑1 (MQ=255)
                  cGGGTTCAGCTTCAAGTTGTACTGCACCCATTttgt                   >  1:3448277/1‑36 (MQ=255)
                   gggTTCAGCTTCAAGTTGTACTGCACCCATTttgtt                  >  1:2150398/1‑36 (MQ=255)
                   gggTTCAGCTTCAAGTTGTACTGCACCCATTttgtt                  <  1:2566923/36‑1 (MQ=255)
                    ggTTCAGCTTCAAGTTGTACTGCACCCATTttgttg                 >  1:6153228/1‑36 (MQ=255)
                     gTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTgg                >  1:6697971/1‑36 (MQ=255)
                     gTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTgg                <  1:5710602/36‑1 (MQ=255)
                     gTTCAGCTTCAAGTTGTACTGCACCAATTtttttat                >  1:4195875/1‑34 (MQ=25)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           <  1:1465385/36‑1 (MQ=255)
                          gCTTCAAGTTGTACTGCACCCATTTTGTTGGACGAt           <  1:5966054/36‑1 (MQ=255)
                            ttCAAGTTGTACTGCACCCATTTTGTTGGACGATGa         <  1:7196062/36‑1 (MQ=255)
                            ttCAAGTTGTACTGCACCCATTTTGTTGGACGATGa         <  1:6801628/36‑1 (MQ=255)
                            ttCAAGTTGTACTGCACCCATTTTGTTGGACGATGa         <  1:384068/36‑1 (MQ=255)
                            ttCAAGTTGTACTGCACCCATTTTGTTGGACGATGa         <  1:812874/36‑1 (MQ=255)
                             tCAAGTTGTACTGCACCCATTTTGTTGGACGATGaa        >  1:7050029/1‑36 (MQ=255)
                              cAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa       >  1:4837487/1‑36 (MQ=40)
                              cAAGTTGTACTGCACCCATTTTGTTGGACGATGaaa       <  1:7504958/36‑1 (MQ=40)
                               aaGTTGTACTGCACCCATTTTGTTGGACGATGAAAt      >  1:6916749/1‑36 (MQ=38)
                               aaGTTGTACTGCACCCATTTTGTTGGACGATGAAAt      >  1:1404889/1‑36 (MQ=38)
                               aaGTTGTACTGCACCCATCTTGTTGGACGATGAAAt      >  1:3417385/1‑36 (MQ=255)
                                tgtTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:1508495‑M1/33‑1 (MQ=255)
                                cgtTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:4396939‑M1/4‑36 (MQ=255)
                                cagTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:1824451‑M1/33‑1 (MQ=255)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:4095914/36‑1 (MQ=35)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     <  1:3862668/36‑1 (MQ=35)
                                aGTTGTACTGCACCCATTTTGTTGGACGATGAAATg     >  1:7450188/1‑36 (MQ=35)
                                 gtTGTACTGCACCCATTTTGTTGGACGATGAAATgg    <  1:1640207‑M1/34‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:4135358‑M1/35‑1 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   >  1:6026420‑M1/2‑36 (MQ=255)
                                  tTGTACTGCACCCATTTTGTTGGACGATGAAATGGa   <  1:4243758‑M1/35‑1 (MQ=255)
                                   TGTACTGCATCCATTTTGTTGGACGATGAAATGGaa  <  1:1517688‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGTACGATGAAATGGaa  >  1:1809295‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:455503‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  <  1:4810476‑M1/36‑1 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:5234743‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:2746268‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:183081‑M1/1‑36 (MQ=255)
                                   TGTACTGCACCCATTTTGTTGGACGATGAAATGGaa  >  1:1670066‑M1/1‑36 (MQ=255)
                                   |                                   
GAGTCATGCAGAAAAAAACGGGTTCAGCTTCAAGTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAA  >  REL606/588460‑588530

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG