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 =4 (0.130)29 (1.010) 26/146 0.2 92.2% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 1 (0.030)coding (331/1413 nt) pykF pyruvate kinase

AAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTT  >  REL606/1733221‑1733366
                                                                          |                                                                       
nagCTGGTCAGACCTTTACTTNCACCACNGATAATTCTGTTATCGGGAACAGCGAAATGGTTGCGGTAACGTATGa                                                                        >  1:805847‑M2/2‑75 (MQ=255)
cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                        <  1:483844‑M1/9‑1 (MQ=255)
 agacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatgaa                                                                       <  1:649075‑M1/10‑1 (MQ=255)
 aGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                                                                       >  1:398517‑M2/1‑74 (MQ=255)
   cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                     <  1:443191‑M2/76‑5 (MQ=255)
      tttactttcaccactgataaatctgatatcggcancagcgaaatggttgcgntgacgtatgAACGTATNAAGGttt                                                                  <  2:856292‑M1/15‑1 (MQ=255)
      gTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                                                                  >  2:363692‑M2/1‑69 (MQ=255)
       tCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatn                                                                 <  1:845701‑M2/76‑9 (MQ=255)
        tactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCa                                                                <  2:719861‑M1/17‑1 (MQ=255)
        cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                                <  1:483844‑M2/76‑10 (MQ=255)
         aGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaa                                                               <  1:649075‑M2/76‑11 (MQ=255)
            ttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactac                                                            >  1:719861‑M1/56‑76 (MQ=255)
            ttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactac                                                            <  2:704827‑M1/21‑1 (MQ=255)
             tcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactact                                                           >  2:18841‑M1/55‑76 (MQ=255)
              tttACTTTCACCACTGATAAATCTGATATCGGCANCAGCGAAATGGTTGCGNTGACGTATGaacgtatnaaggttt                                                          <  2:856292‑M2/76‑16 (MQ=255)
                tACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttca                                                        <  2:719861‑M2/76‑18 (MQ=255)
                ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAc                                                        <  1:227837‑M1/25‑1 (MQ=255)
                    ttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactac                                                    >  1:719861‑M2/1‑55 (MQ=255)
                    ttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactac                                                    <  2:704827‑M2/76‑22 (MQ=255)
                     tCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactact                                                   >  2:18841‑M2/1‑54 (MQ=255)
                      ataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtct                                                  >  2:467239‑M1/46‑76 (MQ=255)
                        ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgac                                                <  1:227837‑M2/76‑26 (MQ=255)
                           tctgttatcggcaacagcgaaatggttgaggtaacgtatgAACGTATGAAGGTTTCACTACAGACCTGTCTGTTgg                                             <  2:937548‑M1/36‑1 (MQ=255)
                             tgttatcggcaacagcgnaatggttgcggtaacgtatgAACGTATGAAGGTNTCACTAATGCCATGTCTGTTGGCa                                           <  2:793795‑M1/38‑1 (MQ=255)
                              gttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCaa                                          <  2:131528‑M1/39‑1 (MQ=255)
                              aTAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtct                                          >  2:467239‑M2/1‑45 (MQ=255)
                                  tcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAANTTTGACTACTGACCTGTCTGTTGGCAACAcc                                      >  1:78884‑M1/34‑76 (MQ=255)
                                   tCTGTTATCGGCAACAGCGAAATGGTTGAGGTAACGTATGaacgtatgaaggtttcactacagacctgtctgttgg                                     <  2:937548‑M2/76‑37 (MQ=255)
                                    ngcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGt                                    >  1:842273‑M1/32‑76 (MQ=255)
                                    ggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGt                                    >  2:122622‑M1/32‑76 (MQ=255)
                                     tGTTATCGGCAACAGCGNAATGGTTGCGGTAACGTATGaacgtatgaaggtntcactaatgccatgtctgttggca                                   <  2:793795‑M2/76‑39 (MQ=255)
                                     gcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTa                                   <  1:587952‑M1/46‑1 (MQ=255)
                                      gTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaa                                  <  2:131528‑M2/76‑40 (MQ=255)
                                      caacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTAc                                  <  1:43275‑M1/47‑1 (MQ=255)
                                       aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                                 <  1:797508‑M1/48‑1 (MQ=255)
                                       aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                                 <  2:410433‑M1/48‑1 (MQ=255)
                                       aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACATGTCTGTTGGCAACACCGTACt                                 >  1:509557‑M1/29‑76 (MQ=255)
                                          tCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaantttgactactgacctgtctgttggcaacacc                              >  1:78884‑M2/1‑33 (MQ=255)
                                            ngcAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgt                            >  1:842273‑M2/2‑31 (MQ=255)
                                            ggCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgt                            >  2:122622‑M2/1‑31 (MQ=255)
                                            cgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTg                            <  2:89474‑M1/53‑1 (MQ=255)
                                             gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgta                           <  1:587952‑M2/76‑47 (MQ=255)
                                              cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtac                          <  1:43275‑M2/76‑48 (MQ=255)
                                              aaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGAc                          <  2:966488‑M1/55‑1 (MQ=255)
                                              aaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGAc                          <  2:204652‑M1/55‑1 (MQ=255)
                                               aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                         <  2:410433‑M2/76‑49 (MQ=255)
                                               aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                         <  1:797508‑M2/76‑49 (MQ=255)
                                               aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacatgtctgttggcaacaccgtact                         >  1:509557‑M2/1‑28 (MQ=255)
                                                  ggttgcggtaacgtatgAACGTATGAGGGTTTAACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATg                      <  1:850270‑M1/59‑1 (MQ=255)
                                                  gcttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATg                      >  2:102717‑M1/18‑76 (MQ=255)
                                                    cGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttg                    <  2:89474‑M2/76‑54 (MQ=255)
                                                      aaaTGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgac                  <  2:966488‑M2/76‑56 (MQ=255)
                                                      aaaTGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgac                  <  2:204652‑M2/76‑56 (MQ=255)
                                                          ggTTGCGGTAACGTATGaacgtatgagggtttaactactgacctgtctgttggcaacaccgtactggttgacgatg              <  1:850270‑M2/76‑60 (MQ=255)
                                                          gcttgcGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatg              >  2:102717‑M2/3‑17 (MQ=255)
                                                            acgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACATCGTACTGGTTGACGATGGTCTGATCgg            >  1:89474‑M1/8‑76 (MQ=255)
                                                                   AACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGaa     >  2:194355‑M1/1‑76 (MQ=255)
                                                                   AACGTATGAAGGTTTCACTACTAACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGaa     >  1:10268‑M1/1‑76 (MQ=255)
                                                                    acgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacatcgtactggttgacgatggtctgatcgg    >  1:89474‑M2/1‑7 (MQ=255)
                                                                    ACGTATGAAGGTTTCACTACTGACCTTTCTGTTNGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAg    >  1:778529‑M1/1‑76 (MQ=255)
                                                                      gTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGtt  <  1:23341/76‑1 (MQ=255)
                                                                          |                                                                       
AAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTT  >  REL606/1733221‑1733366

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