Method: Exodb::GeneLocationField#dl_seq!

Defined in:
lib/exodb/datamodel/genelocfield.rb

#dl_seq!Object

Download sequence from web service please use by caution. NCBI will block scamming sequest



66
67
68
69
70
71
72
73
74
75
76
77
# File 'lib/exodb/datamodel/genelocfield.rb', line 66

def dl_seq!
  
  case self.chrrefseq
  when /\Aurn:miriam:refseq:/
    self.sequence = Bio::FastaFormat.new(Bio::NCBI::REST.efetch(self.chrrefseq.split(':', 4), {"db"=>"nucleotide", "rettype"=>"fasta", "retmode"=>"text", "seq_start"=>self.start, "seq_stop"=>self.end})).seq
  else 
    self.sequence = Exodb::Ensembl::REST.sequence_region()
  end
  
  self.save!
  
end