[Bioperl-l] Assignment problem :(

Kurt Gobain enrique_rulz at yahoo.com
Thu Jan 4 22:15:16 UTC 2007

Thanx for the replys...As suggested b Florin...I m posting the prog which I
had done so far...
use strict;
use warnings;

my $dir= "$ARGV[0]" or die "inadequate argument at command line use 
my $database="$ARGV[1]" or die "inadequate argument at commandline use 
/databases/swissprot/swissprot\n" ;
my $file;
my ($i,$seq);
my $length=0;
open (OUTPUT,">result.fasta");
opendir (DIR,"$dir") or die "cannot open the specified directory\n";
foreach $file (readdir DIR)
	if ($file =~ /\.embl/)
		open(FILE,"$dir/$file") or die "cannot open the specified file\n";
		#print ">$file\n";
		#print OUTPUT ">$file\n";
		while (<FILE>)
				if (/^\s/)
					my @array = split ('',$sequence);
					foreach my $fasta_seq (@array)
						if ( $fasta_seq =~ /[A-Za-z]/)
							print "$fasta_seq";
							print OUTPUT "$fasta_seq";
		print "length of fasta sequence is $length\n";
		my $fasta_seq=@_;
		for ($i=0,$i<6,$i++)
		my $result="blastall -F T -p blastp -d $database -i result.fasta -o 
	print "\n";
	print OUTPUT"\n";

close (DIR);
close (FILE);

I am stuck on bit no 5..

""for each sequence print the swissprot name and code of its close
homologues to an HTML formatted file that can be viewed in a web browser.Is
it possible to extract the evalues for each hit & record them in the HTML

I am told we need to use sum Bioperl Packages. like Bio::DB::SwissProt,
Bio::SeqIO;....It wud be great if any 1 can shed lighte on this....

Thanx in advance..
View this message in context: http://www.nabble.com/Assignment-problem-%3A%28-tf2913859.html#a8168772
Sent from the Perl - Bioperl-L mailing list archive at Nabble.com.

More information about the Bioperl-l mailing list