码迷,mamicode.com
首页 > 其他好文 > 详细

使用bedtools提取vcf多个位置的变异(extract multi-region of genotypes by bedtools)

时间:2018-07-25 23:09:11      阅读:397      评论:0      收藏:0      [点我收藏+]

标签:div   class   信息   标准   int   染色体   下载安装   reg   head   

1、下载安装bedtools

2、生成bed文件;标准的bed文件格式如下:

chr7    127471196  127472363  Pos1  0  +  127471196  127472363  255,0,0
chr7    127472363  127473530  Pos2  0  +  127472363  127473530  255,0,0
chr7    127473530  127474697  Pos3  0  +  127473530  127474697  255,0,0
chr7    127474697  127475864  Pos4  0  +  127474697  127475864  255,0,0

如果你只有染色体、起始位置和终止位置信息的话,也无大碍。不大标准但是不伤大雅的bed文件格式如下:

chr7    127471196  127472363 
chr7    127472363  127473530  
chr7    127473530  127474697  
chr7    127474697  127475864  

3、提取多个位置的vcf文件;

bedtools intersect -a myfile.vcf.gz -b mutil-region.bed -header > output.vcf

 






使用bedtools提取vcf多个位置的变异(extract multi-region of genotypes by bedtools)

标签:div   class   信息   标准   int   染色体   下载安装   reg   head   

原文地址:https://www.cnblogs.com/chenwenyan/p/9368833.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!