[Bioperl-l] Very basic Perl/BioPerl Help

Terry Jones tcj25 at cam.ac.uk
Thu Apr 14 12:12:59 EDT 2005


| my $f1=shift;
| my $f2=shift;
| open (F1,$f1)||die;
| open (F2,$f2)||die;
| my @accn1=<F1>;
| my @accn2=<F2>;
| my @unique1;
| foreach my $accn (@accn1) {
| push @unique1,$accn unless (grep(/$accn/, at accn2));
| }

This also just does one half (F1-F2). You need to add F2-F1 to it.

T



More information about the Bioperl-l mailing list