# this works as a bash include
#
#COMM SeqLab::Align - split alignment (&AlnSplit)
#
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);
    my ($pAln5,$pAln3) = &AlnSplit($pAln,15);
    print &AlnClustalSprint($pAln5);
    print &AlnClustalSprint($pAln3);
  ' -- "$dircode" "$dirtest"/${setstamp}_data.aln  \
| flt_clwhead -
