def main():
""" prints the description of ##FORMAT metalines in a VCF """
usage = "usage: %prog [options] file.vcf"
parser = OptionParser(usage)
parser.add_option("--formatag", type="string", dest="formatag", help="prints the description for the FORMAT id formatag")
parser.add_option("--all", action="store_true", dest="all", help="prints the description for *every* FORMAT tag in VCF")
parser.add_option
(options, args)=parser.parse_args()
vcfilename=args[0]
vcfh=open(vcfilename, 'r')
#instantiate a VcfFile object
vcfobj=VcfFile(vcfilename)
#parse its metainfo lines (ones that begin with ##)
vcfobj.parseMetaLines(vcfh)
descriptors = vcfobj.getMetaFormatDescription()
if len(descriptors) == 0:
sys.stderr.write("No ##FORMAT lines in vcf!\n")
return
found_tag=0
for (id, description) in descriptors:
if options.all==True:
print id, description
found_tag=1
continue
if id == options.formatag:
print id, "\t", description
found_tag=1
if found_tag ==0 : sys.stderr.write(options.formatag + " not in ##FORMAT headers\n")
