Re: Downloading genbank (full) format
Peter Cock <[email protected]> Wed, 23 Sep 2015 08:36:56 +0100
| Newsgroups | gmane.comp.lang.perl.bio.general |
|---|---|
| Message-ID | <CAKVJ-_7XsKmReG3wX_WvNkMnN3HEjc29vgRhi=m2ieVqBJ-H+Q@mail.gmail.com> |
Hi Cacau, It seems "gb" is just an alias for "genbank", but as you've noticed some records are more like manifests for how to build up the full record from its parts using a CONTIG entry. There is also "gbwithparts" which in the NCBI web interface is "Genbank (full)", but there can be subtle errors with the processing done at the NCBI end. http://blastedbio.blogspot.co.uk/2012/03/missing-external-exons-in-genbank-with.html http://blastedbio.blogspot.co.uk/2012/04/missing-feature-locations-in-genbank.html Another possible but not yet confirmed issue: http://lists.open-bio.org/pipermail/biopython/2015-September/015746.html Peter On Wed, Sep 23, 2015 at 5:57 AM, Cacau Centurion <[email protected]> wrote: > Hi, > > I would like to download sequences in genbank (full) format in batch. What > should the format be (see the codes)? > > I tried 'gb' but got sequences in genbank format. Sometimes the full > genome sequences might not be downloaded. > > > #code > $seqin = Bio::SeqIO->new(-file => $out, > -format => $format, > ); > > > Yours, > Cacau > > _______________________________________________ > Bioperl-l mailing list > [email protected] > http://mailman.open-bio.org/mailman/listinfo/bioperl-l > _______________________________________________ Bioperl-l mailing list [email protected] http://mailman.open-bio.org/mailman/listinfo/bioperl-l