<div dir="ltr"><div><div><div>Hello:<br><br>I want to assign gene function to the predicted genes. I collected UniProt/SwissProt human, mouse and rat protein sequences (including both canonical and isoforms). Then use "makeblastdb" to build the database. And then run "blastp -db ... -max_hsps_per_subject 1" following the example in the protocol. But it returns me an error: Unknown argument: "max_hsps_per_subject".  Why this happens? Is it because I am using a different version of blastp?<br><br>USAGE<br>  blastp [-h] [-help] [-import_search_strategy filename]<br>    [-export_search_strategy filename] [-task task_name] [-db database_name]<br>    [-dbsize num_letters] [-gilist filename] [-seqidlist filename]<br>    [-negative_gilist filename] [-entrez_query entrez_query]<br>    [-db_soft_mask filtering_algorithm] [-db_hard_mask filtering_algorithm]<br>    [-subject subject_input_file] [-subject_loc range] [-query input_file]<br>    [-out output_file] [-evalue evalue] [-word_size int_value]<br>    [-gapopen open_penalty] [-gapextend extend_penalty]<br>    [-qcov_hsp_perc float_value] [-max_hsps int_value]<br>    [-xdrop_ungap float_value] [-xdrop_gap float_value]<br>    [-xdrop_gap_final float_value] [-searchsp int_value]<br>    [-sum_stats bool_value] [-seg SEG_options] [-soft_masking soft_masking]<br>    [-matrix matrix_name] [-threshold float_value] [-culling_limit int_value]<br>    [-best_hit_overhang float_value] [-best_hit_score_edge float_value]<br>    [-window_size int_value] [-lcase_masking] [-query_loc range]<br>    [-parse_deflines] [-outfmt format] [-show_gis]<br>    [-num_descriptions int_value] [-num_alignments int_value]<br>    [-line_length line_length] [-html] [-max_target_seqs num_sequences]<br>    [-num_threads int_value] [-ungapped] [-remote] [-comp_based_stats compo]<br>    [-use_sw_tback] [-version]<br><br></div> Thanks<br></div>Best<br></div>Quanwei<br></div>