使用bcftools提取指定样本的vcf文件（extract specified samples in vcf format）

1、下载安装bcftools。

2、准备样本ID文件，这里命名为samplelistname.txt，一个样本一行，如下所示：

sample1

sample2

sample3

3、输入命令：

bcftools view -S samplelistname.txt /1000genomes/ALL.chr16.phase3_shapeit2_mvncall_integrated_v5a.20130502.genotypes.vcf.gz -Ov > samplelist_1000Genomes.vcf

参考链接：

https://www.biostars.org/p/184950/

https://samtools.github.io/bcftools/bcftools.html#view

转载于:https://www.cnblogs.com/chenwenyan/p/9212170.html

使用bcftools提取指定样本的vcf文件（extract specified samples in vcf format）相关推荐

bat 批量提取指定目录下的文件
如下是批量获取指定目录下的文件名. @echo off echo text inputset input= set /p input=: echo %input% is input cd %input ...
R语言丨根据VCF文件设计引物，自动识别两样本差异SNP位点，调用samtools获取上下游参考序列
根据变异位点设计引物序列今天碰到一个新问题:假如有一个vcf文件储存了两个样品的变异位点基因型数据,每行代表一个位点,我现在想找出两样本差异的SNP位点,再把差异位点用[REF/ALT]的形式表示, ...
python提取数据包中的文件_Python-对Pcap文件进行处理，获取指定TCP流
通过对TCP/IP协议的学习,本人写了一个可以实现对PCAP文件中的IPV4下的TCP流提取,以及提取指定的TCP流,鉴于为了学习,没有采用第三方包解析pcap,而是对bytes流进行解析,其核心思想 ...
java批量提取文件夹名称_bat 批量提取指定目录下的文件名
bat 批量提取指定目录下的文件名下面是批量获取指定目录下的文件名的核心代码 @echo off echo text input set input= set /p input=: echo %in ...
2021.06.08|提取、比较各样品vcf文件中snp突变频率
目录摘要环境与方法使用代码分析结果总结摘要接到一个wgs项目,要帮助客户统计vcf文件中snp突变频率,比较两个样品的突变位点.这个工作在上一个项目中是手动处理的,当时参考序列短,突变位 ...
python csv文件和xlsx文件混杂时，提取指定列数据并合并
这篇文章是R语言如何合并csv文件(批量读取csv文件)的姊妹篇.提供更为强大的功能. 解决的问题是提取csv文件和xlsx文件混杂时,文件合并问题.具体来说,下面的代码是提取指定列文本,简单清洗后 ...
python实现批量提取指定文件夹下同类型文件
python实现批量提取指定文件夹下同类型文件讲解略代码 import os import shutildef take_samefile(or_path, tar_path, tar_type ...
Excel 从多个文件提取指定的多行（固定行）进行汇总的工具
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档文章目录前言一.工具位置二.使用步骤总结前言 Excel 如果从多个文件提取指定的行汇总到一个文件?在办公过程中我们往往需要 ...
python提取文本中的手机号_Python从vcf文件中读取手机号并进行去重操作
文章目录 1. Python代码 file = open('test.vcf', 'r', encoding='utf-8') tels = [] for line in file: line = l ...

使用bcftools提取指定样本的vcf文件（extract specified samples in vcf format）

使用bcftools提取指定样本的vcf文件（extract specified samples in vcf format）相关推荐

最新文章

热门文章