Sha256: 152b707c16d2259955941a5f036f31c0a8057b82382085730a1f698fb1ac3daf
Contents?: true
Size: 545 Bytes
Versions: 7
Compression:
Stored size: 545 Bytes
Contents
#!/usr/local/bin/ruby-1.9
#require 'bio-ucsc'
require '../lib/bio-ucsc'
include Bio
Ucsc::Hg19::DBConnection.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 = Ucsc::Hg19::RefGene.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
7 entries across 7 versions & 1 rubygems