Skip to main content
Join the official 2020 Python Developers SurveyStart the survey!

probability of mendelian error in trios

Project description


This package attempts to assign probability to a mendelian error event in a trio.
For example, given a mother with genotype 'C/C' and a father with genotype 'C/C'
a child with genotype 'C/T' will be a "mendelian error", in this case, a candidate
*de novo* mutation.

We can filter candidates based on their genotype likelihoods. For example if the
child had a genotype likelihood of -2,-1,-20 then we are likely to consider this
a genotyping error because the homozygous reference ("C/C" with GL -2) is fairly
close the the GL for het (with GL -1). If the genotype likelihood field was
-20,0,-20, then the call is confidently het. We can use the genotype likelihoods
to assign a probability:


>>> from mendelianerror import mendelian_error
# everyone is homref. this should have a low probability of an error:
>>> father = mother = child = [-0.1, -8.0, -8.0]
>>> mendelian_error(mother, father, child)

# parents are hom, child is het. this is a likely mendelian error:
>>> father = mother = [-0.6, -2.5, -2.5]
>>> child = [-2.5, -0.6, -2.5]
>>> mendelian_error(mother, father, child)


So the input is the 3 GL numbers for each of the father, mother, child.


pip install mendelianerror


After installation, one can use this on a multi-sample VCF file like:
mendelianerror $input.vcf father_id mother_id child_id > $new.vcf

to get only the putative mendelian errors, set a high cutoff:

mendelianerror --cutoff 0.999 $input.vcf father_id mother_id child_id > $new.vcf


+ Only make sense for autosomal variants.
+ Only works on trios (doesn't consider extended pedigrees or siblings).

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for mendelianerror, version 0.0.3
Filename, size File type Python version Upload date Hashes
Filename, size mendelianerror-0.0.3.tar.gz (4.5 kB) File type Source Python version None Upload date Hashes View

Supported by

Pingdom Pingdom Monitoring Google Google Object Storage and Download Analytics Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page