#!/bin/bash


### cat Dry2014.faa Wet2014.faa > Dry_Wet.faa
date
#mmseqs easy-cluster --threads 32 -c 0.95 --min-seq-id 0.95 --cov-mode 2 Dry_Wet.faa Dry_Wet_cluster tmp 
#### ~1min
date

date
#bwa index Dry_Wet_cluster.fnn_rep_seq.fasta -p Dry_Wet
date

date

#run_dbcan Dry_Wet.cluster.faa protein --tools hmmer --stp_cpu 32 -c Dry_Wet.gff --cgc_substrate --out_dir Dry_Wet.PUL.Sub --dia_cpu 32 --hmm_cpu 32 --tf_cpu 32 --dbcan_thread 32
run_dbcan Dry_Wet.cluster.faa protein --out_dir Dry_Wet.HMM --dia_cpu 32 --hmm_cpu 32 
date