Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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="dna.fasta";
- my $pepfile="pep.fasta";
- &TranslateDNAFile($DNAfile,$pepfile);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement