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.
Diese Übung ist Teil des Kurses
<Kurs>ChIP-seq with Bioconductor in R</Kurs>Übungsanweisungen
- 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.
Interaktive praktische Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# Load peak calls from chr20_peaks
peaks <- ___(chr20_peaks)
# Create a BamViews object
bam_views <- ___(___, bamRanges=peaks)
# Load the reads
reads <- ___(___)