Get cells in a region

Extract cell names containing reads mapped within a given genomic region

Usage

GetCellsInRegion(tabix, region, cells = NULL)

Arguments

tabix

Tabix object

region

A string giving the region to extract from the fragments file

cells

Vector of cells to include in output. If NULL, include all cells

Value

Returns a list

Examples

fpath <- system.file("extdata", "fragments.tsv.gz", package = "Signac")
GetCellsInRegion(tabix = fpath, region = "chr1:10245-762629")
#> $`chr1:10245-762629`
#>  [1] "AAAGATGCACGCTGTG-1" "AAACTGCCAGTAGGCA-1" "AAACTGCAGGCAAGCT-1"
#>  [4] "AAAGATGAGACTAGGC-1" "AAACTCGAGAGGCCTA-1" "AAAGGATCAATACTGC-1"
#>  [7] "AAACGAATCTCGTAGA-1" "AAACTGCAGGCAAGCT-1" "AAAGATGGTGCAAGAC-1"
#> [10] "AAACGAACACCCTTTG-1" "AAAGATGAGATAGGTT-1" "AAACTCGTCTTCCACG-1"
#> [13] "AAACTCGAGGTCTTTG-1" "AAACTGCGTCCGCTTT-1" "AAAGATGCACCGTTGG-1"
#> [16] "AAACTGCAGATCTAAG-1" "AAACGAACACCCTTTG-1"
#>