NAME

Bio::MUST::Core::Ali::Temporary - Thin wrapper for a temporary mapped Ali written on disk

VERSION

version 0.180140

SYNOPSIS

    #!/usr/bin/env perl

    use Modern::Perl '2011';
    # same as:
    # use strict;
    # use warnings;
    # use feature qw(say);

    use Bio::MUST::Core;
    use aliased 'Bio::MUST::Core::Ali::Temporary';

    # build Ali::Temporary object from existing ALI file
    my $temp_db = Temporary->new( seqs => 'database.ali');

    # get properties
    my $db = $temp_db->filename;
    my $dbtype = $temp_db->type;

    # pass it to external program
    system("makeblastdb -in $db -dbtype $dbtype");

    # alternative constructor call
    # build Ali::Temporary object from existing Ali object
    use aliased 'Bio::MUST::Core::Ali';
    my $ali = Ali->load('queries.ali');
    my $temp_qu = Temporary->new( seqs => $ali );

    # pass it to external program
    use File::Temp;
    my $query = $temp_qu->filename;
    my $out = File::Temp->new( UNLINK => 0, SUFFIX => '.blastp' );
    system("blastp -query $query -db $db -out $out");
    say "report: $out";

    # later... when parsing the BLAST report
    # let's say $id is a BLAST hit in database.ali
    my $id = 'seq2';
    my $long_id = $temp_db->long_id_for($id);
    say "hit id: $long_id";
    # ...

    # more alternative constructor calls
    # build Ali::Temporary object from list of Seq objects
    my @seqs = $ali->filter_seqs( sub { $_->seq_len >= 500 } );
    my $temp_ls = Temporary->new( seqs => \@seqs );

    # build Ali::Temporary object preserving gaps in Seq objects
    my $temp_gp = Temporary->new( seqs => \@seqs, args => { degap => 0 } );

DESCRIPTION

This module implements a class representing a temporary FASTA file where sequence ids are automatically abbreviated (seq1, seq2...) for maximum compatibility with external programs. To this end, it combines an internal Bio::MUST::Core::Ali object and a Bio::MUST::Core::IdMapper object.

An Ali::Temporary can be built from an existing ALI (or FASTA) file or on-the-fly from a list (ArrayRef) of Bio::MUST::Core::Seq objects (see the SYNOPSIS for examples).

Its sequences can be aligned or not but by default sequences are degapped before writing the associated temporary FASTA file. If gaps are to be preserved, this behavior can be altered via the optional args attribute.

ATTRIBUTES

seqs

Bio::MUST::Core::Ali object (required)

This required attribute contains the Bio::MUST::Core::Seq objects that are written in the associated temporary FASTA file. It can be specified either as a path to an ALI/FASTA file or as an Ali object or as an ArrayRef of Seq objects (see the SYNOPSIS for examples).

It provides the following methods: all_seq_ids, has_uniq_ids, is_protein, is_aligned, get_seq, get_seq_with_id, all_seqs, filter_seqs and count_seqs (see Bio::MUST::Core::Ali).

args

HashRef (optional)

When specified this optional attribute is passed as is to the temp_fasta method of the internal Ali object. Its purpose is to allow the fine-tuning of the format of the associated temporary FASTA file.

By default, its contents is clean = 1> and degap = 1>, so as to generate a FASTA file of deggaped sequences where ambiguous and missing states are replaced by X.

file

Path::Class::File object (auto)

This attribute is automatically initialized with the path of the associated temporary FASTA file. Thus, it cannot be user-specified.

It provides the following methods: remove and filename (see below).

mapper

Bio::MUST::Core::IdMapper object (auto)

This attribute is automatically initialized with the mapper associating the long ids of the internal Ali object to the abbreviated ids used in the associated temporary FASTA file. Thus, it cannot be user-specified.

It provides the following methods: all_long_ids, all_abbr_ids, long_id_for and abbr_id_for (see Bio::MUST::Core::IdMapper).

ACCESSORS

filename

Returns the stringified filename of the associated temporary FASTA file.

This method does not accept any arguments.

type

Returns the type of the sequences in the internal Ali object using BLAST denomination (prot or nucl). See Bio::MUST::Core::Seq::is_protein for the exact test performed.

This method does not accept any arguments.

MISC METHODS

remove

Remove (unlink) the associated temporary FASTA file. Since this method is automatically invoked at object destruction, users should not need it.

AUTHOR

Denis BAURAIN <denis.baurain@uliege.be>

COPYRIGHT AND LICENSE

This software is copyright (c) 2013 by University of Liege / Unit of Eukaryotic Phylogenomics / Denis BAURAIN.

This is free software; you can redistribute it and/or modify it under the same terms as the Perl 5 programming language system itself.

cpanm Bio::MUST::Core

perl -MCPAN -e shell
install Bio::MUST::Core