基因组功能注释uniprot报错

运行下面这行代码时报错Not a valid version 4 database

blastp -query pep.fa -db ${db} \

-num_threads 10  -evalue 0.000001 -seg yes -soft_masking true -lcase_masking -max_hsps 1 -max_target_seqs 1\

 -outfmt "6 qseqid sseqid qlen slen pident length  mismatch gapopen qstart qend sstart send evalue bitscore qcovs stitle"\

  -out uniprot.blast.out.tsv

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

我们已经重新建立数据库,上传了文件,你重新下载数据即可

请先 登录 后评论
  • 1 关注
  • 0 收藏,400 浏览
  • 橙子 提出于 2023-09-22 09:59

相似问题