Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,420,709 +GA :: IS3 (+) +4 bp 14.6% intergenic (+21/‑9) ynaK → / → ECB_01341 conserved protein/hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1420709 =60 (0.780)11 (0.120) 10/44 NT 13.0% intergenic (+21/‑12) ynaK/ECB_01341 conserved protein/hypothetical protein
?REL606 = 1443893 NA (NA)noncoding (1255/1255 nt) IS3 repeat region
* ? REL606 = 142071269 (0.810)11 (0.110) 11/46 NT 12.2% intergenic (+24/‑9) ynaK/ECB_01341 conserved protein/hypothetical protein
?REL606 1607918 = NA (NA)coding (15/1320 nt) ECB_01510 putative tail component of prophage

Only 100 of 183 total aligned reads displayed.

CAGGGATTTAAATAAATCATTTATTACAACAGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGATCCG  >  REL606/1420678‑1420746
                                  |                                  
cAGGGATTTAAATAAATCATTTATTACAACAGGCT                                    <  1:13249009‑M1/35‑1 (MQ=255)
cAGGGATTTAAATAAATCATTTATTACAACAGGCT                                    >  1:15606382‑M1/1‑35 (MQ=255)
 aGGGATTTAAATAAATCATTTATTACAACAGGCTa                                   >  1:10071638/1‑35 (MQ=255)
 aGGGATTTAAATAAATCATTTATTACAACAGGCTa                                   >  1:11023111/1‑35 (MQ=255)
  gggATTTAAATAAATCATTTATTACAACAGGCTAc                                  <  1:14419043/35‑1 (MQ=255)
  gggATTTAAATAAATCATTTATTACAACAGGCTAc                                  >  1:7364088/1‑35 (MQ=255)
  gggATTTAAATAAATCATTTATTACAACAGGCTAc                                  <  1:18530205/35‑1 (MQ=255)
   ggATTTAAATAAATCATTTATTGCAACAGGCTAcc                                 <  1:8388073/35‑1 (MQ=255)
   ggATTTAAATAAATCATTTATTACAACAGGCTAcc                                 >  1:1370189/1‑35 (MQ=255)
    gATTTAAATAAATCATTTATTACAACAGGCTAcat                                <  1:12261458/35‑3 (MQ=255)
    gATTTAAATAAATCATTTATTACAACAGGCTAcat                                <  1:1257396/35‑3 (MQ=255)
    gATTTAAATAAATCATTTATTACAACAGGCTACCt                                >  1:6956318/1‑35 (MQ=255)
    gATTTAAATAAATCATTTATTACAACAGGCTACCt                                <  1:5961773/35‑1 (MQ=255)
    gATTTAAATAAATCATTTATTACAACAGGCTACCt                                <  1:8616683/35‑1 (MQ=255)
     gtgtccatattacgtgggtaggatcaGGCTACCtt                               >  1:16358749‑M1/27‑35 (MQ=255)
     aTTTAAATAAATCATTTATTACAACAGGCTACCtt                               <  1:1625811/35‑1 (MQ=255)
     aTTTAAATAAATCATTTATTACAACAGGCTACCtt                               <  1:18934868/35‑1 (MQ=255)
      tttAAATAAATCATTTATTACAACATGCTACCTTc                              >  1:15862088/1‑35 (MQ=255)
      tttAAATAAATCATTTATTACAACAGGCTACCTTc                              <  1:14060952/35‑1 (MQ=255)
       ttAAATAAATCATTTATTACAACAGGCTACCTTCg                             >  1:4566623/1‑35 (MQ=255)
       ttAAATAAATCATTTATTACAACAGGCTACCTTCg                             >  1:7195448/1‑35 (MQ=255)
        tccatattacgtgggtaggatcaGGCTACCTTCgg                            >  1:8243402‑M1/24‑35 (MQ=255)
        taaataaaTCATTTATTACAACAGGCTgatcttac                            <  1:20108480‑M1/35‑9 (MQ=255)
        taaataaaTCATTTATTACAACAGGCTACCTTCgg                            >  1:16050821/1‑35 (MQ=255)
        taaataaaTCATTTATTACAACAGGCTACCTTCgg                            <  1:16355890/35‑1 (MQ=255)
        taaataaaTCATTTATTACAACAGGCTACCTTCgg                            <  1:17070938/35‑1 (MQ=255)
         aaataaatCATTTATTACAACAGGCTgatcttacc                           >  1:12954301‑M1/1‑26 (MQ=255)
          aataaatCATTTATTACAACAGGCTACCTTCGGGt                          >  1:13123100/1‑35 (MQ=255)
           ataaatCATTTATTACAACAGGCTACCTTCGGGTg                         >  1:10486392/1‑35 (MQ=255)
           ataaatCATTTATTACAACAGGCTACCTTCGGGTg                         <  1:9152614/35‑1 (MQ=255)
            taaatCATTTATTACAACAGGCTACCTTCGGGTgg                        >  1:6973893/1‑35 (MQ=255)
             attacgtgggtaggatcaGGCTACCTTCGGGTGGc                       >  1:3548724‑M1/19‑35 (MQ=255)
             aaatCATTTATTACAACAGGCTACCTTCGGGTGGc                       <  1:4328128/35‑1 (MQ=255)
             aaatCATTTATTACAACAGGCTACCTTCGGGTGGc                       <  1:7424869/35‑1 (MQ=255)
              aatCATTTATTACAACAGGCTgatcttacccagca                      >  1:11546969‑M1/1‑21 (MQ=255)
              aatCATTTATTACAACAGGCTACCTTCGGGTGGcc                      <  1:13659660/35‑1 (MQ=255)
               atCATTTATTACAACAGGCTACCTTCGGGGGGCCt                     >  1:14422836/1‑35 (MQ=255)
                tCATTTATTACAACAGGCTACCTTCGGGTGGCCtt                    <  1:20533476/35‑1 (MQ=255)
                tCATTTATTACAACAGGCTACCTTCGGGTGGCCtt                    >  1:6755049/1‑35 (MQ=255)
                 cATTTATTACAACAGGCTgatcttacccagcaata                   <  1:6636569‑M1/35‑18 (MQ=255)
                 cATTTATTACAACAGGCTgatcttacccagcaata                   >  1:1871483‑M1/1‑18 (MQ=255)
                  gtgggtaggatcaGGCTACCTTCGGGTGGCCtttt                  >  1:4195193‑M1/14‑35 (MQ=255)
                  gtgggtaggatcaGGCTACCTTCGGGGGGCCtttt                  <  1:6506581‑M1/22‑1 (MQ=255)
                  aTTTATTACAACAGGCTACCTTCGGGTGGCCtttt                  <  1:13980635/35‑1 (MQ=255)
                  aTTTATTACAACAGGCTACCTTCGGGTGGCCtttt                  <  1:12898822/35‑1 (MQ=255)
                  aTTTATTACAACAGGCTACCTTCGGGTGGCCtttt                  >  1:4614020/1‑35 (MQ=255)
                  aTTTATTACAACAGGCTACCTTCGGGTGGCCtttt                  >  1:20519891/1‑35 (MQ=255)
                   tttATTACAACAGGCTgatcttacccagcaatagt                 <  1:6059269‑M1/35‑20 (MQ=255)
                   tttATTACAACAGGCTACCTTAGGGTGGCCttttt                 <  1:3359008/35‑1 (MQ=255)
                   tgggtaggatcaGGCTACCTTCGGGTGGCCttttt                 <  1:13378089‑M1/23‑1 (MQ=255)
                    ttattaCAACAGGCTgatcttacccagcaatagtg                <  1:4018553‑M1/35‑21 (MQ=255)
                    ttattaCAACAGGCTACCTTCGGGTGGCCtttttt                >  1:16642782/1‑35 (MQ=255)
                    ttattaCAACAGGCTACCTTCGGGTGGCCtttttt                >  1:5077988/1‑35 (MQ=255)
                    ttattaCAACAGGCTACCTTCGGGTGGCCtttttt                >  1:14131479/1‑35 (MQ=255)
                    ttattaCAACAGGCTACCTTCGGGTGGCCtttttt                >  1:19833314/1‑35 (MQ=255)
                     tattaCAACAGGCTACCTTCGGGTGGCCTTTTTTa               <  1:10344306/35‑1 (MQ=255)
                     tattaCAACAGGCTACCTTCGGGTGGCCTTTTTTa               <  1:1726656/35‑1 (MQ=255)
                     tattaCAACAGGCTACCTTCGGGTGGCCTTTTTTa               <  1:5998228/35‑1 (MQ=255)
                      gtaggatcaGGCTACCTTCGGGTGGCCTTTTTTAt              <  1:10082081‑M1/26‑1 (MQ=255)
                      attaCAACAGGCTACCTTCGGGTGGCCTTTTTTAt              <  1:15420536/35‑1 (MQ=255)
                       ttaCAACAGGCTACCTTCGGGTGGCCTTTTTTAtt             <  1:19300463/35‑1 (MQ=255)
                        taCAACAGGCTACCTTCGGGTGGCCTTTTTTATtg            >  1:1485437/1‑35 (MQ=255)
                        taCAACAGGCTACCTTCGGGTGGCCTTTTTTATtg            <  1:15935798/35‑1 (MQ=255)
                         ggatcaGGCTACCTTCGGGTGGCCTTTTTTATtgc           <  1:14869056‑M1/29‑1 (MQ=255)
                         acaacaGGCTgatcttacccagcaatagtggacac           >  1:18241532‑M1/1‑10 (MQ=255)
                         acaacaGGCTACCTTCGGGTGGCCTTTTTTATtgc           <  1:15485120/35‑1 (MQ=255)
                         acaacaGGCTACCTTCGGGTGGCCTTTTTTATtgc           >  1:11792965/1‑35 (MQ=255)
                          caacaGGCTACCTTCTGGTGGCCTTTTTTATtgct          >  1:18889876/1‑35 (MQ=255)
                          caacaGGCTACCTTCGGGTGGCCTTTTTTATtgct          <  1:15369131/35‑1 (MQ=255)
                          caacaGGCTACCTTCGGGTGGCCTTTTTTATtgct          >  1:17411519/1‑35 (MQ=255)
                           aacaGGCTACCTTCGGGTGGCCTTTTTTATtgctg         <  1:3719674/35‑1 (MQ=255)
                           aacaGGCTACCTTCGGGTGGCCTTTTTTATtgctg         <  1:185686/35‑1 (MQ=255)
                           aacaGGCTACCTTCGGGTGGCCTTTTTTATtgctg         >  1:17667603/1‑35 (MQ=255)
                            tcaGGCTACCTTCGGGTGGCCTTTTTTATtgctgc        >  1:7127080‑M1/4‑35 (MQ=255)
                             naGGCTACCTTCGGGTGGCCTTTTTTATTTCTGCg       <  1:7947799‑M1/33‑1 (MQ=255)
                              aGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGa      <  1:2303163‑M1/34‑1 (MQ=255)
                              aGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGa      <  1:777406‑M1/34‑1 (MQ=255)
                               ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt     <  1:10683959/35‑1 (MQ=255)
                               ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt     >  1:5188534/1‑35 (MQ=255)
                               ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt     <  1:1934999/35‑1 (MQ=255)
                               ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt     <  1:9739928/35‑1 (MQ=255)
                                gCTACCTTCGGGTGGCCTTTTTTATTGCTGCGATc    >  1:12206909/1‑35 (MQ=255)
                                gCTACCTTCGGGTGGCCTTTTTTATTGCTGCGATc    >  1:7424315/1‑35 (MQ=255)
                                 cTgatcttacccagcaatagtggacacgcggctaa   <  1:5480941‑M1/35‑34 (MQ=255)
                                 cTgatcttacccagcaatagtggacacgcggctaa   >  1:8677302‑M1/1‑2 (MQ=255)
                                 cTgatcttacccagcaatagtggacacgcggctaa   >  1:2872768‑M1/1‑2 (MQ=255)
                                 cTgatcttacccagcaatagtggacacgcggctaa   <  1:18167789‑M1/35‑34 (MQ=255)
                                 cTACCTTCGGGTGGCCTTTTTTATTGCTGCGATcc   <  1:7252606/35‑1 (MQ=255)
                                 cTACCTTCGGGTGGCCTTTTTTATTGCTGCGATcc   <  1:11125791/35‑1 (MQ=255)
                                  tgatcttacccagcaatagtggacacgctgctaag  >  1:5702385‑M1/1‑1 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  <  1:10307756‑M1/35‑35 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  <  1:8821101‑M1/35‑35 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  <  1:9902548‑M1/35‑35 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  >  1:11286236‑M1/1‑1 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  <  1:13225826‑M1/35‑35 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  >  1:2915025‑M1/1‑1 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  <  1:18905772‑M1/35‑35 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  >  1:18444677‑M1/1‑1 (MQ=255)
                                  tgatcttacccagcaatagtggacacgcggctaag  >  1:16883588‑M1/1‑1 (MQ=255)
                                  tACCTTCGGGTGGCCTTTTTTATTGCTGCGATCCg  <  1:2390822/35‑1 (MQ=255)
                                  |                                  
CAGGGATTTAAATAAATCATTTATTACAACAGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGATCCG  >  REL606/1420678‑1420746

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