Alien-SeqAlignment-edlib

 view release on metacpan or  search on metacpan

README  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...

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 )