qq_41856194的博客

我用的python3.6版本，网络上都是从ftp直接下载psiblast的压缩包，然后解压使用，再配置环境变量，和再anaconda中安装是一样的。

大文件：31万个序列（我看的文章里面说一小时处理100条序列，这里我用了）代码如下：把mulfa.fasta文件拆分成单个fasta文件awk '{if($0~/^>/)a=$0; else{system("echo \\"a"\"\n\""$0a)}}' mulfa.fastamulfa.fasta文件如下：生成的单个fasta文件如下：小文件处理及代码解释如这篇文章所说：http://t.csdn.cn/xTZId感谢这位大佬......