시작하기무료로 시작하기

fastq 파일에서 샘플 추출하기

이제 여러 리드로 이루어진 시퀀스에서 샘플 조각을 뽑아 볼 차례예요.

바로 이전 연습 문제에서 사용한 같은 파일을 씁니다. 이 파일에는 50 bp 길이의 리드가 500개 들어 있습니다. 파일 경로는 f라는 객체에 저장되어 있어요.

FastqSampler(con = file_path, n = length), set.seed(), yield()를 사용하면 시퀀스 파일에서 100개의 리드를 부분 추출할 수 있어요.

이 연습은 강의의 일부입니다

R로 시작하는 Bioconductor

강의 보기

연습 안내

  • ShortRead를 불러오세요.
  • set.seed()에 값 1234를 사용하세요.
  • f에 있는 작은 fastq 파일로 FastqSampler()를 사용하고, 100개의 리드를 선택하세요.
  • yield()를 사용해 부분 시퀀스를 생성하세요.

실습형 인터랙티브 연습

이 예제를 이 샘플 코드를 완성하여 풀어보세요.

# 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
___
코드 편집 및 실행