BioPerl to convert between sequence formats from Fasta to Genbank

Anjana — Wed, 11 May 2016 09:49:04 -0500

#!/usr/local/bin/perl -w

# Sequence formats to choose: Fasta, EMBL. GenBank, Swissprot, PIR and GCG

use Bio::SeqIO;

$inFile = "BRCA2.fa";

$in  = Bio::SeqIO->newFh('-file' => "$inFile" ,
                           '-format' => 'Fasta');
$out = Bio::SeqIO->newFh('-format' => 'Genbank');
print $out $_ while <$in>;

BOL: BioPerl to convert between sequence formats from Fasta to Genbank

BioPerl to convert between sequence formats from Fasta to Genbank