Skip to content

Instantly share code, notes, and snippets.

@sirusb

sirusb/ScanParam_exmaple.r

Created February 20, 2015 15:26
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save sirusb/e3922911ff5281deec4f to your computer and use it in GitHub Desktop.
Save sirusb/e3922911ff5281deec4f to your computer and use it in GitHub Desktop.
Download ZIP
Raw
ScanParam_exmaple.r
require(pasillaBamSubset)
require(Rsamtools)
bamFile <- untreated1_chr4()
region <- GRanges("chr4",IRanges(100000L, 110000L))
param <- ScanBamParam(which=region, what=scanBamWhat())
mappedReads <- scanBam(bamFile, param=param)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment