![](https://img-blog.csdnimg.cn/20201014180756923.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Perl
Janton Wang
这个作者很懒,什么都没留下…
展开
-
导入模块实例【PERL】
自定义模块(iprint.pm文件)package iprint;require Exporter;our @ISA = qw(Exporter);our @EXPORT = qw(ihash);sub ihash{ my $h=shift; for my $k(sort keys %$h){ print $k,"\t",$h->{$k},"\n"; }}1;调...原创 2019-11-11 20:09:31 · 493 阅读 · 0 评论 -
统计fasta文件GC含量【Perl】
my $file="myseq.fa";open(FH,$file) or die $!;my $G_num=0;my $C_num=0;my $g_num=0;my $c_num=0;my $bases=0;while(<FH>){ if($_=~/^>/){next;} if(/C/){$C_num+= $_ =~ tr/C/C/;} if(/c/){$c...原创 2019-11-10 19:03:26 · 2379 阅读 · 0 评论 -
统计fasta文件每条序列的长度【Perl】
#!/usr/bin/perl -w use strict;my $file="myseq.fa";open(FH,$file) or die $!;my %hash;my $id=0;while(<FH>){ chomp; my $temp; if($_=~/^>/){ $_=~ s/^>//g; my $temp=$_; next; }...原创 2019-11-10 19:35:09 · 2441 阅读 · 0 评论 -
模拟遗传漂变【PERL】
#!/usr/bin/perl -w use strict;my $DNA="AAAAAAAAAAAAAAAAAAAAAA";my ($i,$mutant);srand(time|$$);$mutant=mutate($DNA);print "\nMutate DNA\n";print "\nThe original DNA:\n$DNA\n";print "\nThe mut...原创 2019-11-10 21:47:39 · 525 阅读 · 0 评论