1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search→Download
before a call is made. Without the --strict option, reads for individuals below the threshold are ignored.
before a call is made. Without the --strict option, reads for individuals below the threshold are ignored.
--minDepth ''threshold'' Positions where the read depth falls below this threshold will be excluded.
--minDepth ''threshold'' Positions where the read depth falls below this threshold will be excluded.
--maxDepth ''threshold'' Positions where the read depth exceeds this threshold will be excluded.
--maxDepth ''threshold'' Positions where the read depth exceeds this threshold will be excluded.
--hardFilter Filtered positions will be completely absent from output. The default is to use a soft filter, where these
positions are included in output but annotated as failing specific filters.
=== VCF Output ===
--glfAliases ''filename'' By default, GLF filenames are used to label each column in the VCF file. This option allows each filename
to be matched to a more specific individual identifier. The aliases file should include two columns per row,
the first specifying the VCF filename, the second specifying a sample name.
<!--
=== X Chromosome Variant Calling ===
=== X Chromosome Variant Calling ===
--xStart ''sexChromosomeStart'' Start of the non-pseudo-autosomal portion of the X (2,709,521 bp in build 36)
--xStart ''sexChromosomeStart'' Start of the non-pseudo-autosomal portion of the X (2,709,521 bp in build 36)
--xStop ''sexChromosomeEnd'' End of the non-pseudo-autosomal portion of the X (154,584,237 bp in build 36)
--xStop ''sexChromosomeEnd'' End of the non-pseudo-autosomal portion of the X (154,584,237 bp in build 36)
-->
== How It Works ==
For each possible position, glfMultiples considers a series of potential polymorphisms, including transitions and transversions from the reference base, but also bi-allelic polymorphisms where neither of the alleles present in the sample is the reference base. For each potential polymorphism type, the likelihood of the observed bases is maximized with respect to allele frequency. Decisions of which sites are polymorphic take into account the maximized likelihood but also an overall prior for each type of polymorphism (for example, transitions are assumed to account for ~2/3 of all variants).
glfMultiples works with log-likelihoods internally to avoid underflows in samples that may include hundreds or thousands of individuals.
== Download ==
== Download ==
The current version is available for download from [http://www.sph.umich.edu/csg/abecasis/downloads/generic-glfMultiples-2010-05-20.tar.gz here].
The current version is available for download from http://www.sph.umich.edu/csg/abecasis/downloads/generic-glfMultiples-2010-06-16.tar.gz.
== TODO ==
== TODO ==
