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

AAGCAGGGATTTAAATAAATCATTTATTACAACAGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAT  >  REL606/1420675‑1420743
                                  |                                  
ttaggcctgtgtccatattacgtgggtaggatcag                                    >  1:4997309‑M1/35‑35 (MQ=255)
ttaggcctgtgtccatattacgtgggtaggatcag                                    >  1:3249249‑M1/35‑35 (MQ=255)
ttaggcctgtgtccatattacgtgggtaggatcag                                    >  1:270751‑M1/35‑35 (MQ=255)
ttaggcctgtgtccatattacgcgggtaggatcag                                    <  1:846702‑M1/1‑1 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    <  1:15181109/35‑1 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    >  1:20310369/1‑35 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    >  1:2149145/1‑35 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    <  1:12368968/35‑1 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    <  1:25073/35‑1 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    >  1:11439608/1‑35 (MQ=255)
aaGCAGGGATTTAAATAAATCATTTATTACAACAg                                    >  1:3651252/1‑35 (MQ=255)
 aGCAGGGATTTAAATAAATCATTTATTACAACAgg                                   <  1:5509660/35‑1 (MQ=255)
 aGCAGGGATTTAAATAAATCATTTATTACAACAgg                                   >  1:20596729/1‑35 (MQ=255)
 aGCAGGGATTTAAATAAATCATTTATTACAACAgg                                   >  1:6038207/1‑35 (MQ=255)
 aGCAGGGATTTAAATAAATCATTTATTACAACAgg                                   >  1:19851397/1‑35 (MQ=255)
  gCAGGGATTTAAATAAATCATTTATTACAACAGGc                                  <  1:13024818/35‑1 (MQ=255)
  gCAGGGATTTAAATAAATCATTTATTACAACAGGc                                  <  1:18935314/35‑1 (MQ=255)
  gCAGGGATTTAAATAAATCATTTATTACAACAGGc                                  <  1:14774652/35‑1 (MQ=255)
   cAGGGATTTAAATAAATCATTTATTACAACAGGCT                                 <  1:13249009‑M1/35‑1 (MQ=255)
   cAGGGATTTAAATAAATCATTTATTACAACAGGCT                                 <  1:6681764‑M1/35‑1 (MQ=255)
    aGGGATTTAAATAAATCATTTATTACAACAGGCTa                                >  1:11023111/1‑35 (MQ=255)
    aGGGATTTAAATAAATCATTTATTACAACAGGCTa                                >  1:10071638/1‑35 (MQ=255)
     gggATTTAAATAAATCATTTATTACAACAGGCTAc                               <  1:14419043/35‑1 (MQ=255)
     gggATTTAAATAAATCATTTATTACAACAGGCTAc                               >  1:15864580/1‑35 (MQ=255)
      ggATTTAAATAAATCATTTATTGCAACAGGCTAcc                              <  1:8388073/35‑1 (MQ=255)
      ggATTTAAATAAATCATTTATTACAACAGGCTAcc                              >  1:1370189/1‑35 (MQ=255)
       gATTTAAATAAATCATTTATTACAACAGGCTAcat                             <  1:1257396/35‑3 (MQ=255)
       gATTTAAATAAATCATTTATTACAACAGGCTACCt                             <  1:5961773/35‑1 (MQ=255)
        gtgtccatattacgtgggtaggatcaGGCTACCtt                            >  1:16358749‑M1/27‑35 (MQ=255)
        aTTTAAATAAATCATTTATTACAACAGGCTACCtt                            <  1:1625811/35‑1 (MQ=255)
        aTTTAAATAAATCATTTATTACAACAGGCTACCtt                            <  1:9555453/35‑1 (MQ=255)
         tttAAATAAATCATTTATTACAACAGGCTACCTTc                           <  1:14864540/35‑1 (MQ=255)
         tttAAATAAATCATTTATTACAACAGGCTACCTTc                           <  1:14060952/35‑1 (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:17070938/35‑1 (MQ=255)
           taaataaaTCATTTATTACAACAGGCTACCTTCgg                         <  1:7392456/35‑1 (MQ=255)
           taaataaaTCATTTATTACAACAGGCTACCTTCgg                         <  1:16355890/35‑1 (MQ=255)
           taaataaaTCATTTATTACAACAGGCTACCTTCgg                         >  1:16050821/1‑35 (MQ=255)
            ccatattacgtgggtaggatcaGGCTACCTTCggg                        >  1:12725740‑M1/23‑35 (MQ=255)
            aaataaatCATTTATTACAACAGGCTgatcttacc                        >  1:12954301‑M1/1‑26 (MQ=255)
             aataaatCATTTATTACAACAGGCTACCTTCgggg                       >  1:3210635/1‑34 (MQ=255)
             aataaatCATTTATTACAACAGGCTACCTTCgggg                       >  1:14138509/1‑34 (MQ=255)
             aataaatCATTTATTACAACAGGCTACCTTCTGGt                       >  1:6000142/1‑35 (MQ=255)
             aataaatCATTTATTACAACAGGCTACCTTCGGGt                       >  1:11231995/1‑35 (MQ=255)
             aataaatCATTTATTACAACAGGCTACCTTCGGGt                       >  1:18201092/1‑35 (MQ=255)
             aataaatCATTTATTACAACAGGCTACCTTCGGGt                       >  1:13123100/1‑35 (MQ=255)
              ataaatCATTTATTACAACAGGCTACCTTCGGGTg                      <  1:9152614/35‑1 (MQ=255)
              ataaatCATTTATTACAACAGGCTACCTTCGGGTg                      >  1:10486392/1‑35 (MQ=255)
               taaatCATTTATTACAACAGGCTACCTTCGGGTgg                     >  1:6973893/1‑35 (MQ=255)
                aaatCATTTATTACAACAGGCTACCTTCGGGTGGc                    <  1:4328128/35‑1 (MQ=255)
                aaatCATTTATTACAACAGGCTACCTTCGGGTGGc                    >  1:6321106/1‑35 (MQ=255)
                aaatCATTTATTACAACAGGCTACCTTCGGGTGGc                    <  1:7424869/35‑1 (MQ=255)
                aaatCATTTATTACAACAGGCTACCTTCGGGTGGc                    >  1:8760034/1‑35 (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:6755049/1‑35 (MQ=255)
                    cgtgggtaggatcaGGCTACCTTCGGGGGGCCttt                <  1:19310640‑M1/21‑1 (MQ=255)
                     gtgggtaggatcaGGCTACCTTCGGGTGGCCtttt               >  1:4195193‑M1/14‑35 (MQ=255)
                     aTTTATTACAACAGGCTACCTTCGGGTGGCCtttt               >  1:20519891/1‑35 (MQ=255)
                     aTTTATTACAACAGGCTACCTTCGGGTGGCCtttt               >  1:4614020/1‑35 (MQ=255)
                      tttATTACAACAGGCTACCTTCGGGTGGCCttttt              >  1:9764168/1‑35 (MQ=255)
                      tttATTACAACAGGCTACCTTCGGGTGGCCttttt              <  1:14536891/35‑1 (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:15571742/1‑35 (MQ=255)
                       ttattaCAACAGGCTACCTTCGGGTGGCCtttttt             >  1:18615128/1‑35 (MQ=255)
                        tattaCAACAGGCTACCTTCGGGTGGCCTTTTTTa            <  1:1726656/35‑1 (MQ=255)
                        tattaCAACAGGCTACCTTCGGGTGGCCTTTTTTa            <  1:6570056/35‑1 (MQ=255)
                        tattaCAACAGGCTACCTTCGGGGGGCCTTTTTTa            >  1:7134170/1‑35 (MQ=255)
                         gtaggatcaGGCTACCTTCGGGTGGCCTTTTTTAt           <  1:10082081‑M1/26‑1 (MQ=255)
                         attaCAACAGGCTACCTTCGGGTGGCCTTTTTTAt           <  1:15420536/35‑1 (MQ=255)
                          ttaCAACAGGCTACCTTCGGGTGGCCTTTTTTAtt          >  1:4516084/1‑35 (MQ=255)
                           taCAACAGGCTACCTTCGGGTGGCCTTTTTTATtg         <  1:15935798/35‑1 (MQ=255)
                            acaacaGGCTgatcttacccagcaatagtggacac        >  1:18241532‑M1/1‑10 (MQ=255)
                            acaacaGGCTACCTTCGGGTGGCCTTTTTTATtgc        >  1:8229198/1‑35 (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:17667603/1‑35 (MQ=255)
                               tcaGGCTACCTTCGGGTGGCCTTTTTTATtgctgc     >  1:7127080‑M1/4‑35 (MQ=255)
                               tcaGGCTACCTTCGGGTGGCCTTTTTTATtgctgc     >  1:18539655‑M1/4‑35 (MQ=255)
                               acaGGCTgatcttacccagcaatagtggacacgcg     <  1:13884621‑M1/35‑29 (MQ=255)
                               acaGGCTACCTTCGGGTGGCCTTTTTTATtgctgc     <  1:18815966/35‑1 (MQ=255)
                                naGGCTACCTTCGGGTGGCCTTTTTTATTTCTGCg    <  1:7947799‑M1/33‑1 (MQ=255)
                                caGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCg    <  1:7326418‑M1/33‑1 (MQ=255)
                                 aGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGa   <  1:777406‑M1/34‑1 (MQ=255)
                                 aGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGa   <  1:2303163‑M1/34‑1 (MQ=255)
                                 aGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGa   >  1:16723613‑M1/2‑35 (MQ=255)
                                 aGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGa   <  1:10334803‑M1/34‑1 (MQ=255)
                                  ggCTACCTTCGGGTTGCCTTTTTTATTGCTGCGAt  >  1:3710844/1‑35 (MQ=255)
                                  ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt  <  1:1934999/35‑1 (MQ=255)
                                  ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt  <  1:10683959/35‑1 (MQ=255)
                                  ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt  <  1:158583/35‑1 (MQ=255)
                                  ggCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAt  >  1:5188534/1‑35 (MQ=255)
                                  |                                  
AAGCAGGGATTTAAATAAATCATTTATTACAACAGGCTACCTTCGGGTGGCCTTTTTTATTGCTGCGAT  >  REL606/1420675‑1420743

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.