cf37e57c66f855418bc05aff039ab846ad317354 hiram Mon Mar 2 13:20:42 2020 -0800 correct one line header message in usage() refs #25071 diff --git src/hg/utils/vai.pl src/hg/utils/vai.pl index abfb991..cdab417 100755 --- src/hg/utils/vai.pl +++ src/hg/utils/vai.pl @@ -65,32 +65,34 @@ 'Include "variants" with no observed variation in output'], variantLimit => [10000, '=N', 'Maximum number of variants to process'], ); # CGI params and values to be passed to hgVai are collected here: my %hgVaiParams = ( $doQueryParam => 'go', ); sub usage($) { # Show usage and exit with the given status. my ($status) = @_; print STDERR <<EOF +vai.pl - Invokes hgVai (Variant Annotation Integrator) on a set of variant calls to add functional effect predictions and other data relevant to function. -usage: vai.pl [options] db input.(vcf|pgsnp|pgSnp|txt)[.gz] > output.tab +usage: + vai.pl [options] db input.(vcf|pgsnp|pgSnp|txt)[.gz] > output.tab Invokes hgVai (Variant Annotation Integrator) on a set of variant calls to add functional effect predictions (e.g. does the variant fall within a regulatory region or part of a gene) and other data relevant to function. input.(...) must be a file or URL containing either variants formatted as VCF or pgSnp, or a sequence of dbSNP rs# IDs, optionally compressed by gzip. Output is printed to stdout. options: --hgVai=/path/to/hgVai Path to hgVai executable (default: $hgVai) --position=chrX:N-M Sequence name, start and end of range to query (default: genome-wide query) --rsId Attempt to match dbSNP rs# ID with variant