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

登录

该问题已被锁定！

2: 关注

840: 浏览

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

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

相关问题

如何获得cas9蛋白的高效低效的sgRNA: 713 浏览 2 关注 1 回答 0 评论

如何在Ubuntu 14.04 server版上安装Rstudio，并在网页上实现绘图功能？: 1046 浏览 4 关注 3 回答 0 评论

如何使用wget命令下载“清华云盘”内的文件: 2289 浏览 2 关注 2 回答 0 评论

chipseq原始文件fastq如何截取DNA长度小于等于260bp进行后续mapping: 852 浏览 2 关注 2 回答 0 评论

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

有台Dell工作站，如何配置成服务器: 726 浏览 4 关注 2 回答 0 评论

如何结合ChatGPT和其他相关软件查询文献: 761 浏览 3 关注 1 回答 0 评论

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

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

seurat中如何通过坐标/细胞名删掉细胞: 1468 浏览 2 关注 1 回答 0 评论

推荐内容

问信号通路特征基因集的查询？用于作ssGSEA分析。: 743 浏览 2 关注 1 回答 0 评论

问 RNA Seq中什么样的reads是chimeric reads？是不是多数的chimeric reads都是circRNA的？: 686 浏览 2 关注 1 回答 0 评论

问 logistic回归、异常值: 769 浏览 2 关注 1 回答 0 评论

问 MCPcounter输入TCGA矩阵的要求？: 715 浏览 2 关注 1 回答 0 评论

问尿代谢组正负离子数据标准化是否可以均用正离子检测出来的肌酐峰: 662 浏览 2 关注 1 回答 0 评论

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

问关于ceRNA网络构建的后续分析有哪些？: 609 浏览 2 关注 1 回答 0 评论

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

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

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024