SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bioperl installation chandrugcg Bioinformatics 2 01-28-2013 02:45 PM
does bioperl::seqIO:: Kojo Bioinformatics 2 10-08-2012 04:49 AM
Simulations in Bioperl? A_Morozov Bioinformatics 0 04-11-2012 11:10 PM
bioperl for samtools m_elena_bioinfo Bioinformatics 4 09-07-2010 05:39 AM
Bioperl jsun529 Bioinformatics 7 09-01-2009 07:27 AM

Reply
 
Thread Tools
Old 02-19-2013, 07:36 PM   #1
noobie
Junior Member
 
Location: USA

Join Date: Jun 2012
Posts: 7
Default BioPerl Issues

I am having issues with a BioPerl script. I have a blastxml file from a blastx blast and the original multifasta file containing the original nucleotides sequences.

I want to take the blast result (ie. the blast description) and annotated my multifasta file.

I have written 2 while loops that extract the blast descriptions as well as the nucleotide sequence from the multifasta file.

My problem is that I cannot incorporate one of the while loops into the other without loosing the loop property of one of the loops.

See script below:


use warnings;
use strict;
use Bio::SearchIO;
use Bio::SeqIO;


my $search_in = Bio::SearchIO->new(-format => 'blastxml', -file => "$ARGV[0]");
while (my $result = $search_in->next_result) {
while (my $hit = $result->next_hit) {
while (my $hsp = $hit->next_hsp) {
my $qd = $hit->description;
print $qd, "\n";
}
}
}

my $seqio = Bio::SeqIO->new(-format => 'fasta', -file => "$ARGV[1]");
while (my $seqobj = $seqio->next_seq) {
my $nuc = $seqobj->seq();
print $nuc, "\n";
}
noobie is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:46 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO