Skip to content

Instantly share code, notes, and snippets.

@helios

helios/gist:2658172

Created May 11, 2012
Embed
What would you like to do?
Parsing Cufflinks Gtf file/data, in an ActiveRecord style
g=Bio::Ngs::Cufflinks::Gtf.new "./spec/fixture/transcripts.gtf"
g.set_lazy
a=g.multi_exon.coverage_gt(3).length_gt(200)
#not yet parsed until you iterate of the transcripts
a.multi_exons.each_transcript {|t| puts t; break}
#accessing to an exon using the indexin:
g.read_transcript(10) #return a string of the block
g.get_tranascript(10) #return an object transcript.
@helios

This comment has been minimized.

Copy link
Owner Author

@helios helios commented May 11, 2012

also g[] works.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment