Lecture de fichiers BED
Dans cet exercice, vous allez charger des appels de pics à partir d'un fichier BED et utiliser l'information sur l'emplacement des pics pour extraire, d'un fichier BAM, les lectures qui chevauchent ces pics.
Cette activité fait partie du cours
ChIP-seq avec Bioconductor en R
Instructions de l’exercice
- Utilisez la fonction
import.bedpour charger les appels de pics depuis chr20_peaks. - Utilisez ces appels de pics pour créer un objet
BamViewspour chr20_bam. - Chargez les lectures qui chevauchent les appels de pics.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Load peak calls from chr20_peaks
peaks <- ___(chr20_peaks)
# Create a BamViews object
bam_views <- ___(___, bamRanges=peaks)
# Load the reads
reads <- ___(___)