shell脚本提取txt文件中关键词并统计出现次数

有一个test.txt文件，里面想提取特定关键词，并对各项关键词统计出现过几次，使用shell脚本。

关键词的特征是：**SRA.关键词.**，开头“SRA.”+关键词+结尾“.”

SRA:SRR10168379.12205864.1关键词是：SRR10168379
SRA:SRR10168392.8392060.2关键词是：SRR10168392

希望运行完成的结果res.txt，关键词——号后面加上次数，按照次数多少排序：
SRR10168379——7次
SRR10168392——2次

test.txt文件内容是

# tblastn
# Iteration: 0
# Query: 
# RID: ZRD35BAF013
# Database: SRR10168375 SRR10168376 SRR10168377 SRR10168378 SRR10168379 SRR10168381 SRR10168392 SRR10168393 SRR13285085 SRR13285570
# Fields: query acc.ver, subject acc.ver, % identity, alignment length, mismatches, gap opens, q. start, q. end, s. start, s. end, evalue, bit score, % positives, query/sbjct frames
# 100 hits found
Query_45991    SRA:SRR10168379.12205864.1    34.694    49    32    0    527    575    148    2    4.6    36.6    53.06    0/-3
Query_45991    SRA:SRR10168379.10841544.1    41.667    48    20    1    187    226    144    1    5.7    36.2    58.33    0/-1
Query_45991    SRA:SRR10168392.8392060.2    47.059    34    13    1    194    222    27    128    11    35.4    61.76    0/3
Query_45991    SRA:SRR10168393.9810230.1    41.304    46    19    1    187    224    1    138    15    35.0    58.70    0/1
Query_45991    SRA:SRR10168379.2460949.2    41.304    46    19    1    187    224    1    138    18    34.7    58.70    0/1
Query_45991    SRA:SRR10168393.20965295.2    42.222    45    18    1    188    224    1    135    20    34.7    57.78    0/1
Query_45991    SRA:SRR10168376.8708660.2    43.902    41    15    1    192    224    1    123    28    34.3    58.54    0/1
Query_45991    SRA:SRR10168379.12533534.1    40.000    50    22    1    187    228    150    1    31    34.3    56.00    0/-1
Query_45991    SRA:SRR10168379.6639135.2    41.304    46    19    1    187    224    141    4    34    33.9    58.70    0/-1
Query_45991    SRA:SRR10168379.13010027.2    39.583    48    21    1    187    226    1    144    41    33.9    58.33    0/1
Query_45991    SRA:SRR10168381.1806861.1    40.816    49    21    1    188    228    150    4    41    33.9    55.10    0/-1
Query_45991    SRA:SRR10168379.3721520.2    40.816    49    21    1    188    228    150    4    41    33.9    55.10    0/-1
Query_45991    SRA:SRR10168378.17299083.1    41.304    46    19    1    187    224    139    2    42    33.9    56.52    0/-3
Query_45991    SRA:SRR10168393.17810994.2    39.583    48    21    1    187    226    3    146    46    33.5    58.33    0/3
Query_45991    SRA:SRR10168379.2656880.1    41.304    46    19    1    187    224    144    7    53    33.5    58.70    0/-1
Query_45991    SRA:SRR10168379.1997738.2    41.304    46    19    1    187    224    146    9    53    33.5    58.70    0/-2
Query_45991    SRA:SRR10168379.11604415.1    41.304    46    19    1    187    224    149    12    55    33.5    56.52    0/-2
Query_45991    SRA:SRR10168379.11899618.1    39.583    48    21    1    187    226    147    4    57    33.5    56.25    0/-1
Query_45991    SRA:SRR10168379.4610022.2    39.583    48    21    1    187    226    147    4    57    33.5    56.25    0/-1

写回答
好问题 0 提建议
追加酬金
关注问题
分享
邀请回答
编辑收藏删除
收藏举报

4条回答默认最新

_GX_ 2022-02-03 22:48

关注

