view tools/stats/count_gff_features.py @ 2:c2a356708570

Uploaded
author xuebing
date Fri, 09 Mar 2012 19:45:42 -0500
parents 9071e359b9a3
children
line wrap: on
line source

#!/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