-
Notifications
You must be signed in to change notification settings - Fork 30
/
Copy pathtest_AnnotationSet.py
53 lines (48 loc) · 1.99 KB
/
test_AnnotationSet.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
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")