Predicted mutation
evidence seq id position mutation annotation gene description
JC REL606 1,733,297 (CGTATGAA)1→2 coding (333/1413 nt) pykF → pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1733288 =3 (0.140)19 (0.900) 18/146 0.3 88.5% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 2 (0.100)coding (331/1413 nt) pykF pyruvate kinase

TTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGA  >  REL606/1733213‑1733362
                                                                           |                                                                          
ttCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTa                                                                            >  1:220448/1‑76 (MQ=255)
   tctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACg                                                                         <  1:634792/76‑1 (MQ=255)
    ctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGt                                                                        >  1:53344/1‑76 (MQ=255)
        cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                    <  2:305771‑M1/9‑1 (MQ=255)
        cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                    >  1:95211‑M1/68‑76 (MQ=255)
          gacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAg                                                                  <  1:32947‑M1/11‑1 (MQ=255)
              tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                              <  1:157211‑M1/15‑1 (MQ=255)
              tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                              <  1:157213‑M1/15‑1 (MQ=255)
                cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                            <  2:305771‑M2/76‑10 (MQ=255)
                cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                            >  1:95211‑M2/1‑67 (MQ=255)
                  gACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                                                          <  1:32947‑M2/76‑12 (MQ=255)
                  ctttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCact                                                          <  2:628610‑M1/19‑1 (MQ=255)
                   tttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCacta                                                         >  2:6130‑M1/57‑76 (MQ=255)
                      tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                      <  1:157211‑M2/76‑16 (MQ=255)
                      tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                      <  1:157213‑M2/76‑16 (MQ=255)
                          cTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcact                                                  <  2:628610‑M2/76‑20 (MQ=255)
                           tttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcacta                                                 >  2:6130‑M2/1‑56 (MQ=255)
                                 aatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGtt                                           <  1:172718‑M1/34‑1 (MQ=255)
                                        tatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAaca                                    <  2:606336‑M1/41‑1 (MQ=255)
                                         aaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgtt                                   <  1:172718‑M2/76‑35 (MQ=255)
                                            ggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGt                                >  2:149409‑M1/32‑76 (MQ=255)
                                              caacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTAc                              >  1:286512‑M1/30‑76 (MQ=255)
                                              caacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTAc                              <  1:85329‑M1/47‑1 (MQ=255)
                                                tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaaca                            <  2:606336‑M2/76‑42 (MQ=255)
                                                 cagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTgg                           >  1:438223‑M1/27‑76 (MQ=255)
                                                   gcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGtt                         >  2:526990‑M1/25‑76 (MQ=255)
                                                    ggCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgt                        >  2:149409‑M2/1‑31 (MQ=255)
                                                    cgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTg                        <  1:170140‑M1/53‑1 (MQ=255)
                                                    cgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTg                        >  1:284956‑M1/24‑76 (MQ=255)
                                                      cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtac                      >  1:286512‑M2/1‑29 (MQ=255)
                                                      cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtac                      <  1:85329‑M2/76‑48 (MQ=255)
                                                        atggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGa                    <  1:182326‑M1/57‑1 (MQ=255)
                                                         cAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactgg                   >  1:438223‑M2/1‑26 (MQ=255)
                                                           gCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggtt                 >  2:526990‑M2/1‑24 (MQ=255)
                                                            cGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttg                <  1:170140‑M2/76‑54 (MQ=255)
                                                            cGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttg                >  1:284956‑M2/1‑23 (MQ=255)
                                                             tgcggtaacgtatgAACGTATGAAGGCTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTc               <  2:571187‑M1/62‑1 (MQ=255)
                                                                aTGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacga            <  1:182326‑M2/76‑58 (MQ=255)
                                                                  taacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATc          <  2:166283‑M1/67‑1 (MQ=255)
                                                                     tGCGGTAACGTATGaacgtatgaaggcttcactactgacctgtctgttggcaacaccgtactggttgacgatggtc       <  2:571187‑M2/76‑63 (MQ=255)
                                                                          tAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatc  <  2:166283‑M2/76‑68 (MQ=255)
                                                                           |                                                                          
TTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGA  >  REL606/1733213‑1733362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.