728x90
반응형
1. unique하게 mapping된 read만 뽑기
samtools view -S -h -q 1 -F 4 ~~.sam -o ~~.sam
-F : 특정 flag의 read를 제거함
-f : 특정 flag의 read를 가짐
-q : 최소 minimum mapping quality score (MAPQ) 요구량
https://www.biostars.org/p/56246/
https://broadinstitute.github.io/picard/explain-flags.html # flag에 대한 값들 설명
728x90
반응형
'Bioinformatics(생정보학)' 카테고리의 다른 글
| GATK-somaticindeldetector (0) | 2017.03.31 |
|---|---|
| GATK-base quality score recalibration (0) | 2017.03.30 |
| gatk bundle설치 (0) | 2017.03.21 |
| Mutect 설치 (0) | 2017.03.20 |
| GATK(Genome Analysis Toolkit) 설치 (0) | 2017.03.20 |