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

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1733579 =33 (0.570)27 (0.470) 22/70 NT 49.5% coding (615/1413 nt) pykF pyruvate kinase
?REL606 = 1733585 22 (0.380)coding (621/1413 nt) pykF pyruvate kinase

TGTGGATGTTTTCGCCGCCGTGCGCTTTCAGGTGC...................................  <  REL606/1733613‑1733579
...................................CAGGTGCTCACGGATTTCGATAACGTCAGAACGCT  <  REL606/1733585‑1733551
                                                                      
TGTGGATGTTTTCGCCGCCGTGCGCTTTCAGGTGCC                                    <  1:5522938/36‑1 (MQ=255)
 GTGGATGTTTTCGCCGCCGTGCGCTTTCAGGTGCCA                                   >  1:3014077/1‑36 (MQ=255)
 GTGGATGTTTTCGCCGCCGTGCGCTTTCAGGTGCCA                                   >  1:6303443/1‑36 (MQ=255)
    GATGTTTTCGCCGCCGTGCGCTTTCAGGTGCCAGGT                                <  1:6052798/36‑1 (MQ=255)
    GATGTTTTCGCCGCCGTGCGCTTTCAGGTGCCAGGT                                >  1:5350008/1‑36 (MQ=255)
        TTTTCGCCGCCGTGCGCTTTCAGGTGCCAGGTGCTC                            <  1:3566314/36‑1 (MQ=255)
          TTCGCCGCCGTGCGCTTTCAGGTGCCAGGTGCTCAC                          <  1:3710981/36‑1 (MQ=255)
           TCGCCGCCGTGCGCTTTCAGGTGCCAGGTGCTCACG                         >  1:6648471/1‑36 (MQ=255)
            CGCCGCCGTGCGCTTTCAGGTGCCAGGTGCTCACGG                        <  1:196453/36‑1 (MQ=255)
             GCCGCCGTGCGCTTTCAGGTGCCAGGTGCTCACGGA                       <  1:2829465/36‑1 (MQ=255)
              CCGCCGTGCGCTTTCAGGTGCCAGGTGCTCACGGAT                      >  1:964955/1‑36 (MQ=255)
              CCGCCGTGCGCTTTCAGGTGCCAGGTGCTCACGGAT                      >  1:1249450/1‑36 (MQ=255)
               CGCCGTGCGCTTTCAGGTGCCAGGTGCTCACGGATT                     >  1:3544126/1‑36 (MQ=255)
                GCCGTGCGCTTTCAGGTGCCAGGTGCTCACGGATTT                    <  1:3797703/36‑1 (MQ=255)
                GCCGTGCGCTTTCAGGTGCCAGGTGCTCACGGATTT                    <  1:5895799/36‑1 (MQ=255)
                  CGTGCGCTTTCAGGTGCCAGGTGCTCACGGATTTCG                  <  1:6030752/36‑1 (MQ=255)
                   GTGCGCTTTCAGGTGCCAGGTGCTCACGGATTTCGA                 >  1:5842308/1‑36 (MQ=40)
                    TGCGCTTTCAGGTGCCAGGTGCTCACGGATTTCGAT                <  1:4681567/36‑1 (MQ=38)
                    TGCGCTTTCAGGTGCCAGGTGCTCACGGATTTCGAT                <  1:5012000/36‑1 (MQ=38)
                    TGCGCTTTCAGGTGCCAGGTGCTCACGGATTTCGAT                >  1:6730858/1‑36 (MQ=38)
                     GCGCTTTCAGGTGCCAGGTGCTCACGGATTTCGATA               >  1:7467176/1‑36 (MQ=38)
                        aTTTCAGGTGCCAGGTGCTCACGGATTTCGATAACG            >  1:4585445/2‑36 (MQ=1)
                          TTCAGGTGCCAGGTGCTCACGGATTTCGATAACGTC          <  1:933320/36‑1 (MQ=1)
                             AGGTGCCAGGTGCTCACGGATTTCGATAACGTCAGA       <  1:3607571/36‑1 (MQ=1)
                              GGTGCCAGGTGCTCACGGATTTCGATAACGTCAGAA      >  1:5750833/1‑36 (MQ=1)
                                 GCCAGGTGCTCACGGATTTCGATAACGTCAGAACGC   <  1:3185511/36‑1 (MQ=1)
                                  CCAGGTGCTCACGGAGTTCGATAACGTCAGAACGCT  <  1:4607761/36‑1 (MQ=0)
                                                                      
TGTGGATGTTTTCGCCGCCGTGCGCTTTCAGGTGC...................................  <  REL606/1733613‑1733579
...................................CAGGTGCTCACGGATTTCGATAACGTCAGAACGCT  <  REL606/1733585‑1733551

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG