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.050)43 (0.790) 37/146 0.3 92.5% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 4 (0.070)coding (331/1413 nt) pykF pyruvate kinase

AAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACC  >  REL606/1733220‑1733369
                                                                           |                                                                          
tcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatg                                                                            >  1:910796‑M1/69‑76 (MQ=255)
 cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                           >  2:756354‑M1/68‑76 (MQ=255)
 cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                           <  2:1098264‑M1/9‑1 (MQ=255)
 aaGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGa                                                                           >  2:1422802‑M2/1‑75 (MQ=255)
  aGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                                                                          <  1:638791‑M2/76‑3 (MQ=255)
    cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                        <  1:1398568‑M2/76‑5 (MQ=255)
      ggTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgta                                                                      <  2:1226724‑M2/76‑7 (MQ=255)
      ggTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgta                                                                      <  1:381892‑M2/76‑7 (MQ=255)
       gTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                                                                     >  1:724642‑M2/1‑69 (MQ=255)
       gTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtat                                                                     >  1:696073‑M2/1‑69 (MQ=255)
        ttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTc                                                                    <  2:1110206‑M1/16‑1 (MQ=255)
        ttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTc                                                                    <  2:1429132‑M1/16‑1 (MQ=255)
        tCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                                                                    >  1:910796‑M2/1‑68 (MQ=255)
         cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                                   <  2:1098264‑M2/76‑10 (MQ=255)
         cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                                   >  2:756354‑M2/1‑67 (MQ=255)
          actttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCac                                                                  <  2:416258‑M1/18‑1 (MQ=255)
             ttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactac                                                               <  2:1157034‑M1/21‑1 (MQ=255)
                ttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttc                                                            <  2:1110206‑M2/76‑17 (MQ=255)
                ttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttc                                                            <  2:1429132‑M2/76‑17 (MQ=255)
                accactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGa                                                            >  1:576497‑M1/53‑76 (MQ=255)
                 ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAc                                                           >  2:1549782‑M1/52‑76 (MQ=255)
                  aCTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcac                                                          <  2:416258‑M2/76‑19 (MQ=255)
                   actgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACct                                                         >  2:1553946‑M1/50‑76 (MQ=255)
                    ctgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctg                                                        <  2:170084‑M1/28‑1 (MQ=255)
                     ttCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactac                                                       <  2:1157034‑M2/76‑22 (MQ=255)
                       ataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtct                                                     <  1:676262‑M1/31‑1 (MQ=255)
                        accacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactga                                                    >  1:576497‑M2/1‑52 (MQ=255)
                         ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgac                                                   >  2:1549782‑M2/1‑51 (MQ=255)
                           acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacct                                                 >  2:1553946‑M2/1‑49 (MQ=255)
                            tctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTgg                                                <  2:754968‑M1/36‑1 (MQ=255)
                            tctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTgg                                                <  2:1237264‑M1/36‑1 (MQ=255)
                            cTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctg                                                <  2:170084‑M2/76‑29 (MQ=255)
                              tgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCa                                              >  1:1504218‑M1/39‑76 (MQ=255)
                               gttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCaa                                             >  1:1364205‑M1/38‑76 (MQ=255)
                               gttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCaa                                             <  1:1724592‑M1/39‑1 (MQ=255)
                               aTAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtct                                             <  1:676262‑M2/76‑32 (MQ=255)
                                ttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAac                                            <  1:191549‑M1/40‑1 (MQ=255)
                                ttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAac                                            <  2:914623‑M1/40‑1 (MQ=255)
                                ttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAac                                            <  2:965275‑M1/40‑1 (MQ=255)
                                  atcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAacac                                          <  1:1184003‑M1/42‑1 (MQ=255)
                                   tcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACAcc                                         <  2:201962‑M1/43‑1 (MQ=255)
                                    tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttgg                                        <  2:1237264‑M2/76‑37 (MQ=255)
                                    tCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttgg                                        <  2:754968‑M2/76‑37 (MQ=255)
                                    cggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCg                                        <  2:1738075‑M1/44‑1 (MQ=255)
                                     ggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGt                                       >  1:668817‑M1/32‑76 (MQ=255)
                                      tGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggca                                      >  1:1504218‑M2/1‑38 (MQ=255)
                                      gcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTa                                      <  2:287380‑M1/46‑1 (MQ=255)
                                       gTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaa                                     <  1:1724592‑M2/76‑40 (MQ=255)
                                       gTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaa                                     >  1:1364205‑M2/1‑37 (MQ=255)
                                       caacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTAc                                     <  2:618624‑M1/47‑1 (MQ=255)
                                        ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaac                                    <  2:965275‑M2/76‑41 (MQ=255)
                                        ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaac                                    <  2:914623‑M2/76‑41 (MQ=255)
                                        ttATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaac                                    <  1:191549‑M2/76‑41 (MQ=255)
                                        aacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACt                                    <  2:1547184‑M1/48‑1 (MQ=255)
                                          cagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTgg                                  >  1:1707557‑M1/27‑76 (MQ=255)
                                          aTCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacac                                  <  1:1184003‑M2/76‑43 (MQ=255)
                                           tCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacacc                                 <  2:201962‑M2/76‑44 (MQ=255)
                                            cGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccg                                <  2:1738075‑M2/76‑45 (MQ=255)
                                             ggCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgt                               >  1:668817‑M2/1‑31 (MQ=255)
                                              gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgta                              <  2:287380‑M2/76‑47 (MQ=255)
                                               cAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtac                             <  2:618624‑M2/76‑48 (MQ=255)
                                               aaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGAc                             >  1:424595‑M1/22‑76 (MQ=255)
                                                aaCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtact                            <  2:1547184‑M2/76‑49 (MQ=255)
                                                  cAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactgg                          >  1:1707557‑M2/1‑26 (MQ=255)
                                                    gttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATgg                        <  1:1560817‑M1/60‑1 (MQ=255)
                                                       aaaTGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgac                     >  1:424595‑M2/1‑21 (MQ=255)
                                                         ggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGa                   >  2:465079‑M1/12‑76 (MQ=255)
                                                           taacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATc                 >  2:15552‑M1/10‑76 (MQ=255)
                                                            gTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatgg                <  1:1560817‑M2/76‑61 (MQ=255)
                                                            aacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCg                <  2:671199‑M1/68‑1 (MQ=255)
                                                             acgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCgg               <  2:450443‑M1/69‑1 (MQ=255)
                                                              cgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGt              <  1:311989‑M1/70‑1 (MQ=255)
                                                               gtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTa             >  2:827985‑M1/6‑76 (MQ=255)
                                                                tatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTAt            >  2:1240329‑M1/5‑76 (MQ=255)
                                                                 ggTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctga           >  2:465079‑M2/1‑11 (MQ=255)
                                                                  tgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATgg          <  2:1666689‑M1/74‑1 (MQ=255)
                                                                   tAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatc         >  2:15552‑M2/1‑9 (MQ=255)
                                                                    aacgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcg        <  2:671199‑M2/76‑69 (MQ=255)
                                                                    AACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGaa        >  1:817256‑M1/1‑76 (MQ=255)
                                                                     acgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcgg       <  2:450443‑M2/76‑70 (MQ=255)
                                                                      cgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggt      <  1:311989‑M2/76‑71 (MQ=255)
                                                                       gtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggta     >  2:827985‑M2/1‑5 (MQ=255)
                                                                       gTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGtt     <  2:634280/76‑1 (MQ=255)
                                                                       gTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGtt     <  1:6844/76‑1 (MQ=255)
                                                                        tatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtat    >  2:1240329‑M2/1‑4 (MQ=255)
                                                                        tATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTa    <  1:950212/76‑1 (MQ=255)
                                                                        tATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTa    <  1:1489030/76‑1 (MQ=255)
                                                                          tgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatgg  <  2:1666689‑M2/76‑75 (MQ=255)
                                                                           |                                                                          
AAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTACC  >  REL606/1733220‑1733369

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