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 =2 (0.050)33 (0.900) 28/146 0.4 91.7% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 4 (0.110)coding (331/1413 nt) pykF pyruvate kinase

TCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733216‑1733363
                                                                        |                                                                           
ctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacg                                                                          <  1:813556‑M1/4‑1 (MQ=255)
 ctGAAAGCTGGTCAGGCCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTNCGGTAACGt                                                                         <  1:34509/76‑1 (MQ=255)
    aaaGCCTGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATg                                                                      <  1:268523/76‑1 (MQ=255)
     cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                     >  2:562345‑M1/68‑76 (MQ=255)
     cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                     <  1:326830‑M1/9‑1 (MQ=255)
        cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                  <  1:813556‑M2/76‑5 (MQ=255)
           tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                               >  1:856687‑M1/62‑76 (MQ=255)
           tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                               >  2:762910‑M1/62‑76 (MQ=255)
             cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                             >  2:562345‑M2/1‑67 (MQ=255)
             cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                             <  1:326830‑M2/76‑10 (MQ=255)
              actttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCac                                                            <  1:960052‑M1/18‑1 (MQ=255)
                   tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                       >  2:762910‑M2/1‑61 (MQ=255)
                   tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                       >  1:856687‑M2/1‑61 (MQ=255)
                   caccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTg                                                       <  1:696910‑M1/23‑1 (MQ=255)
                   caccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTg                                                       <  1:656765‑M1/23‑1 (MQ=255)
                     ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAc                                                     >  1:91241‑M1/52‑76 (MQ=255)
                      cactgataaatctgttatcggcaacagcgacatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAcc                                                    <  1:1095871‑M1/26‑1 (MQ=255)
                      aCTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcac                                                    <  1:960052‑M2/76‑19 (MQ=255)
                        ctgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctg                                                  <  1:810431‑M1/28‑1 (MQ=255)
                         tgataaatctgttatcggcaacatcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgt                                                 <  1:721358‑M1/29‑1 (MQ=255)
                         tgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgt                                                 >  1:518396‑M1/48‑76 (MQ=255)
                          nataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtc                                                <  1:1024103‑M1/30‑1 (MQ=255)
                           caccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactg                                               <  1:656765‑M2/76‑24 (MQ=255)
                           caccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactg                                               <  1:696910‑M2/76‑24 (MQ=255)
                           ataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtct                                               >  2:563124‑M1/46‑76 (MQ=255)
                             ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgac                                             >  1:91241‑M2/1‑51 (MQ=255)
                             aaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctgt                                             <  1:1070606‑M1/33‑1 (MQ=255)
                              cacTGATAAATCTGTTATCGGCAACAGCGACATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacc                                            <  1:1095871‑M2/76‑27 (MQ=255)
                              aatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGtt                                            >  1:467945‑M1/43‑76 (MQ=255)
                               atctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTg                                           <  1:238603‑M1/35‑1 (MQ=255)
                                tctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTgg                                          <  1:1011293‑M1/36‑1 (MQ=255)
                                cTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctg                                          <  1:810431‑M2/76‑29 (MQ=255)
                                 tGATAAATCTGTTATCGGCAACATCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgt                                         <  1:721358‑M2/76‑30 (MQ=255)
                                 tGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgt                                         >  1:518396‑M2/1‑47 (MQ=255)
                                  tgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCa                                        >  2:78046‑M1/39‑76 (MQ=255)
                                  natAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtc                                        <  1:1024103‑M2/75‑31 (MQ=255)
                                   ttcatcggcaacancgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCaa                                       >  1:163634‑M1/38‑76 (MQ=255)
                                   aTAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtct                                       >  2:563124‑M2/1‑45 (MQ=255)
                                    ntatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACNTGTCTGTTGGCAac                                      <  2:169820‑M1/40‑1 (MQ=255)
                                     tatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAaca                                     >  1:356033‑M1/36‑76 (MQ=255)
                                     aaaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgt                                     <  1:1070606‑M2/76‑34 (MQ=255)
                                      aaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgtt                                    >  1:467945‑M2/1‑42 (MQ=255)
                                       aTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttg                                   <  1:238603‑M2/76‑36 (MQ=255)
                                        tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttgg                                  <  1:1011293‑M2/76‑37 (MQ=255)
                                          tGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggca                                >  2:78046‑M2/1‑38 (MQ=255)
                                           ttcatcgGCAACANCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaa                               >  1:163634‑M2/4‑37 (MQ=255)
                                            ntaTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacntgtctgttggcaac                              <  2:169820‑M2/75‑41 (MQ=255)
                                            aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                              <  2:322325‑M1/48‑1 (MQ=255)
                                             tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaaca                             >  1:356033‑M2/1‑35 (MQ=255)
                                               agcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGt                           <  2:1176271‑M1/51‑1 (MQ=255)
                                                    aatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACg                      >  1:663096‑M1/21‑76 (MQ=255)
                                                    aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                      <  2:322325‑M2/76‑49 (MQ=255)
                                                       ggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATg                   <  2:342019‑M1/59‑1 (MQ=255)
                                                       aGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggt                   <  2:1176271‑M2/76‑52 (MQ=255)
                                                           gcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCt               <  1:942686‑M1/63‑1 (MQ=255)
                                                            aaTGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacg              >  1:663096‑M2/1‑20 (MQ=255)
                                                               taacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATc           >  2:832183‑M1/10‑76 (MQ=255)
                                                               ggTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatg           <  2:342019‑M2/76‑60 (MQ=255)
                                                                aacgtatgAACGTATGAAGGTTTCACTACTGTCCTGTCTCNNGGCAACACCGNNCTGGTTGACGATGGTCTGATCg          <  1:1087750‑M1/68‑1 (MQ=255)
                                                                  cgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGt        <  2:702037‑M1/70‑1 (MQ=255)
                                                                   gtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTTGTTGACGATGGTCTGATCGGTa       >  2:530823‑M1/6‑76 (MQ=255)
                                                                   gCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtct       <  1:942686‑M2/76‑64 (MQ=255)
                                                                     atgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATg     <  1:947641‑M1/73‑1 (MQ=255)
                                                                       tAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatc   >  2:832183‑M2/1‑9 (MQ=255)
                                                                        aacgtatgaacgtatgaaggtttcactactgtcctgtctcnnggcaacaccgnnctggttgacgatggtctgatcg  <  1:1087750‑M2/76‑69 (MQ=255)
                                                                        |                                                                           
TCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733216‑1733363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ 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.