J00153:42:HC5NCBBXX:6:1101:10896:14959
J00153:42:HC5NCBBXX:6:1101:10896:14959
J00153:42:HC5NCBBXX:6:1101:26616:20709
J00153:42:HC5NCBBXX:6:1101:27549:19935
####
J00153:42:HC5NCBBXX:6:1101:10896:14959 99 gnl|Btau_4.6.1|chr16 72729218 1 12M
J00153:42:HC5NCBBXX:6:1101:27549:19935 83 gnl|Btau_4.6.1|chr8 49556412 1 7M
####
# Open query file and read into array
$queryfile = $ARGV[0];
open (QUERYFILE, $queryfile) or die "Cannot open query file\n";
@queries = ;
close QUERYFILE;
# Open main file
$mainfile = $ARGV[1];
open (MAINFILE, $mainfile) or die "Cannot open searchable file\n";
# Search through main file
while ($inline = ) {
@split = split /\t/, $inline;
$ID = $split[0];
if (grep /$ID/, @queries) {
print $inline;
}
else {
}
}
exit;
####
# Open query file and read into array
$queryfile = $ARGV[0];
open (QUERYFILE, $queryfile) or die "Cannot open query file\n";
@queries = ;
close QUERYFILE;
# Convert array to hash
%hash = map {$_ => 1} @queries;
# Open main file
$mainfile = $ARGV[1];
open (MAINFILE, $mainfile) or die "Cannot open searchable file\n";
# Search through main file
while ($inline = ) {
@split = split /\t/, $inline;
$ID = $split[0];
if (defined $hash{$ID}) {
print $inline;
}
else {
}
}
exit;