ChIP-seq with Bioconductor in R
Peter Humburg
Statistician, Macquarie University
Load peak calls
peaks <- import.bed(peak_file, genome="hg19")
Create a BamViews object
bam_views <- BamViews(bam_file, bamRanges=peaks)
Load the reads
reads <- readGAlignments(bam_views)
ChIP-seq with Bioconductor in R