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

登录

该问题已被锁定！

2: 关注

944: 浏览

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

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

相关问题

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

在riboseq分析中，如何对没有起始密码子和终止密码子的转录本进行三碱基准确性分析？: 785 浏览 2 关注 1 回答 0 评论

用scanpy做单细胞分析，如何去除周期细胞的影响: 1082 浏览 2 关注 1 回答 0 评论

如何理解单细胞测序中用fastQC对低质量细胞进行过滤: 1132 浏览 2 关注 1 回答 0 评论

统计学，如何从入门到放弃？: 849 浏览 3 关注 1 回答 0 评论

Affy包中的MAplot函数的作用是什么？画出的MAplot图如何解读？为何要进行MAplot检查？: 1050 浏览 3 关注 1 回答 0 评论

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

菌侵染植物的转录组分析中，若菌的一个基因在当前时间没有引起植物基因的变化，过了三个小时后才引起植物基因的变化，如何分析这种现象或如何将这一类基因关联起来，有没有类似的文献推荐？: 843 浏览 2 关注 1 回答 0 评论

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

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

推荐内容

问 pheatmap画图的数据导入问题: 1095 浏览 3 关注 2 回答 0 评论

问 cellranger使用问题: 1322 浏览 2 关注 1 回答 0 评论

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

问使用Tracking Tumor ImmunoPhenotype(TIP)网站分析TCGA的BLCA_tpm数据: 943 浏览 2 关注 2 回答 0 评论

问 infercnv运行报错: 1097 浏览 2 关注 1 回答 0 评论

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

问基因和染色体，lncRNA，mRNA之间的关系是怎样的？: 542 浏览 1 关注 0 回答 1 评论

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

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024