a亚洲精品_精品国产91乱码一区二区三区_亚洲精品在线免费观看视频_欧美日韩亚洲国产综合_久久久久久久久久久成人_在线区

首頁(yè) > 編程 > Perl > 正文

perl生成特定堿基比例的隨機(jī)序列的代碼

2020-06-04 20:29:56
字體:
來(lái)源:轉(zhuǎn)載
供稿:網(wǎng)友

方法一(不使用模塊,by agonyr)

代碼如下:
#!/usr/bin/perl -w

use strict;

my @seq = ( "A", "T", "C", "G" );
my $length = 10000;

undef my %hash;
$hash{"A"} = int( $length * 0.3 );
$hash{"C"} = int( $length * 0.3 );
$hash{"G"} = int( $length * 0.2 );
$hash{"T"} = int( $length * 0.2 );

my $i = 0;
while ( $i 《 $length ) {
    my $word = $seq[ rand(@seq) ];
    if ( $hash{$word} ) {
        print "$word";
        $i++;
    }
    $hash{$word}--;
}
print "n";

方法二(使用模塊,by yixf)

代碼如下:
#!/usr/bin/perl

use strict;
use warnings;

use BioX::SeqUtils::RandomSequence;

my $randomizer = BioX::SeqUtils::RandomSequence-》new(
    {
        l =》 10000,
        s =》 1,
        y =》 "dna",
        a =》 3,
        c =》 3,
        g =》 2,
        t =》 2
    }
);
print $randomizer-》rand_seq(), "n";

兩種方法比較

設(shè)定長(zhǎng)度為10000,ACGT的比例為3:3:2:2。
代碼如下:
withoutModule Length=10000 GC=49.42% A=2558,C=2503,G=2439,T=2500,Others=0
withModule Length=10000 GC=50.00% A=3000,C=3000,G=2000,T=2000,Others=0

發(fā)表評(píng)論 共有條評(píng)論
用戶名: 密碼:
驗(yàn)證碼: 匿名發(fā)表
主站蜘蛛池模板: 男人天堂中文字幕 | 久久久高清 | 精品欧美一区二区三区久久久 | 狠狠躁夜夜躁人人爽天天高潮 | 精品日韩欧美一区二区三区在线播放 | 99国产精品| 超碰日本| 成人精品一区 | 久久综合一区二区三区 | 日韩欧美一级精品久久 | 91国高清视频 | 男女做爰高清无遮挡免费视频 | 少妇一区二区三区毛片免费下载看 | 久草久| 亚洲 中文 欧美 日韩 在线观看 | 国产精品美女久久久久aⅴ国产馆 | 国产美女在线免费 | 欧美成人一区二区三区片免费 | 黄色福利 | 成人欧美一区二区三区在线播放 | 国产美女久久久 | 免费人成在线观看网站 | 久久电影国产免费久久电影 | 久久久久黄色 | 一区| 国产精品久久久久久久久久 | 亚洲一区亚洲二区 | 好看的一级毛片 | 中文一区二区 | 国产精品美女一区二区三区四区 | 羞羞视频官网 | 精品乱码一区二区 | 国产午夜精品久久久久久久 | theporn国产在线精品 | 成人情趣视频 | 日韩中文字幕一区二区 | 91在线播 | 久久国产精品成人免费观看的软件 | 亚洲一区在线视频 | 最新中文字幕在线观看 | 精品在线免费播放 |