根据染色体提取vcf中的信息：
bcftools view -r chr20,chr21,chr22 NA12878.vcf.gz -o NA12878.vcf.chr20_22.gz -O z
解释：
提取金标准变异NA12878.vcf.gz中20、21和22染色体的变异，保存到NA12878.vcf.chr20_22.gz文件中
-O z：表示使用压缩形式保存文件
提取之后的文件作为其它软件输入时经常需要索引，可按如下建立索引：
bcftools index -t NA12878.vcf.chr20_22.gz
解释：
-t：产生TBI格式的索引（如果没有-t选项，默认会产生csi的索引）
命令结束后会产生NA12878.vcf.chr20_22.gz.tbi的索引文件
————————————————
版权声明：本文为CSDN博主「bigcindy」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/Jwenxue/article/details/109330562
根据染色体提取vcf中的信息： bcftools view -r chr20,chr21,chr22 NA1287...
根据id提取vcf（提取重测序数据不会有日志等啰嗦行） 根据染色体提取vcf 根据位置提取vcf（pos文件含染色...
提取染色体片段 提取文件中的某几列 根据位置提取vcf文件对应位点的信息 提取某一列数值满足条件的列 提取某些样本...
建立索引，最新版的bcftools 1.11 支持不压缩排序，不过vcf 中染色体顺序要与header中的一致Ne...
1.提取个体list bcftools query -l SNP.filtered.vcf > ind.txt 2...
1、提取3k材料的名称 ''' bcftools query -l 3k.snp.vcf > 3k.list ''...
tabix index 提取指定染色体 BCFTOOLS 设置 统计有多少个SNP： 过滤 F_MISSING M...
bcftools view sample.vcf -Oz -o sample.vcf.gzbcftools ind...
文件1：需要提出的样品的id，文件名DB_sample.list 文件2：含有不同样本的变异位点的vcf文件 利用...
问题描述 在生物信息分析过程中，常常需要修改 vcf 染色体名称，即删除染色体名中的chr字符或添加chr字符到v...
本文标题：bcftools：根据染色体提取vcf中的信息
本文链接：https://www.haomeiwen.com/subject/waspzltx.html
网友评论