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 138 total aligned reads displayed.

AGGCCTGTGTCCATATTACGTGGGTAGGATCAATCTCGGGCTGGACGTGCTGGATGTGGTTAACCT  >  REL606/1443862‑1443927
                               |                                  
aGGCCTGTGTCCATATTACGTGGGTAGGATCAAtc                                 <  1:17591161/35‑1 (MQ=32)
aGGCCTGTGTCCATATTACGTGGGTAGGATCAAtc                                 >  1:18239303/1‑35 (MQ=32)
aGGCCTGTGTCCATATTACGTGGGTAGGATCAAtc                                 <  1:10600630/35‑1 (MQ=32)
 nnccTGNGTCCATATTACGTGGGTAGGATCAAtcn                                <  1:19468137/33‑1 (MQ=255)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                <  1:17915480/35‑1 (MQ=35)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                >  1:18576722/1‑35 (MQ=35)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                >  1:11988478/1‑35 (MQ=35)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                >  1:9249800/1‑35 (MQ=35)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                >  1:16838799/1‑35 (MQ=35)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                <  1:19994543/35‑1 (MQ=35)
 ggCCTGTGTCCATATTACGTGGGTAGGATCAAtct                                >  1:16837490/1‑35 (MQ=35)
  gCCTGTGTCCATATTACGTGGGTAGGATCAAtctc                               >  1:2473481/1‑35 (MQ=38)
   ccTGTGTCCATATTACGTGGGTAGGATCAATCTCg                              <  1:9455450/35‑1 (MQ=40)
   ccTGTGTCCATATTACGTGGGTAGGATCAATCTCg                              >  1:10896264/1‑35 (MQ=40)
   ccTGTGTCCATATTACGTGGGTAGGATCAATCTCg                              <  1:3144205/35‑1 (MQ=40)
    cTGTGTCCATATTACGTGGGTAGGATCAATCTCgg                             <  1:9247396/35‑1 (MQ=255)
    cTGTGTCCATATTACGTGGGTAGGATCAATCTCgg                             >  1:8043761/1‑35 (MQ=255)
    cTGTGTCCATATTACGTGGGTAGGATCAATCTCgg                             >  1:4824284/1‑35 (MQ=255)
    cTGTGTCCATATGACGTGGGTAGGATCAATCTCgg                             <  1:7125043/35‑1 (MQ=255)
     tttgtCCATATTACGTGGGTAGGATCAATCTCggg                            >  1:13152507/3‑35 (MQ=255)
     tgtgTCCATATTACGTGGGTAGGATCAATCTCggg                            <  1:1753630/35‑1 (MQ=255)
     tgtgTCCATATTACGTGGGTAGGATCAATCTCggg                            <  1:16478141/35‑1 (MQ=255)
     tgtgTCCATATTACGTGGGTAGGATCAATCTCggg                            >  1:11833125/1‑35 (MQ=255)
     tgtgTCCATATTACGTGGGTAGGATCAATCTCggg                            <  1:20175076/35‑1 (MQ=255)
     tgtgTCCATATCACGTGGGTAGGATCAATCTCggg                            <  1:7473228/35‑1 (MQ=255)
      gtgtCCATATTACGTGGGTAGGATCAATCTCGGGc                           <  1:6649322/35‑1 (MQ=255)
      gtgtCCATATTACGTGGGTAGGATCAATCTCGGGc                           <  1:835076/35‑1 (MQ=255)
      gtgtCCATATTACGTGGGTAGGATCAATCTCGGGc                           >  1:14879779/1‑35 (MQ=255)
      gtgtCCATATTACGTGGGTAGGATCAATCTCGGGc                           >  1:6328617/1‑35 (MQ=255)
       tgtCCATATTACGTGGGTAGGATCAATCTCTGGCt                          >  1:5408798/1‑35 (MQ=255)
       tgtCCATATTACGTGGGTAGGATCAATCTCGGGCt                          <  1:19747283/35‑1 (MQ=255)
       tgtCCATATTACGTGGGTAGGATCAATCTCGGGCt                          >  1:10520497/1‑35 (MQ=255)
       tgtCCATATTACGTGGGTAGGATCAATCTCGGGCt                          <  1:6889593/35‑1 (MQ=255)
       tgtCCATATTACGTGGGTAGGATCAATCTCGGGCt                          >  1:17757923/1‑35 (MQ=255)
       tgtCCATATTACGTGGGTAGGATCAATCTCGGGCt                          >  1:13776481/1‑35 (MQ=255)
        gtCCATATTACGTGGGTAGGATCAATCTCTGGCTg                         >  1:3859125/1‑35 (MQ=255)
         tCCATATTACGTGGGTAGGATCAATCTCGGGCTgg                        >  1:17084239/1‑35 (MQ=255)
         tCCATATTACGTGGGTAGGATCAATCTCGGGCTgg                        <  1:12879745/35‑1 (MQ=255)
         tCCATATTACGTGGGTAGGATCAATCTCGGGCTgg                        <  1:4754160/35‑1 (MQ=255)
          ccATATTACGTGGGTAGGATCAATCTCGGGCTGGa                       >  1:20004894/1‑35 (MQ=255)
          ccATATTACGTGGGTAGGATCAATCTCGGGCTGGa                       <  1:12463821/35‑1 (MQ=255)
           cATATTACGTGGGTAGGATCAATCTCGGGCTGGAc                      >  1:9355968/1‑35 (MQ=255)
           cATATTACGTGGGTAGGATCAATCTCGGGCTGGAc                      <  1:3924863/35‑1 (MQ=255)
           cATATTACGTGGGTAGGATCAATCTCGGGCTGGAc                      >  1:13852241/1‑35 (MQ=255)
           cATATTACGTGGGTAGGATCAATCTCGGGCTGGAc                      <  1:15294822/35‑1 (MQ=255)
           cATATTACGTGGGTAGGATCAATCTCGGGCTGGAc                      >  1:16415078/1‑35 (MQ=255)
            ntatTACGTGGGTAGGATCAATCTCGGGCTGGACg                     >  1:8511458/2‑35 (MQ=255)
            atatTACGTGGGTAGGATCAATCTCGGGCTGGACg                     >  1:4323320/1‑35 (MQ=255)
            atatTACGTGGGTAGGATCAATCTCGGGCTGGACg                     >  1:14191643/1‑35 (MQ=255)
            atatTACGTGGGTAGGATCAATCTCGGGCTGGACg                     <  1:17207555/35‑1 (MQ=255)
            atatTACGTGGGTAGGATCAATCTCGGGCTGGACg                     <  1:14614165/35‑1 (MQ=255)
             tatTACGTGGGTAGGATCAATCTCGGGCTGGACGt                    >  1:946508/1‑35 (MQ=255)
             tatTACGTGGGTAGGATCAATCTCGGGCTGGACGt                    <  1:4890425/35‑1 (MQ=255)
             tatTACGTGGGTAGGATCAATCTCGGGCTGGACGt                    <  1:19699521/35‑1 (MQ=255)
             tatTACGTGGGTAGGATCAATCTCGGGCTGGACGt                    <  1:12384467/35‑1 (MQ=255)
              atTACGTGGGTAGGATCAATCTCGGGCTGGACggg                   >  1:13693369/1‑33 (MQ=255)
              atTACGTGGGTAGGATCAATCTCGGGCTGGACggg                   >  1:19536590/1‑33 (MQ=255)
              atTACGTGGGTAGGATCAATCTCGGGCTGGACGTg                   <  1:1956721/35‑1 (MQ=255)
              atTACGTGGGTAGGATCAATCTCGGGCTGGACGTg                   <  1:12830803/35‑1 (MQ=255)
              atTACGTGGGTAGGATCAATCTCGGGCTGGACGTg                   <  1:16855609/35‑1 (MQ=255)
              atTACGTGGGTAAGATCAATCTCGGGCTGGACGTg                   <  1:18461488/35‑1 (MQ=255)
               ttgcGTGAGTAGGATCAATCTCGGGCTGGACGTGc                  <  1:5612892/32‑1 (MQ=37)
                tACGTGGGTAGGATCAATCTCGGGCTGGACGTGCt                 >  1:10272710/1‑35 (MQ=255)
                 aCGTGGGTAGGATCAATCTCGGGCTGGACGTGCTg                <  1:18026182/35‑1 (MQ=255)
                 aCGTGGGTAGGATCAATCTCGGGCTGGACGTGCTg                <  1:1794711/35‑1 (MQ=255)
                 aCGTGGGTAGGATCAATCTCGGGCTGGACGGGCTg                >  1:10585167/1‑35 (MQ=255)
                  cGTGGGTAGGATCAATCTCGGGCTGGACGTGCTgg               <  1:3413481/35‑1 (MQ=255)
                    tGGGTAGGATCAATCTCGGGCTGGACGTGCTggct             >  1:15141814/1‑33 (MQ=255)
                    tGGGTAGGATCAATCTCGGGCTGGACGTGCTGGAt             <  1:2869184/35‑1 (MQ=255)
                    tGGGTAGGATCAATCTCGGGCTGGACGTGCTGGAt             >  1:6630698/1‑35 (MQ=255)
                     gggTAGGATCAATCTCGGGCTGGACGTGCTGGAtg            >  1:4153358/1‑35 (MQ=255)
                      ggTAGGATCAATCTCGGGCTGGACGTGGCGGAtag           >  1:679655/1‑33 (MQ=25)
                      ggTAGGATCAATCTCGGGCTGGACGTGCTGGAtgg           >  1:6191855/1‑34 (MQ=255)
                      ggTAGGATCAATCTCGGGCTGGACGTGCTGGAtgg           >  1:10718686/1‑34 (MQ=255)
                       gTAGGATCAATCTCGGGCTGGACGTGCTGGAtgtg          >  1:9800120/1‑35 (MQ=255)
                       gTAGGATCAATCTCGGGCTGGACGTGCTGGAtgtg          <  1:19823962/35‑1 (MQ=255)
                       gTAGGATCAATCTCGGGCTGGACGTGCTGGAtgtg          >  1:16532044/1‑35 (MQ=255)
                       gTAGGATCAATCTCGGGCTGGACGTGCTGGAtggg          >  1:11614773/1‑33 (MQ=255)
                        tAGGATCAATCTCGGGCTGGACGTGCTGGATGTgg         >  1:8321965/1‑35 (MQ=255)
                        tAGGATCAATCTCGGGCTGGACGTGCTGGATGGgg         >  1:11707597/1‑35 (MQ=255)
                         aGGATCAATCTCGGGCTGGACGTGCTGGATGTGGt        <  1:11967829/35‑1 (MQ=255)
                         aGGATCAATCTCGGGCTGGACGTGCTGGATGTGGt        <  1:18571648/35‑1 (MQ=255)
                         aGGATCAATCTCGGGCTGGACGTGCTGGATGTGGt        <  1:18421294/35‑1 (MQ=255)
                          ggATCAATCTCGGGCTGGACGTGCTGGATGTGGtt       <  1:18066338/35‑1 (MQ=255)
                          ggATCAATCTCGGGCTGGACGTGCTGGATGTGGtt       <  1:7646567/35‑1 (MQ=255)
                           gATCAATCTCGGGCTGGACGTGCTGGATGTGGTTa      <  1:1021462/35‑1 (MQ=255)
                           gATCAATCTCGGGCTGGACGTGCTGGATGTGGTTa      <  1:4643439/35‑1 (MQ=255)
                           gATCAATCTCGGGCTGGACGTGCTGGATGGGGttt      >  1:6868297/1‑34 (MQ=255)
                           gATCAATCTCGGGCTGGACGTGCTGGATGGGGGTa      >  1:10173737/1‑35 (MQ=37)
                            aTCAATCTCGGGCTGGACGTGCTGGATGTGGTTaa     <  1:14296505/35‑1 (MQ=255)
                            aTCAATCTCGGGCTGGACGTGCTGGATGTGGTTaa     <  1:4021836/35‑1 (MQ=255)
                            aTCAATCTCGGGCTGGACGTGCTGGATGGGGTTaa     >  1:8208636/1‑35 (MQ=255)
                            aTCAATCTCGGGCTGGACGTGCTGGATGGGGGTaa     >  1:5167586/1‑35 (MQ=37)
                             tCAATCTCGGGCTGGACGTGCTGGATGTGGTTAAc    <  1:4810966/35‑1 (MQ=255)
                             tCAATCTCGGGCTGGACGTGCTGGATGTGGGTaaa    >  1:9646521/1‑34 (MQ=255)
                              cAATCTCGGGCTGGACGTGCTGGATGTGGTTAAcc   <  1:5533750/35‑1 (MQ=255)
                               aaTCTCGGGCTGGACGTGCTGGATGTGGTTAACCt  <  1:17906483/35‑1 (MQ=255)
                               aaTCTCGGGCTGGACGTGCTGGATGTGGTTAACCt  <  1:14215035/35‑1 (MQ=255)
                               aaTCTCGGGCTGGACGTGCTGGATGTGGTTAACCt  <  1:12512319/35‑1 (MQ=255)
                               aaTCTCGGGCTGGACGTGCTGGATGTGGTTAACCt  <  1:4894308/35‑1 (MQ=255)
                               |                                  
AGGCCTGTGTCCATATTACGTGGGTAGGATCAATCTCGGGCTGGACGTGCTGGATGTGGTTAACCT  >  REL606/1443862‑1443927

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