小编Man*_*uel的帖子

对多个文件运行命令

我有一个文件夹,里面有很多文件(xyz1、xyz2,一直到 xyz5025),我需要在每个文件上运行一个脚本,得到 xyz1.faa、xyz2.faa 等等作为输出。

单个文件的命令是:

./transeq xyz1 xyz1.faa -table 11
Run Code Online (Sandbox Code Playgroud)

有没有办法自动做到这一点?也许是一个for-do组合?

scripting shell-script text-processing files bioinformatics

20
推荐指数
3
解决办法
5870
查看次数

将命令应用于目录中的所有文件

我想将以下命令应用于当前目录中的所有文件

clustalw -align -infile=*here goes the input filename* -outfile=*here goes the output filename*
Run Code Online (Sandbox Code Playgroud)

我还希望输出文件名与输入加上“.aln”相同。例如:

clustalw -align -infile=nexus0000.fna -outfile=nexus000.fna.aln
Run Code Online (Sandbox Code Playgroud)

谢谢您的帮助!

scripting shell-script files

2
推荐指数
1
解决办法
6580
查看次数

提取包含模式的行

我想提取包含以下模式的文件中的所有行:“#1:”和“tree length for”。

输入:

#1: nexus0002_Pseudomonas_10M     

 branch          t       N       S   dN/dS      dN      dS  N*dN  S*dS

   6..5      0.000   390.0   195.0  0.0668  0.0000  0.0000   0.0   0.0
   6..7      0.013   390.0   195.0  0.0668  0.0008  0.0114   0.3   2.2
   7..1      0.000   390.0   195.0  0.0668  0.0000  0.0000   0.0   0.0
   7..4      0.000   390.0   195.0  0.0668  0.0000  0.0000   0.0   0.0
   6..8      0.000   390.0   195.0  0.0668  0.0000  0.0000   0.0   0.0
   8..2      0.013   390.0   195.0  0.0668  0.0008  0.0114   0.3   2.2
   8..3      0.013   390.0   195.0  0.0668  0.0008  0.0114   0.3   2.2

tree length …
Run Code Online (Sandbox Code Playgroud)

sed text-processing bioinformatics

2
推荐指数
1
解决办法
8729
查看次数