[Bioperl-l] using Bio::SeqFeature::Tools::Unflattener

[Jason Stajich]

you can try this but it is a bit of an involved script because it is 
setup for dealing with multiple genomes in multiple folders so you might 
want to simplify it.
  http://github.com/hyphaltip/genome-scripts/blob/master/data_format/genbank_gbk2gff3_unflatten.pl

But I thought the perldoc was a good starting point - have you tried it

Generally I do:
  GENBANK -> GFF3  --> genbank_gbk2gff3_unflatten.pl
  GFF3 -> {CDS,PEP,GENE}  --> 
http://github.com/hyphaltip/genome-scripts/blob/master/gene_prediction/gff3_to_cdspep.pl 
(or equivalent)

-jason

Jessica Sun wrote, On 3/22/10 11:48 AM:
> Does any know how to get CDS of the corresponding mRNA accession(NM_) using
> this function?
> *Bio::SeqFeature::Tools::Unflattener
>
> many thanks in advance.
>
> *
>