Given a protein sequence query, the script performs a blast search against a nucleotide subject, filters hits to keep only hits with >30% sequence identity and >90% match length (i.e., 90% of the query sequence length) and prints the number of matches identified.
conda create -n new_env
conda activate new_env
conda install -c bioconda blast
./find_homologs.sh $QUERY_SEQ $SUBJECT_SEQ $OUTPUT