bioruby

convert psl to gff

Bio::FlatFile.open(Bio::Blat::Report,f).each{|r| r.hits.each{|h| puts [ h.target.name , source , feature , h.target.start , h.target.end , "." , h.strand , ".", "Sequence " + h.query.name ].join("\t") } }

parse GFF files

Bio::FlatFile.open(Bio::GFF,f).each{|gff| gff.records.each{|r| p r } }

parse BLAT results (psl file)

Bio::FlatFile.open(Bio::Blat::Report,f).each{|r| r.hits.each{|h| p h } }

parse FASTA formatted sequence file

Bio::FlatFile.open(Bio::FastaFormat,f).each{|e| p e }