Predicted mutation
evidence seq id position mutation freq annotation gene description
JC REL606 1,733,288 8 bp x 2 100% duplication pykF → pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1733288 =60 (1.110)14 (0.260) 11/66 NT 17.8% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 69 (1.200)coding (331/1413 nt) pykF pyruvate kinase

GCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT...................................  <  REL606/1733324‑1733288
...................................TTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA  <  REL606/1733295‑1733261
                                                                        
  CAACAGACAGGTCAGTAGTGAAACCTTCATACGTTC                                    >  1:45104/1‑36 (MQ=255)
  CAACAGACAGGTCAGTAGTGAAACCCTCATACGTTC                                    >  1:7050129/1‑36 (MQ=255)
      AGACAGGTCAGTAGTGAAACCTTCATACGTTCATAC                                <  1:4302961/36‑1 (MQ=255)
       GACAGGTCAGTAGTGAAACCTTCATACGTTCATACG                               >  1:996711/1‑36 (MQ=255)
        ACAGGTCAGTAGTGAAACCTTCATACGTTCATACGT                              >  1:5666042/1‑36 (MQ=255)
                GTAGTGAAACCTTCATACGTTCATACGTTACCGCAA                      <  1:6802879/36‑1 (MQ=255)
                       AACCTTCATACGTTCATACGTTACCGCAACCATTTC               <  1:4320091/36‑1 (MQ=255)
                       AACCTTCATACGTTCATACGTTACCGCAACCATTTC               <  1:6638639/36‑1 (MQ=255)
                       AACCTTCATACGTTCATACGTTACCGCAACCATTTC               <  1:7401000/36‑1 (MQ=255)
                        ACCTTCATACGTTCATACGTTACCGCAACCATTTCG              <  1:6673123/36‑1 (MQ=255)
                          CTTCATACGTTCATACGTTACCGCAACCATTTCGCT            <  1:2711219/36‑1 (MQ=255)
                               TACGTTCATACGTTACCGCAACCATTTCGCTGTTGC       >  1:7108537/1‑36 (MQ=255)
                                ACGTTCATACGTTACCGCAACCATTTCGCTGTTGCC      >  1:2416581/1‑36 (MQ=255)
                                  GTTCATACGTTACCGCAACCATTTCGCTGTTGCCGA    >  1:4758208/1‑36 (MQ=255)
                                                                        
GCCAACAGACAGGTCAGTAGTGAAACCTTCATACGTT...................................  <  REL606/1733324‑1733288
...................................TTCATACGTTACCGCAACCATTTCGCTGTTGCCGATA  <  REL606/1733295‑1733261

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