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的主要内容,如果未能解决你的问题,请参考以下文章

从细菌GFF文件提取CDS序列并转换为氨基酸序列

python版的MCScan绘图

利用R语言及perl语言提取基因CDS区实战

perl脚本:读取fasta序列

实用案例精讲!如何用perl写一个截序列的脚本?

sh 从CDS #fasta中提取fasta序列