view getdata/get_1_gb.pl @ 0:5b9a38ec4a39 draft default tip

First commit of old repositories
author osiris_phylogenetics <ucsb_phylogenetics@lifesci.ucsb.edu>
date Tue, 11 Mar 2014 12:19:13 -0700
parents
children
line wrap: on
line source

#!/usr/bin/perl -w
use strict;

#use FindBin;
#use lib "$FindBin::Bin/lib";
use Bio::DB::GenBank;
use Bio::SeqIO;

my $accession = $ARGV[0];
my $datatype = $ARGV[1];
my $outtype = $ARGV[2];
my $outfile = $ARGV[3];


        my $qry_string .= $accession."[accession]";
        my $fh = Bio::SeqIO->newFh(-format=>$outtype, -file=>">$outfile");

        my $GBseq;
        my $gb = new Bio::DB::GenBank;
        my $query = Bio::DB::Query::GenBank->new
                (-query   =>$qry_string,
                 -db      =>$datatype);

        my $count;
        my $species;
        my $seqio = $gb->get_Stream_by_query($query);
        while( defined ($GBseq = $seqio->next_seq )) {
                my $sequence = $GBseq;   # read a sequence object
                print $fh $sequence; # write a sequence object
        }

exit;