我有这个小脚本名为my_script.sh
#!/bin/bash
size=$(infoseq -auto -nocolumns -delimiter ',' -only -noheading -length lambda_8.fa)
size=$(( $size - 200 ))
pos=$(($RANDOM % $size))
size2=$(( $pos + 200 ))
extractseq lambda_8.fa -regions $pos-$size2
我想要的是能够将它与不同的文件一起使用,并根据输入文件将输出保存在另一个文件中,这样我就可以像这样使用它:
bash my_script.sh file.fa
输出可能如下所示200_pb_file.fa
。