Sha256: a1b47185b92d28904a2b82f12186846b7ba7facdb2ce6c60622ce215aebea1e7
Contents?: true
Size: 508 Bytes
Versions: 6
Compression:
Stored size: 508 Bytes
Contents
#!/usr/bin/env ruby
#require 'bio-ucsc'
require '../lib/bio-ucsc'
DB = Bio::Ucsc::Hg19
DB.connect
genes = Array.new
ARGF.each_line do |row|
row.chomp!
next if row.empty?
chr, chr_start, chr_end = row.split("\t")
gi = GenomicInterval.zero_based(chr,
Integer(chr_start),
Integer(chr_end),)
results = DB.with_interval(gi).select(:name2).find(:all)
genes.concat(results.map{|e|e.name2})
end
genes.uniq.each{|e|puts e} if genes
Version data entries
6 entries across 6 versions & 1 rubygems