Reading BED files
In this exercise, you will load peak calls from a BED file and use the information about peak locations to extract reads that overlap with the peaks from a BAM file.
Deze oefening maakt deel uit van de cursus
ChIP-seq with Bioconductor in R
Oefeninstructies
- Use the
import.bedfunction to load peak calls from chr20_peaks. - Use these peak calls to create a
BamViewsobject for chr20_bam. - Load the reads overlapping with the peak calls.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Load peak calls from chr20_peaks
peaks <- ___(chr20_peaks)
# Create a BamViews object
bam_views <- ___(___, bamRanges=peaks)
# Load the reads
reads <- ___(___)