<div dir="ltr">Hello!<div><br></div><div>I am trying to perform online a PSI-BLAST. In the tutorial there isn't any example so looking on the help I found what I thought it was the way:</div><div><br></div><div><div>record = SeqIO.read("protein.fasta", format="fasta")</div><div>result_handle = NCBIWWW.qblast('blastp', 'nr', record.format("fasta"),</div><div>                                    service='psi')</div></div><div>blast_record = NCBIXML.read(result_handle)</div><div><br></div><div>However it returns a normal Blast record. But the object I get is not the "<span style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium">PSIBlast record</span><span style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"> " </span>described on the tutorial with information about. I haven't found the way to set the number of round of the PSI-BLAST (neither in the NCBI help at <a href="http://www.ncbi.nlm.nih.gov/BLAST/Doc/urlapi.html">http://www.ncbi.nlm.nih.gov/BLAST/Doc/urlapi.html</a> ). </div><div><br></div><div><br></div><div>How can I do a PSI-BLAST? </div><div><br></div><div>Many thanks,</div><div><br></div><div>Lluís</div></div>