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
54
55
56
57
|
from .common import eq_
from .common import run_multiple_genomes
from .data import TP53_gene_id
@run_multiple_genomes()
def test_TP53_gene_object_by_id(genome):
# when we look up TP53 by its gene ID, we should get the
# correct gene back
gene = genome.gene_by_id(TP53_gene_id)
assert gene.name == "TP53", "Incorrect gene name %s for gene ID %s in %s" % (
gene.name,
gene.id,
genome,
)
assert gene.contig == "17", "Incorrect gene contig %s for gene ID %s in %s" % (
gene.contig,
gene.id,
genome,
)
@run_multiple_genomes()
def test_TP53_gene_object_by_name(genome):
genes = genome.genes_by_name("TP53")
# we should only have one TP53 gene (there aren't any copies)
assert len(genes) == 1, "Expected only one gene with name TP53, got %s" % (genes,)
# make sure it has the correct gene ID
assert genes[0].id == TP53_gene_id, "Expected gene to have ID %s, got %s" % (
TP53_gene_id,
genes[0].id,
)
@run_multiple_genomes()
def test_equal_genes(genome):
gene1 = genome.genes_by_name("TP53")[0]
# get an identical gene
gene2 = genome.gene_by_id(gene1.id)
assert hash(gene1) == hash(gene2)
assert gene1 == gene2
@run_multiple_genomes()
def test_not_equal_genes(genome):
gene1 = genome.genes_by_name("MUC1")[0]
gene2 = genome.genes_by_name("BRCA1")[0]
assert hash(gene1) != hash(gene2)
assert gene1 != gene2
@run_multiple_genomes()
def test_BRCA1_protein_coding_biotype(genome):
gene = genome.genes_by_name("BRCA1")[0]
assert gene.is_protein_coding
eq_(gene.biotype, "protein_coding")
|
