티스토리 뷰
Mothur SOP 를 따라하려면 Stability파일을 만들어야 하는데, 이게 생각보다 쉽지 않다. 아래는 이 파일을 쉽게 만들 수 있는 방법을 소개한다.
시퀀싱 파일이 각 샘플별로 이미 demultiplex 되어 있고, 파일 이름이 samplename_R1_001.fastq, samplename_R2_001.fastq 와 같은 식으로 되어 있다고 하자. 시퀀스 파일이 있는 폴더 안에서 아래와 같은 방법으로 stability파일을 만든다.
ls *R1_001.fastq | cut -f1 -d "_" > groupnames.txt
ls *R1_001.fastq > forward.txt
ls *R2_001.fastq > reverse.txt
paste groupnames.txt forward.txt > temp.txt
paste temp.txt reverse.txt > stability.files
이렇게 하면 statility.files 안에 샘플이름, Forward시퀀스, Reverse시퀀스 의 순서로 파일이 만들어 진다.
'앰플리콘 시퀀싱' 카테고리의 다른 글
Mothur 파이프라인 (pipeline) (17) | 2017.01.19 |
---|---|
R1, R2, I1, I2, 도대체 어떤게 시퀀스 파일이가? (4) | 2016.12.31 |
OTU 테이블 표준화 (normalization) (0) | 2016.12.30 |
QIIME을 이용한 키메라 제거 (8) | 2016.12.07 |
QIIME 사용 방법 (파이프라인) (6) | 2016.11.26 |
댓글