Linux下根据ID文件筛选另一个文件的列

最新推荐文章于 2023-12-06 10:12:56 发布
最新推荐文章于 2023-12-06 10:12:56 发布
阅读量1.2k 收藏 1
点赞数
文章标签： linux python shell
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_46605479/article/details/112511181
版权

linux下筛选需要的ID列，这个方法仅适用于在vcf文件中。
vcf文件：
在这里插入图片描述
前8列是固定的，后面的列是个体ID，现在我要根据需要的ID筛选这个vcf文件。需要的ID放到yes.txt文件里
在这里插入图片描述
使用的命令为：

vcftools --vcf in.vcf --recode --recode-INFO-all --stdout  --keep id.txt   > out.vcf

如果是要去掉yes.txt文件ID列，那么就把keep改成remove。

vcftools --vcf in.vcf --recode --recode-INFO-all --stdout  --remove id.txt   > out.vcf

参考链接：https://www.omicsclass.com/article/647

小八四爱吃甜食
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Linux/Shell中根据一列信息批量提取文件中匹配的特定行
Timo_CSDN的博客
03-28 1万+
1. 非压缩文件提取 看过很多代码，我觉得这个最好用 awk 'NR==FNR{a[$1]++}NR>FNR&&($2 in a)' fileA.txt <(cat fileB) > ../out.txt fileA.txt是一列你用来匹配的信息，fileB是你要从中提取数据的文件，导出到out.txt 之前也尝试过用 #!/bin/bash cat fileA.txt | while read id; do real_id=`echo $id| cut -d
linux tail 查看日志文件方法
09-15
Linux中的`tail`命令是一个非常实用的工具，主要用于查看文件的末尾内容，尤其是在系统日志监控、程序调试和故障排查中发挥着重要作用。通过这个命令，你可以实时追踪日志文件的变化，无需手动刷新，从而高效地获取...
【linux基础】指定分隔符分割文件 cut 、awk(指定多个分隔符分割文件提取列，如提取注释文件gff3内的基因ID)
jxm的博客
04-19 1681
如：你想提取QTL区间的基因，那你得先准备一个参考基因组上基因ID 、基因Name、起始、终止的文件。2. 思路：先提取第三列是基因---再提取有Name 的基因--再分隔。匹配也可以借助 awk '/Name/' file.txt。示例：如果你想从GFF3里提取出基因怎么做？2. awk ## 同时指定多个分割符。如果指定多个分隔符，使用中括号。
根据一列数据 抽取另一个文件的行数据
weixin_34015860的博客
07-29 2143
方法一： 使用awk处理，先读入文件中的一列数据，然后在第二文件中做判断 awk -F'\t' 'FILENAME=="commUsers_Hotel"{F[$0]=1}FILENAME=="Hotel3"{if($0 in F){print}}' commUsers_Hotel Hotel3 &gt; fugai  ...
python dict批量选择_根据ID从FASTA文件中批量提取序列【Python脚本】
weixin_39640573的博客
11-27 467
博主是一个刚刚接触生信的新手，正在学习Linux和Python，偶尔会在该博客上面发布自己练习编程写的脚本，用来记录自己的学习之路。介绍根据序列的ID号从FASTA文件中批量提取序列是在平时常常要做的工作，Linux当中grep和awk工具、Perl语言和Python语言，以及samtools等都可以实现，以下是博主用Python实现的从FASTA文件中批量提取序列的脚本。说明需要用到fasta...
使用linux命令 根据序列ID快速提取fastq序列
dujidan的博客
04-02 6643
做生信的人经常会接触到一些fastq数据，有时想要从fastq文件中提取出某些序列来查看。 一般情况可以使用grep命令，就可以实现 grep -A 3 seq_id fastq.file # 匹配seq_id 并向下去3行 但是，如果需要的序列很多，在使用grep就会很慢了，所以这里给出了 awk 的命令，速度简直快的飞起。 awk -F ' ' 'NR==FNR {a[$1]=1; next} { if (a[substr($1,2)]) {print $0; getline b; print b.
linux 文件夹属性 id,Linuxchmod命令修改文件与文件夹权限方法
weixin_42297714的博客
05-11 162
在linux中要修改一个文件夹或文件的权限我们需要用到linux chmod命令来做，接下来是小编为大家收集的Linuxchmod命令修改文件与文件夹权限方法，欢迎大家阅读：Linuxchmod命令修改文件与文件夹权限方法语法：chmod [who] [+ | - | =] [mode] 文件名命令中各选项的含义为u 表示“用户(user)”，即文件或目录的所有者。g 表示“同组(group)用户...
Linux下重启多个 tomcat 服务的脚本(推荐)
09-15
在实际操作中，我们通常需要获取进程ID，这可以通过`awk`命令实现，它可以根据特定列（如第二列，即PID）进行处理。 在脚本中，我们首先创建一个名为`restart.sh`的文件，并设置其为可执行。在文件中，我们定义了...
常用查看日志的linux命令
01-20
在Linux系统中，日志文件是诊断和解决各种问题的重要工具。...通过以上方法，您可以根据需求灵活地在Linux日志文件中查找、筛选和分析信息。在日常的系统管理和故障排查中，熟练运用`grep`命令将极大地提高工作效率。
Linux下Qt实现任务管理器
03-12
在Linux环境下，使用Qt框架开发任务管理器是一个有趣且实用的项目。Qt是一个跨平台的C++图形用户界面应用程序开发框架，广泛应用于桌面、移动和嵌入式系统。结合Linux的强大功能，我们可以创建一个功能丰富的任务...
linux下查看系统进程占用的句柄数方法
09-15
因此，建议根据实际需求谨慎设置这个值，通常不超过4096是一个合理的上限。 总的来说，理解如何查看和管理Linux进程的句柄数对于优化系统性能、排查问题以及确保服务稳定运行至关重要。通过熟练掌握`ulimit`、`lsof...
Linux命令awk，逐行按列提取文件【超级简单】
SunJHe的博客
03-13 1万+
例：cat 文件名 | awk '{print "*"$2"\n"$5}' &gt; 路径/新文件名 &amp;解：cat：打开文件命令       文件名：打开的文件名字       |：管道      awk：awk文件分析命令      ‘{}’：文件操作      print：在新文件中输出      “*”：在新文件中输出的内容，*可以是任何内容      $2：在新文件中输出，旧文件的...
linux两个字符串拼接输出,linux join命令参数及用法详解---连接两个文件中的指定内容...
weixin_39620578的博客
04-28 1117
功能说明：将两个文件中，指定栏位内容相同的行连接起来。语　　法：~~~~~~~~~~~~~~~~~~~经常使用这个命令。但老是搞忘一些参数，特此记录join类似 db里面的join方法，同样有left join right join inner join等指定参数-a 可以指定join的方式。 -a1表示 显示第一个文件中不匹配的行，即为left join-a2 表示 显示第二个文件中不匹配的行，...
linux 筛选文件,「Linux基础知识」grep文件内容筛选命令的使用
weixin_32534073的博客
04-29 1094
原标题：「Linux基础知识」grep文件内容筛选命令的使用grep命令用于从文档中抓取显示包含指定字符的行，grep命令的使用格式如下：grep [选项] 匹配模式 文件1 文件2 ......grep常见的选项有：-n 显示匹配文档行的行号-i 忽略大小写按照匹配模式匹配文档的内容-v 显示复合匹配模式之外的文档行的内容(反选) grep命令的使用方式举例如下：①grep "name" f1：...
VET：基因变异VCF数据集便捷提取工具
青笋的博客
08-18 663
VET是一个基于R语言开发的变异位点信息批量提取工具，主要功能是根据VCF数据集，按照基因ID、样品ID、变异位点ID等参数，实现批量提取，同时支持变异位点结构注释，一步搞定变异数据的快速提取。，则执行下面的步骤，主要调用Python程序进行信息检索，并由bcftools工具批量提取变异信息，若需要根据指定样品进行过滤，则利用view功能对样品进行筛选，最后生成结果压缩文件。支持流式操作，提取后筛选指定样品并对每个变异位点进行结构注释（判断错义突变、移码突变等），最终将结果文件打包生成压缩包。
Linux如何查找指定内容所在的行?
最新发布
oldboyedu1的博客
12-06 2565
其中，pattern是要查找的字符串，action是匹配到pattern后要执行的操作，如print $0、print $1等。其中，pattern是要查找的字符串，files为要查找的文件。其中，pattern是要查找的字符串，files为要查找的文件。在Linux中，查找指定内容所在的行是非常重要的存在，其使用频繁非常高，因为在多种场景下都会运用到它，比如日志分析、系统维护等，那么Linux如何查找指定内容所在的行?以下是常用的方法介绍。grep命令可以在文件中查找指定的字符串，并输出包含该字符串行。
Linux 中找出两个文件的不同
不辣红烧肉的博客
11-15 1万+
今天遇到一个问题就是两个大文件，其中一个文件是另一个文件的子集，需要找出不同的部分，上网查询后得到如下方法： 方法一：comm命令实现 首先看看comm命令的帮助 comm --help Usage: comm [OPTION]... FILE1 FILE2 Compare sorted files FILE1 and FILE2 line by line. -1
字节跳动面试题汇总 -- C++后端（含答案）
热门推荐
qq_40989769的博客
09-05 1万+
malloc和new的区别排序算法时间复杂度 如何用Linux shell命令统计一个文本中各个单词的个数more log.txt | tr ' ' '\n' | sort | uniq -cLinux下需要打开或者查看大文件linux 怎么查看进程，怎么结束进程？原理是什么？Linux怎么查看当前的负载情况Http CodeMySQL的底层索引结构，InnoDB里面的B+Tree不同引擎索引的区别i++是否原子操作锁的底层实现B Tree 和 B+ Tree的区别b树和b+树MySQL索引的发展过程？是一
vcftools的使用小结
小八四爱吃甜食的博客
04-08 4803
vcftools的使用小结计算等位基因频率计算Fst转换为plink格式比较两个vcf文件根据id保留或去除vcf文件的样本基因型数据转换为012格式参考链接 计算等位基因频率 用--freq计算等位基因频率 vcftools --vcf test.vcf --freq --out output 对来自chr1的每一个位点统计其基因频率 vcftools --gzvcf combined200.vcf.gz --freq --chr chr1 --out chr1_analysis 计算Fst 使用--
Linux综合实践试卷生成及评阅程序是设计与分析详细设计说明用文件实现
06-04
Linux综合实践试卷生成及评阅程序的详细设计说明包括以下几个方面： 1. 数据库设计 试题库和答卷库是程序的核心数据库，需要设计合理的数据库结构和表关系，确保数据的完整性和一致性。试题库包括试题ID、题目、难度系数、类型和答案等信息，答卷库包括学生ID、学生姓名、考试时间、试卷ID和答案等信息。可以使用MySQL数据库来实现。 2. 试卷生成算法 试卷生成算法是根据试题库的信息，以一定的难度系数和数量，自动组合试题，生成一份试卷。可以使用Python编程语言来实现。具体实现过程如下： （1）从试题库中根据难度系数和类型筛选出若干试题。 （2）根据试卷要求的试题数量和各类型试题的比例，计算每个类型所需的试题数量。 （3）从筛选出的试题中，按照类型分别选取所需数量的试题。 （4）将选取的试题按照一定的顺序组合成试卷。 （5）将试卷保存到文件或数据库中。 3. 评分算法 评分算法是根据试题的难度系数和答案，以及学生的答案，自动评估学生的答卷，计算学生的得分。可以使用Python编程语言来实现。具体实现过程如下： （1）从答卷库中根据试卷ID和学生ID获取学生的答卷信息。 （2）从试题库中根据试题ID获取每个试题的难度系数和答案。 （3）将学生的答案和试题的答案进行比对，计算每个试题的得分。 （4）根据试题的难度系数和得分，计算学生的总得分。 （5）将学生的得分保存到答卷库中。 4. 文件实现 试题库、答卷库和试卷可以使用文件来实现。可以使用CSV或JSON格式的文件来存储数据。具体实现过程如下： （1）试题库和答卷库：使用CSV或JSON格式的文件来存储数据，每行代表一个试题或答卷，每列代表一个数据项。 （2）试卷：使用文本文件来存储试卷，每行代表一道试题，按照一定的顺序组合成试卷。 文件实现的优点是简单易用，适合小规模数据存储和读取。缺点是不适合大规模数据存储和复杂查询，需要额外的查询和索引机制来提高效率。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值