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 =7 (0.320)24 (1.110) 21/146 0.2 84.5% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 2 (0.090)coding (331/1413 nt) pykF pyruvate kinase

AAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAC  >  REL606/1733220‑1733368
                                                                           |                                                                         
aaaGCTGGTCAGACCTTTACTTTCACCACTGATAACTCTGTTATCGGCAACAGCGAAATGGTCGCGGTAACGTATg                                                                           >  2:641007/1‑76 (MQ=255)
  agacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatgaa                                                                         >  2:435041‑M1/67‑76 (MQ=255)
   gacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAg                                                                        >  1:431408‑M1/66‑76 (MQ=255)
    cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                       >  1:205727‑M2/1‑72 (MQ=255)
    cTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacg                                                                       <  2:544453‑M2/76‑5 (MQ=255)
    acctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAgg                                                                       <  1:623485‑M1/12‑1 (MQ=255)
         tactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCa                                                                  >  2:657525‑M1/60‑76 (MQ=255)
          aGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaa                                                                 >  2:435041‑M2/1‑66 (MQ=255)
           gACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaag                                                                >  1:431408‑M2/1‑65 (MQ=255)
           ctttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCact                                                                >  2:97170‑M1/58‑76 (MQ=255)
            aCCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaagg                                                               <  1:623485‑M2/76‑13 (MQ=255)
              tcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCactact                                                             >  2:695864‑M1/55‑76 (MQ=255)
                 tACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttca                                                          >  2:657525‑M2/1‑59 (MQ=255)
                 ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAc                                                          >  2:337924‑M1/52‑76 (MQ=255)
                   cTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcact                                                        >  2:97170‑M2/1‑57 (MQ=255)
                      tCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactact                                                     >  2:695864‑M2/1‑54 (MQ=255)
                         ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgac                                                  >  2:337924‑M2/1‑51 (MQ=255)
                         aaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctgt                                                  <  1:233552‑M1/33‑1 (MQ=255)
                         aaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACctgtctgt                                                  <  1:603714‑M1/33‑1 (MQ=255)
                           atctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTg                                                <  1:471717‑M1/35‑1 (MQ=255)
                                 aaaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgt                                          <  1:233552‑M2/76‑34 (MQ=255)
                                 aaaTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgt                                          <  1:603714‑M2/76‑34 (MQ=255)
                                   aTCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttg                                        <  1:471717‑M2/76‑36 (MQ=255)
                                    cggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGGTTCACTACTGACCTGTCTGTTGGCAACACCg                                       <  2:313228‑M1/44‑1 (MQ=255)
                                      gcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTa                                     >  1:522822‑M1/31‑76 (MQ=255)
                                           agcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGt                                <  2:6429‑M1/51‑1 (MQ=255)
                                            gcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGtt                               >  1:47893‑M1/25‑76 (MQ=255)
                                            cGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaagggttcactactgacctgtctgttggcaacaccg                               <  2:313228‑M2/76‑45 (MQ=255)
                                             cgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTg                              >  1:60595‑M1/24‑76 (MQ=255)
                                              gaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGa                             <  2:267876‑M1/54‑1 (MQ=255)
                                              gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgta                             >  1:522822‑M2/1‑30 (MQ=255)
                                                   ggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATg                        <  1:715080‑M1/59‑1 (MQ=255)
                                                   gggtgcggtaacgtatgAACGTATGCAGGTTTCACTACTGACCTGTCTGTTGGCACCACCGTACTGGTTGACGATg                        <  2:234991‑M1/59‑1 (MQ=255)
                                                   aGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggt                        <  2:6429‑M2/76‑52 (MQ=255)
                                                    gCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggtt                       >  1:47893‑M2/1‑24 (MQ=255)
                                                     cGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttg                      >  1:60595‑M2/1‑23 (MQ=255)
                                                      gAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttga                     <  2:267876‑M2/76‑55 (MQ=255)
                                                       gcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCt                    >  2:541698‑M1/14‑76 (MQ=255)
                                                       gcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCt                    >  2:631340‑M1/14‑76 (MQ=255)
                                                       gcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCt                    <  1:57292‑M1/63‑1 (MQ=255)
                                                           gggTGCGGTAACGTATGaacgtatgcaggtttcactactgacctgtctgttggcaccaccgtactggttgacgatg                <  2:234991‑M2/76‑60 (MQ=255)
                                                           ggTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatg                <  1:715080‑M2/76‑60 (MQ=255)
                                                               gCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtct            >  2:541698‑M2/1‑13 (MQ=255)
                                                               gCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtct            >  2:631340‑M2/1‑13 (MQ=255)
                                                               gCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtct            <  1:57292‑M2/76‑64 (MQ=255)
                                                                 atgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATg          <  2:171136‑M1/73‑1 (MQ=255)
                                                                        tATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTa   >  2:217302/1‑76 (MQ=255)
                                                                         atgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcggtatg  <  2:171136‑M2/76‑74 (MQ=255)
                                                                         aTGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAc  >  1:315809/1‑76 (MQ=255)
                                                                           |                                                                         
AAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAAGTTAC  >  REL606/1733220‑1733368

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