Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,733,647 IS150 (–) +3 bp 98.2% 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)52 (0.970) 35/62 NT 98.1% noncoding (1/1977 nt) IS150 repeat region
?REL606 1733647 = 1 (0.070)coding (683/1413 nt) pykF pyruvate kinase
* ? REL606 = 666130NA (NA)55 (0.930) 40/68 NT 98.2% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 1733649 1 (0.020)coding (685/1413 nt) pykF pyruvate kinase

Only 100 of 109 total aligned reads displayed.

CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666096‑666165
                                  |                                   
cGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAt                                    >  1:1263016‑M1/1‑35 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGTGGTCAGTACAtt                                   <  1:3182884/36‑1 (MQ=255)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtt                                   <  1:5147933/36‑1 (MQ=35)
 gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtt                                   >  1:5419213/1‑36 (MQ=35)
  tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAttt                                  <  1:6139893/36‑1 (MQ=38)
  tgtTTAACTGTCCAACTTTTTGGGGTCAGTACAtgt                                  <  1:5977080‑M1/36‑4 (MQ=255)
  tgtTTAACTGTCCAACTTTTCGGGGTCAGTACAttt                                  <  1:3418354/36‑1 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAtttg                                 >  1:1544219/1‑35 (MQ=38)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAtgtt                                 >  1:1457406‑M1/1‑32 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAggtt                                 >  1:1375041‑M1/1‑32 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACAggtt                                 <  1:589376‑M1/36‑5 (MQ=255)
   gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTc                                 >  1:2870494/1‑36 (MQ=40)
    tttAACTGTCCAACTTTTTGGGGTCAGTACAtgttg                                <  1:2354904‑M1/36‑6 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACAtgttg                                <  1:7625522‑M1/36‑6 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACAtgttg                                >  1:6911413‑M1/1‑31 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                                <  1:3536179/36‑1 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                                <  1:5770817/36‑1 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                                <  1:6019451/36‑1 (MQ=255)
    tttAACTGTCCAACTTTTTGGGGTCAGTACATTTCa                                <  1:1319348/36‑1 (MQ=255)
     ttAACTGTCCACCTTTTTGGGGTCAGTACAtgttga                               <  1:1244634‑M1/36‑7 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACAtgttgag                              >  1:6590643‑M1/1‑29 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACAtgttgag                              >  1:4623959‑M1/1‑29 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACAtgttgag                              <  1:4061294‑M1/36‑8 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACAtgttgag                              >  1:5940827‑M1/1‑29 (MQ=255)
      tAACTGTCCAACTTTTTGGGGTCAGTACATTTCAgg                              <  1:6436053/36‑1 (MQ=255)
       aaCTTTCCAACTTTTTGGGGTCAGTACAtgttgagg                             <  1:4471424‑M1/36‑9 (MQ=255)
       aaCTGTCCAACTTTTTGGGGTCAGTACAtgttgagg                             >  1:1352476‑M1/1‑28 (MQ=255)
       aaCTGTCCAACTTTTTGGGGTCAGTACATTTCAggg                             <  1:944358/36‑1 (MQ=255)
        aCTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                            <  1:7206924/36‑1 (MQ=255)
        aCTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGc                            <  1:2578358/36‑1 (MQ=255)
         cTGTCCAACTTTTTGGGGTCAGTACAtgttgaggcc                           <  1:6963972‑M1/36‑11 (MQ=255)
         cTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                           >  1:493323/1‑36 (MQ=255)
         cTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCt                           >  1:3316015/1‑36 (MQ=255)
         cTGTCCAACTTTTTGGGGTCAGTACAAAACAGGGcg                           >  1:7305737/1‑35 (MQ=11)
          tGTCCAACTTTTTGGGGTCAGTACAtgttgaggcct                          >  1:3144897‑M1/1‑25 (MQ=255)
          tGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCtt                          >  1:7371294/1‑36 (MQ=255)
          tGTCCAACTTTTTGCGGTCAGTACAtgttgaggcct                          >  1:6345976‑M1/1‑25 (MQ=255)
          tGTCCAACTTTTGGGGGTCAGTACATTTCAGGGCtt                          <  1:6279/36‑1 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         >  1:4922092‑M1/1‑24 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACAtgttgaggcctt                         >  1:1370259‑M1/1‑24 (MQ=255)
           gTCCAACTTTTTGGGGTCAGTACATTTCAGGGCttt                         >  1:3928618/1‑36 (MQ=255)
             ccAACTTTTTGGGGTCAGTACAtgttgaggccttcc                       <  1:2873443‑M1/36‑15 (MQ=255)
             ccAACTTTTTGGGGTCAGTACAtgttgaggccttcc                       <  1:6720214‑M1/36‑15 (MQ=255)
             ccAACTTTTTGGGGTCAGTACATTTCAGGGCttttt                       >  1:4968353/1‑36 (MQ=255)
             ccAACTTTTTGGGGTCAGTACATTTCAGGGCttttt                       <  1:2054839/36‑1 (MQ=255)
               aaCTTTTTGGGGTCAGTACAtgttgaggccttcctg                     <  1:6183802‑M1/36‑17 (MQ=255)
               aaCTTTTTGGGGTCAGTACAtgttgaggccttcctg                     <  1:4310337‑M1/36‑17 (MQ=255)
               aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCa                     >  1:6621435/1‑36 (MQ=255)
               aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCa                     <  1:482264/36‑1 (MQ=255)
               aaCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCa                     <  1:2100496/36‑1 (MQ=255)
                aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:2244270/36‑1 (MQ=255)
                aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:5823215/36‑1 (MQ=255)
                aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:599043/36‑1 (MQ=255)
                aCTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCAc                    <  1:6169371/36‑1 (MQ=255)
                aCTTTTTGGGGTCAGTACATTTCAGGGCTGTTTccc                    <  1:3410190/36‑3 (MQ=25)
                 cTTTTTGGGGTCAGTACAtgttgaggccttcctggt                   <  1:369749‑M1/36‑19 (MQ=255)
                 cTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACg                   >  1:5130405/1‑36 (MQ=255)
                  tttttGGGGTCAGTACAtgttgaggccttcctggtt                  >  1:5942946‑M1/1‑17 (MQ=255)
                  tttttGGGGTCAGTACAtgttgaggccttcctggtt                  <  1:4653535‑M1/36‑20 (MQ=255)
                   ttttGGGGTCAGTACAtgttgaggccttcctggttt                 <  1:604363‑M1/36‑21 (MQ=255)
                   ttttGGGGTCAGTACATTTCAGGGCTTTTTCACGAg                 >  1:940716/1‑36 (MQ=255)
                    tttGGGGTCAGTACAtgttgaggccttcctggtttt                >  1:6904157‑M1/1‑15 (MQ=255)
                    tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                <  1:5351359/36‑1 (MQ=255)
                    tttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGc                <  1:2724144/36‑1 (MQ=255)
                     ttGGGGTCAGTACAtgttgaggccttcctggttttc               <  1:2555365‑M1/36‑23 (MQ=255)
                     ttGGGGTCAGTACATTTCAGGGCTTTTTCACGAGcc               >  1:3302300/1‑36 (MQ=255)
                     ttGGGGTCAATACAtgttgaggccttcctggttttc               >  1:4305331‑M1/1‑14 (MQ=255)
                      tGGGGTCAGTACAtgttgaggccttcctggttttcg              <  1:4811789‑M1/36‑24 (MQ=255)
                      tGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCg              >  1:420660/1‑36 (MQ=255)
                       ggggTCAGTACAtgttgaggccttcctggttttcga             >  1:1843478‑M1/1‑12 (MQ=255)
                       ggggTCAGTACAtgttgaggccttcctggttttcga             >  1:1213674‑M1/1‑12 (MQ=255)
                       ggggTCAGTACATTTCAGGGCTTTTTCACGAGCCgg             >  1:3599682/1‑36 (MQ=255)
                        gggTCAGTACAtgttgaggccttcctggttttcgat            <  1:6054337‑M1/36‑26 (MQ=255)
                        gggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgc            >  1:3435129/1‑36 (MQ=255)
                        gggTCAGTACATTTCAGGGCTTTTTCACGAGCCGgc            >  1:3358650/1‑36 (MQ=255)
                         ggTCAGTACAtgttgaggccttcctggttttcgatt           <  1:1055340‑M1/36‑27 (MQ=255)
                          gTCAGTACAtgttgaggccttcctggttttcgattt          >  1:5670162‑M1/1‑9 (MQ=255)
                          gTCAGTACAtgttgaggccttcctggttttcgattt          <  1:4683240‑M1/36‑28 (MQ=255)
                          gTCAGTACATTTCAGGGCTTTTTCACGAGCCGgcgc          <  1:4397692/36‑1 (MQ=255)
                          gTCAGTACATTTCAGGGCTTTTTCACGAGACGgcgc          >  1:684059/1‑36 (MQ=255)
                           tCAGTACAtgttgaggccttcctggttttcgatttt         <  1:7858885‑M1/36‑29 (MQ=255)
                           tCAGTACAtgttgaggccttcctggttttcgatttt         <  1:3909576‑M1/36‑29 (MQ=255)
                           tCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCt         <  1:1670776/36‑1 (MQ=255)
                            cAGTACAtgttgaggccttcctggttttcgattttg        >  1:6852691‑M1/1‑7 (MQ=255)
                             aGTACAtgttgaggccttcctggttttcgattttgg       <  1:7771319‑M1/36‑31 (MQ=255)
                             aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt       >  1:3821315/1‑36 (MQ=255)
                             aGTACATTTCAGGGCTTTTTCACGAGCCGGCGCttt       >  1:5988230/1‑36 (MQ=255)
                              gTACAtgttgaggccttcctggttttcgattttgga      <  1:6826474‑M1/36‑32 (MQ=255)
                              gTACAtgttgaggccttcctggttttcgattttgga      >  1:3944634‑M1/1‑5 (MQ=255)
                              gTACAtgttgaggccttcctggttttcgattttgga      >  1:6411285‑M1/1‑5 (MQ=255)
                              gTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTa      >  1:1770850/1‑36 (MQ=255)
                               tACAtgttgaggccttcctggttttcgattttggag     >  1:7837541‑M1/1‑4 (MQ=255)
                               tACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTAt     <  1:4655474/36‑1 (MQ=255)
                                aCAtgttgaggccttcctggttttcgattttggaga    >  1:3188733‑M1/1‑3 (MQ=255)
                                aCATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc    <  1:143382/36‑1 (MQ=255)
                                aCATTTCAGGGCTGTTTCACGAGCCGGCGCTTTATc    >  1:6563456/1‑36 (MQ=255)
                                 cAtgttgaggccttcctggttttcgattttggagat   >  1:4309542‑M1/1‑2 (MQ=255)
                                 cAtgttgaggccttcctggttttcgattttggagat   <  1:7331300‑M1/36‑35 (MQ=255)
                                  atgttgaggccttcctggttttcgattttggagatg  <  1:4221973‑M1/36‑36 (MQ=255)
                                  atgttgaggccttcctggttttcgattttggagatg  >  1:5005022‑M1/1‑1 (MQ=255)
                                  |                                   
CGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGC  >  REL606/666096‑666165

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG