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 =1 (0.040)22 (0.950) 18/146 0.4 93.7% coding (324/1413 nt) pykF pyruvate kinase
?REL606 = 1733295 2 (0.090)coding (331/1413 nt) pykF pyruvate kinase

TTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733213‑1733363
                                                                           |                                                                           
aagctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatga                                                                             <  2:725548‑M1/1‑1 (MQ=255)
 agctggtcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaa                                                                            >  1:397491‑M1/75‑76 (MQ=255)
  ctctGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAAc                                                                           <  2:193164/76‑1 (MQ=255)
       tcagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatg                                                                      <  2:462127‑M1/8‑1 (MQ=255)
        cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                     >  2:343484‑M1/68‑76 (MQ=255)
        cagacctttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgaacgtatga                                                                     >  2:130396‑M1/68‑76 (MQ=255)
        aaGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGa                                                                     <  2:725548‑M2/76‑2 (MQ=255)
         aGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaa                                                                    >  1:397491‑M2/1‑74 (MQ=255)
              tttactttcaccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGttt                                                               >  2:597096‑M1/62‑76 (MQ=255)
               tCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatg                                                              <  2:462127‑M2/76‑9 (MQ=255)
                cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                             >  2:343484‑M2/1‑67 (MQ=255)
                cAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatga                                                             >  2:130396‑M2/1‑67 (MQ=255)
                 actttcacgactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCac                                                            <  2:484709‑M1/18‑1 (MQ=255)
                      tttACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggttt                                                       >  2:597096‑M2/1‑61 (MQ=255)
                        ccactgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGAc                                                     <  2:698252‑M1/25‑1 (MQ=255)
                         aCTTTCACGACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcac                                                    <  2:484709‑M2/76‑19 (MQ=255)
                          actgataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACct                                                   >  2:639748‑M1/50‑76 (MQ=255)
                              ataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGGTTCACTACTGACctgtct                                               <  2:583704‑M1/31‑1 (MQ=255)
                              ataaatctgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGGTTCAATACTGACctgtct                                               <  2:583700‑M1/31‑1 (MQ=255)
                                ccacTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgac                                             <  2:698252‑M2/76‑26 (MQ=255)
                                  acTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacct                                           >  2:639748‑M2/1‑49 (MQ=255)
                                    ntgttatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGc                                         <  1:45859‑M1/37‑1 (MQ=255)
                                      aTAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaagggttcactactgacctgtct                                       <  2:583704‑M2/76‑32 (MQ=255)
                                      aTAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaagggttcaatactgacctgtct                                       <  2:583700‑M2/76‑32 (MQ=255)
                                        tatcggcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAaca                                     >  1:107315‑M1/36‑76 (MQ=255)
                                            ntgTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggc                                 <  1:45859‑M2/75‑38 (MQ=255)
                                             gcaacagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTa                                >  2:750491‑M1/31‑76 (MQ=255)
                                                tATCGGCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaaca                             >  1:107315‑M2/1‑35 (MQ=255)
                                                acagcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTg                             >  1:365564‑M1/28‑76 (MQ=255)
                                                   gcgaaatggttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGtt                          >  1:415701‑M1/25‑76 (MQ=255)
                                                     gCAACAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgta                        >  2:750491‑M2/1‑30 (MQ=255)
                                                        aCAGCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactg                     >  1:365564‑M2/1‑27 (MQ=255)
                                                           gttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATgg                  >  2:539238‑M1/17‑76 (MQ=255)
                                                           gCGAAATGGTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggtt                  >  1:415701‑M2/1‑24 (MQ=255)
                                                            ttgcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGt                 <  1:219576‑M1/61‑1 (MQ=255)
                                                             ttcgataacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTc                <  2:664438‑M1/62‑1 (MQ=255)
                                                              gcggtaacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCt               <  2:199253‑M1/63‑1 (MQ=255)
                                                                   gTTGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatgg          >  2:539238‑M2/1‑16 (MQ=255)
                                                                   aacgtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCg          <  2:277210‑M1/68‑1 (MQ=255)
                                                                    ttGCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggt         <  1:219576‑M2/76‑62 (MQ=255)
                                                                     ttcgaTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtc        <  2:664438‑M2/74‑63 (MQ=255)
                                                                      gtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTa       <  1:534082‑M1/71‑1 (MQ=255)
                                                                      gtatgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTa       >  1:512691‑M1/6‑76 (MQ=255)
                                                                      gCGGTAACGTATGaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtct       <  2:199253‑M2/76‑64 (MQ=255)
                                                                         tgAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATgg    <  1:602629‑M1/74‑1 (MQ=255)
                                                                           aacgtatgaacgtatgaaggtttcactactgacctgtctgttggcaacaccgtactggttgacgatggtctgatcg  <  2:277210‑M2/76‑69 (MQ=255)
                                                                           |                                                                           
TTCTCTGAAAGCTGGTCAGACCTTTACTTTCACCACTGATAAATCTGTTATCGGCAACAGCGAAATGGTTGCGGTAACGTATGAAGGTTTCACTACTGACCTGTCTGTTGGCAACACCGTACTGGTTGACGATGGTCTGATCGGTATGGAA  >  REL606/1733213‑1733363

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