Hello,
Using standalone Blast, I used BlastP to align the human proteome against a bacterial proteome, this gave me a number of hits where the human proteome has been found within the bacterial proteome, I need to delete these hits from the original bacterial proteome, in order to gain non-homologous genes, I was wondering if anyone knew a method that would automatically purge the original sequence of the hits instead of doing it manually. The other option is getting blastp to show the non-hits, however I cannot work out the command for this.
Many thanks
Daniella
•
written
by
daniellaplummer •
1 answer
The easiest would be to filter for significant hits and then collect the IDs of query sequences with a significant hit. After you have this list of IDs you could write a simple script to handle this or find one of the many tools online (e.g. http://qiime.org/scripts/filter_fasta.html).
Log in to answer this question.
-
Running blastp with BLAST+ 2.15.0 against custom database; need to identify hits
written by rebecca.calvo •
-
Assigning GO ids to blast hits
written by Mousumi Akter •
-
subtractive Proteomics
written by MISGANA •
-
conducting individual blastp for multiple protein sequences
written by rb77 •
-
When will the blast program stops giving out hits if there are too many (for web server and standal…
written by johnnytam100
-
How to retrieve non homologous protein sequences after blasting - Blastp
written by nombaliphumelele •
-
To find missing proteins at blast output
written by Talip Zengin •
-
blastp in linux
written by joepiwity •
-
Best strategy to get GO terms for a proteome?
written by biotech
-
Low Mge Identities Vs Aclame, The Norm?
written by biotech
how can i filter for significant hits? and how can i collect the ids of query sequences with a significant hits? can you share grep command for it?