$ cat test.awk
!/^[ \t]*#/ {
    keyword = substr($2, 5, 11);
    count[keyword]++;
} END {
    for (keyword in count)
        print keyword "---" count[keyword];
}
$ cat test.txt
# tblastn
# Iteration: 0
# Query: 
# RID: ZRD35BAF013
# Database: SRR10168375 SRR10168376 SRR10168377 SRR10168378 SRR10168379 SRR10168381 SRR10168392 SRR10168393 SRR13285085 SRR13285570
# Fields: query acc.ver, subject acc.ver, % identity, alignment length, mismatches, gap opens, q. start, q. end, s. start, s. end, evalue, bit score, % positives, query/sbjct frames
# 100 hits found
Query_45991    SRA:SRR10168379.12205864.1    34.694    49    32    0    527    575    148    2    4.6    36.6    53.06    0/-3
Query_45991    SRA:SRR10168379.10841544.1    41.667    48    20    1    187    226    144    1    5.7    36.2    58.33    0/-1
Query_45991    SRA:SRR10168392.8392060.2    47.059    34    13    1    194    222    27    128    11    35.4    61.76    0/3
Query_45991    SRA:SRR10168393.9810230.1    41.304    46    19    1    187    224    1    138    15    35.0    58.70    0/1
Query_45991    SRA:SRR10168379.2460949.2    41.304    46    19    1    187    224    1    138    18    34.7    58.70    0/1
Query_45991    SRA:SRR10168393.20965295.2    42.222    45    18    1    188    224    1    135    20    34.7    57.78    0/1
Query_45991    SRA:SRR10168376.8708660.2    43.902    41    15    1    192    224    1    123    28    34.3    58.54    0/1
Query_45991    SRA:SRR10168379.12533534.1    40.000    50    22    1    187    228    150    1    31    34.3    56.00    0/-1
Query_45991    SRA:SRR10168379.6639135.2    41.304    46    19    1    187    224    141    4    34    33.9    58.70    0/-1
Query_45991    SRA:SRR10168379.13010027.2    39.583    48    21    1    187    226    1    144    41    33.9    58.33    0/1
Query_45991    SRA:SRR10168381.1806861.1    40.816    49    21    1    188    228    150    4    41    33.9    55.10    0/-1
Query_45991    SRA:SRR10168379.3721520.2    40.816    49    21    1    188    228    150    4    41    33.9    55.10    0/-1
Query_45991    SRA:SRR10168378.17299083.1    41.304    46    19    1    187    224    139    2    42    33.9    56.52    0/-3
Query_45991    SRA:SRR10168393.17810994.2    39.583    48    21    1    187    226    3    146    46    33.5    58.33    0/3
Query_45991    SRA:SRR10168379.2656880.1    41.304    46    19    1    187    224    144    7    53    33.5    58.70    0/-1
Query_45991    SRA:SRR10168379.1997738.2    41.304    46    19    1    187    224    146    9    53    33.5    58.70    0/-2
Query_45991    SRA:SRR10168379.11604415.1    41.304    46    19    1    187    224    149    12    55    33.5    56.52    0/-2
Query_45991    SRA:SRR10168379.11899618.1    39.583    48    21    1    187    226    147    4    57    33.5    56.25    0/-1
Query_45991    SRA:SRR10168379.4610022.2    39.583    48    21    1    187    226    147    4    57    33.5    56.25    0/-1
$ cat test.txt | awk -f test.awk > res.txt
$ cat res.txt
SRR10168392---1
SRR10168376---1
SRR10168381---1
SRR10168393---3
SRR10168378---1
SRR10168379---12

本回答被题主选为最佳回答 , 对您是否有帮助呢?

查看更多回答(3条)

报告相同问题？

关注问题

shell脚本提取txt文件中关键词并统计出现次数 bash 全文检索有问必答正则表达式
2022-02-03 21:50

回答 4 已采纳 $ cat test.awk !/^[ \t]*#/ { keyword = substr($2, 5, 11); count[keyword]++; } END { for
shell或者python脚本提取xml文件中信息 bash xml 有问必答
2022-03-05 19:19

回答 3 已采纳将test.xml文件放在xml-info.py的同一目录下，运行如下代码即可： import argparse from lxml import etree parser=argparse.Argu
shell脚本提取关键字 linux ssh
2023-02-01 14:21

回答 2 已采纳先用下面的第一句确认一下是不是正确，然后用第二句输出结果 grep -P "(?<= as )[a-z_]+" grep -Po "(?<= as )[a-z_]+"
常用的Shell脚本集合
2022-04-30 18:10

头发莫的了呀的博客 Linux 系统发送告警脚本3. MySQL 数据库备份单循环4. MySQL 数据库备份多循环5. Nginx 访问日志按天切割6. Nginx 访问日志分析脚本7. 查看网卡实时流量脚本8. 服务器系统配置初始化脚本9. 监控 100 台服务器磁盘...
shell脚本调用可执行文件？ bash c++ linux
2021-07-16 13:35

