Pseudocode for Generating Synthetic Tags with Error Rate

[foolishbrat]

Dear all,

I am aware of Metasim program that generates synthetic tags,
but it doesn't seem to give me much control for what I want to do.

What I intend to do is this:

Given these inputs:

 my $init_seq = "AAAAAAAAAA" #length 10 bp 
    my $sub_rate = 0.003;
    my $nof_tags = 1000;
    my @dna = qw( A C G T );

I want to generate:

1. 1. One thousand length-10 tags
2. 2. Substitution rate for each position in a tag is 0.003

Yielding output like:

   AAAAAAAAAA
    AATAACAAAA
    .....
    AAGGAAAAGA # 1000th tags

Is there a compact way to do it in Perl?

I am stuck with this script logic as core:

#!/usr/bin/perl
    
    my $init_seq = "AAAAAAAAAA" #length 10 bp 
    my $sub_rate = 0.003;
    my $nof_tags = 1000;
    my @dna = qw( A C G T );
    
        $i = 0;
        while ($i < length($init_seq)) {
        	$roll = int(rand 4) + 1;       # $roll is now an integer between 1 and 4
        	
        	if ($roll == 1) {$base = A;}
        	elsif ($roll == 2) {$base = T;}
        	elsif ($roll == 3) {$base = C;}
        	elsif ($roll == 4) {$base = G;};
        	
        	print $base;
        }
        continue {
        	$i++;
        }

[andpet]

One problem could be that $i is not increased in the while - loop ? Something like this would make more sense:

while ($i < length($init_seq)) {
...
print $base;
$i++;
}

Maybe your "continue" statement is supposed to do this ?

Andreas