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

CATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAT  >  REL606/2775845‑2775912
                                |                                   
cATTTCATCGTCCAACAAAATGGTAGCAGTACAGGt                                  >  1:3045433/1‑36 (MQ=11)
      aTCGTGCGACAAAATGGGTGCAGTACAGGTTGATAg                            <  2:833501/36‑1 (MQ=18)
      aTCGTCCAACAAAATGGGGGTAGTACAGGTCGAttg                            >  2:485289/1‑34 (MQ=14)
       tCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGa                           <  1:3100740/36‑1 (MQ=255)
        cGTCCATCATGATGGGTGCCGTGCAGGTTGATAgta                          <  1:1463930/36‑3 (MQ=9)
        cGTCCAACAAAATGGGTGCAGTACAGGTTGATAGaa                          >  2:1264660/1‑36 (MQ=255)
         gTCCAACAAAATGGGTGCAGTACAGGTTGATAGaaa                         <  1:156173/36‑1 (MQ=255)
         gTCCAACAAAATGGGTGCAGTACAGGTTGATAGaaa                         <  2:2079666/36‑1 (MQ=255)
          tCCTACAGAATGGGTGCAGTACAGGTTGATAGAAat                        <  2:3580019/36‑1 (MQ=25)
          tCCAACAAAATGGGTGCAGTACAGGTTGATAGAAat                        <  1:3333768/36‑1 (MQ=255)
            caacaaAATGGGTGCAGTACAGGTTGATAGAAatat                      <  2:1707059/36‑1 (MQ=255)
            caacaaAATGGGTGCAGTACAGGTTGATAGAAatat                      <  2:2629300/36‑1 (MQ=255)
             aacaaAATGGGTGCAGTACAGGTTGATAGAAATATc                     >  1:1843326/1‑36 (MQ=255)
             aacaaAATGGGTGCAGTACAGGTTGATAGAAATATc                     >  1:894172/1‑36 (MQ=255)
             aacaaAATGGGTGCAGTACAGGTTGATAGAAACATc                     >  2:578100/1‑36 (MQ=255)
               caaAATGGGTGCAGTACAGGTTGATAGAAATATCGc                   <  2:1648107/36‑1 (MQ=255)
                aaaaTGGGTGCAGTACAGGTTGATAGAAATATCGcc                  >  1:1833756/1‑36 (MQ=255)
                aaaaTGGGTGCAGTACAGGTTGATAGAAATATCAcc                  >  1:101053/1‑36 (MQ=255)
                  aaTGGGTGCAGTACAGGTTGATAGAAATATCGCCTg                <  1:3359118/36‑1 (MQ=255)
                  aaTGGGTGCAGTACAGGTTGATAGAAATATCGCCTg                >  2:2549125/1‑36 (MQ=255)
                    tGGGTGCAGTACAGGTTGATAGAAATATCGCCTggg              >  1:1873813/1‑36 (MQ=255)
                     gggTGCAGTACAGGTTGATAGAAATATCGCCTgggg             <  2:3509995/36‑1 (MQ=255)
                        tGCAGTACAGGTTGATAGAAATATCGCCTGGGGCtt          <  2:3177660/36‑1 (MQ=255)
                         gCAGTGCTGTTTGATTGAAATATCGCCTGGGGCttt         <  2:3049427/36‑1 (MQ=19)
                         gCAGTACAGGTTGATAGAAATATCGCCTGGGGCttt         <  1:848925/36‑1 (MQ=255)
                         gCAGTACAGGTTGATAGAAATATCGCCTGGGGCttt         <  1:2221371/36‑1 (MQ=255)
                           aGTACAGGTTGATAGAAATATCGCCTGGGGGTTTTg       >  2:258623/1‑36 (MQ=255)
                           aGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTg       >  1:1149831/1‑36 (MQ=255)
                            gTACATGTTGATAGAAATATCGCCTGGGGCTTTTGt      <  1:374849/36‑1 (MQ=255)
                            gTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGt      <  1:3434025/36‑1 (MQ=255)
                               cAGGTTGATAGAAATATCGCCTGGGGCTTTTTTCCa   >  1:2855589/1‑36 (MQ=255)
                               cAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa   >  1:1452388/1‑36 (MQ=255)
                                aGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt  >  2:1205291/1‑36 (MQ=255)
                                |                                   
CATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAT  >  REL606/2775845‑2775912

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG