perl脚本:cds序列转换为pep
Posted erinaceinae
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了perl脚本:cds序列转换为pep相关的知识,希望对你有一定的参考价值。
#!/usr/bin/env perl
use strict;
use warnings;
sub TranslateDNAFile(){
use Bio::SeqIO;
(my $infile,my $outfile)[email protected]_;
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);
以上是关于perl脚本:cds序列转换为pep的主要内容,如果未能解决你的问题,请参考以下文章