Mercurial > repos > xuebing > sharplabtool
diff tools/stats/count_gff_features.py @ 0:9071e359b9a3
Uploaded
author | xuebing |
---|---|
date | Fri, 09 Mar 2012 19:37:19 -0500 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tools/stats/count_gff_features.py Fri Mar 09 19:37:19 2012 -0500 @@ -0,0 +1,18 @@ +#!/usr/bin/env python +# This tool takes a gff file as input and counts the number of features in it. + +import sys, fileinput +from galaxy import eggs +from galaxy.datatypes.util.gff_util import GFFReaderWrapper +from bx.intervals.io import GenomicInterval + +# Get args. +input_file = sys.argv[1:] + +# Count features. +count = 0 +for feature in GFFReaderWrapper( fileinput.FileInput( input_file ), fix_strand=True ): + if isinstance( feature, GenomicInterval ): + count += 1 + +print count \ No newline at end of file