Perl 编程实例(原创)
- 格式:rtf
- 大小:2.78 KB
- 文档页数:3
} } my $GC = $n/$length; if($GC>=0.4 && $GC<=0.6){
$primer{$key}=$value; print OUT2 ">$key\n$primer{$key}\n"; } } } } close IN;
open IN2, "$file2"; open OUT, ">$name";
my %double; while(<IN2>){
chomp; my @output; @output = split; if($output[7]=22 && $double{$output[0]}){
print OUT2 "$output[0]\n"; delete $primer{$output[0]}; }elsif($output[7]=22 && (!($double{$output[0]}))){ $double{$output[0]=1}; } }
close IN; close OUT;
#!/usr/bin/perl -w use strict;
@ARGV != 4 && die "perl program.pl <file1> <file2> <name> <name2>\n"; my $file1=shift @ARGV; my $file2=shift @ARGV; my $name =shift @ARGV; my $name2=shift @ARGV;
}
foreach(sort keys %primer){ print "$primer{$_}\n"; print OUT ">$_\n", "$primer{$_}\n" ;
}
close Il -w use strict;
@ARGV !=2 && die "perl program.pl <> <>\n";
my $file = shift @ARGV; my $name = shift @ARGV;
open IN, "$file" or die; open OUT, ">$name";
while(<IN>){ chomp; s/>//; my $key=$_; my $a= <IN>; my $power=<IN>; chomp $power; if($power =~ m/((..\d\.\d\d)\)$/){ next if (abs($1)>=2); delete $primer{$key}; print OUT ">$key\n"; }
my %primer; my $key; open IN, "$file1"; open OUT2, ">$name2";
while(<IN>){ chomp; if(s/>//){ my $key = $_; my $value = <IN>; chomp $value; unless ($value =~ /A$/i || $value =~ /(G|C)(G|C)(G|C)$/i){ my @value = split(//, $value); my $n; my $length = scalar @value; while($_=shift(@value)){ if(/G|C/i){ $n++;