R/get_gene.R

#' Get gene in region
#' 
#' Get gene using \code{query_genes}; see \code{qtl2db}.
#' 
#' @param chr_id chromosome identifier
#' @param start start position in Mbp
#' @param stop  stop position in Mbp
#' 
#' @export
#' @importFrom dplyr filter
#' 
get_gene <- function(chr_id, start, stop) {
  out <- dplyr::filter(
    query_genes(chr_id, start, stop),
    !is.na(Name))
  class(out) <- c("feature_tbl", class(out))
  out
}
byandell/CCSanger documentation built on May 13, 2019, 9:26 a.m.