org.bdgenomics.avocado

genotyping

package genotyping

Visibility
  1. Public
  2. All

Type Members

  1. case class DiscoveredVariant(contigName: String, start: Int, referenceAllele: String, alternateAllele: Option[String]) extends Product with Serializable

    A variant site and alleles.

    A variant site and alleles.

    contigName

    The contig this variant is on.

    start

    The position this variant starts at.

    referenceAllele

    The reference allele this variant varies from.

    alternateAllele

    The substituted allele.

Value Members

  1. object DiscoverVariants extends Serializable with Logging

    Discovers the variants present in a set of aligned reads.

    Discovers the variants present in a set of aligned reads.

    Useful for force-calling variants.

  2. object JointAnnotatorCaller extends Serializable

    Jointly calls variants and computes variant annotations.

  3. object SquareOffReferenceModel

    Squares off a set of genotypes with reference models.

    Squares off a set of genotypes with reference models.

    Many joint genotyping workflows use a "Genome VCF" (gVCF) based approach to incrementally compute genotype likelihoods across their dataset. In this methodology, we generate genotype likelihoods at all positions in all samples. For sites where we do not see evidence of a variant, we compute a "reference model", which is a set of genotype likelihoods assuming that we saw an unknown alternate allele. These likelihoods are then used in a joint genotyping step.

    The alternative to this approach is to discover variants across all samples simultaneously, and to then score these variants. This approach is generally considered too computationally expensive for large cohorts.

    This singleton object "squares off" the reference model by discovering all sites where we called a variant in at least one sample, joining these discovered variants back against the input genotypes, and then excising the genotype likelihoods from the reference models.

  4. object TrioCaller extends Serializable

    Variant caller for genotyping the offspring of two parents.

Ungrouped