Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,517,306 IS186 (+) +8 bp :: Δ2 56.2% coding (187‑194/1743 nt) ggt ← gamma‑glutamyltranspeptidase periplasmic precursor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 590715 =NA (NA)34 (0.380) 22/42 NT 41.2% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 3517313 30 (0.540)coding (187/1743 nt) ggt gamma‑glutamyltranspeptidase periplasmic precursor
* ? REL606 = 2774194NA (NA)52 (0.530) 28/46 NT 47.8% noncoding (1341/1343 nt) IS186 repeat region
?REL606 3517306 = 37 (0.580)coding (194/1743 nt) ggt gamma‑glutamyltranspeptidase periplasmic precursor

Only 100 of 233 total aligned reads displayed.

TCAACGGCATTCCCGCCCTCCTTGAGGATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCT  >  REL606/3517279‑3517347
                                  |                                  
tCAACGGCATTCCCGCCCTCCTTGAGGATATCCAc                                    <  1:18779442/35‑1 (MQ=255)
tCAACGGCATTCCCGCCCTCCTTGAGGATATCCAc                                    <  1:10559467/35‑1 (MQ=255)
tCAACGGCATTCCCGCCCTCCTTGAGGATATCCAc                                    <  1:1258034/35‑1 (MQ=255)
tCAACGGCATTCCCGCCCTCCTTGAGGATATCCAc                                    <  1:7666085/35‑1 (MQ=255)
tCAACGGCATTCCCGCCCTCCTTGAGGATATCCAc                                    >  1:20529180/1‑35 (MQ=255)
ctaaCGGCATTCCCGCCCTCCTTGAGGATATCCAc                                    <  1:3880129/33‑1 (MQ=255)
 cAACGGCATTCCCGCCCTCCTTGAGGATATCCACc                                   >  1:3579304‑M2/1‑34 (MQ=255)
  aaCGGCATTCCCGCCCTCCTTGAGGATATCCACcc                                  >  1:2106486‑M2/1‑33 (MQ=255)
   aCGGCATTCCCGCCCTCCTTGAGGATATCCACccc                                 >  1:19239467‑M2/1‑32 (MQ=255)
   aCGGCATTCCCGCCCTCCTTGAGGATATCCACccc                                 >  1:1642104‑M2/1‑32 (MQ=255)
    tggcATTCCCGCCCTCCTTGAGGATATCCACccca                                <  1:7353114‑M2/34‑5 (MQ=255)
    caacccttaagttagcgcttatgATATCCAccccc                                <  1:7207245‑M2/12‑1 (MQ=255)
    cGGCATTCCCGCCCTCCTTGAGGATATCCACccca                                >  1:16231546‑M2/1‑31 (MQ=255)
     ggCATTCCCGCCCTCCTTGAGGATATCCACCCCCa                               <  1:8907556/35‑1 (MQ=255)
     aacccttaagttagcgcttatgATATCCACCCCCa                               <  1:8719247‑M2/13‑1 (MQ=255)
      gCATTCCCGCCCTCCTTGAGGATATCCAAcccata                              >  1:13961622‑M2/1‑29 (MQ=255)
      acccttaagttagcgcttatgATATCCACCCCCAc                              <  1:15650415‑M2/14‑1 (MQ=255)
       cATTCCCGCCCTCCTTGAGGATATCCACcccataa                             <  1:1994146‑M2/35‑8 (MQ=255)
       cATTCCCGCCCTCCTTGAGGATATCCACCCCCAcc                             <  1:14072894/35‑1 (MQ=255)
        ccttaagttagcgcttatgATATCCACCCCCACCt                            <  1:19681203‑M2/16‑1 (MQ=255)
        ccttaagttagcgcttatgATATCCACCCCCACCt                            <  1:14950497‑M2/16‑1 (MQ=255)
        ccttaagttagcgcttatgATATCCACCCCCACCt                            <  1:5330067‑M2/16‑1 (MQ=255)
        aTTCCCGCCCTCCTTGAGGATATCCACcccataag                            <  1:9568272‑M2/35‑9 (MQ=255)
        aTTCCCGCCCTCCTTGAGGATATCCACCCCCACCt                            <  1:18363625/35‑1 (MQ=255)
        aTTCCCGCCCTCCTTGAGGATATCCACCCCCACCt                            <  1:7613652/35‑1 (MQ=255)
         ttCCCGCCCTCCTTGAGGATATCCACcccataagc                           >  1:6922249‑M2/1‑26 (MQ=255)
         ttCCCGCCCTCCTTGAGGATATCCACCCCCACCTg                           >  1:13176073/1‑35 (MQ=255)
         ttCCCGCCCTCCTTGAGGATATCCACCCCCACCTg                           >  1:17102818/1‑35 (MQ=255)
         cttaagttagcgcttatgATATCCACCCCCACCTg                           <  1:15986211‑M2/17‑1 (MQ=255)
         cttaagttagcgcttatgATATCCACCCCCACCTg                           <  1:2684643‑M2/17‑1 (MQ=255)
          tCCCGCCCTCCTTGAGGATATCCACCCCCACCTGa                          <  1:14394177/35‑1 (MQ=255)
          tCCCGCCCTCCTTGAGGATATCCACCCCCACCTGa                          >  1:5953009/1‑35 (MQ=255)
          tCCCGCCCTCCTTGAGGATATCCACCCCCACCTGa                          <  1:12423896/35‑1 (MQ=255)
           taagttagcgcttatgATATCCACCCCCACCTGag                         >  1:12997294‑M2/17‑35 (MQ=255)
           taagttagcgcttatgATATCCACCCCCACATGag                         <  1:20602169‑M2/19‑1 (MQ=255)
           cccGCCCTCCTTGAGGATATCCACcccataagcgc                         <  1:11711293‑M2/35‑12 (MQ=255)
            ccGCCCTCCTTGAGGATATCCACcccataagcgct                        <  1:18074974‑M2/35‑13 (MQ=255)
            ccGCCCTCCTTGAGGATATCCACcccataagcgct                        <  1:12713030‑M2/35‑13 (MQ=255)
            ccGCCCTCCTTGAGGATATCCACCCCCACCTGagt                        >  1:9652288/1‑35 (MQ=255)
             cGCCCTCCTTGAGGATATCCACcccataagcgcta                       >  1:17045165‑M2/1‑22 (MQ=255)
             cGCCCTCCTTGAGGATATCCACCCCCACCTGagtg                       <  1:20466816/35‑1 (MQ=255)
             agttagcgcttatgATATCCACCCCCACCTGagtg                       <  1:12458464‑M2/21‑1 (MQ=255)
             agttagcgcttatgATATCCACCCCCACCTGagtg                       >  1:9234131‑M2/15‑35 (MQ=255)
              gttagcgcttatgATATCCACCCCCACCTGagtgg                      <  1:2305173‑M2/22‑1 (MQ=255)
              gCCCTCCTTGAGCATATCCACCCCCACCTGagtgg                      >  1:784272/1‑35 (MQ=255)
              gCCCTCCNTGAGGATATCCACCCCCNCCTGagtgg                      >  1:5204290/1‑35 (MQ=255)
              ctcctcctTGAGGATATCCACCCCCACCTGagtgg                      <  1:10018741/33‑1 (MQ=255)
               ttagcgcttatgATATCCACCCCCACCTGagtggc                     <  1:15535501‑M2/23‑1 (MQ=255)
               cccTCCTTGAGGATATCCACCCCCACCTGagtggc                     <  1:4187498/35‑1 (MQ=255)
                tagcgcttatgATATCCACCCCCACCTGagtggca                    >  1:5400960‑M2/12‑35 (MQ=255)
                cctcctTGAGGATATCCACCCCCACCTGagtggca                    >  1:10761003/1‑35 (MQ=255)
                 ctcctTGAGGATATCCACcccataagcgctaactt                   >  1:4709469‑M2/1‑18 (MQ=255)
                 ctcctTGAGGATATCCACcccataagcgctaactt                   >  1:10404178‑M2/1‑18 (MQ=255)
                  tcctTGAGGATATCCACCCCCACCTGagtggcagt                  <  1:11540476/35‑1 (MQ=255)
                  gcgcttatgATATCCACCCCCACCTgagtggcagg                  >  1:576251‑M2/10‑34 (MQ=255)
                    gcttatgATATCCACCCCCACCTGagtggcagtgg                <  1:2754140‑M2/28‑1 (MQ=255)
                    ctTGAGGATATCCACcccataagcgctaacttaag                <  1:19629749‑M2/35‑21 (MQ=255)
                     ttGAGGATATCCACcccataagcgctaacttaagg               <  1:14508016‑M2/35‑22 (MQ=255)
                     ttGAGGATATCCACCCCCACCTGagtggcagtggc               <  1:11470381/35‑1 (MQ=255)
                     ttGAGGATATCCACCCCCACCTGagtggcagtggc               <  1:19376479/35‑1 (MQ=255)
                     ttGAGGATATCCACCCCCACCTGagtggcaggggg               >  1:7043363/1‑34 (MQ=255)
                     cttatgATATCCACCCACACCTGagtggcagtggc               <  1:10935883‑M2/29‑1 (MQ=255)
                      ttatgATATCCACCCCCACCTGAGTGGCAGTGGCg              <  1:4275925‑M2/30‑1 (MQ=255)
                      ttatgATATCCACCCCCACCTGAGTGGCAGTGGCg              <  1:18621781‑M2/30‑1 (MQ=255)
                      tGAGGATATCCACcccataagcgctaacttaaggg              <  1:17019764‑M2/35‑23 (MQ=255)
                      tGAGGATATCCACcccataagcgctaacttaaggg              <  1:3060388‑M2/35‑23 (MQ=255)
                      tGAGGATATCCACCCCCACCTGagtggcagggggg              >  1:8549444/1‑33 (MQ=38)
                       gAGGATATCCACaccataagcgctaacttaagggt             >  1:4357556‑M2/1‑12 (MQ=255)
                       gAGGATATCCACCCCCACCTGAGTGGCAGTGGCGt             <  1:3987416/35‑1 (MQ=255)
                       gAGGATATCCACCCCCACCTGAGTGGCAGTGGCGt             <  1:18000105/35‑1 (MQ=255)
                        aGGATATCCACcccataagcgctaacttaagggtt            <  1:15805813‑M2/35‑25 (MQ=255)
                        aGGATATCCACCCCCACCTGAGTGGCAGTGGCGTc            <  1:9729486/35‑1 (MQ=255)
                          gtgatCCCCCCCCGCCTGAGTGGCAGTGGCGTCca          <  1:20202456/32‑1 (MQ=25)
                          gATATCCACCCCCACCTGAGTGGCAGTGGCGTCCa          <  1:10788600‑M2/34‑1 (MQ=255)
                           atatCCACCCCCACCTGAGTGGCAGTGGCGTCcac         <  1:11874262/35‑1 (MQ=255)
                           atatCAACCCCCACCTGAGTGGCAGTGGCGTCcac         <  1:6788009/35‑1 (MQ=255)
                            tatCCCCCCCCACCTGAGTGGCAGTGGCGTCcaca        <  1:13559207/35‑1 (MQ=255)
                            tatCCACcccataagcgctaacttaagggttgtgg        >  1:18920703‑M2/1‑7 (MQ=255)
                             atCCACcccataagcgctaacttaagggttgtggt       <  1:10642343‑M2/35‑30 (MQ=255)
                             atCCACCCCCACCTGAGTGGCAGTGGCGTCCACAg       <  1:16074050/35‑1 (MQ=255)
                              tCCACcccataagcgctaacttaagggttgtggta      <  1:5067338‑M2/35‑31 (MQ=255)
                               ccACcccataagcgctaacttaagggttgtggtat     <  1:2817891‑M2/35‑32 (MQ=255)
                               ccACCCCCACCTGAGTGGCAGTGGCGTCCACAGAc     <  1:2949106/35‑1 (MQ=255)
                               ccACACCCCCCTGAGTGGCAGTGGCGTCCACAGAc     <  1:6835948/35‑1 (MQ=38)
                                 aCcccataagcgctaacttaagggttgtggtatta   <  1:18115747‑M2/35‑34 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:82570‑M2/1‑2 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:17412133‑M2/1‑2 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:16691976‑M2/1‑2 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:12033125‑M2/1‑2 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:12128879‑M2/1‑2 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:6204484‑M2/1‑2 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   <  1:344730‑M2/35‑34 (MQ=255)
                                 aCcccataagcgctaacttaagggttgtggtatta   >  1:14664863‑M2/1‑2 (MQ=255)
                                  ccccataatcgctaacttaagggttgtggtattac  <  1:4533329‑M2/35‑35 (MQ=255)
                                  ccccataagggctaacttaagggttgtggtattac  >  1:3793453‑M2/1‑1 (MQ=255)
                                  ccccataagcgctaacttaaggtttgtggtattac  <  1:17773658‑M2/35‑35 (MQ=255)
                                  ccccataagcgctaacttaagggttgtggtattac  >  1:7573233‑M2/1‑1 (MQ=255)
                                  ccccataagcgctaacttaagggttgtggtattac  >  1:6283330‑M2/1‑1 (MQ=255)
                                  ccccataagcgctaacttaagggttgtggtattac  <  1:3072927‑M2/35‑35 (MQ=255)
                                  ccccataagcgctaacttaagggttgtggtattac  >  1:18454627‑M2/1‑1 (MQ=255)
                                  |                                  
TCAACGGCATTCCCGCCCTCCTTGAGGATATCCACCCCCACCTGAGTGGCAGTGGCGTCCACAGACGCT  >  REL606/3517279‑3517347

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.