Beefy Boxes and Bandwidth Generously Provided by pair Networks
Do you know where your variables are?
 
PerlMonks  

Re: Extracting BLAST hits from a list of sequences

by no_slogan (Deacon)
on Jan 20, 2014 at 00:38 UTC ( [id://1071249]=note: print w/replies, xml ) Need Help??


in reply to Extracting BLAST hits from a list of sequences

I know nothing about BLAST files, but it's probably because

my $blastcut = substr $blastline, 7;

is sometimes an empty string. Maybe do something like this:

if ($blastline =~ /$ender/ && $headerline =~ /\Q$blastcut\E/) {

or this:

my $printing = 0; while (<BLASTFILE>) { if (/Query=\s*(.*)/) { $printing = (index($headerline,$1) >= 0); } print if $printing; }

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://1071249]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others browsing the Monastery: (8)
As of 2024-04-19 07:24 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found