编写Bash脚本定位核苷酸序列中T碱基位置的问题求助

原脚本问题分析 #

你的脚本存在三个核心问题：

1. 多余的$@行：会尝试把传入的文件当作命令执行，完全无效，必须删除。
2. grep在脚本（非交互式终端）环境下默认关闭颜色输出，所以看不到标红效果。
3. 原命令仅输出含T的行，没有标记每个T的具体位置，不符合需求。

修正后的脚本方案 #

根据你的需求，提供两种常用实现：

方案1：输出每个T碱基的具体位置（行号+字节偏移） #

#!/bin/bash
# 匹配大小写的T/t，输出行号、字节位置并高亮显示
grep -n -b --color=always -i -o 'T' "$@"

参数说明：

• -n：显示匹配行的行号
• -b：显示T在该行的字节偏移位置（从0开始计数）
• -o：仅输出匹配的T/t字符，而非整行
• --color=always：强制开启颜色高亮，解决脚本中无颜色的问题
• -i：忽略大小写，同时匹配小写t和大写T
• "$@"：安全处理传入的文件名参数，支持多文件输入

方案2：显示含T的整行，高亮T并标注行号 #

如果你需要保留整行内容，同时高亮T并显示行号，用这个版本：

#!/bin/bash
# 显示含T/t的整行，高亮T并标注行号
grep -n --color=always -i 'T' "$@"

运行方式 #

执行脚本时直接传入核苷酸序列文件即可：

bash nucleotide.sh new.dat

内容的提问来源于stack exchange，提问作者Eman