Alien-SeqAlignment-edlib
view release on metacpan or search on metacpan
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...
lib/Alien/SeqAlignment/edlib.pm view on Meta::CPAN
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...
( run in 1.566 second using v1.01-cache-2.11-cpan-483215c6ad5 )