junk/findcontigs.R

findcontigs <- function( bedfile, mergewithin = 50, minsize = 500 ){
  mergebed <- bed.merge(bedfile, flank = mergewithin )
  mergebed <- bed.parselengths( mergebed, brks=c( minsize, Inf ) )
  return( mergebed )
}
dvera/travis documentation built on June 5, 2019, 5:12 a.m.