annotate variant_effect_predictor/Bio/SearchIO/exonerate.pm @ 1:d6778b5d8382 draft default tip

Deleted selected files
author willmclaren
date Fri, 03 Aug 2012 10:05:43 -0400
parents 21066c0abaf5
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
1 # $Id: exonerate.pm,v 1.3.2.3 2003/03/29 20:30:54 jason Exp $
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
2 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
3 # BioPerl module for Bio::SearchIO::exonerate
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
4 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
5 # Cared for by Jason Stajich <jason@bioperl.org>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
6 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
7 # Copyright Jason Stajich
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
8 #
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
9 # You may distribute this module under the same terms as perl itself
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
10
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
11 # POD documentation - main docs before the code
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
12
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
13 =head1 NAME
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
14
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
15 Bio::SearchIO::exonerate - parser for Exonerate
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
16
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
17 =head1 SYNOPSIS
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
18
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
19 # do not use this module directly, it is a driver for SearchIO
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
20
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
21 use Bio::SearchIO;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
22 my $searchio = new Bio::SearchIO(-file => 'file.exonerate',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
23 -format => 'exonerate');
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
24
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
25
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
26 while( my $r = $searchio->next_result ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
27 print $r->query_name, "\n";
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
28 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
29
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
30 =head1 DESCRIPTION
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
31
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
32 This is a driver for the SearchIO system for parsing Exonerate (Guy
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
33 Slater) output. You can get Exonerate at
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
34 http://cvsweb.sanger.ac.uk/cgi-bin/cvsweb.cgi/exonerate/?cvsroot=Ensembl
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
35 [until Guy puts up a Web reference,publication for it.]).
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
36
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
37 An optional parameter -min_intron is supported by the L<new>
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
38 initialization method. This is if you run Exonerate with a different
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
39 minimum intron length (default is 30) the parser will be able to
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
40 detect the difference between standard deletions and an intron. Still
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
41 some room to play with there that might cause this to get
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
42 misinterpreted that has not been fully tested or explored.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
43
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
44 =head1 FEEDBACK
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
45
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
46 =head2 Mailing Lists
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
47
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
48 User feedback is an integral part of the evolution of this and other
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
49 Bioperl modules. Send your comments and suggestions preferably to
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
50 the Bioperl mailing list. Your participation is much appreciated.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
51
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
52 bioperl-l@bioperl.org - General discussion
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
53 http://bioperl.org/MailList.shtml - About the mailing lists
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
54
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
55 =head2 Reporting Bugs
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
56
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
57 Report bugs to the Bioperl bug tracking system to help us keep track
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
58 of the bugs and their resolution. Bug reports can be submitted via
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
59 email or the web:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
60
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
61 bioperl-bugs@bioperl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
62 http://bioperl.org/bioperl-bugs/
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
63
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
64 =head1 AUTHOR - Jason Stajich
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
65
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
66 Email jason@bioperl.org
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
67
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
68 Describe contact details here
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
69
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
70 =head1 CONTRIBUTORS
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
71
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
72 Additional contributors names and emails here
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
73
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
74 =head1 APPENDIX
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
75
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
76 The rest of the documentation details each of the object methods.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
77 Internal methods are usually preceded with a _
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
78
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
79 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
80
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
81
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
82 # Let the code begin...
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
83
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
84
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
85 package Bio::SearchIO::exonerate;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
86 use strict;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
87 use vars qw(@ISA @STATES %MAPPING %MODEMAP $DEFAULT_WRITER_CLASS $MIN_INTRON);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
88 use Bio::SearchIO;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
89
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
90 @ISA = qw(Bio::SearchIO );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
91
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
92 use POSIX;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
93
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
94
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
95 %MODEMAP = ('ExonerateOutput' => 'result',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
96 'Hit' => 'hit',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
97 'Hsp' => 'hsp'
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
98 );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
99 %MAPPING =
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
100 (
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
101 'Hsp_query-from'=> 'HSP-query_start',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
102 'Hsp_query-to' => 'HSP-query_end',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
103 'Hsp_hit-from' => 'HSP-hit_start',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
104 'Hsp_hit-to' => 'HSP-hit_end',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
105 'Hsp_qseq' => 'HSP-query_seq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
106 'Hsp_hseq' => 'HSP-hit_seq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
107 'Hsp_midline' => 'HSP-homology_seq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
108 'Hsp_score' => 'HSP-score',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
109 'Hsp_qlength' => 'HSP-query_length',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
110 'Hsp_hlength' => 'HSP-hit_length',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
111 'Hsp_align-len' => 'HSP-hsp_length',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
112 'Hsp_identity' => 'HSP-identical',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
113 'Hsp_gaps' => 'HSP-hsp_gaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
114 'Hsp_hitgaps' => 'HSP-hit_gaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
115 'Hsp_querygaps' => 'HSP-query_gaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
116
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
117 'Hit_id' => 'HIT-name',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
118 'Hit_desc' => 'HIT-description',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
119 'Hit_len' => 'HIT-length',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
120 'Hit_score' => 'HIT-score',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
121
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
122 'ExonerateOutput_program' => 'RESULT-algorithm_name',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
123 'ExonerateOutput_query-def' => 'RESULT-query_name',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
124 'ExonerateOutput_query-desc'=> 'RESULT-query_description',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
125 'ExonerateOutput_query-len' => 'RESULT-query_length',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
126 );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
127
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
128 $DEFAULT_WRITER_CLASS = 'Bio::Search::Writer::HitTableWriter';
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
129
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
130 $MIN_INTRON=30; # This is the minimum intron size
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
131
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
132 =head2 new
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
133
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
134 Title : new
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
135 Usage : my $obj = new Bio::SearchIO::exonerate();
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
136 Function: Builds a new Bio::SearchIO::exonerate object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
137 Returns : an instance of Bio::SearchIO::exonerate
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
138 Args :
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
139
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
140
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
141 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
142
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
143 sub new {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
144 my ($class) = shift;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
145 my $self = $class->SUPER::new(@_);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
146
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
147 my ($min_intron) = $self->_rearrange([qw(MIN_INTRON)], @_);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
148 if( $min_intron ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
149 $MIN_INTRON = $min_intron;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
150 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
151 $self;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
152 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
153
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
154 =head2 next_result
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
155
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
156 Title : next_result
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
157 Usage : my $hit = $searchio->next_result;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
158 Function: Returns the next Result from a search
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
159 Returns : Bio::Search::Result::ResultI object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
160 Args : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
161
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
162 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
163
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
164 sub next_result{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
165 my ($self) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
166 $self->{'_last_data'} = '';
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
167 my ($reporttype,$seenquery,$reportline);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
168 $self->start_document();
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
169 my @hit_signifs;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
170 my $seentop;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
171 my (@q_ex, @m_ex, @h_ex); ## gc addition
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
172 while( defined($_ = $self->_readline) ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
173 if( /^Query:\s+(\S+)(\s+(.+))?/ ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
174 if( $seentop ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
175 $self->end_element({'Name' => 'ExonerateOutput'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
176 $self->_pushback($_);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
177 return $self->end_document();
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
178 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
179 $seentop = 1;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
180 my ($nm,$desc) = ($1,$2);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
181 chomp($desc) if defined $desc;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
182 $self->{'_result_count'}++;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
183 $self->start_element({'Name' => 'ExonerateOutput'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
184 $self->element({'Name' => 'ExonerateOutput_query-def',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
185 'Data' => $nm });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
186 $self->element({'Name' => 'ExonerateOutput_query-desc',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
187 'Data' => $desc });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
188 $self->element({'Name' => 'ExonerateOutput_program',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
189 'Data' => 'Exonerate' });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
190
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
191 } elsif ( /^Target:\s+(\S+)(\s+(.+))?/ ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
192 my ($nm,$desc) = ($1,$2);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
193 chomp($desc) if defined $desc;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
194 $self->start_element({'Name' => 'Hit'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
195 $self->element({'Name' => 'Hit_id',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
196 'Data' => $nm});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
197 $self->element({'Name' => 'Hit_desc',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
198 'Data' => $desc});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
199 } elsif( s/^cigar:\s+(\S+)\s+ # query sequence id
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
200 (\d+)\s+(\d+)\s+([\-\+])\s+ # query start-end-strand
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
201 (\S+)\s+ # target sequence id
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
202 (\d+)\s+(\d+)\s+([\-\+])\s+ # target start-end-strand
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
203 (\d+)\s+ # score
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
204 //ox ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
205
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
206 ## gc note:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
207 ## $qe and $he are no longer used for calculating the ends,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
208 ## just the $qs and $hs values and the alignment and insert lenghts
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
209 my ($qs,$qe,$qstrand) = ($2,$3,$4);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
210 my ($hs,$he,$hstrand) = ($6,$7,$8);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
211 my $score = $9;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
212 # $self->element({'Name' => 'ExonerateOutput_query-len',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
213 # 'Data' => $qe});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
214 # $self->element({'Name' => 'Hit_len',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
215 # 'Data' => $he});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
216
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
217 my @rest = split;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
218 if( $qstrand eq '-' ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
219 $qstrand = -1;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
220 ($qs,$qe) = ($qe,$qs); # flip-flop if we're on opp strand
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
221 $qs--; $qe++;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
222 } else { $qstrand = 1; }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
223 if( $hstrand eq '-' ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
224 $hstrand = -1;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
225 ($hs,$he) = ($he,$hs); # flip-flop if we're on opp strand
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
226 $hs--; $he++;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
227 } else { $hstrand = 1; }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
228 # okay let's do this right and generate a set of HSPs
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
229 # from the cigar line
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
230
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
231 ## gc note:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
232 ## add one because these values are zero-based
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
233 ## this calculation was originally done lower in the code,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
234 ## but it's clearer to do it just once at the start
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
235 $qs++;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
236 $hs++;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
237
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
238 my ($aln_len,$inserts,$deletes) = (0,0,0);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
239 while( @rest >= 2 ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
240 my ($state,$len) = (shift @rest, shift @rest);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
241 if( $state eq 'I' ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
242 $inserts+=$len;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
243 } elsif( $state eq 'D' ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
244 if( $len >= $MIN_INTRON ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
245 $self->start_element({'Name' => 'Hsp'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
246
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
247 $self->element({'Name' => 'Hsp_score',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
248 'Data' => $score});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
249 $self->element({'Name' => 'Hsp_align-len',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
250 'Data' => $aln_len});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
251 $self->element({'Name' => 'Hsp_identity',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
252 'Data' => $aln_len -
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
253 ($inserts + $deletes)});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
254
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
255 # HSP ends where the other begins
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
256 $self->element({'Name' => 'Hsp_query-from',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
257 'Data' => $qs});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
258 ## gc note:
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
259 ## $qs is now the start of the next hsp
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
260 ## the end of this hsp is 1 before this position
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
261 ## (or 1 after in case of reverse strand)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
262 $qs += $aln_len*$qstrand;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
263 $self->element({'Name' => 'Hsp_query-to',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
264 'Data' => $qs - ($qstrand*1)});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
265
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
266 $hs += $deletes*$hstrand;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
267 $self->element({'Name' => 'Hsp_hit-from',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
268 'Data' => $hs});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
269 $hs += $aln_len*$hstrand;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
270 $self->element({'Name' => 'Hsp_hit-to',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
271 'Data' => $hs-($hstrand*1)});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
272
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
273 $self->element({'Name' => 'Hsp_align-len',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
274 'Data' => $aln_len + $inserts
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
275 + $deletes});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
276 $self->element({'Name' => 'Hsp_identity',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
277 'Data' => $aln_len });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
278
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
279 $self->element({'Name' => 'Hsp_gaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
280 'Data' => $inserts + $deletes});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
281 $self->element({'Name' => 'Hsp_querygaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
282 'Data' => $inserts});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
283 $self->element({'Name' => 'Hsp_hitgaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
284 'Data' => $deletes});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
285
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
286 ## gc addition start
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
287
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
288 $self->element({'Name' => 'Hsp_qseq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
289 'Data' => shift @q_ex,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
290 });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
291 $self->element({'Name' => 'Hsp_hseq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
292 'Data' => shift @h_ex,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
293 });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
294 $self->element({'Name' => 'Hsp_midline',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
295 'Data' => shift @m_ex,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
296 });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
297 ## gc addition end
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
298 $self->end_element({'Name' => 'Hsp'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
299
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
300 $aln_len = $inserts = $deletes = 0;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
301 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
302 $deletes+=$len;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
303 } else {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
304 $aln_len += $len;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
305 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
306 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
307 $self->start_element({'Name' => 'Hsp'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
308
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
309 ## gc addition start
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
310
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
311 $self->element({'Name' => 'Hsp_qseq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
312 'Data' => shift @q_ex,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
313 });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
314 $self->element({'Name' => 'Hsp_hseq',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
315 'Data' => shift @h_ex,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
316 });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
317 $self->element({'Name' => 'Hsp_midline',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
318 'Data' => shift @m_ex,
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
319 });
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
320 ## gc addition end
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
321
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
322 $self->element({'Name' => 'Hsp_score',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
323 'Data' => $score});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
324
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
325 $self->element({'Name' => 'Hsp_query-from',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
326 'Data' => $qs});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
327
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
328 $qs += $aln_len*$qstrand;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
329 $self->element({'Name' => 'Hsp_query-to',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
330 'Data' => $qs - ($qstrand*1)});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
331
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
332 $hs += $deletes*$hstrand;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
333 $self->element({'Name' => 'Hsp_hit-from',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
334 'Data' => $hs});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
335 $hs += $aln_len*$hstrand;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
336 $self->element({'Name' => 'Hsp_hit-to',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
337 'Data' => $hs -($hstrand*1)});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
338
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
339 $self->element({'Name' => 'Hsp_align-len',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
340 'Data' => $aln_len});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
341
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
342 $self->element({'Name' => 'Hsp_identity',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
343 'Data' => $aln_len - ($inserts + $deletes)});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
344
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
345 $self->element({'Name' => 'Hsp_gaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
346 'Data' => $inserts + $deletes});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
347
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
348 $self->element({'Name' => 'Hsp_querygaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
349 'Data' => $inserts});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
350 $self->element({'Name' => 'Hsp_hitgaps',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
351 'Data' => $deletes});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
352 $self->end_element({'Name' => 'Hsp'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
353 $self->element({'Name' => 'Hit_score',
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
354 'Data' => $score});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
355 $self->end_element({'Name' => 'Hit'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
356 $self->end_element({'Name' => 'ExonerateOutput'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
357
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
358 return $self->end_document();
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
359 } else {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
360 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
361 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
362 return $self->end_document() if( $seentop );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
363 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
364
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
365 =head2 start_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
366
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
367 Title : start_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
368 Usage : $eventgenerator->start_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
369 Function: Handles a start element event
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
370 Returns : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
371 Args : hashref with at least 2 keys 'Data' and 'Name'
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
372
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
373
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
374 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
375
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
376 sub start_element{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
377 my ($self,$data) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
378 # we currently don't care about attributes
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
379 my $nm = $data->{'Name'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
380 my $type = $MODEMAP{$nm};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
381
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
382 if( $type ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
383 if( $self->_eventHandler->will_handle($type) ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
384 my $func = sprintf("start_%s",lc $type);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
385 $self->_eventHandler->$func($data->{'Attributes'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
386 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
387 unshift @{$self->{'_elements'}}, $type;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
388
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
389 if($type eq 'result') {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
390 $self->{'_values'} = {};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
391 $self->{'_result'}= undef;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
392 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
393 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
394
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
395 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
396
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
397 =head2 end_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
398
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
399 Title : start_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
400 Usage : $eventgenerator->end_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
401 Function: Handles an end element event
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
402 Returns : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
403 Args : hashref with at least 2 keys 'Data' and 'Name'
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
404
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
405
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
406 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
407
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
408 sub end_element {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
409 my ($self,$data) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
410 my $nm = $data->{'Name'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
411 my $type = $MODEMAP{$nm};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
412 my $rc;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
413
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
414 if( $type = $MODEMAP{$nm} ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
415 if( $self->_eventHandler->will_handle($type) ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
416 my $func = sprintf("end_%s",lc $type);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
417 $rc = $self->_eventHandler->$func($self->{'_reporttype'},
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
418 $self->{'_values'});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
419 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
420 shift @{$self->{'_elements'}};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
421
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
422 } elsif( $MAPPING{$nm} ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
423
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
424 if ( ref($MAPPING{$nm}) =~ /hash/i ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
425 my $key = (keys %{$MAPPING{$nm}})[0];
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
426 $self->{'_values'}->{$key}->{$MAPPING{$nm}->{$key}} = $self->{'_last_data'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
427 } else {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
428 $self->{'_values'}->{$MAPPING{$nm}} = $self->{'_last_data'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
429 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
430 } else {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
431 $self->debug( "unknown nm $nm, ignoring\n");
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
432 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
433 $self->{'_last_data'} = ''; # remove read data if we are at
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
434 # end of an element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
435 $self->{'_result'} = $rc if( defined $type && $type eq 'result' );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
436 return $rc;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
437 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
438
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
439 =head2 element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
440
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
441 Title : element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
442 Usage : $eventhandler->element({'Name' => $name, 'Data' => $str});
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
443 Function: Convience method that calls start_element, characters, end_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
444 Returns : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
445 Args : Hash ref with the keys 'Name' and 'Data'
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
446
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
447
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
448 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
449
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
450 sub element{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
451 my ($self,$data) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
452 $self->start_element($data);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
453 $self->characters($data);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
454 $self->end_element($data);
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
455 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
456
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
457 =head2 characters
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
458
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
459 Title : characters
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
460 Usage : $eventgenerator->characters($str)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
461 Function: Send a character events
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
462 Returns : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
463 Args : string
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
464
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
465
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
466 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
467
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
468 sub characters{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
469 my ($self,$data) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
470
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
471 return unless ( defined $data->{'Data'} && $data->{'Data'} !~ /^\s+$/ );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
472
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
473 $self->{'_last_data'} = $data->{'Data'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
474 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
475
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
476 =head2 within_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
477
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
478 Title : within_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
479 Usage : if( $eventgenerator->within_element($element) ) {}
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
480 Function: Test if we are within a particular element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
481 This is different than 'in' because within can be tested
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
482 for a whole block.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
483 Returns : boolean
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
484 Args : string element name
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
485
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
486
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
487 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
488
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
489 sub within_element{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
490 my ($self,$name) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
491 return 0 if ( ! defined $name &&
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
492 ! defined $self->{'_elements'} ||
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
493 scalar @{$self->{'_elements'}} == 0) ;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
494 foreach ( @{$self->{'_elements'}} ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
495 if( $_ eq $name ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
496 return 1;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
497 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
498 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
499 return 0;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
500 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
501
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
502
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
503 =head2 in_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
504
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
505 Title : in_element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
506 Usage : if( $eventgenerator->in_element($element) ) {}
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
507 Function: Test if we are in a particular element
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
508 This is different than 'in' because within can be tested
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
509 for a whole block.
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
510 Returns : boolean
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
511 Args : string element name
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
512
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
513
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
514 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
515
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
516 sub in_element{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
517 my ($self,$name) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
518 return 0 if ! defined $self->{'_elements'}->[0];
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
519 return ( $self->{'_elements'}->[0] eq $name)
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
520 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
521
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
522 =head2 start_document
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
523
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
524 Title : start_document
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
525 Usage : $eventgenerator->start_document
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
526 Function: Handle a start document event
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
527 Returns : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
528 Args : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
529
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
530
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
531 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
532
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
533 sub start_document{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
534 my ($self) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
535 $self->{'_lasttype'} = '';
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
536 $self->{'_values'} = {};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
537 $self->{'_result'}= undef;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
538 $self->{'_elements'} = [];
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
539 $self->{'_reporttype'} = 'exonerate';
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
540 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
541
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
542
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
543 =head2 end_document
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
544
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
545 Title : end_document
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
546 Usage : $eventgenerator->end_document
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
547 Function: Handles an end document event
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
548 Returns : Bio::Search::Result::ResultI object
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
549 Args : none
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
550
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
551
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
552 =cut
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
553
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
554 sub end_document{
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
555 my ($self,@args) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
556 return $self->{'_result'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
557 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
558
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
559
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
560 sub write_result {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
561 my ($self, $blast, @args) = @_;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
562
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
563 if( not defined($self->writer) ) {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
564 $self->warn("Writer not defined. Using a $DEFAULT_WRITER_CLASS");
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
565 $self->writer( $DEFAULT_WRITER_CLASS->new() );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
566 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
567 $self->SUPER::write_result( $blast, @args );
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
568 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
569
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
570 sub result_count {
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
571 my $self = shift;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
572 return $self->{'_result_count'};
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
573 }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
574
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
575 sub report_count { shift->result_count }
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
576
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
577 1;
21066c0abaf5 Uploaded
willmclaren
parents:
diff changeset
578