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

BioScripts
Anjana
BioPerl to convert between sequence formats from Fasta to Genbank

BioPerl to convert between sequence formats from Fasta to Genbank

By Anjana 3121 days ago

#!/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>;