2019-11-12 21:17:52 +00:00
|
|
|
# bedtools
|
|
|
|
|
|
|
|
|
|
> 유전자 분석 작업을 위한 도구의 swiss-army knife. BAM, BED, GFF/GTF, VCF 형식으로 데이터를 교차, 그룹화, 변환 및 카운트하는 데 사용.
|
|
|
|
|
> 자세한 정보: <https://bedtools.readthedocs.io/en/latest/>.
|
|
|
|
|
|
|
|
|
|
- sequence의 strand를 기준으로 두개의 파일을 교차하고 결과를 {{path/to/output_file}}의 경로에 저장:
|
|
|
|
|
|
|
|
|
|
`bedtools intersect -a {{path/to/file_1}} -b {{path/to/file_2}} -s > {{path/to/output_file}}`
|
|
|
|
|
|
|
|
|
|
- 외부 조인이 왼쪽인 두개의 파일을 교차, 예시. {{file_1}}에서 각 기능을 보고하고 {{file_2}}와 겹치지 않으면 NULL:
|
|
|
|
|
|
|
|
|
|
`bedtools intersect -a {{path/to/file_1}} -b {{path/to/file_2}} -lof > {{path/to/output_file}}`
|
|
|
|
|
|
|
|
|
|
- 더 효율적인 알고리즘을 사용하여 두개의 사전 정렬된 파일을 교차:
|
|
|
|
|
|
|
|
|
|
`bedtools intersect -a {{path/to/file_1}} -b {{path/to/file_2}} -sorted > {{path/to/output_file}}`
|
|
|
|
|
|
|
|
|
|
- 첫 3열과 5열을 기준으로 {{path/to/file}}을 그룹화하여 6열을 요약:
|
|
|
|
|
|
|
|
|
|
`bedtools groupby -i {{path/to/file}} -c 1-3,5 -g 6 -o sum`
|
|
|
|
|
|
|
|
|
|
- bam-formated 파일을 bed-formated 파일로 변환:
|
|
|
|
|
|
|
|
|
|
`bedtools bamtobed -i {{path/to/file}}.bam > {{path/to/file}}.bed`
|
|
|
|
|
|
|
|
|
|
- {{file_2}}.bed와 가장 가까운 {{file_1}}.bed에서의 모든 기능을 찾고,그들의 거리와 추가 열을 기록 (입력 파일 정렬 필요):
|
|
|
|
|
|
2019-11-15 15:59:10 +00:00
|
|
|
`bedtools closest -a {{path/to/file_1}}.bed -b {{path/to/file_2}}.bed -d`
|
