Question: How to Delete Hits from an original sequence
0
gravatar for daniellaplummer
5.2 years ago by
European Union
daniellaplummer0 wrote:

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

 

blast alignment genome • 1.2k views
ADD COMMENTlink modified 5.2 years ago by pld4.8k • written 5.2 years ago by daniellaplummer0
0
gravatar for pld
5.2 years ago by
pld4.8k
United States
pld4.8k wrote:

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).

ADD COMMENTlink modified 5.2 years ago • written 5.2 years ago by pld4.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1730 users visited in the last hour