Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 41826

Select Sequence From A Psiblast, And Launch A New Iteration From It

$
0
0

I would like on my local machine or remotely if it is possible too, launch a psiblast from a sequence. Then, like on the website of ncbi, select my own sequences from the result file. And launch another iteration from these sequences.

When I launch psiblast on local I obtain a pairwise alignment in xml or other format. How can I select my sequences to go to run another iteration? If I select sequences, I need to make a MSA so, I destroy the blast result if evalue if very high. I would like to do exactly the same choice of psiblast but in a script (which launch blast localy or remotely)

for instance :

blastpgp -i mysequence -j 1 -o mysequence.blast

--> I create a file just with the sequences that I want to keep (not only significative sequence), -> mynewsequenceTo2d_iteration

blastpgp -i mynewsequenceTo2d_iteration -j 1 -o 2diteration.blast

etc...

Is it possible? I did not find any way to do it... I make a multiple alignment between each iteration but Psiblast in ncbi do not that...

I take every way to do it, even with Biopython or CGI


Viewing all articles
Browse latest Browse all 41826

Trending Articles