perl生成特定堿基比例的隨機(jī)序列的代碼
更新時(shí)間：2013年03月01日 11:24:56 作者：
怎么用perl程序，隨機(jī)生成一條序列，使ACGT四種堿基的含量分別為0.3，0.3，0.2，0.2！
方法一（不使用模塊，by agonyr）
復(fù)制代碼 代碼如下:
#!/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）
復(fù)制代碼 代碼如下:
#!/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è)定長度為10000，ACGT的比例為3:3:2:2。
復(fù)制代碼 代碼如下:
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
相關(guān)文章
perl特殊符號(hào)及默認(rèn)的內(nèi)部變量
perl特殊符號(hào)及默認(rèn)的內(nèi)部變量，有需要的朋友不妨參考下2013-02-02
這篇文章主要介紹了Perl學(xué)習(xí)筆記之文件操作,本文分別給出了打開文件、讀取文件、寫入文件代碼實(shí)例,需要的朋友可以參考下2015-06-06
這篇文章主要介紹了Perl腳本實(shí)現(xiàn)遞歸遍歷目錄下的文件,本文直接給出實(shí)現(xiàn)代碼,代碼中包含明細(xì)注釋,需要的朋友可以參考下2015-04-04
在EditPlus中配置Perl開發(fā)編譯環(huán)境
這篇文章主要介紹了在EditPlus中配置Perl開發(fā)編譯環(huán)境,配置起來蠻簡單,適合小型開發(fā)環(huán)境,需要的朋友可以參考下2015-06-06
Perl字符串處理學(xué)習(xí)備忘錄，包括chop和chomp函數(shù)、getc函數(shù)、拼接字符串、分割字符串、重復(fù)拼接字符串、字串替換等。有需要的朋友，可以參考下2013-02-02
強(qiáng)大的Perl正則表達(dá)式實(shí)例詳解
正則表達(dá)式是Perl語言的特色，基本的語法不是很難，但是編寫一個(gè)符合需求、高效的正則表達(dá)式，還是有一些挑戰(zhàn)的。下面這篇文章主要給大家介紹了關(guān)于強(qiáng)大的Perl正則表達(dá)式的相關(guān)資料，需要的朋友可以參考借鑒，下面來一起看看吧。2017-11-11
perl中的字符串操作函數(shù)chomp與chop介紹
這篇文章主要介紹了perl中的字符串操作函數(shù)chomp與chop介紹,chomp與chop都是用于去除字符串變量尾部的字符，但它們有各自的區(qū)別,需要的朋友可以參考下2015-02-02
這篇文章主要介紹了Perl讀寫文件簡單示例,本文直接給出實(shí)現(xiàn)代碼,需要的朋友可以參考下2015-04-04
這篇文章主要介紹了python urllib中的編碼處理示例,本文直接給出幾個(gè)處理實(shí)例,需要的朋友可以參考下2015-04-04
在Perl排名持續(xù)下降的情況下學(xué)Perl，似乎是有點(diǎn)不明智。但是，工作需要，不得不學(xué)啊。再說，Perl現(xiàn)在在測試領(lǐng)域還是用得非常多的。Phython雖然也在測試領(lǐng)域開始活躍起來，不過我始終還是不太喜歡Phython的語法。2009-04-04