perl程序实现对fasta序列反向互补
发布时间
阅读量：
阅读量
#!/usr/bin/perl
use strict;
use warnings;
open IN,"$ARGV[0]" or die "Can't open fasta file!";
while (<IN>) {
chomp;
if (/>\S+/) {
my $id = $_;
print "$id\n";
}else{
my $seq = $_;
my $rseq = rev_and_com($seq);
print "$rseq";
}
}
sub rev_and_com {
my $s ="";
my $a ="";
$s = shift;
$a = $s;
$a =~ tr/atcgATCG/tagcTAGC/;
return reverse($a);
}
全部评论 (0)
还没有任何评论哟~