LoslegenKostenlos loslegen

Eine Stichprobe aus einer fastq-Datei ziehen

Jetzt bist du dran, aus einer Sequenz mit vielen Reads ein Beispielstück zu ziehen.

Du verwendest dieselbe Datei wie in der vorherigen Übung. Diese Datei enthält 500 Reads, jeweils 50 bp lang. Der Dateipfad ist in einem Objekt namens f gespeichert.

Mit FastqSampler(con = file_path, n = length), set.seed() und yield() kannst du 100 Reads aus deiner Sequenzdatei auswählen.

Diese Übung ist Teil des Kurses

Einführung in Bioconductor in R

Kurs anzeigen

Anleitung zur Übung

  • Lade ShortRead.
  • Verwende set.seed() mit dem Wert 1234.
  • Nutze FastqSampler() mit der kleinen fastq-Datei in f und wähle 100 Reads aus.
  • Verwende yield(), um die Teilsequenz zu erzeugen.

Interaktive Übung

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

# Load ShortRead
library(ShortRead)

# Set a seed for sampling
___

# Use FastqSampler with f and select 100 reads
fs <- ___(con = ___, ___ = ___)

# Generate new sample yield
my_sample <- ___

# Print my_sample
___
Code bearbeiten und ausführen