my @locus_small = (); my $line; foreach $line (@locus) { my @tokens = split(/\t+/, $line); unless(scalar @tokens < 6) { push(@locus_small, "$tokens[0]\t$tokens[1]\t"); } } #### foreach my $molecule (@molecules) { my @found = grep /\Q$molecule\E/i, @locus_small; if (@found) { print OUTDATA ($molecule, ": \n\t", join "\t", @found); } }