Parser for the Genscan report output.
Returns
Returns C+G content of the query sequence.
Returns Array of Bio::Genscan::Report::Gene.
Returns Genscan version.
Returns
Returns Length of the query sequence.
Returns
Returns Name of query sequence.
Returns Array of Bio::Genscan::Report::Gene.
Returns Array of Bio::Genscan::Report::Gene.
Returns Name of query sequence.
Returns Name of query sequence.
Returns
Parse a Genscan report output string.
# File lib/bio/appl/genscan/report.rb, line 67 def initialize(report) @predictions = [] @genscan_version = nil @date_run = nil @time = nil @query_name = nil @length = nil @gccontent = nil @isochore = nil @matrix = nil report.each_line("\n") do |line| case line when %r^GENSCAN/ parse_headline(line) when %r^Sequence/ parse_sequence(line) when %r^Parameter/ parse_parameter(line) when %r^Predicted genes/ break end end # rests i = report.index(%r^Predicted gene/) j = report.index(%r^Predicted peptide sequence/) # genes/exons genes_region = report[i...j] genes_region.each_line("\n") do |line| if %rInit|Intr|Term|PlyA|Prom|Sngl/ =~ line gn, en = line.strip.split(" +")[0].split(%r\./).map {|i| i.to_i } add_exon(gn, en, line) end end # sequences (peptide|CDS) sequence_region = report[j...report.size] sequence_region.gsub!(%r^Predicted .+?:/, '') sequence_region.gsub!(%r^\s*$/, '') sequence_region.split(Bio::FastaFormat::RS).each do |ff| add_seq(Bio::FastaFormat.new(ff)) end end