prodigal+cd-hit建立基因集

最新推荐文章于 2021-05-16 16:42:05 发布

马志远的生信笔记

最新推荐文章于 2021-05-16 16:42:05 发布

阅读量474

点赞数 1

分类专栏：序列数据处理

本文链接：https://blog.csdn.net/weixin_42480153/article/details/109077255

版权

序列数据处理专栏收录该内容

35 篇文章 8 订阅

订阅专栏

cd /mnt/mzy/dairycow/24samples/mag/genset/raw
for i in *fa
do
prodigal -i $i -a /mnt/mzy/dairycow/24samples/mag/genset/$i.faatmp -d /mnt/mzy/dairycow/24samples/mag/genset/$i.fatmp;
done
cd /mnt/mzy/dairycow/24samples/mag/genset/
cat *.fatmp >all.fa
cat *.faatmp>all.faa && rm *tmp
cd-hit-est -c 0.95 -G 0 -M 0 -d 0 -aS 0.9 -r 1 -T 32 -n 9 -i all.fa -o unique.fa