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

登录

该问题已被锁定！

2: 关注

3264: 浏览

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

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

相关问题

基因symbol注释GEO芯片时，如何把mRNA和lncRNA分别标注出来？: 2341 浏览 2 关注 1 回答 0 评论

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

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

你好ChatGPT，我身高1.75m，体重50kg，应该如何减肥: 1946 浏览 1 关注 1 回答 0 评论

请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？: 4486 浏览 2 关注 2 回答 0 评论

如何在20个蛋白序列里找3个motif: 1306 浏览 2 关注 0 回答 0 评论

如何针对affy表达谱芯片计算结果分析不同样本之间是否表达差异: 2060 浏览 2 关注 1 回答 0 评论

请问有谁知道这张图是如何绘制的吗，在NCBI的Conserved Domians模块没找到可以直接导出的地方: 2445 浏览 2 关注 1 回答 0 评论

rMATS是如何得到外显子的两种isform的？: 2402 浏览 2 关注 1 回答 0 评论

请教一个问题，一个转录因子TF只会调控位于同一条链的基因，还是有可能调控反义链的基因？: 2908 浏览 2 关注 2 回答 0 评论

推荐内容

问群体内同源基因的所有变异: 2023 浏览 2 关注 1 回答 0 评论

问去除批次效应选定靶标细胞群后如何进行亚群分析？: 3813 浏览 2 关注 2 回答 0 评论

问 Ancestry_hmm: 1475 浏览 1 关注 0 回答 0 评论

问 prokka数据库更新: 2718 浏览 2 关注 1 回答 0 评论

问 GAPIT包FarmCPU和Blink模型进行GWAS分析报错: 3551 浏览 2 关注 1 回答 0 评论

问 linux下非root用户设置所运行任务的CPU占用率和线程数: 3601 浏览 2 关注 3 回答 0 评论

问 bulk-RNAseq数据集整合: 3182 浏览 2 关注 1 回答 0 评论

问关于GO分析的问题: 2026 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2026