perl脚本:cds序列转换为pep

#!/usr/bin/env perl
use strict;
use warnings;

sub TranslateDNAFile(){
use Bio::SeqIO;
(my $infile,my $outfile)=@_;
  my $in=Bio::SeqIO->new(-file=>"$infile",-format=>"fasta");
  my $out=Bio::SeqIO->new(-file=>">$outfile", -format=>"fasta");
  while (my $seq=$in->next_seq()){
  $out->write_seq($seq->translate);
}
}
my $DNAfile="in_cds";
my $pepfile="out_pep";
&TranslateDNAFile($DNAfile,$pepfile);

posted @ 2018-08-12 09:28  萃英山下  阅读(598)  评论(0编辑  收藏  举报