Add a CpG density track by using bsseq to find all CpG sites in the genome and adding up CpG sites at each distance from the repeat start and end, plot as tile – light blue = high CpG content, dark blue = low CpG content

Removed 18 row(s) containing missing values (geom_path).Removed 18 rows containing missing values (geom_point).

LS0tCnRpdGxlOiAiQ29tcG9zaXRlIEVsZW1lbnQgTWV0aHlsYXRpb24gVEVMT2NsdXN0ZXJzIgphdXRob3I6ICJBcmllbCIKZGF0ZTogImByIGZvcm1hdChTeXMudGltZSgpLCAnJUIgJWQsICVZJylgIgpvdXRwdXQ6CiAgcGRmX2RvY3VtZW50OiBkZWZhdWx0CiAgaHRtbF9kb2N1bWVudDoKICAgIGRmX3ByaW50OiBwYWdlZAogIGh0bWxfbm90ZWJvb2s6IGRlZmF1bHQKbWFpbmZvbnQ6IERlamFWdSBTYW5zCnRpdGxlZm9udDogRGVqYVZ1IFNhbnMKLS0tCgpgYGB7ciBpbmNsdWRlPUZBTFNFfQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IEZBTFNFLCB3YXJuaW5nPUZBTFNFLCAgbWVzc2FnZT1GQUxTRSkKCiMgdGhlc2UgYXJlIG15IHBlcnNvbmFsIHV0aWxzIC0tIGZpbmQgdGhlbSBpbiB0aGUgdXRpbHMgZGlyIG9mIHRoZSBUMlQtRXBpZ2VuZXRpY3MgZ2l0IHJlcG8Kc291cmNlKCJ+L1QyVC1FcGlnZW5ldGljcy91dGlscy9pbGVlX3Bsb3RfdXRpbHMuUiIpCnNvdXJjZSgifi9UMlQtRXBpZ2VuZXRpY3MvdXRpbHMvbWV0aHlsYXRpb25fUl91dGlscy5SIikKbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkoY293cGxvdCkKbGlicmFyeSh6b28pCmxpYnJhcnkoYnNzZXEpCiMgQlNnZW5vbWUgcGFja2FnZSBmb3IgVDJUIGdlbm9tZSAtLSBmaW5kIG9uIGdsb2J1cyBhbmQgaW5zdGFsbApsaWJyYXJ5KEJTZ2Vub21lLnQydC52MS4wLnJlbGVhc2UpCgpmaWdzPSIva3liZXIvRGF0YS9OYW5vcG9yZS9BbmFseXNpcy9nbW9uZXkvQ0hNMTMvdjEuMF9maW5hbF9hc3NlbWJseS9maWd1cmVzIgpkYXQ9Ii9reWJlci9EYXRhL05hbm9wb3JlL0FuYWx5c2lzL2dtb25leS9DSE0xMy92MS4wX2ZpbmFsX2Fzc2VtYmx5IgpgYGAKCmBgYHtyIGxvYWREYXRhLCBpbmNsdWRlPUZ9CmRhdD0iL2t5YmVyL0RhdGEvTmFub3BvcmUvQW5hbHlzaXMvZ21vbmV5L0NITTEzL3YxLjBfZmluYWxfYXNzZW1ibHkiCmNobTEzX21ldGggPC0gcmVhZF90c3YocGFzdGUwKGRhdCwgIi9tZXRoeWxhdGlvbl9jYWxscy9tZXRoeWxhdGlvbl9mcmVxdWVuY3lfNTBrYl9zcGxpdC50c3YiKSkgJT4lCiAgR1JhbmdlcygpCgpyZWdzIDwtIHJlYWRfdHN2KHBhc3RlMChkYXQsICIvVEUvVGVsb0NvbXBvc2l0ZXMvVEVMTy1jb21wb3NpdGUtdW5pdEFfUk0tcG9saXNoX0FQUjIxMjFfU0hveXQudHh0IiksIGNvbF9uYW1lcyA9IFQpICU+JQogIG11dGF0ZShsZW49ZW5kLXN0YXJ0KSAlPiUKICBmaWx0ZXIobGVuPjApICU+JQogIGZpbHRlcihsZW48MzAwMCkgJT4lCiAgbXV0YXRlKElEPW4oKSkKCmBgYAoKYGBge3IgNWtiLCBlY2hvPUZ9CgojIHNpemUgb2YgZmxhbmtpbmcgcmVnaW9uCmZsYW5rbiA8LSAyMDAwMAojIHNpemUgb2YgcmVwZWF0IHJlbGF0aXZlIHRvIGZsYW5raW5nIHJlZ2lvbgpib2R5bGVuIDwtIDEvNQoKIyBhZGQgZmxhbmtzCmwxX3JlZ2lvbnMgPC0gcmVncyAlPiUKICBtdXRhdGUoc3RhcnQgPSBzdGFydCAtIGZsYW5rbiwgZW5kID0gZW5kICsgZmxhbmtuKSAlPiUKICBHUmFuZ2VzKCkKCiMgYWRkIG1ldGh5bGF0aW9uIGRhdGEKb3ZsIDwtIGZpbmRPdmVybGFwcyhjaG0xM19tZXRoLCBsMV9yZWdpb25zKQpnZW5lcy5vdmwgPC0gYXMuZGF0YS5mcmFtZShyZWdzKVtzdWJqZWN0SGl0cyhvdmwpLF0gJT4lCiAgZHBseXI6Om11dGF0ZShnZW5ld2lkdGggPSBlbmQgLSBzdGFydCkgJT4lCiAgZHBseXI6OnJlbmFtZShnZW5lX3N0YXJ0ID0gc3RhcnQsIGdlbmVfZW5kID0gZW5kKSAKCiMgY2FsY3VsYXRlIGRpc3RhbmNlIGZyb20gc3RhcnQgYW5kIGVuZCwgZGVwZW5kaW5nIG9uIHN0cmFuZCAKY2htMTMub3ZsIDwtIGFzLmRhdGEuZnJhbWUoY2htMTNfbWV0aFtxdWVyeUhpdHMob3ZsKSxdKSAlPiUKICBiaW5kX2NvbHMoZ2VuZXMub3ZsKSAlPiUKICBkcGx5cjo6cmVuYW1lKHNlcW5hbWVzID0gMSkgJT4lCiAgbXV0YXRlKGRpc3QgPSBpZmVsc2UoT3JpZW50YXRpb249PSJwb3MiLCBzdGFydCAtIGdlbmVfc3RhcnQsIGdlbmVfZW5kLXN0YXJ0KSkgJT4lCiAgbXV0YXRlKGRpc3QgPSBpZmVsc2UoZGlzdCA8IDAsIGRpc3QvZmxhbmtuLAogICAgICBpZmVsc2UoZGlzdCA8IGdlbmV3aWR0aCwKICAgICAgICBib2R5bGVuICogZGlzdCAvIGdlbmV3aWR0aCwKICAgICAgICBib2R5bGVuICsgKGRpc3QgLSBnZW5ld2lkdGgpL2ZsYW5rbikpLCAKICAgIGRpc3QgPSByb3VuZChkaXN0LDIpCiAgICApCgojIGJpbiBkYXRhIHRvIG1ha2UgcGxvdCBwcmV0dHkKbl93aW5kb3dzPTEwMDAKY2htMTMub3ZsJGN1dCA9IGN1dChjaG0xMy5vdmwkZGlzdCwgYnJlYWtzPW5fd2luZG93cykKCmNobTEzLm92bC5sYWJzIDwtIGNobTEzLm92bCAlPiUKICBncm91cF9ieShjdXQsYDQ5YnBfeWVzL25vYCkgJT4lCiAgc3VtbWFyaXNlKG1lZCA9IG1lZGlhbihtZXRoeWxhdGVkX2ZyZXF1ZW5jeSksIHRvcCA9IHF1YW50aWxlKG1ldGh5bGF0ZWRfZnJlcXVlbmN5LCAwLjc1KSwgYm90ID0gcXVhbnRpbGUobWV0aHlsYXRlZF9mcmVxdWVuY3ksIDAuMjUpLCBuX2dlbmVzID0gbGVuZ3RoKG1ldGh5bGF0ZWRfZnJlcXVlbmN5KSkgJT4lCiAgbXV0YXRlKHhfdG1wID0gc3RyX3N1YihjdXQsIDIsIC0yKSkgJT4lIAogIHNlcGFyYXRlKHhfdG1wLCBjKCJtaW4iLCAibWF4IiksIHNlcCA9ICIsIikgJT4lIAogIG11dGF0ZV9hdChjKCJtaW4iLCAibWF4IiksIGFzLmRvdWJsZSkgJT4lCiAgdW5ncm91cCgpICU+JQogIGdyb3VwX2J5KGA0OWJwX3llcy9ub2ApICU+JQogIGFycmFuZ2UobWluKSAlPiUKICBtdXRhdGUobWVkX3Ntb290aCA9IHJvbGxtZWFuKG1lZCwgMTAsIE5BKSkKCiMgcGxvdCAtLSBpbmNsdWRlIGFubm90YXRpb24gZm9yIG51bWJlciBvZiByZXBlYXRzIGluIGVhY2ggZ3JvdXAKcCA8LSBnZ3Bsb3QoY2htMTMub3ZsLmxhYnMsYWVzKCB4ID0gbWluLCB5ID0gbWVkX3Ntb290aCwgY29sb3I9YDQ5YnBfeWVzL25vYCwgZmlsbD1gNDlicF95ZXMvbm9gKSwgYWxwaGE9LjUpKwogIGdlb21fcmliYm9uKGFlcyh5bWluPWJvdCwgeW1heD10b3ApLCBhbHBoYT0wLjIsIHNpemU9MCkrCiAgZ2VvbV9saW5lKGFlcyh5PW1lZF9zbW9vdGgpLCBhbHBoYT0uNSwgc2l6ZT0xKSArCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMCkgKwogIGdlb21fdmxpbmUoeGludGVyY2VwdCA9IGJvZHlsZW4pKwogIGdlb21fcG9pbnQoYWxwaGE9LjUpKwogIGxhYnMoIHggPSAiR2Vub21pYyBQb3NpdGlvbiIsIHkgPSAiQWdncmVnYXRlZCBNZXRoeWxhdGlvbiBGcmVxdWVuY3kiKSArCiAgdGhlbWVfY2xhc3NpYygpK3NjYWxlX3hfY29udGludW91cyhicmVha3M9IGMoLTEsMCxib2R5bGVuLGJvZHlsZW4gKyAxKSwgbGFiZWxzID0gYyhwYXN0ZTAoIi0iLGZsYW5rbi8xZTMsImtiIiksIlN0YXJ0IiwiRW5kIixwYXN0ZTAoIisiLGZsYW5rbi8xZTMsImtiIikpKSsKZ2VvbV90ZXh0KGRhdGEgPSBjaG0xMy5vdmwgJT4lIAogICAgICAgICAgICAgIGdyb3VwX2J5KGA0OWJwX3llcy9ub2ApICU+JQogICAgICAgICAgICAgIHNlbGVjdChgNDlicF95ZXMvbm9gLGdlbmVfc3RhcnQpICU+JQogICAgICAgICAgICAgIGRpc3RpbmN0KCkgJT4lCiAgICAgICAgICAgICAgc3VtbWFyaXNlKG5fZ2VuZXM9bigpKSAlPiUKICAgICAgICAgICAgICBkaXN0aW5jdCgpICU+JQogICAgICAgICAgICAgIG11dGF0ZShoZWlnaHQ9Y2FzZV93aGVuKGA0OWJwX3llcy9ub2AgPT0gInllcyIgfiAwLCAKICAgICAgICAgICAgICAgICAgICAgICAgICBgNDlicF95ZXMvbm9gID09ICJubyIgfiAuMSwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgVFJVRSB+IC4wMikpLAogICAgICAgICAgICBhZXMoeD0xLAogICAgICAgICAgICAgICAgeT0uODUraGVpZ2h0LAogICAgICAgICAgICAgICAgbGFiZWw9cGFzdGUoIm4gPSIsY29tbWEobl9nZW5lcykpCiAgICAgICAgICAgICAgICApLCBoanVzdD0xKSsKdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gInRvcCIpCgoKcAoKCmdnc2F2ZSgKICBwYXN0ZTAoZmlncywgIi9ldm9sX21ldGgvVEUvIiwiQ0hNMTNfc3ViVEVMT19tZXRhcGxvdC5wZGYiKSwKICBwbG90ID0gcCwKICBzY2FsZSA9IDEsCiAgd2lkdGggPSAxMiwKICBoZWlnaHQgPSA1LAopCmBgYAoKCgpBZGQgYSBDcEcgZGVuc2l0eSB0cmFjayBieSB1c2luZyBic3NlcSB0byBmaW5kIGFsbCBDcEcgc2l0ZXMgaW4gdGhlIGdlbm9tZSBhbmQgYWRkaW5nIHVwIENwRyBzaXRlcyBhdCBlYWNoIGRpc3RhbmNlIGZyb20gdGhlIHJlcGVhdCBzdGFydCBhbmQgZW5kLCBwbG90IGFzIHRpbGUgLS0gbGlnaHQgYmx1ZSA9IGhpZ2ggQ3BHIGNvbnRlbnQsIGRhcmsgYmx1ZSA9IGxvdyBDcEcgY29udGVudCAKCmBgYHtyfQoKIyBjb3VudCBDRyBzaXRlcyBmcm9tIGJzc2VxIHBhY2thZ2UgLS0gdGFrZXMgYSBob3Qgc2VjLCBzYXZlIGFzIHJicyBmaWxlIGZvciBmdXR1cmUgdXNlCnQydC5jcGcubG9jaSA8LSBmaW5kTG9jaShwYXR0ZXJuID0gIkNHIiwKICAgICAgICAgICAgICAgICAgICAgICAgIHN1YmplY3QgPSBCU2dlbm9tZS50MnQudjEuMC5yZWxlYXNlOjpCU2dlbm9tZS50MnQudjEuMC5yZWxlYXNlLCAKICAgICAgICAgICAgICAgICAgICAgICAgIHN0cmFuZCA9ICIrIikKIyB5b3Ugc2hvdWxkIGNvdW50IH4zMk0gQ3BHcwpsZW5ndGgodDJ0LmNwZy5sb2NpKQoKIyBub3cgZmluZCBvdmVybGFwcyBiZXR3ZWVuIENwRyBzaXRlcyBhbmQgcmVwZWF0cwpsMV9yZWdpb25zIDwtIHJlZ3MgJT4lCiAgbXV0YXRlKHN0YXJ0ID0gc3RhcnQgLSBmbGFua24sIGVuZCA9IGVuZCArIGZsYW5rbikgJT4lCiAgR1JhbmdlcygpCgpvdmwgPC0gZmluZE92ZXJsYXBzKHQydC5jcGcubG9jaSwgbDFfcmVnaW9ucykKZ2VuZXMub3ZsIDwtIGFzLmRhdGEuZnJhbWUocmVncylbc3ViamVjdEhpdHMob3ZsKSxdICU+JQogIGRwbHlyOjptdXRhdGUoZ2VuZXdpZHRoID0gZW5kIC0gc3RhcnQpICU+JQogIGRwbHlyOjpyZW5hbWUoZ2VuZV9zdGFydCA9IHN0YXJ0LCBnZW5lX2VuZCA9IGVuZCkgCgpjaG0xMy5vdmwgPC0gYXMuZGF0YS5mcmFtZSh0MnQuY3BnLmxvY2lbcXVlcnlIaXRzKG92bCksXSkgJT4lCiAgYmluZF9jb2xzKGdlbmVzLm92bCkgJT4lCiAgZHBseXI6OnJlbmFtZShzZXFuYW1lcyA9IDEpICU+JQogIGRwbHlyOjptdXRhdGUoZGlzdCA9IGlmZWxzZShPcmllbnRhdGlvbiA9PSAicG9zIixzdGFydCAtIGdlbmVfc3RhcnQsIGdlbmVfZW5kIC0gc3RhcnQpLAogICAgZGlzdCA9IGlmZWxzZShkaXN0IDwgMCwgZGlzdC9mbGFua24sCiAgICAgIGlmZWxzZShkaXN0IDwgZ2VuZXdpZHRoLAogICAgICAgIGJvZHlsZW4gKiBkaXN0IC8gZ2VuZXdpZHRoLAogICAgICAgIGJvZHlsZW4gKyAoZGlzdCAtIGdlbmV3aWR0aCkvZmxhbmtuKSksIAogICAgZGlzdCA9IHJvdW5kKGRpc3QsMikKICAgICkKCiMgcGxvdCBDcEcgZGVuc2l0eSBhcyBoZWF0bWFwIC0tIHRoZSBhY3R1YWwgbnVtYmVyIGRvZXNudCBtZWFuIGFueXRoaW5nLCBpdHMganVzdCBmb3IgdmlzdWFsaXphdGlvbgpjaG0xMy5vdmwubGFicyA8LSBjaG0xMy5vdmwgJT4lCiAgZ3JvdXBfYnkoSUQpICU+JQogIG11dGF0ZSh0b3RhbD1uKCkpICU+JQogIHVuZ3JvdXAoKSAlPiUKICBncm91cF9ieShkaXN0KSAlPiUKICBzdW1tYXJpc2UobWVkID0gc3VtKHdpZHRoKSwgY2c9bWVkL3RvdGFsKSAlPiUKICBkaXN0aW5jdCgpCiAgCnEgPC0gZ2dwbG90KGNobTEzLm92bC5sYWJzLCBhZXMoeD1kaXN0LHk9MSwgZmlsbD1jZykpK2dlb21fdGlsZSgpK3NjYWxlX3hfY29udGludW91cyhicmVha3M9IGMoLTEsMCxib2R5bGVuLGJvZHlsZW4gKyAxKSwgbGFiZWxzID0gYyhwYXN0ZTAoIi0iLGZsYW5rbi8xZTMsImtiIiksIlN0YXJ0IiwiRW5kIixwYXN0ZTAoIisiLGZsYW5rbi8xZTMsImtiIikpKSt0aGVtZV92b2lkKCkrdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKQoKIyB1c2UgY293cGxvdCB0byBjb21iaW5lIHBsb3RzIG9uIHRoZSBzYW1lIFgtYXhpcwpwbG90X2dyaWQocCxxLCBuY29sPTEsIGFsaWduID0gInYiLHJlbF9oZWlnaHRzID0gYygzLC41KSkKYGBgCgoKCgoKCgo=