This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Local BLAST search

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
• 1,393 views • link
written by v.shapovalova1 • v.shapovalova1

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.

• 0 views • link
written by Joe Joe

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

• 0 views • link
written by v.shapovalova1 • v.shapovalova1

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

• 0 views • link
written by young_bioinformatician young_bioinformatician

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

• 0 views • link
written by v.shapovalova1 • v.shapovalova1

0 answers

No answers yet.

Log in to answer this question.

More posts like this