广州明生医药有限公司


perl 指定长度并生成一个随机的DNA序列的脚本代码

网络编程 perl 指定长度并生成一个随机的DNA序列的脚本代码 06-21


#!/bin/perl

use strict;
use warnings;

#进行定义
my @dna;
my $dna_length;
my $newbase;
my $i=0;

print "please input the DNA lengthn";
chomp($dna_length=<>);

while($i<$dna_length)
{
#从四个碱基中随机选取一个
my(@nucleotides)=qw/A T G C/;
$newbase=$nucleotides[rand @nucleotides];

#将随机生成的序列添加到@dna的数组中
push(@dna,$newbase);
++$i;
}
print "@dna";

perl生成特定碱基比例的随机序列的代码
方法一(不使用模块,byagonyr)#!/usr/bin/perl-wusestrict;my@seq=("A","T","C","G");my$length=10000;undefmy%hash;$hash{"A"}=int($length*0.3);$hash{"C"}=int($length*0.3);$hash{"G"}=int($length

求婚示爱的Perl代码之注释篇
原始代码代码出处:http://www.perlmonks.org/index.pl?node_id=384100#!/usr/bin/perl-wusestrict;my$f=$[;my$ch=0;subl{length}subr{join"",reversesplit("",$_[$[])}subss{substr($_[0],$_[1],$_[2])}su

求婚示爱的Perl代码之改写篇
改写代码#!/usr/bin/perl-wusestrict;my$f=$[;my$ch=0;subl{length}subr{join"",reversesplit("",$_[$[])}subss{substr($_[0],$_[1],$_[2])}subbe{$_=$_[0];p(ss($_,$f,1));$f+=l()/2;$f%=l();$f++if$ch%2;$ch++


编辑:广州明生医药有限公司

标签:代码,碱基,序列,数组,注释