if (check.binary("bedtools")) {
index <- get.example.regions();
a <- index[[1]];
a <- bedr(engine = "bedtools", input = list(i = a), method = "sort", params = "");
b <- flank.region(a, n.add = 20);
}
Run the code above in your browser using DataLab