def test_Xlinked_Norine_BMX(self):
"""This gene has no X-linked variant"""
db = 'lgueneau'
ss = samples_selection_factory(db=db,
groups = {'affected': ['09818','09819'], 'not_affected':['09960','09961']})
qs = Variant.objects.using(db).filter(gene_symbol__in=['BMX'], chrom='chrX')
# Using GenotypesFilterXLinked.apply()
variants = variants_collection_factory(db=db, qs=qs)
var = GenotypesFilterXLinked(ss).apply(variants).variants
self.assertEqual(len(var), 0)
# Using FilterCollection.apply()
fc = FiltersCollection([GenotypesFilterXLinked(ss)])
var = fc.apply(initqs=qs, db=db).variants
self.assertEqual(len(var), 0)
############################
# COMPOUND HET #
############################
#@unittest.skip('')
评论列表
文章目录