Hi

I am using the following Centipede tutorial for ATAC foot printing and wondering if I could get some advice:

http://slowkow.com/CENTIPEDE.tutorial/#determine_if_tf_sites_are_bound

I have successfully gotten some number of sites with a posterior probability of 1 out of my total set and now have cen and fit value lists in R. So for instance out of 150K peaks I have something like 25K sites with a PostPr of 1 and I would like to pull out those specific sites in a BED file but not exactly sure the best way to do it from the fit values. If anyone has any ideas that would be very helpful, thank you.

Rob.