From Genome Analysis Wiki
Jump to navigationJump to search
453 bytes added
, 12:06, 8 October 2012
Line 98: |
Line 98: |
| == Other options == | | == Other options == |
| | | |
− | Some of command line options are explained below and others are self-explanatory. | + | Some of command line options are explained below and others are self-explanatory. |
| | | |
| -c : minimum cutoff of posterior probability to output a variant [''Default: 0.5''] | | -c : minimum cutoff of posterior probability to output a variant [''Default: 0.5''] |
− | --theta : scaled mutation rate per site [''Default: 0.001''] | + | --theta : scaled mutation rate per site for single nucleotide variants [''Default: 0.001''] |
| + | --indel_theta : scaled mutation rate per site for Indels (works only for VCF input) [''Default: 0.0001''] |
| --tstv: prior of ts:tv ratio [''Default: 2.0''] | | --tstv: prior of ts:tv ratio [''Default: 2.0''] |
| --nthreads : number of threads to run and it is recommended to use 4 threads for small number of input files [''Default: 1'] | | --nthreads : number of threads to run and it is recommended to use 4 threads for small number of input files [''Default: 1'] |
| | | |
| + | * The following applies only to GLF input |
| --denovo : a boolean flag to turn on ''de novo'' mutation detection. The following options take effect only when this flag is ON | | --denovo : a boolean flag to turn on ''de novo'' mutation detection. The following options take effect only when this flag is ON |
| --rate_denovo : mutation rate per haplotype per generation. [''Default: 1.5e-08''] | | --rate_denovo : mutation rate per haplotype per generation. [''Default: 1.5e-08''] |
Line 110: |
Line 112: |
| --minLLR_denovo : minimum value of log10 likelihood ratio of allowing vs. disallowing ''de novo'' mutations in the data to output [''Default: 1.0''] | | --minLLR_denovo : minimum value of log10 likelihood ratio of allowing vs. disallowing ''de novo'' mutations in the data to output [''Default: 1.0''] |
| | | |
− | --chr2process: the chromosome names to process. Default is empty and is to process all chromosomes in the input. | + | --pos : a file with two columns (chr pos) to output genotypes of all individuals, even if the sites are monomorphic |
− | If multiple chromosomes are provided, they should be separated by comma, e.g. --chr2process 2,10 or --chr2process chr2,chr10
| + | --all_sites : If turned on, all sites with at least one read coverage will be output. |
− | | + | --gl_off : If turned on, not to output genotype likelihood values for each individual. Default is to output 3 GLs for polymorphisms and 10 GLs for de novo mutations |
− | --gl_off: not to output genotype likelihood values for each individual. Default is to output 3 GLs for polymorphisms and 10 GLs for de novo mutations | + | --quick_call : If turned on, it will perform variant calling assuming that all individuals are unrelated, and if a site is detected as a variant site then the family-aware variant calling will be performed. This will be beneficial for complex pedigrees for which the likelihood calculation may be demanding. |
| | | |
| == Output files == | | == Output files == |