[Bioperl-l] bl2seq parsing

Emilie Lalonde emilie.lalonde at gmail.com
Thu Aug 1 13:21:58 UTC 2013


Hi all,

I'm trying to parse bl2seq output but am not getting any data in the 
returned object. My bl2seq output file contains an alignment, but it isn't 
being loaded into an object, which is undefined.

I've tried by simply parsing the bl2seq report:

# make Seq objects for forward and reverse reads
my $rev_hsp = Bio::Seq->new(-seq =>$proper_pairs{$ampl}{R}->query_string, 
-alphabet => 'dna');
my $fwd_hsp = Bio::Seq->new(-seq =>$proper_pairs{$ampl}{F}->query_string, 
-alphabet => 'dna');

# make Blast factory and run bl2seq
my $paired_blast_obj = Bio::Tools::Run::StandAloneBlast->new(-program => 
'blastn', -outfile => 'blast/paired_alignment/bl2seq.out');
my $bl2seq_report = $paired_blast_obj->bl2seq($rev_hsp, $fwd_hsp);

# parse bl2seq report
if(my $r = $bl2seq_report->next_result){
   while( my $hit = $r->next_hit ) {
     while( my $hsp = $hit->next_hsp ) {
         print $hsp->get_aln->consensus_iupac()."\n";
       }
     }
 }


And I've also tried with AlignIO:
# read and parse bl2seq output file
my $str = Bio::AlignIO->new(-file => 'blast/paired_alignment/bl2seq.out', 
'-format' => 'bl2seq'); 
while(my $aln = $str->next_aln() ){
     print $aln->length . "\n";
}

I've attached my bl2seq report (output file).



Any help would be greatly appreciated! 


Emilie
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CPCG0003F_PIK3CA_chr3_exon13_1_bl2seq.out
Type: application/octet-stream
Size: 2316 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20130801/24aa342f/attachment-0004.obj>


More information about the Bioperl-l mailing list