# this works as a bash include
#
#COMM SeqLab::Align - clustalw format I/O
#
function flt_clwhead {
  perl -ne 'if(m/^CLUSTAL W \(/){ scalar <> } else{ print }' "$@"
}

perl -e '
    use lib shift();  # give include precedence to code directory, perl -I would not work!
    use SeqLab::Align;
    my ($fin) = @ARGV;
    open (IN,$fin); my @d; while(<IN>){ push @d,$_ }; my $d=join(q(),@d);
    my $pAln = &AlnClustalRead(\$d);
    print &AlnClustalSprint($pAln);
  ' -- "$dircode" "$dirtest"/${setstamp}_data.aln  \
| flt_clwhead -
