in reply to how do i obtain blast result from the given file
---------------------------------------------------------------------
#!usr/bin/perl -w
use strict;
use warnings;
open(IN,"/home/maize/sequence-ID.txt") || die $!;
open(OUT,">abc.txt") || die $!;
open(OUT1,">blastresults.txt") || die $!;
while(my $id=<IN>)
{
$id=~s/\n|\r//g;
open(BLAST,"output/$id") || die $!;
while(my $file=<BLAST>)
{
if ($file=~ m/No hits found/)
{
print OUT "$id \n";
}
}
close(BLAST) || die $!;
}
close(IN) || die $!;
close(OUT) || die $!;
close(OUT1) || die $!;
---------------------------------------------------------------------
If you prefer you can use normal Perl variables instead of file handles
---------------------------------------------------------------------
open(my $blast,"output/$id") || die $!;
while(my $file=<$blast>)
{
if ($file=~ m/No hits found/)
{
print OUT "$id \n";
}
}
close($blast) || die $!;
---------------------------------------------------------------------
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^2: how do i obtain blast result from the given file
by bingalee (Acolyte) on Jun 18, 2013 at 19:56 UTC |