Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 100% 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 588495 =NA (NA)40 (1.010) 26/62 NT 100% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 0 (0.040)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)56 (1.280) 35/68 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 0 (0.020)coding (685/1413 nt) pykF pyruvate kinase

Only 100 of 213 total aligned reads displayed.

TCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666095‑666165
                                   |                                   
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     <  1:1060891‑M2/36‑1 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:1098186‑M2/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:5719968‑M2/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     >  1:1178651‑M2/1‑36 (MQ=255)
tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA                                     <  1:5302311‑M2/36‑1 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    <  1:2958898‑M1/36‑2 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    <  1:2978426‑M1/36‑2 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    >  1:3066004‑M1/1‑35 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    >  1:5308791‑M1/1‑35 (MQ=255)
 cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    <  1:5589793‑M1/36‑2 (MQ=255)
  gtgtTTAACTGTCCAACTTTTTTGGGTCAGTACAta                                   >  1:5093032‑M1/1‑34 (MQ=255)
  gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAta                                   <  1:4248795‑M1/36‑3 (MQ=255)
   tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtac                                  >  1:5638412‑M1/1‑33 (MQ=255)
    gtTTAACTGTCCAACTTTTTGGGGTCAGTACAggtt                                 <  1:3560027‑M1/36‑5 (MQ=255)
    gtTTAACTGTCCAACTTTTTGGGGTCAGTACAggtt                                 <  1:775052‑M1/36‑5 (MQ=255)
    gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTc                                 >  1:727802/1‑36 (MQ=40)
     tttAACTGTCCAACTTTTTTGGGTCAGTACATTTct                                >  1:3890815/1‑35 (MQ=255)
     tttAACTGTCCAACTTTTTGGTGTCAGTACAtgttg                                <  1:1540506‑M1/36‑6 (MQ=255)
     tttAACTGTCCAACTTTTTGGGGTCAGTACAtgttg                                >  1:1897360‑M1/1‑31 (MQ=255)
     tttAACTGTCCAACTTTTTGGGGTCAGTACAtacgc                                <  1:4686943‑M2/36‑6 (MQ=255)
     tttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                                <  1:37358/36‑1 (MQ=255)
      ttAACTGTGTAGCTTTTTGGGGTCAGTACATTTCAg                               <  1:3659586/36‑1 (MQ=17)
      ttAACTGTCTAACTTTTTGGGGTCAGTACATTTCAg                               <  1:3006734/36‑1 (MQ=255)
      ttAACTGTCCAACTTTTTGGGGTCAGTACAtgttga                               <  1:5311735‑M1/36‑7 (MQ=255)
      ttAACTGTCCAACTTTTTGGGGTCAGTACAtacgcg                               <  1:3464020‑M2/36‑7 (MQ=255)
      ttAACTGTCCAACTTTTTGGGGTCAGTACAggttgc                               <  1:3929976‑M1/36‑7 (MQ=255)
      ttAACTGTCCAACTTTTTGGGGTCAGTACATTTCAg                               >  1:3756559/1‑36 (MQ=255)
       tAACTGTCCAACTTTTTGGGGTCAGTACAtgttgag                              >  1:414481‑M1/1‑29 (MQ=255)
        aaCTGTCCAACTTTTTGGGGTCAGTACAtgttgagg                             >  1:3983917‑M1/1‑28 (MQ=255)
        aaCTGTCCAACTTTTTGGGGTCAGTACAtacgcgta                             <  1:504952‑M2/36‑9 (MQ=255)
        aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                             >  1:3189216/1‑36 (MQ=255)
         aCTGTCCAACTTTTTGGGGTCAGTACAtgttgaggc                            <  1:2055804‑M1/36‑10 (MQ=255)
         aCTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                            >  1:1254028/1‑36 (MQ=255)
         aCTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                            <  1:2372166/36‑1 (MQ=255)
          cTTTCCAACTTTTTGGGGTCAGTACAtacgcgtaaa                           >  1:1592056‑M2/1‑26 (MQ=255)
          cTGTCCAACTTTTTGGGGTCAGTACAtgttgaggcc                           >  1:1127081‑M1/1‑26 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAtgttgaggcct                          <  1:905139‑M1/36‑12 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAtacgcgtaaaa                          >  1:196983‑M2/1‑25 (MQ=255)
           tGTCCAACTTTTTGGGGTCAGTACAtacgcgtaaaa                          >  1:589968‑M2/1‑25 (MQ=255)
           tGTCCAACTGTTTGGGGTCAGTACAtacgcgtaaaa                          >  1:762742‑M2/1‑25 (MQ=255)
            gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         >  1:1755164‑M1/1‑24 (MQ=255)
            gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         <  1:5642823‑M1/36‑13 (MQ=255)
            gTCCAACTTTTTGGGGTCAGTACAtacgcgtaaaac                         >  1:4154537‑M2/1‑24 (MQ=255)
             tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        <  1:2556547/36‑1 (MQ=255)
             tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        <  1:4652287/36‑1 (MQ=255)
             tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        >  1:4498702/1‑36 (MQ=255)
             tCCAACTTTTTGGGGTCAGTACATTTCAGGGCtttt                        <  1:5851747/36‑1 (MQ=255)
              ccAACTTTTTGGGGTCAGTACAtgttgaggccttcc                       <  1:3436889‑M1/36‑15 (MQ=255)
              ccAACTTTTTGGGGTCAGTACAtgttgaggccttcc                       <  1:4480476‑M1/36‑15 (MQ=255)
              ccAACTTTTTGGGGTCAGTACATTTCAGGGCttttt                       <  1:4282657/36‑1 (MQ=255)
              ccAACTTTTTGGGGTCAGTACATTTCAGGGCttttt                       <  1:1268447/36‑1 (MQ=255)
               cAACTTTTTGGGGTCAGTACAtgttgaggccttcct                      <  1:5912946‑M1/36‑16 (MQ=255)
               cAACTTTTTGGGGTCAGTACAtgttgaggccttcct                      <  1:291479‑M1/36‑16 (MQ=255)
               cAACTTTTTGGGGTCAGTACAtacgcgtaaaacagc                      >  1:1571398‑M2/1‑21 (MQ=255)
               cAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTc                      <  1:2431258/36‑1 (MQ=255)
                  cTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACg                   <  1:4080319/36‑1 (MQ=255)
                     tttGGGGTCAGTACAtgttgaggccttcctggtttt                <  1:5221215‑M1/36‑22 (MQ=255)
                     tttGGGGTCAGTACAtacgcgtaaaacagccagcgc                >  1:435081‑M2/1‑15 (MQ=255)
                     tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                <  1:2255716/36‑1 (MQ=255)
                         gggTCAGTACAtgttgaggccttcctggttttcgat            <  1:1447377‑M1/36‑26 (MQ=255)
                         gggTCAGTACAtgttgaggccttcctggttttcgat            >  1:1886920‑M1/1‑11 (MQ=255)
                         gggTCAGTAAAtacgcgtaaaacagccagcgctggc            >  1:1141244‑M2/1‑11 (MQ=255)
                          ggTCAGTACAtgttgaggccttcctggttttcgatt           <  1:3381031‑M1/36‑27 (MQ=255)
                          ggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgcg           >  1:678565/1‑36 (MQ=255)
                             aagTAAAtacgcgtaaaacagccagcgctggcgcga        >  1:38402‑M2/2‑7 (MQ=255)
                              aGTACAtgttgaggccttcctggttttcgattttgg       >  1:2339346‑M1/1‑6 (MQ=255)
                              aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt       <  1:314631/36‑1 (MQ=255)
                               gTACAtacgcgtaaaacagccagcgctggcgcgatt      <  1:1769336‑M2/36‑32 (MQ=255)
                               gTACATTTCAGGGCTTTTTCACGAGCCGGCGAttat      >  1:3480276/1‑34 (MQ=38)
                                gtaatacgcgtaaaacagccagcgctggcgcgattt     >  1:5403394‑M2/4‑4 (MQ=255)
                                 taatacgcgtaaaacagccagcgctggcgcgattta    >  1:4797163‑M2/3‑3 (MQ=255)
                                 caatacgcgtaaaacagccagcgctggcgcgattta    >  1:4920368‑M2/3‑3 (MQ=255)
                                 aCATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:1950657/36‑1 (MQ=255)
                                  cAtacgcgtaaaacatccagcgctggcgcgatttag   <  1:2094446‑M2/36‑35 (MQ=255)
                                  cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:4914467‑M2/1‑2 (MQ=255)
                                  cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:4571308‑M2/1‑2 (MQ=255)
                                  cAtacgcgtaaaacagccagcgctggcgcgatttag   <  1:881767‑M2/36‑35 (MQ=255)
                                  cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:3969159‑M2/1‑2 (MQ=255)
                                  cAtacgcgtaaaacagccagcgctggcgcgatttag   <  1:5688398‑M2/36‑35 (MQ=255)
                                  cAtacgcgtaaaacagccagcgctggcgcgatttag   >  1:258382‑M2/1‑2 (MQ=255)
                                  cAtacgagtaaaacagccagcgctggcgcgatttag   >  1:651271‑M2/1‑2 (MQ=255)
                                  cATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg   <  1:944761/36‑1 (MQ=255)
                                  cATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg   <  1:3338359/36‑1 (MQ=255)
                                  aatacgcgtaaaacagccagcgctggcgcgatttag   >  1:2684438‑M2/2‑2 (MQ=255)
                                  aatacgcataaaacagccagcgctggcgcgatttag   >  1:5028387‑M2/2‑2 (MQ=255)
                                   atacgcgtaaaccagccagcgctggcgcgatttagc  >  1:464779‑M2/1‑1 (MQ=255)
                                   atacgcgtaaaacagncagcgctggcgcgatttagc  <  1:5161319‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcggtttagc  <  1:2723267‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:5273721‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  >  1:589280‑M2/1‑1 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:2517830‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:453842‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:2124706‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:2564539‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:1105263‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:852518‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:4436916‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctggcgcgatttagc  <  1:4223455‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgctagcgcgatttagc  <  1:1386363‑M2/36‑36 (MQ=255)
                                   atacgcgtaaaacagccagcgccggcgcgatttagc  <  1:1637082‑M2/36‑36 (MQ=255)
                                   |                                   
TCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666095‑666165

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG