Local BLAST search
0
0
Entering edit mode
3.7 years ago

Hello,

I plan to make blast db and find sequences with two genes. So the sequences must have both two genes.(AND AND) Not only the first or only the second.(OR) There is a flag -query in local blast where I'll put genes in fasta format like this:

1 ATGGGGCC.. 2 ATGTTGC..

When I tried to do this I had sequences with the only first gene, the only second gene and the first and the second. How could I filter the variants with both two genes?

A lot of thanks, Valery

blast • 685 views
ADD COMMENT
1
Entering edit mode

You can't do that with blast.

You would need to run 2 separate blast searches, one for each input sequence, and find the commonalities in the outputs.

ADD REPLY
0
Entering edit mode

Seems that it could be so easy) Thank you, I'll stop searching for a magic flag.

ADD REPLY
0
Entering edit mode

First, I do not know if it works because I have not tried it before. The way is that;

  • concatenate the genes - do not write any ids between them
  • run blast using -qcov_hsp_perc parameter, at least over 70%, 80%
  • check query coverage and/or match position in the results If you know the length of gene you have

If it works, let me know :)

ADD REPLY
0
Entering edit mode

Thank you for your reply, but i need different thersholds for identity for different genes, so i've done a lot of repeats.

ADD REPLY

Login before adding your answer.

Traffic: 1646 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6