[Bioperl-guts-l] bioperl-live/Bio/Tools RNAMotif.pm,1.10,1.10.4.1
Senduran Balasubramaniam
sendu at dev.open-bio.org
Sun Aug 20 18:03:32 EDT 2006
Update of /home/repository/bioperl/bioperl-live/Bio/Tools
In directory dev.open-bio.org:/tmp/cvs-serv10682/Bio/Tools
Modified Files:
Tag: branch-1-5-2
RNAMotif.pm
Log Message:
merge from HEAD, bugfixes and new PullParserI and hmmer parser
Index: RNAMotif.pm
===================================================================
RCS file: /home/repository/bioperl/bioperl-live/Bio/Tools/RNAMotif.pm,v
retrieving revision 1.10
retrieving revision 1.10.4.1
diff -C2 -d -r1.10 -r1.10.4.1
*** RNAMotif.pm 7 Aug 2006 17:24:09 -0000 1.10
--- RNAMotif.pm 20 Aug 2006 22:03:30 -0000 1.10.4.1
***************
*** 316,336 ****
if(/^#RM\sdescr\s(.*)$/) { # contains sec structure
$self->{'_sec_structure'}=$1;
- $self->debug("Found descriptor: $self->{'_sec_structure'}\n");
}
if(/^#RM\sdfile\s(.*)$/) { # contains dfile
$self->{'_dfile'}=$1;
- $self->debug("Found dfile : $self->{'_dfile'}\n");
}
$_ = $self->_readline;
}
if(m/^>((\S*)\s.*)$/) {
- $self->debug("Found a description: $1\n");
$seqid = $2;
$description = $1; # contains entire description line if needed
- # try to get accession from FASTA-like header (if from Genbank) otherwise punt
- # Thanks to J. Stajich for this idea, which is part of SeqIO now
if($seqid =~ /(gb|emb|dbj|sp|pdb|bbs|ref|lcl)\|(.*)\|/) {
$seqid = $2; # pulls out gid
- $self->debug("Genbank gid: $seqid\n");
}
}
--- 316,330 ----
***************
*** 353,360 ****
($start, $end, $strand) = ($end, $start, -1);
}
- if($score =~ /[a-zA-Z\s]/) { #probably containing extra data from sprintf()
- $self->debug
- ("packing sprintf() data into score!\n");
- }
my $gene = Bio::SeqFeature::Generic->new(-seq_id => $seqid,
-start => $start,
--- 347,350 ----
More information about the Bioperl-guts-l
mailing list