Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp coding (683‑685/1413 nt) pykF → pyruvate kinase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1733647 =0 (0.040)66 (1.270) 36/66 0.1 100% coding (683/1413 nt) pykF pyruvate kinase
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region
* ? REL606 = 17336490 (0.050)47 (0.850) 29/70 0.5 100% coding (685/1413 nt) pykF pyruvate kinase
?REL606 2774435 = NA (NA)noncoding (1443/1443 nt) IS150 repeat region

ATCTCCAAAATCGAAAACCAGGAAGGCCTCAACA...................................  >  REL606/1733616‑1733649
..................................TGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  REL606/2774435‑2774469
                                                                     
ATCTCCAAAATCGAAAACCAGGAAGGCCTCAACATG                                   >  2:3050603/1‑36 (MQ=255)
  CTCCAAAATCGAAAACCAGGAAGGCCTCAACATGTA                                 >  1:2616775/1‑36 (MQ=255)
   TCCAAAATCGAAAACCAGGAAGGCCTCAACATGTAC                                >  1:802586/1‑36 (MQ=255)
    CCAAGATCGAAAACCAGGAAGGCCTCAACATGTACT                               <  1:1588985/36‑1 (MQ=255)
    CCAAAATCGAAAACCAGGAAGGCCTCAACATGTACT                               >  2:988009/1‑36 (MQ=255)
    CCAAAATCGAAAACCAGGAAGGCCTCAACATGTACT                               <  2:2234905/36‑1 (MQ=255)
      AAAATCGAAAACCAGGAAGGCCTCAACATGTACTGA                             <  2:2740799/36‑1 (MQ=255)
        AATCGAAAACCAGGAAGGCCTCAACATGTACTGACC                           >  1:1371748/1‑36 (MQ=255)
        AATCGAAAACCAGGAAGGCCTCAACATGTACTGACC                           >  2:429035/1‑36 (MQ=255)
         ATCGAAAACCAGGAAGGCCTCAACATGTACTGACCC                          >  2:2006914/1‑36 (MQ=255)
         ATCGAAAACCAGGAAGGCCTCAACATGTACTGACCC                          <  1:691462/36‑1 (MQ=255)
          TCGAAAACCAGGAAGGCCTCAACATGTACTGACCCC                         <  1:2723664/36‑1 (MQ=255)
            GAAAACCAGGAAGGCCTCAACATGTACTGACCCCAA                       >  2:3340523/1‑36 (MQ=255)
             AAAACCAGGAAGGCCTCAACATGTACTGACCCCAAA                      >  2:1760710/1‑36 (MQ=255)
             AAAACCAGGAAGGCCTCAACATGTACTGACCCCAAA                      <  1:3726680/36‑1 (MQ=255)
             AAAACCAGGAAGGCCTCAACATGTACTGACCCCAAA                      <  2:323199/36‑1 (MQ=255)
               AACCAGGAAGGCCTCAACATGTACTGACCCCAAAAA                    <  1:2471673/36‑1 (MQ=255)
                ACCAGGAAGGCCTCAACATGTACTGACCCCAAAAAG                   >  2:1039357/1‑36 (MQ=255)
                  CAGGAAGGCCTCAACATGTACTGACCCCAAAAAGTT                 <  1:3356538/36‑1 (MQ=255)
                   AGGAGGGCCTCAACATGTACTGACCCCAAAAAGTTG                <  1:589780/36‑1 (MQ=255)
                   AGGAAGGCCTCAACATGTACTGACCCCAAAAAGTTG                >  1:3398128/1‑36 (MQ=255)
                   AGGAAGGCCTCAACATGTACTGACCCCAAAAAGTTG                <  2:3221748/36‑1 (MQ=255)
                    GGAAGGCCTCAACATGTACTGACCCCAAAAAGTTGG               <  1:651316/36‑1 (MQ=255)
                     GAAGGCCTCAACATGTACTGACCCCAAAAAGTTGGA              >  1:3644252/1‑36 (MQ=255)
                          CCTCAACATGTACTGACCCCAAAAAGTTGGACAGTT         >  2:2046875/1‑36 (MQ=255)
                          CCTCAACATGTACTGACCCCAAAAAGTTGGACAGTT         <  2:1176510/36‑1 (MQ=255)
                           CTCAACATGTACTGACCCCAAAAAGTTGGACAGTTA        >  1:2387214/1‑36 (MQ=255)
                            TCAACATGTACTGACCCCAAAAAGTTGGACAGTTAA       >  2:2521362/1‑36 (MQ=40)
                            TCAACATGTACTGACCCCAAAAAGTTGGACAGTTAA       <  2:2754198/36‑1 (MQ=40)
                             CAACCTGTACTGACCCCAAAAAGTTGGACAGTTAAA      >  1:40355/1‑36 (MQ=32)
                             CAACATGTACTGACCGCAAAAAGTTGGACAGTTAAA      <  2:2258718/36‑1 (MQ=255)
                             CAACATCTACTGACCCCAAAAAGTTGGACAGTTAAA      <  2:2762656/36‑1 (MQ=255)
                              cgtATGTACTGACCCCAAAAAGTTGGACAGTTAAAC     >  1:2886302/4‑36 (MQ=1)
                              AACCTGTACTGACCCCAAAAAGTTGGACAGTTAAAC     <  2:682036/36‑1 (MQ=1)
                              AACCTGTACTGACCCCAAAAAGTTGGACAGTTAAAC     >  2:2986398/1‑36 (MQ=1)
                              AACATGTACTGACCCCAAAAAGTTGGACAGTTAAAC     >  1:2608350/1‑36 (MQ=38)
                               gtATGTACTGACCCCAAAAAGTTGGACAGTTAAACA    <  1:635345/34‑1 (MQ=1)
                               gtATGTACTGACCCCAAAAAGTTGGACAGTTAAACA    <  2:2592512/34‑1 (MQ=1)
                               gtATGTACTGACCCCAAAAAGTTGGACAGTTAAACA    >  2:591026/3‑36 (MQ=1)
                               ACATGTACTGACCCCAAAAAGTTGGACAGTTAAACA    >  1:1265789/1‑36 (MQ=35)
                                tATGTACTGACCCCAAAAAGTTGGACAGTTAAACAC   >  1:2026751/2‑36 (MQ=1)
                                tATGTACTGACCCCAAAAAGTTGGACAGTTAAACAC   >  1:3440988/2‑36 (MQ=1)
                                CATGTACTGACCCCAAAAAGTTGGACAGTTAAACAC   >  2:3314998/1‑36 (MQ=32)
                                 ATGTGCTGACCCCAAAAAGTTGGACAGTTAAACACG  <  2:1264826/36‑1 (MQ=0)
                                 ATGTACTGACCCCAAAAAGTTGGACAGTTATACACG  >  2:2847805/1‑36 (MQ=0)
                                 ATGTACTGACCCCAAAAAGTTGGACAGTTAAACAag  >  1:2474577/1‑34 (MQ=1)
                                 ATGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  2:632853/1‑36 (MQ=1)
                                                                     
ATCTCCAAAATCGAAAACCAGGAAGGCCTCAACA...................................  >  REL606/1733616‑1733649
..................................TGTACTGACCCCAAAAAGTTGGACAGTTAAACACG  >  REL606/2774435‑2774469

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG