File Coverage

blib/lib/TAIR/Blast.pm

Criterion	Covered	Total	%
statement	12	60	20.0
branch	0	16	0.0
condition			n/a
subroutine	4	7	57.1
pod	3	3	100.0
total	19	86	22.0

package TAIR::Blast;

use warnings;
use strict;

use LWP::UserAgent;
use HTTP::Request::Common;

=head1 NAME

TAIR::Blast - A module to gather automated BLAST result from TAIR (http://www.arabidopsis.org/Blast/index.jsp)

=head1 VERSION

Version 1.01

=cut

our $VERSION = '1.01';


=head1 SYNOPSIS

This module simply automatically! BLAST any type of sequences (nucleotide, protein) with using different type of algorithm (blastp, blastn, tblastx etc.) by using TAIR Blast engine.

	use TAIR::Blast;
	use Bio::DB::Fasta;

	my $TB = TAIR::Blast->new();

	my $fasta_file = "seqs.fasta";
	my $algorithm = "blastn"; #other set : blastp, blastx, tblastx,tblastn
	my $maxscore;
	my $blast_target_set;
	my $verbose = 1;

	my $stream  = Bio::DB::Fasta->new($fasta_file)->get_PrimarySeq_stream;
	  while (my $query = $stream->next_seq) {
		my $query_seq = $query->seq;
		my $result = $TB->connect($query,$query_seq,$algorithm,$maxscore,$blast_target_set,$verbose);
	  }
	my $output_file = "output_file.txt";
	$TB->write($output_file);

=head1 SUBROUTINES/METHODS

=head2 new

Object-oriented master-hash!

=cut

sub new
  {
    my %master_hash;
    return(bless(\%master_hash, __PACKAGE__));
  }

=head2 connect

Blast parameters, do not mess with them unless you know what you are doing.

=cut

sub connect
 {
	my $tair_master = shift;
	my $query = shift;
	my $query_seq = shift;
	my $algorithm = shift;
	#my $out_type = shift; #will be used in v2
	my $maxscore = shift;
	my $blast_target_set = shift;
	my $verbose = shift;

	if (!defined($blast_target_set)) {
		if ($algorithm =~ /^blast[px]$/i){
			$blast_target_set = "ATH1_pep";
		}elsif ($algorithm =~ /^blastn|tblastx|tblastn$/i){
			$blast_target_set = "At_transcripts";
		}
	}

	my $URLtoPostTo = "http://www.arabidopsis.org/cgi-bin/Blast/TAIRblast.pl";
	my %Fields = (
   "Algorithm" => $algorithm,
   "default_db" => "At_transcripts",
   "BlastTargetSet" => $blast_target_set,
   "textbox" => "seq",
   "QueryText" => $query_seq,
	"QueryFilter" => "T",
	"Matrix" => "Blosum62",
	"MaxScore" => "3",
	"Expectation" => "10",
	"MaxAlignments" => $maxscore,
	"NucleicMismatch" => "-3",
	"GappedAlignment" => "T",
	"NucleicMatch" => "1",
	"OpenPenalty" => "0",
	"ExtensionThreshold" => "0",
	"ExtendPenalty" => "0",
	"WordSize" => "0",
	"QueryGeneticCode" => "1",
	"Comment" => "optional, will be added to output for your use",
	"ReplyTo" => "",
	"ReplyVia" => "BROWSER",
	#"ReplyFormat" => $out_type, #will be used in v2
	"ReplyFormat" => "TABULAR",
	"PageType" => "JavaScr"
	);

	my $BrowserName = "TAIR";

	my $Browser = new LWP::UserAgent;

	if($BrowserName) { $Browser->agent($BrowserName); }

	my $Page = $Browser->request(POST $URLtoPostTo,\%Fields);

	if ($Page->is_success) {
	my $output_raw = $Page->content;
	$output_raw =~ s/TAIR Blast Job Pending//g;
	my @output_array = split("\n", $output_raw); #[0]query/[1]hit/[10]e-value/[11]Score(bits)
	my %output_hsps;
	my @output_hsps;
	my $i = 1;
		foreach (@output_array){
			my @details = split("\t",$_);
				if (defined($details[1])) {
					my %details = (
					"hsp" => $i,
					"hits" => $details[1],
					"e-value" => $details[10],
					"score" => $details[11]
					);
					shift(@details);
					push (@output_hsps, \%details);
					$i++;
				}
		}
	$tair_master->{'output'}->{$query}=\@output_hsps;
	if ($verbose == 1){
		print $query, " is blasted!\n";
	}
 	}
	else { print $Page->message; }

 }

=head2 write

This subroutine handles output of the program, it writes the blast results into a specified file name.

=cut

sub write
	{
		my $tair_master = shift;
		my $output_file = shift;
		my ($output_file_hn);
		my $output_hash = $tair_master->{'output'};
		open($output_file_hn, ">", $output_file) or die("Cannot open $output_file for writing: $!");
		print $output_file_hn "query\thits\tscore\te-value\n";
			foreach (keys %$output_hash){
			my $next_query = $_;
				foreach(@{$output_hash->{$_}}){
				print $output_file_hn $next_query,"\t",
						$_->{'hits'},"\t",
						$_->{'score'},"\t",
						$_->{'e-value'},"\n";
				}
			}
	}

=head1 AUTHOR

Haktan Suren, <<  >>

=head1 BUGS

Please report any bugs or feature requests to C, or through
the web interface at L.  I will be notified, and then you'll
automatically be notified of progress on your bug as I make changes.




=head1 SUPPORT

You can find documentation for this module with the perldoc command.

    perldoc TAIR::Blast


You can also look for information at:

=over 4

=item * RT: CPAN's request tracker

L

=item * AnnoCPAN: Annotated CPAN documentation

L

=item * CPAN Ratings

L

=item * Search CPAN

L

=back

=head1 DEPENDENCIES

	LWP::UserAgent
	HTTP::Request::Common

=head1 ACKNOWLEDGEMENTS

Special thanks to LC :)

=head1 LICENSE AND COPYRIGHT

Copyright 2011 Haktan Suren.

This program is free software; you can redistribute it and/or modify it
under the terms of either: the GNU General Public License as published
by the Free Software Foundation; or the Artistic License.

See http://dev.perl.org/licenses/ for more information.


=cut

1; # End of TAIR::Blast