(Fyi: your post came through in a tiny (illegible) font.)
> #!/usr/sbin/perl
I recommend you start your scripts:
#!/usr/sbin/perl -w
use strict;
This will help identify a lot of basic mistakes.
> if ($line =~ /^([^ ]+)\s+([-a-zA-Z*.]+) *$/) {
The [^] bit is the
First I have some questions.
Why is ALIGNMENT not a hash ? It seems perfect for what you are need. If
you use ALIGNMENT as a hash, then you can do away with the ID array.
Also you are checking for existence in alignment not ALIGNMENT, this will
be causing problems.
Simplify like this:
mine:
[
Dear ME, sorry I am making myself not clear. Let's try
in another way. Here is the script again
#!/usr/sbin/perl
if (!@ARGV) {
print STDERR "usage: $0 alignment_file
[threshold%]...\n";
print_sets();
exit 0;
}
my $FILE = shift
@ARGV;
my @THRESHOLD;
if (@ARGV) {
@THRESHOLD =
> Hi, I have a question regarding the following script.
> [code]
Did you write the code?
> This should be the output.
> [data1]
> However, the output is the following:
> [data2]
Did you know that [data2] is near enough just [data1], twice?
Which is probably because the somethings in this cod
Hi Paul, thank for your interest. I am talking
about multiple protein sequence alignments generated by the program
clustalW (see http://www.ebi.ac.uk/clustalw/help.html for additional information).
Since the sequences to be aligned can be very long, in the output clustalW
split the sequences
--- Pedro A Reche Gallardo <[EMAIL PROTECTED]> wrote:
> Hi, I have a question regarding the following script.
[snip]
> Heres is the question.
> This script should take an alignment with sequences spread in two
> or more blocks, and print them in one single block. Se below
Um, what does that mean
Hi, I have a question regarding the following script.
#!/usr/sbin/perl
if (!@ARGV) {
print STDERR "usage: $0 alignment_file [threshold%]...\n";
print_sets();
exit 0;
}
my $FILE = shift @ARGV;
my @THRESHOLD;
if (@ARGV) {
@THRESHOLD = @ARGV;
} else {
@THRESHOLD = (90, 80,
