LoslegenKostenlos loslegen

Partitioning the yeast genome

Genomes are often big, but interest usually lies in specific regions of them. Therefore, we need to subset a genome by extracting parts of it. To pick a sequence interval, use getSeq() and specify the name of the chromosome and the start and end of the sequence interval.

The following example will select the bases of "chrI" from 100 to 150.

getSeq(yeastGenome, names = "chrI", start = 100, end = 150)

Note: names is optional; if not specified, it will return all chromosomes. The parameters start and end are also optional and, if not specified, will take the default values 1 and the length of the sequence, respectively.

Diese Übung ist Teil des Kurses

Introduction to Bioconductor in R

Kurs anzeigen

Anleitung zur Übung

  • Use getSeq() to get the first 30 bases of the M chromosome ("chrM") in the yeastGenome object.

Interaktive Übung

Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.

# Load the yeast genome
library(BSgenome.Scerevisiae.UCSC.sacCer3)

# Assign data to the yeastGenome object
yeastGenome <- BSgenome.Scerevisiae.UCSC.sacCer3

# Get the first 30 bases of chrM
___
Code bearbeiten und ausführen