Entering edit mode
16 months ago
Riyad
•
0
Hi, I am new at bioinforamtic world. I am looking for homology of my target sequences through blastp with Swissport database. My subjects and querry bothfiles are large and expecting a large number of hits as well. Currently, I am running the code-
blastp -query pep.fasta -db /data/swissprot -out blstp.txt -evalue 1e-10 -outfmt "6 qseqid evalue pident ssciname sskingdom stitle" -max_target_seqs 1
I want to have high confidence hits/match to select as target seqs.
What is the best strategy to do it?