序列生成_隨機生成序列
阿新 • • 發佈：2021-01-17
技術標籤：序列生成
之前推送過一篇關於用R語言生成的隨機序列，最近在用perl分析巨集基因組的資料，想著是不是可以寫一個隨機生成一段鹼基序列的程式碼。所以就有了下面的程式碼。
#!/usr/bin/perl-w
usestrict;
usewarnings;
useList::Utilqw/shuffle/;
my@seq;
my@newseq;
@seq=qw/A T C G/; ##生成一個含有ATGC四個鹼基的陣列。
my$num="$ARGV[0]"; ##設定一個數字引數，數字代表生成的鹼基數目。
for(my$i=1;$i<=$num;$i++){
@newseq=shuffle@seq;
print"$sequnce";
}
print"\n";
執行perl $0 50 後，將會得到以下含有50個鹼基的序列：
GCTCCTCGGCACCTAAAAGGGGCACGCTGTAGAATCAGACTTCTGGATTA
是不是很不務正業？整天搞這些無聊又沒用的東西。