from rgt.AnnotationSet import *
from rgt.Util import GenomeData

# genome = "hg19"
# print("Checking " + genome)
# annot = AnnotationSet(genome,filter_havana=True,protein_coding=True,known_only=True)
# print("\tloading AnnotationSet... succeeds")
# promoters = annot.get_promoters()
# print("\tPromoters "+str(len(promoters)))
# gd = GenomeData(organism=genome)
# print("\t"+gd.get_annotation())
# print("\tloading GenomeData... succeeds")

# genome = "hg38"
# print("Checking " + genome)
# annot = AnnotationSet(genome, filter_havana=False, protein_coding=True, known_only=False)
# # annot = AnnotationSet(genome,filter_havana=True,protein_coding=True,known_only=True)
# print("\tloading AnnotationSet... succeeds")
# promoters = annot.get_promoters()
# print("\tPromoters " + str(len(promoters)))
# gd = GenomeData(organism=genome)
# print("\t" + gd.get_annotation())
# print("\tloading GenomeData... succeeds")

# genome = "mm9"
# print("Checking " + genome)
# annot = AnnotationSet(genome,filter_havana=True,protein_coding=True,known_only=True)
# print("\tloading AnnotationSet... succeeds")
# promoters = annot.get_promoters()
# print("\tPromoters "+str(len(promoters)))
# gd = GenomeData(organism=genome)
# print("\t"+gd.get_annotation())
# print("\tloading GenomeData... succeeds")

# genome = "zv9"
# print("Checking " + genome)
# annot = AnnotationSet(genome,filter_havana=True,protein_coding=True,known_only=True)
# print("\tloading AnnotationSet... succeeds")
# promoters = annot.get_promoters()
# print("\tPromoters "+str(len(promoters)))
# gd = GenomeData(organism=genome)
# print("\t"+gd.get_annotation())
# print("\tloading GenomeData... succeeds")

# genome = "zv10"
# print("Checking " + genome)
# annot = AnnotationSet(genome,filter_havana=True,protein_coding=True,known_only=True)
# print("\tloading AnnotationSet... succeeds")
# promoters = annot.get_promoters()
# print("\tPromoters "+str(len(promoters)))
# gd = GenomeData(organism=genome)
# print("\t"+gd.get_annotation())
# print("\tloading GenomeData... succeeds")
