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

AAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAC  >  REL606/1733220‑1733368
                                                                           |                                                                         
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                                                            >  1:856687‑M2/1‑61 (MQ=255)
               tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                            >  2:762910‑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:696910‑M2/76‑24 (MQ=255)
                       caccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactg                                                    <  1:656765‑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)
                                                                      cgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggt     <  2:702037‑M2/76‑71 (MQ=255)
                                                                      cGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGt     <  2:81793/76‑1 (MQ=255)
                                                                      cGTATGAAGGTTTCACTACTGACCTGTATGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTAAGCAAGt     >  1:1188695/1‑76 (MQ=255)
                                                                       gtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtacttgttgacgatggtctgatcggta    >  2:530823‑M2/1‑5 (MQ=255)
                                                                         atgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatg  <  1:947641‑M2/76‑74 (MQ=255)
                                                                         aTGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAc  <  2:721358/76‑1 (MQ=255)
                                                                         aTGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAc  >  2:858409/1‑76 (MQ=255)
                                                                           |                                                                         
AAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAC  >  REL606/1733220‑1733368

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.