prodigal+cd-hit建立基因集

 

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