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 annotation gene product
* ? REL606 588495 =NA (NA)57 (1.170) 35/62 0.1 noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 37 (0.680)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)37 (0.690) 29/68 0.5 noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 65 (1.190)coding (685/1413 nt) pykF pyruvate kinase

GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666097‑666165
                                 |                                   
gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtt                                   >  1:2834650/1‑36 (MQ=35)
 tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  >  1:4510187/1‑36 (MQ=38)
  gtTTATCTGTCCAACTTTTTGGGGTCAGTACATTTc                                 <  1:6116441/36‑1 (MQ=255)
  gtTTAACTGTCCAACTTTTTGGGGTCAGTATAtaag                                 >  1:4276424/1‑33 (MQ=1)
  gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTc                                 <  1:5495358/36‑1 (MQ=40)
    ttAACTGTCCAACTTTTTGGGGTCAGTACATTTCAg                               <  1:4400187/36‑1 (MQ=255)
     tAACTGTCCAACTTTTTGGGGTCAGTACATTTCAgg                              <  1:734120/36‑1 (MQ=255)
      aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                             <  1:4975521/36‑1 (MQ=255)
       aCTGTCCAACTTTTTGGGGTCAGTACATTTCAgggt                            >  1:3117196/1‑35 (MQ=255)
       aCTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                            >  1:4791000/1‑36 (MQ=255)
        cTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                           <  1:2362447/36‑1 (MQ=255)
        cTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                           <  1:268639/36‑1 (MQ=255)
        cTGTCCAACTTTTTGGAGTCAGTACATTTCAGGGCt                           <  1:2356072/36‑1 (MQ=255)
         tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                          <  1:7393464/36‑1 (MQ=255)
          gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                         <  1:2609805/36‑1 (MQ=255)
           tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        <  1:25394/36‑1 (MQ=255)
           tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        >  1:5723053/1‑36 (MQ=255)
            ccAGCTTTTTGGGGTCAGTACATTTCAGGGCttttt                       <  1:4238117/36‑1 (MQ=255)
             cAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                      <  1:47539/36‑1 (MQ=255)
              aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCa                     <  1:3534359/36‑1 (MQ=255)
               aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTctt                    >  1:1011769/1‑34 (MQ=255)
               aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTctc                    <  1:5638076/36‑3 (MQ=255)
               aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:434520/36‑1 (MQ=255)
               aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:3730938/36‑1 (MQ=255)
               aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:165044/36‑1 (MQ=255)
               aCTGTTTGGTGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:4211610/36‑1 (MQ=25)
                cTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACg                   <  1:2935771/36‑1 (MQ=255)
                cTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACg                   >  1:6355774/1‑36 (MQ=255)
                 tttttGGGGTCAGTACATTTCAGGGCTTTTTCACGa                  >  1:6103096/1‑36 (MQ=255)
                  ttttTGGGTCAGTACATTTCAGGGCTTTTTCACGAg                 >  1:3331400/1‑36 (MQ=255)
                  ttttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                 <  1:1048408/36‑1 (MQ=255)
                  ttttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                 >  1:5797432/1‑36 (MQ=255)
                    ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc               <  1:596373/36‑1 (MQ=255)
                    ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc               >  1:3716214/1‑36 (MQ=255)
                     tGGGGTCAGTACATTTCACGCCTTTTTCACGAGCCg              >  1:6568890/1‑36 (MQ=25)
                     gggggTCAGTACATTTCGGGGTTTTTTCACGAGCCg              <  1:3728967/35‑1 (MQ=25)
                       gggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgc            <  1:6000522/36‑1 (MQ=255)
                        ggTCAGTACATTTCAGGGCTTTTTCACGAGCTgggg           >  1:2871147/1‑34 (MQ=25)
                         gccAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc          <  1:6613486/34‑1 (MQ=255)
                            aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt       >  1:2068973/1‑36 (MQ=255)
                             gTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTa      <  1:6362899/36‑1 (MQ=255)
                              tACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTAt     <  1:1063360/36‑1 (MQ=255)
                              tACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTAt     <  1:6416888/36‑1 (MQ=255)
                               aCATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:2377041/36‑1 (MQ=255)
                                 aTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc  <  1:3104670/36‑1 (MQ=255)
                                 |                                   
GTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666097‑666165

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG