my @seq1; my @seq2; load_seqs(\@seq1,\@seq2); sub load_seqs{ open(my $fh, "<", "q1.fa") or die "cannot open sequence file: $!"; while(my $line=<$fh>){ if($line =~ m/^(?!\>).*\n$/){push (@{$_[0]},$line)}; } open(my $fh, "<", "q2.fa") or die "cannot open sequence file: $!"; while(my $line=<$fh>){ if($line =~ m/^(?!\>).*\n$/){push (@{$_[1]},$line)}; } }