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

登录

该问题已被锁定！

2: 关注

3419: 浏览

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

现在有一个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 人关注

相关问题

二代测序得到测序数据组装完成后如何进行丰度计算: 2991 浏览 1 关注 1 回答 0 评论

如何判断一篇CHIPseq的文章用的是hg19还是hg38: 3476 浏览 5 关注 3 回答 0 评论

stringtie 得到的gtf 通过DESeq2分析后stringtieID 如何转换成esmbleID: 2721 浏览 3 关注 2 回答 0 评论

请问一下，如何根据基因名批量下载CDS区的序列？: 2843 浏览 2 关注 1 回答 0 评论

如何知道自己测序的adapter序列？: 3137 浏览 3 关注 1 回答 0 评论

linux 中less -S 如何查看过长被遮盖的内容: 2720 浏览 2 关注 1 回答 0 评论

两次运行同一个FQ文件产生的bam文件分析结果存在差异: 3002 浏览 2 关注 2 回答 0 评论

HHblits使用相同的数据库，得出的结果明显差于HHpred网页版的结果。该如何使用HHsuite命令行得到与网页端类似的效果: 4381 浏览 2 关注 2 回答 0 评论

【求助】如何确定时期特定表达基因: 2364 浏览 2 关注 1 回答 0 评论

x <- x[keep.exprs,, keep.lib.sizes=FALSE] 请问一下这条命令该如何解读呢？: 1605 浏览 1 关注 0 回答 2 评论

推荐内容

问染色体号是罗马数字怎么写sh循环: 2163 浏览 1 关注 1 回答 0 评论

问 3D基因组里compartment里一般是包含好多TAD的，但这图为什么compartment数量比TAD多这么多呢: 3161 浏览 2 关注 2 回答 0 评论

问请问cytoscape是怎么构建网络的？: 2225 浏览 3 关注 1 回答 0 评论

问如何对特征数量少的空间蛋白组数据进行细胞聚类？: 3070 浏览 2 关注 1 回答 0 评论

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

问 ASR祖先序列重建，最后一步使用PAML时出现一些问题: 5373 浏览 1 关注 3 回答 0 评论

问 kraken2软件运行时内存分配的问题: 5112 浏览 1 关注 2 回答 0 评论

问使用bedtools提取序列时的deyond length问题；bed负值: 3810 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2026