query a bedbase bigbed file

query_bb(md5sum, ...)

Arguments

md5sum	character(1) md5sum used as key into bedbase resources
...	passed to rtracklayer::import.bb, which will have format="bigBed" set

Value

GRanges

Examples

b1 = query_bb(md5sum="78c0e4753d04b238fc07e4ebe5a02984", 
    which=GenomicRanges::GRanges("chr17:37000000-39000000"))
b1
#> GRanges object with 70 ranges and 6 metadata columns:
#>        seqnames            ranges strand |        name     score signalValue
#>           <Rle>         <IRanges>  <Rle> | <character> <integer>   <numeric>
#>    [1]    chr17 37288824-37289289      * |           .      1000     71.8859
#>    [2]    chr17 37288824-37289289      * |           .      1000    288.8652
#>    [3]    chr17 37474227-37474736      * |           .       812    122.0867
#>    [4]    chr17 37475042-37475678      * |           .       667     64.2821
#>    [5]    chr17 37477039-37477947      * |           .      1000    125.7743
#>    ...      ...               ...    ... .         ...       ...         ...
#>   [66]    chr17 38867635-38868271      * |           .       538     54.7645
#>   [67]    chr17 38873198-38873637      * |           .      1000    370.8526
#>   [68]    chr17 38991100-38991632      * |           .      1000     23.2802
#>   [69]    chr17 38991100-38991632      * |           .      1000    257.4294
#>   [70]    chr17 38992296-38992932      * |           .       543     72.8908
#>           pValue    qValue      peak
#>        <numeric> <numeric> <integer>
#>    [1]        -1   3.64371       407
#>    [2]        -1   4.66339       205
#>    [3]        -1   4.27231       204
#>    [4]        -1   3.53542       318
#>    [5]        -1   4.25559       679
#>    ...       ...       ...       ...
#>   [66]        -1   3.28290       318
#>   [67]        -1   4.66339       215
#>   [68]        -1   1.13672       474
#>   [69]        -1   4.66339       275
#>   [70]        -1   3.72042       318
#>   -------
#>   seqinfo: 34 sequences from an unspecified genome