in reply to Re: Can you make this code shorter and/or quicker as well?
in thread Can you make this code shorter and/or quicker as well?

How can I put it in a sub-routine that I could use inline my script?
Something like:
sub extract_query_data($file);

[download]

Replies are listed 'Best First'.
Re^3: Can you make this code shorter and/or quicker as well?
by Anonymous Monk on Feb 25, 2014 at 00:27 UTC
    Sorry, I meant: 
    extract_query_data($file);

    [download]
[reply]
[d/l]

      Here's one option:

      use strict;
use warnings;

for my $fastaRec ( extract_query_data('fastaFile.fa') ) {
    print $fastaRec, "\n";
}

sub extract_query_data {
    my ($file) = @_;
    my @arr;
    local $/ = '>';

    open my $fh, '<', $file or die $!;

    while (<$fh>) {
        chomp;
        my ( $id, $seq ) = /(.+?\n)(.+)/s or next;
        $seq =~ s/\s+//g;
        push @arr, ">$id$seq";
    }

    close $fh;

    return @arr;
}

      [download]

      You'll note that the subroutine returns an array of fasta records.

[reply]
[d/l]
        Aha,
        I see! Thank you very much for your time!
[reply]

In Section Seekers of Perl Wisdom