[Bioperl-l] Blast to Clustalw Format
Mark A. Jensen
maj at fortinbras.us
Thu Jan 14 13:54:31 UTC 2010
Thanks Marcelo-- code snips always appreciated! MAJ
----- Original Message -----
From: "Marcelo Iwata" <marcelo011982 at gmail.com>
To: <bioperl-l at lists.open-bio.org>
Sent: Thursday, January 14, 2010 8:46 AM
Subject: Re: [Bioperl-l] Blast to Clustalw Format
> Sorry , the correct code is:
>
>
>
> #!/usr/bin/perl -w
>
> use strict;
> use Bio::SearchIO;
> use Bio::AlignIO;
>
> my $in = new Bio::SearchIO(-format => 'blast',
> -file => '
> ../../fontes/exemplos/blat/teste2/output.blast ');
> my $aln;
> my $alnIO;
> $alnIO = Bio::AlignIO->new(-format =>"clustalw", -file => ">hsp.aln");
> while ( my $result = $in->next_result ) {
> ## $result is a Bio::Search::Result::ResultI compliant object
> while ( my $hit = $result->next_hit ) {
> ## $hit is a Bio::Search::Hit::HitI compliant object
> while ( my $hsp = $hit->next_hsp ) {
> ## $hsp is a Bio::Search::HSP::HSPI compliant object
> $aln = $hsp->get_aln;
> $alnIO->write_aln($aln);
>
> }
> }
> }
>
>
> On Thu, Jan 14, 2010 at 11:44 AM, Marcelo Iwata
> <marcelo011982 at gmail.com>wrote:
>
>> Thanks Mark.
>> I think that most of you already know it.
>> But , i'll put it for new users:
>>
>>
>> #!/usr/bin/perl -w
>>
>> use strict;
>> use Bio::SearchIO;
>> use Bio::AlignIO;
>>
>> my $in = new Bio::SearchIO(-format => 'blast',
>> -file => '
>> ../../fontes/exemplos/blat/teste2/output.blast ');
>> my $aln;
>> my $alnIO;
>> $alnIO = Bio::AlignIO->new(-format =>"clustalw", -file => ">hsp.aln");
>> while ( my $result = $in->next_result ) {
>> ## $result is a Bio::Search::Result::ResultI compliant object
>> while ( my $hit = $result->next_hit ) {
>> ## $hit is a Bio::Search::Hit::HitI compliant object
>> while ( my $hsp = $hit->next_hsp ) {
>> ## $hsp is a Bio::Search::HSP::HSPI compliant object
>> $aln = $hsp->get_aln;
>> $alnIO->write_aln($aln);
>>
>>
>> }
>> }
>> }
>>
>>
>> On Wed, Jan 13, 2010 at 4:44 PM, Mark A. Jensen <maj at fortinbras.us> wrote:
>>
>>> Marcelo-
>>> Yes-- look at the code snip at
>>> http://www.bioperl.org/wiki/HOWTO:SearchIO#Using_SearchIO
>>> combined with the snip at
>>> http://www.bioperl.org/wiki/HOWTO:SearchIO#Using_the_methods
>>> (using -format => 'clustalw')
>>> cheers MAJ
>>> ----- Original Message ----- From: "Marcelo Iwata" <
>>> marcelo011982 at gmail.com>
>>> To: <bioperl-l at lists.open-bio.org>
>>> Sent: Wednesday, January 13, 2010 1:12 PM
>>> Subject: [Bioperl-l] Blast to Clustalw Format
>>>
>>>
>>> Hi..
>>>> I have an simple Blast result, such as blastn.
>>>> Is there an scrip to transform such result to Clustalw format in
>>>> Bioperl
>>>> ?(.aln)
>>>>
>>>> Thanx for any help.
>>>> _______________________________________________
>>>> Bioperl-l mailing list
>>>> Bioperl-l at lists.open-bio.org
>>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>>
>>>>
>>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
More information about the Bioperl-l
mailing list