\n Here’s a handy one(ish)-liner to mutate an input sequence using Perl’s RegEx engine:\n<\/p>\n
epiphyte:~ rmp$ perl -e '$seq=\"ACTAGCTACGACTAGCATCGACT\"; $mutants = [qw(A C T G)];\r\n print \"$seq\\n\";\r\n $seq =~ s{([ATGC])}{ rand() < 0.5 ? $mutants->[int rand 4] : $1 }smiexg;\r\n print \"$seq\\n\";'\r\nACTAGCTACGACTAGCATCGACT\r\nACAATCGCGGACCAGAATCTCTT<\/code><\/pre>\n\nThis gives each base in $seq a 50% chance (rand() < 0.5) of mutating to something, but as the original base is in the available $mutants array it has a further 25% chance of changing to itself. If you wanted to improve it by excluding the original base for each mutation you might do something like:\n<\/p>\n
epiphyte:~ rmp$ perl -e '$seq=\"ACTAGCTACGACTAGCATCGACT\"; $mutants = [qw(A C T G)];\r\n $mutsize=scalar @{$mutants}; print \"$seq\\n\";\r\n $seq =~ s{([ATGC])}{ rand() < 0.5 ? [grep { $_ ne $1 } @{$mutants}]->[int rand $mutsize-1] : $1 }smiexg;\r\n print \"$seq\\n\";'\r\nACTAGCTACGACTAGCATCGACT\r\nTGTAGATAATGTGATACGAGACT<\/code><\/pre>\n<\/p>\n\nThis (quite inefficiently) builds an array of all available options from $mutants, excluding $1 the matched base at each position.<\/p>\n
Unrolling it and tidying it up a little for readability looks like this:\n<\/p>\n
my $seq = 'ACTAGCTACGACTAGCATCGACT';\r\nmy $mutants = [qw(A C T G)];\r\nmy $mutsize = scalar @{$mutants};\r\n\r\nprint \"$seq\\n\";\r\n\r\n$seq =~ s{([ATGC])}{\r\n rand() < 0.5\r\n ?\r\n [grep { $_ ne $1 } @{$mutants}]->[int rand $mutsize-1]\r\n :\r\n $1\r\n }smiexg;\r\n\r\nprint \"$seq\\n\";'<\/code><\/pre>\n","protected":false},"excerpt":{"rendered":"Here’s a handy one(ish)-liner to mutate an input sequence using Perl’s RegEx engine: epiphyte:~ rmp$ perl -e ‘$seq=”ACTAGCTACGACTAGCATCGACT”; $mutants = [qw(A C T G)]; print “$seq\\n”; $seq =~ s{([ATGC])}{ rand() < 0.5 ? $mutants->[int rand 4] : $1 }smiexg; print “$seq\\n”;’ ACTAGCTACGACTAGCATCGACT ACAATCGCGGACCAGAATCTCTT This gives each base in $seq a 50% chance (rand() < 0.5) … Continue reading “Random Sequence Mutator”<\/span><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"om_disable_all_campaigns":false,"_monsterinsights_skip_tracking":false,"_uf_show_specific_survey":0,"_uf_disable_surveys":false,"footnotes":""},"categories":[11],"tags":[934,935,778],"class_list":["post-744","post","type-post","status-publish","format-standard","hentry","category-programming","tag-mutation","tag-random","tag-sequence"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t