Beefy Boxes and Bandwidth Generously Provided by pair Networks
Welcome to the Monastery
 
PerlMonks  

Re: use of diamond operator for multiple files

by hdb (Prior)
on Apr 24, 2014 at 12:16 UTC ( #1083573=note: print w/ replies, xml ) Need Help??


in reply to use of diamond operator for multiple files

Two comments from my side. First, your subroutine does the same thing twice. I would change it such that it reads the filtered contents of a given file into a given array and call it twice. Second, your while loop is essential equivalent to grep, so in summary, I would rewrite your code as:

my @seq1 = load_seqs( "q1.fa" ); my @seq2 = load_seqs( "q2.fa" ); sub load_seqs{ open my $fh, "<", shift or die "cannot open sequence file: $!"; return grep { /^(?!\>).*\n$/ } <$fh>; }


Comment on Re: use of diamond operator for multiple files
Download Code

Log In?
Username:
Password:

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

How do I use this? | Other CB clients
Other Users?
Others rifling through the Monastery: (10)
As of 2014-09-16 21:03 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    How do you remember the number of days in each month?











    Results (48 votes), past polls