Blastp command line to get only full matches of query sequence
Hi, I have a downloaded blast database already. Is it possible to blast this database and only get all the sequences in the database that contain EXACTLY this sequence ADARS ?
In other words, if any protein sequence contains ADARS EXACTLY, then I would you in my XML output/, this is my code so far. Thanks.
blastp -outfmt 5 -query /usr/local/projects/heptamer/peptide.fa -word_size 2 -matrix=PAM30 -db /usr/local/projects/heptamer/blastDB/specialBlastDB/BacteriaNR/all_bacteria_nonredundant_protein.faa -out /usr/local/projects/heptamer/blastResults/blastp/Bacteria/BacteriaFTP.xml -evalue 1000000000000000000000000000 -max_target_seqs 5000000
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
BLAST Database error: No alias or index file found for protein database
written by hrza_mar •
-
BLASTP - Searching for homologous sequences of small peptides in several species
written by 瑞宏 •
-
How to run BLAST locally for multiple fasta files
written by Paula •
-
How can I solve the database error?
written by Kenan.gk •
-
Local blast search against nt db gains less results than online search
written by john
-
blastp NCBI C++ exception error: ThrowNullPointerException()
written by garfield320 •
-
BLAST Database creation error
written by sahilbioinf0 •
-
blast command line and online
written by snehalp •
-
Different results with BLAST depending on if subject is formatted database or FASTA-file
written by driliwyr •
-
Different blast results between CLCBio and local blast
written by biobio •
Have you looked at the
qcovsparameter?What about
grep -B 1 ADARS all_bacteria_nonredundant_protein.faor variants thereof ?This assumes unformatted fasta, though.
Yes, the rest is covered in the 'variants' part :)
Thank you all for your help. The qcov worked!