# 查询序列
querySeq=$1
# 参考数据库
RefSeq=$2.dmnd
# 输出文件
output=$3

./diamond blastp \
--query "$querySeq" \
-d "$RefSeq" \
-o "$output" \
-e 1E-4 \
-p 100 \
-k 10 \
--header simple \
-f 6 qtitle qlen stitle slen length pident mismatch gapopen gaps qstart qend sstart send qcovhsp scovhsp evalue bitscore

# ./diamond makedb --in lucapcycle_extra_p_2_class_v3_case_01_train.fasta --db diamond_db/lucapcycle_extra_p_2_class_v3_case_01_train

# sh run_diamond.sh lucapcycle_extra_p_2_class_v3_case_01_dev.fasta diamond_db/lucapcycle_extra_p_2_class_v3_case_01_train  diamond_against/lucapcycle_extra_p_2_class_v3_case_01_dev_against_lucapcycle_extra_p_2_class_v3_case_01_train.tsv 

# sh run_diamond.sh lucapcycle_extra_p_2_class_v3_case_01_test.fasta diamond_db/lucapcycle_extra_p_2_class_v3_case_01_train  diamond_against/lucapcycle_extra_p_2_class_v3_case_01_test_against_lucapcycle_extra_p_2_class_v3_case_01_train.tsv 


# ./diamond makedb --in lucapcycle_cold_spring_ori_positives.fasta --db diamond_db/lucapcycle_cold_spring_ori_positives

# sh run_diamond.sh lucapcycle_lucapcycle_predicted.fasta diamond_db/lucapcycle_cold_spring_ori_positives  diamond_against/lucapcycle_lucapcycle_predicted_against_lucapcycle_cold_spring_ori_positives.tsv 

# sh run_diamond.sh lucapcycle_cold_spring_ori_positives.fasta diamond_db/lucapcycle_cold_spring_ori_positives  diamond_against/lucapcycle_lucapcycle_cold_spring_ori_positives_against_lucapcycle_cold_spring_ori_positives.tsv 

