linux 对同一类型文件批量循环处理,两种方式

本文详细描述了如何利用for和while循环结构在Linuxshell脚本中遍历目录,对以.id结尾的文件进行内容提取(使用seqkitgrep),并将结果转换为长度信息(seqkitfx2tab)。
摘要由CSDN通过智能技术生成

for循环

RefDir=/mnt/sda2/genome
for i in ./*id; do #这是一个循环结构,遍历当前目录下以".id"结尾的文件。每次循环时,将文件的路径存储在变量"i"中
    base_name=$(basename $i ".id") #将文件名中的".id"后缀去除,得到文件的基本名称,并将其存储在变量"base_name"中
    seqkit grep -f $i $RefDir/genome_pep.fa > ${base_name}_pep.fa #提取匹配"i"文件中的内容,并将结果保存到"${base_name}_pep.fa"文件中。
    seqkit fx2tab --length --name ${base_name}_pep.fa > ${base_name}_pep.len #提取fasta序列的长度并储存在{base_name}_pep.len
done #循环结束

while循环

RefDir=/mnt/sda2/genome
ls ./*id | while read i #这是一个循环结构,遍历当前目录下以".id"结尾的文件。每次循环时,将文件的路径存储在变量"i"中
	do
    base_name=$(basename $i ".id") #将文件名中的".id"后缀去除,得到文件的基本名称,并将其存储在变量"base_name"中
    seqkit grep -f $i $RefDir/genome_pep.fa > ${base_name}_pep.fa #提取匹配"i"文件中的内容,并将结果保存到"${base_name}_pep.fa"文件中。
    seqkit fx2tab --length --name ${base_name}_pep.fa > ${base_name}_pep.len #提取fasta序列的长度并储存在{base_name}_pep.len
done #循环结束
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值