Alien-SeqAlignment-edlib
view release on metacpan or search on metacpan
USAGE
Command line tool
use v5.38;
use Alien::SeqAlignment::edlib;
use Env qw( @PATH );
unshift @PATH, Alien::SeqAlignment::edlib->bin_dir;
my $string1 = "ACGACG";
my $string2 = "CCCCCACGTCG";
# save sequences
open my $fh, '>', 'seq1.fasta';
say $fh ">Seq1\n$string1";
close $fh;
open my $fh, '>', 'seq2.fasta';
say $fh ">Seq2\n$string2";
close $fh;
system Alien::SeqAlignment::edlib->edlib_aligner, '-m', 'HW','-n','0', '-k','-1','-p','-f' ,'NICE','seq1.fasta', 'seq2.fasta';
Output
Using HW alignment mode.
Reading queries...
Read 1 queries, 6 residues total.
Reading target fasta file...
Read target, 11 residues.
Comparing queries to target...
Query #0 (6 residues): score = 1
T: ACGTCG (5 - 10)
||| ||
Q: ACGACG (0 - 5)
SEE ALSO
* edlib <https://github.com/Martinsos/edlib>
Edlib is a lightweight and superfast C/C++ library for sequence
alignment using the edit (Levenshtein) distance between two or more
biological (usually) sequences. It can calculate the edit distance,
find the optimal aligment path and the coordinates (start/end)
locations. It supports multiple alignment modes such as global (NW),
prefix (SHW) and infix (HW). The library does not handle utf8 and
its primary use is to compute edit distances and alignments over
small (255 characters or fewer) alphabets as they occur in
bioinformatic applications.
* Text::Levenshtein::Edlib
<https://metacpan.org/pod/Text::Levenshtein::Edlib>
An XS library that also wraps around the edlib library and returns
edit distances, as well as alignment paths.
* Text::Levenshtein::XS
<https://metacpan.org/pod/Text::Levenshtein::XS>
An XS library that computes edit distances but not alignment paths.
See also its github repository at:
<https://github.com/ugexe/Text--Levenshtein--XS/>)
* Text::LevenshteinXS <https://metacpan.org/pod/Text::LevenshteinXS>
Yet another XS implementation of Levenshtein distance over strings
(no alignment path).
* Alien
Documentation on the Alien concept itself.
* Alien::Base <https://metacpan.org/pod/Alien::Base>
The base class for this Alien. The methods in that class allow you
to use the static and the dynamic edlib library in your code.
* Alien::Build::Manual::AlienUser
<https://metacpan.org/dist/Alien-Build/view/lib/Alien/Build/Manual/A
lienUser.pod>
Detailed manual for users of Alien classes.
* Bio::SeqAlignment <https://metacpan.org/pod/Bio::SeqAlignment>
A collection of tools and libraries for aligning biological
sequences from within Perl.
AUTHOR
Christos Argyropoulos <chrisarg@gmail.com>
COPYRIGHT AND LICENSE
This software is copyright (c) 2023 by Christos Argyropoulos.
This is free software; you can redistribute it and/or modify it under
the same terms as the Perl 5 programming language system itself.
( run in 0.649 second using v1.01-cache-2.11-cpan-df04353d9ac )