回答 2 已采纳 main中用cout<<var;
linux的shell脚本这么读取另一个文件的内容 bash linux
2022-03-08 14:03

回答 1 已采纳 “脚本名文件名”执行即可 tre=`cat $1` for i in $tre; do echo "你好 ${i}！" done
请教shell脚本进行统计目录下的文件数量无法实现问题 linux
2021-11-17 18:20

回答 2 已采纳把命令用$() 包起来你这样赋值肯定不行 eg: w=$(ls)
shell脚本100例
2020-11-18 10:40

超人说的博客判断文件或目录是否存在 #!/bin/bash if [ $# -eq 0 ] ;then echo “未输入任何参数,请输入参数” echo “用法:$0 [文件名|目录名]” fi if [ -f $1 ];then echo “该文件,存在” ls -l $1 else echo “没有该文件...
shell或python脚本批量处理txt文件 bash linux python 有问必答
2022-03-20 04:01

回答 6 已采纳 # -*- utf8 -*- """ pyCode : test.py Created by: Taro Description: 2022/3/20 0020 """ import os impor
Linux中shell脚本怎样实现遍历某个目录，并实现对该目录下文件与目录的分类统计 linux
2017-12-05 11:50

回答 5 已采纳 ![图片说明](https://img-ask.csdn.net/upload/201712/06/1512550787_714026.png)
shell脚本中执行另一个脚本其他
2021-03-22 16:28

回答 1 已采纳 nohup 英文全称 no hang up（不挂起），用于在系统后台不挂断地运行命令，退出终端不会影响程序的运行。 nohup 命令，在默认情况下（非重定向时），会输出一个名叫 nohup.out
统计文件中出现次数最多的前10个单词
2016-08-18 19:17

会飞行的小蜗牛的博客使用linux命令或者shell实现：文件words存放英文单词，格式为每行一个英文单词（单词可以重复），统计这个文件中出现次数最多的前10个单词。 cat words.txt | sort | uniq -c | sort -k1,1nr | head -10 　主要...
shell脚本里的命令加参数参数是从另一个文件中取得怎么设置 linux
2015-05-19 05:43

回答 8 已采纳用awk 可以很简单 #!/bin/bash cat /xxx/arg_conf | while read LINE do eval $(echo $LINE | awk '{prin
Linux脚本能100以内逢三过,shell脚本100例
2021-05-14 02:57

狗哥老司机的博客 /bin/bashweb使用脚本循环建立三位数字的文本文件(111-999 的文件)for i in {1…9}dofor j in {1…9}dofor k in {1…9}dotouch /tmp/iiij$k.txtdonedonedone3二、统计/etc/passwd 中 root 出现的次数bash#!...
100个实用的 Linux Shell 脚本经典案例，收藏备用~
2022-07-15 09:40

Java笔记虾的博客点击关注公众号，利用碎片时间学习1）编写hello world脚本#!/bin/bash #编写helloworld脚本 echo"HelloWorld!"2）通过位置变量创建 Linux 系统账户及密码#!/bin/bash ...#$1是执行脚本的第一个参数,$2是执行脚本...
没有解决我的问题, 去提问

问题事件

关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
系统已结题 2月11日
关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
已采纳回答 2月3日
关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
已采纳回答 2月3日
关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
创建了问题 2月3日

悬赏问题

¥15 metadata提取的PDF元数据，如何转换为一个Excel
¥15 关于arduino编程toCharArray()函数的使用
¥100 vc++混合CEF采用CLR方式编译报错
¥15 coze 的插件输入飞书多维表格 app_token 后一直显示错误，如何解决？
¥15 vite+vue3+plyr播放本地public文件夹下视频无法加载
¥15 c#逐行读取txt文本，但是每一行里面数据之间空格数量不同
¥50 如何openEuler 22.03上安装配置drbd
¥20 ING91680C BLE5.3 芯片怎么实现串口收发数据
¥15 无线连接树莓派，无法执行update，如何解决？（相关搜索：软件下载）
¥15 Windows11, backspace, enter, space键失灵

shell脚本提取txt文件中关键词并统计出现次数

4条回答 默认 最新

问题事件

悬赏问题

4条回答默认最新