Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 3,454,489 IS150 (+) +4 bp 38.1% coding (342‑345/561 nt) yrfE ← ADP‑ribose diphosphatase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3454489 =74 (0.860)46 (0.450) 30/48 NT 34.1% coding (345/561 nt) yrfE ADP‑ribose diphosphatase
?REL606 = 3652533 NA (NA)noncoding (1443/1443 nt) IS150 repeat region
* ? REL606 = 345449274 (0.860)45 (0.440) 28/48 NT 33.6% coding (342/561 nt) yrfE ADP‑ribose diphosphatase
?REL606 3651091 = NA (NA)noncoding (1/1443 nt) IS150 repeat region

GTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCT  >  REL606/3652502‑3652567
                               |                                  
gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                 >  1:14392514/1‑35 (MQ=35)
gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                 <  1:13169817/35‑1 (MQ=35)
gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                 >  1:4686650/1‑35 (MQ=35)
gtTTAACTGTCCAACTTTTTGGGGTAAGTACAAAc                                 <  1:272115/35‑1 (MQ=255)
 tttAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                <  1:4716545/35‑1 (MQ=38)
 tttAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                >  1:3576082/1‑35 (MQ=38)
 tttAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                >  1:9162090/1‑35 (MQ=38)
 tttAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                <  1:7142397/35‑1 (MQ=38)
  ttAACTGTCCAACTTTTTGGGGTCAGTAGAAACtt                               >  1:16172697/1‑35 (MQ=255)
  ttAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                               <  1:18299561/35‑1 (MQ=40)
  ttAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                               <  1:1866610/35‑1 (MQ=40)
  ttAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                               >  1:17973236/1‑35 (MQ=40)
  ttAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                               <  1:3653994/35‑1 (MQ=40)
   tAACTGTCCAACTTTTTGGGGTCAGTACAAACttt                              <  1:12668100/35‑1 (MQ=255)
   tAACTGTCCAACTTTTTGGGGTCAGTACAAACttt                              >  1:18248759/1‑35 (MQ=255)
   tAACTGTCCAACTTTTTGGGGTCAGTACAAACttt                              >  1:5566097/1‑35 (MQ=255)
    aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                             >  1:14438965/1‑35 (MQ=255)
    aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                             <  1:2085632/35‑1 (MQ=255)
    aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                             >  1:3860264/1‑35 (MQ=255)
      cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                           >  1:17396575/1‑35 (MQ=255)
      cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                           >  1:10215734/1‑35 (MQ=255)
      cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                           >  1:13977532/1‑35 (MQ=255)
      cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                           >  1:13535338/1‑35 (MQ=255)
      cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                           >  1:12755630/1‑35 (MQ=255)
       tGTCCAACTTTTTGGGGTCAGTACAAACTTTGAtt                          >  1:9335185/1‑35 (MQ=255)
       tGTCCAACTTTTTGGGGTCAGTACAAACCTTGAtt                          >  1:4523321/1‑35 (MQ=255)
        gTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                         >  1:10439090/1‑35 (MQ=255)
        gTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                         <  1:11457002/35‑1 (MQ=255)
        gTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                         >  1:16042161/1‑35 (MQ=255)
        gTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                         <  1:5508699/35‑1 (MQ=255)
         tCCAACTTTTTGGGGTCAGTACAAACTTTGATTta                        <  1:1091219/35‑1 (MQ=255)
         tCCAACTTTTTGGGGTCAGTACAAACTTTGATTta                        >  1:14579969/1‑35 (MQ=255)
          ccAACTTTTTGGGGTCAGTACAAACTTTGATTtat                       >  1:2881013/1‑35 (MQ=255)
          ccAACTTTTTGGGGTCAGTACAAACTTTGATTtat                       >  1:4693686/1‑35 (MQ=255)
           cAACTTTTTGGGGTCAGTACAAACTTTGATTtata                      >  1:1915204/1‑35 (MQ=255)
            aaCTTTTTGGTGTCAGTACAAACTTTGATTTATAg                     >  1:18026841/1‑35 (MQ=255)
            aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAg                     <  1:18139721/35‑1 (MQ=255)
            aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAg                     <  1:20145590/35‑1 (MQ=255)
            aaCTTTTTGGGGTCAGTACAAACTTTGATTTATAg                     <  1:3517686/35‑1 (MQ=255)
             aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                    <  1:16326421/35‑1 (MQ=255)
             aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                    <  1:2308921/35‑1 (MQ=255)
             aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                    >  1:11772182/1‑35 (MQ=255)
              cTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                   >  1:16051952/1‑35 (MQ=255)
              cTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                   >  1:18902919/1‑35 (MQ=255)
              cTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                   >  1:13120606/1‑35 (MQ=255)
               tttttGGGGTCAGTACAAACTTTGATTTATAGTCa                  <  1:8695499/35‑1 (MQ=255)
                ttttGGGGTCAGTACAAACTTTGATTTATAGTCAg                 <  1:18340144/35‑1 (MQ=255)
                ttttGGGGTCAGTACAAACTTTGATTTATAGTCAg                 <  1:17990599/35‑1 (MQ=255)
                 tttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                <  1:12060688/35‑1 (MQ=255)
                 tttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                <  1:14728234/35‑1 (MQ=255)
                 tttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                >  1:2167411/1‑35 (MQ=255)
                 tttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                <  1:19738911/35‑1 (MQ=255)
                 tttGGGGTCAGTACAAACTTTGATTTATAGTCAgg                <  1:20337575/35‑1 (MQ=255)
                  ttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt               <  1:12962976/35‑1 (MQ=255)
                  ttGGGGTCAGTACAAACTTTGATTTATAGTCAGGt               <  1:4112754/35‑1 (MQ=255)
                   tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTg              <  1:8855925/35‑1 (MQ=255)
                   tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTg              >  1:9392748/1‑35 (MQ=255)
                   tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTg              >  1:10023995/1‑35 (MQ=255)
                    ggggTCAGTACAAACTTTGATTTATAGTCAGGTgg             >  1:12154059/1‑35 (MQ=255)
                     gggTCAGTACAAACTTTGATTTATAGTCAGGTggg            >  1:6161140/1‑35 (MQ=255)
                      ggTCAGTACAAACTTTGATTTATAGTCAGGTgggg           <  1:19392191/35‑1 (MQ=255)
                      ggTCAGTACAAACTTTGATTTATAGTCAGGTgggg           >  1:17359466/1‑35 (MQ=255)
                       gTCAGTACAAACTTTGATTTATAGTCAGGTGGGGc          >  1:20223661/1‑35 (MQ=255)
                       gTCAGTACAAACCTTGATTTATAGTCAGGTGGGGc          <  1:9841823/35‑1 (MQ=255)
                        tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt         >  1:12948715/1‑35 (MQ=255)
                        tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt         <  1:9773238/35‑1 (MQ=255)
                        tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt         <  1:19403271/35‑1 (MQ=255)
                        tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt         <  1:19405325/35‑1 (MQ=255)
                        tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt         <  1:13184045/35‑1 (MQ=255)
                        tCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt         >  1:11412078/1‑35 (MQ=255)
                         cAGTACAAACTTTGATTTATAGTCAGGTGGGGCtt        >  1:5809715/1‑35 (MQ=255)
                         cAGTACAAACTTTGATTTATAGTCAGGTGGGGCtt        >  1:8932020/1‑35 (MQ=255)
                         cAGTACAAACTTTGATTTATAGTCAGGTGGGGCtt        <  1:11192931/35‑1 (MQ=255)
                         cAGTACAAACTTTGATTTATAGTCAGGTGGGGCtt        >  1:16367333/1‑35 (MQ=255)
                          aGTACAAACTTTGATTTATAGTCAGGTGGGTCttt       >  1:18724555/1‑35 (MQ=255)
                          aGTACAAACTTTGATTTATAGTCAGGTGGGGCttt       >  1:15823963/1‑35 (MQ=255)
                          aGTACAAACTTTGATTTATAGTCAGGTGGGGCttt       <  1:16063626/35‑1 (MQ=255)
                          aGTACAAACTTTGATTTATAGTCAGGTGGGGCttt       <  1:19819840/35‑1 (MQ=255)
                           gTACAAACTTTGATTTATAGTCAGGTGGGGCtttt      <  1:14117258/35‑1 (MQ=255)
                           gTACAAACTTTGATTTATAGTCAGGTGGGGCtttt      <  1:16802641/35‑1 (MQ=255)
                            tACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc     <  1:10460608/35‑1 (MQ=255)
                            tACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc     <  1:8341216/35‑1 (MQ=255)
                            tACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc     <  1:5826271/35‑1 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCtttttt    >  1:14996337/1‑33 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtct    <  1:812155/35‑1 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtct    <  1:8672472/35‑1 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtct    >  1:14623740/1‑35 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtct    >  1:171295/1‑35 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtct    >  1:3463600/1‑35 (MQ=255)
                             aCAAACTTTGATTTATAGTCAGGTGGGGCTTTtct    >  1:17655707/1‑35 (MQ=255)
                              cAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc   <  1:2644758/35‑1 (MQ=255)
                              cAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc   >  1:9299535/1‑35 (MQ=255)
                               aaaCTTTGATTTATAGTCAGGTGGGGCTTTTctct  >  1:7552206/1‑35 (MQ=255)
                               aaaCTTTGATTTATAGTCAGGTGGGGCTTTTctct  <  1:5292731/35‑1 (MQ=255)
                               aaaCTTTGATTTATAGTCAGGTGGGGCTTTTctct  <  1:14351732/35‑1 (MQ=255)
                               |                                  
GTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCT  >  REL606/3652502‑3652567

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.