如何按一个列表对基因型文件进行过滤，剔除不需要的样本？ - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

1430: 浏览

如何按一个列表对基因型文件进行过滤，剔除不需要的样本？

现在有一个vcf基因型文件，为了降低后续群体结构计算的内存大小，想要删除没有表型数据的样本列，在linux下如何操作？

好问题 0 评论收藏举报

查看全部 1 个回答

chatGPT机器人机器人用户来自于: 北京市
2023-05-23 11:35

可以使用bcftools工具对vcf基因型文件进行过滤，剔除不需要的样本。以下是具体的操作步骤：

安装bcftools工具，可以使用以下命令进行安装：

sudo apt-get install bcftools

使用bcftools进行过滤，命令如下：

bcftools view -S <样本列表> -Oz -o <输出文件>

其中，-S参数后面跟上样本列表文件的路径，样本列表文件中每行一个样本名；为需要过滤的vcf基因型文件的路径；-Oz和-o参数是输出参数，分别表示输出文件的格式和路径。

例如，如果需要过滤掉样本列表文件中包含的样本，可以使用以下命令：

bcftools view -S samples.txt input.vcf -Oz -o output.vcf.gz

这样就可以生成一个过滤后的vcf基因型文件output.vcf.gz。

赞同 1 0评论

问题动态

发布时间: 2023-05-23 11:15

更新时间: 2023-05-23 15:37

关注人数: 2 人关注

相关问题

如何判断GO号的层级？: 1108 浏览 2 关注 1 回答 0 评论

【求助】如何确定一种基因ID的类型: 1989 浏览 4 关注 2 回答 0 评论

ciriRNA表达量如何计算？（更）: 1088 浏览 2 关注 1 回答 0 评论

在仅有每个样本(有3个repeat)所有转录本的fpkm的情况下，如何利用cuffdiff进行差异基因分析。: 1806 浏览 4 关注 3 回答 0 评论

如何下载指定文献的原始数据？？: 1235 浏览 2 关注 1 回答 0 评论

Ubuntu 14.04 server，如何进行网络配置呢？: 1196 浏览 2 关注 1 回答 0 评论

如何计算碱基偏好性？: 1189 浏览 3 关注 1 回答 0 评论

在N-W算法中，为什么忽略了从x状态道y状态的转化（xi对应到一个空位到yi对应一个空位的转化）: 1413 浏览 2 关注 2 回答 0 评论

生信到底如何入门？？？: 2345 浏览 6 关注 6 回答 0 评论

先埋一个帖子，希望一年之后大家都能收获到自己想要的。: 2870 浏览 11 关注 10 回答 0 评论

推荐内容

问 computeMatrix画图问题: 1809 浏览 1 关注 2 回答 0 评论

问做进化树时，有些分支没有bootstrap值，正常吗？为什么？: 1731 浏览 2 关注 1 回答 0 评论

问请问下这种格式的R语言内容如何选择最小值: 1368 浏览 2 关注 1 回答 0 评论

问 hmmsearch和hmmscan: 1700 浏览 2 关注 1 回答 0 评论

问 juicer_tools.jar hiccups 运行怎么更改物种？: 1968 浏览 2 关注 2 回答 0 评论

问启动子染色质开放性: 1704 浏览 2 关注 2 回答 0 评论

问全基因组关联分析结果与模型选择: 1154 浏览 2 关注 1 回答 0 评论

问如何提取比对上某一个基因或者某一段序列的所有reads: 1082 浏览 2 关注 1 回答 0 评论

问 atac重复样品可视化: 1357 浏览 2 关注 1 回答 0 评论

问使用tophat2和bowtie1寻找环形RNA时报错: 1738 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025