I can't put constraint on fragment length. They could be 1 nt or 2nt ...
• 0 views
•
link
Dear all,
I would like to randomly cleave sequences of length (Ln) until the average length of the resulting fragments is 50nt (+/- 2nt)
I started it with perl, but I've some problems with the average length constraint ...
I wanted to select a random position in [0 ... sequence length] and calculate the length of created segments. But I think it's not the right way to do.
Any suggestions?
One perl solution if there is no overlap:
$Ln=1000;
for ($i=int(rand(50))+1;$i<=$Ln;$i=$i+50) {
printf "%d\n",$i+1-int(rand(3));
}
I'd follow this approach:
I can't put constraint on fragment length. They could be 1 nt or 2nt ...
If you're concerned with performance, use unpack instead of substring.
Log in to answer this question.