Perl script for calculate Levenshtein distance
sub levenshtein_dist { my ($str1, $str2) = @_; my ($len1, $len2) = (length $str1, length $str2); if ($len1 == 0) { return $len2; } if ($len2 == 0) { return $len1; } my %mat; for (my $i = 0; $i2344 days ago
954 days ago
Script to rapid genome clustering based on pairwise ANI
...y -db -outfmt '6 std qlen slen' -max_target_seqs 10000 -o -num_threads 32 Note: using the -perc_identity flag will speed up the search at the cost of sensitivity: blastn -query...639 days ago
Perl script to find edit distance between two sequences !
#!/usr/bin/perl use strict; use warnings; sub edit_distance { my ($s1, $s2) = @_; my $len1 = length($s1); my $len2 = length($s2); my @dp; for (my $i = 0; $i429 days ago