1. 将genome_seq.fasta文件使用makeblastdb构建数据库。成功运行后生成文件makeblastdb.ok。 2. 运行程序同目录下的para_blast程序并行化调用tblastn将homolog_proteins.fasta比对到数据库（阈值为evalue<=1e-9且覆盖率>=0.4）。成功运行后生成文件blast.ok和结果文件out.tab。 3. 运行程序同 ...

如果相等在进行计算，看是否组成的串值最大 python import os pos = [] in_m = [] out_m = [] in_pos = [] out_pos = [] cnt = 0 ans = 0 def cal_ans(a,b): if b==10: return a*100+b else: return a*10+b def check(): global in_m ...