Releases: brentp/slivar
tsv fix
v0.2.1
-
add
INFO.highest_impact_order
so user can do things like:--info "INFO.highest_impact_order == ImpactOrder.missense"
to get only missense mutations. or
--info "INFO.highest_impact_order <= ImpactOrder.missense"
to get variants with an impact of missense or higher. (note that lower values have higher impact)
-
tsv
: bugfix!! previously,slivar tsv
used a< 0
check where it should hvae used<= 0
so
candidates with a parental sample at index 0 in the samples of the vcf would have missing information
for depth, GQ, and AB. (#78 . Thanks @amwenger for finding the problem and its source). -
expr
: always annotate variant and INFO even if no sample (or trio,group,fam) expressions are given so one can use slivar just to get the 'impactful' annotation. -
remove duplicates from default-order.txt
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing. (but you'll likely want the javascript included in this repo and the gnotation files linked below)
users can also use slivar via docker at brentp/slivar:v0.2.1
The pslivar binary allows running slivar expr commands in parallel
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
minor bugfix
v0.2.0
- clear values from expressions for family, group, and sample. was previously only done for
trios. so if you relied on something like:then the expr2 would still see expr1 in INFO even if the preceding expression did not pass but it--group-expr "expr1: ..." \ --group-expr "expr2:('expr1' in INFO) && ..."\
has passed for a previous variant. Thanks Amelia Wallace for finding and reporting a simple test-case. - add /opt/slivar/slivar-functions.js to docker image (#75)
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing. (but you'll likely want the javascript included in this repo and the gnotation files linked below)
users can also use slivar via docker at brentp/slivar:v0.2.0
The pslivar binary allows running slivar expr commands in parallel
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
bugfix for vcfs with >256 info fields
v0.1.13
- fix bug with VCFs with >256 info fields (#74). Thanks @liserjrqlxue for making a great test-case.
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing. (but you'll likely want the javascript included in this repo and the gnotation files linked below)
users can also use slivar via docker at brentp/slivar:v0.1.13
The pslivar binary allows running slivar expr commands in parallel
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
features, bugfixes, compound-het fixes
This adds small improvements except for large bug-fix for compound-hets for families with multiple affected kids. Previously, a shared compound-het would only be reported for one of the kids. This now fixed. See below for full changelog.
v0.1.12
- slivar expr: handle ploidy > 2 (updated hts-nim) (#55 thanks to @markw3lsh for reporting and providing a test-case)
- slivar: optionally support strings in format (sample) fields. default is to ignore these for performance reasons, but users can set the environment variable
SLIVAR_FORMAT_STRINGS
to anything to forceslivar
to populate the javascript objects with any strings in the VCF. e.g.SLIVAR_FORMAT_STRINGS=1 slivar expr ...
and useunset SLIVAR_FORMAT_STRINGS
to remove. (#37) - slivar compound-hets: fix bug with 2 affect sibs sharing same compound-het previously only 1 sib would be reported
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing. (but you'll likely want the javascript included in this repo and the gnotation files linked below)
users can also use slivar via docker at brentp/slivar:v0.1.12
The pslivar binary allows running slivar expr commands in parallel
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
parallel slivar and INFO.genic
This release is the first where pslivar (parallel slivar) is well tested. It handles running slivar expr in parallel (by genomic region) see wiki page for more detail.
This also adds INFO.genic
for VCFs that have been annotated with VEP, snpEff, or bcftools. This complements INFO.impactful
, but mainly adds synonymous
variants to consideration. See the wiki for more detail.
Full changes below:
v0.1.11
- gnotate: ~10-20% speed improvement by inlining cmp in binarysearch
- slivar expr: fix counts in summary table when some families in the cohort have no affected samples
- slivar tsv: improve readability when 1 or both parents are missing
- slivar tsv: fix comma separation with 1 or both parents missing (thanks @brwnj for reporting)
- slivar ddc: when looking at all chromosomes, skip X, Y
- slivar expr: add --exclude option which takes a bed file of exclude regions
- slivar compound-hets: default will now output compound-het pairs that include synonymous variants
- pslivar: fix tsv output from pslivar (which allows running slivar in parallel). more info here
- add
INFO.genic
boolean to complementINFO.impactful
. By default this includes all impacts included byimpactful
along with
synonymous and other exonic, but non-protein-altering variants (but does not includeintronic
. (more info here
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing. (but you'll likely want the javascript included in this repo and the gnotation files linked below)
users can also use slivar via docker at brentp/slivar:v0.1.11
The pslivar binary allows running slivar expr commands in parallel
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
fixes for data-driven cutoffs.
v0.1.10
- [slivar ddc] fix bug when selecting sample filters.
- [slivar ddc] UI fix: remove button when [x] is clicked.
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing.
users can also use slivar via docker at brentp/slivar:v0.1.10
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
data-driven-cutoffs and family expressions for rare disease
v0.1.9
- new command
slivar ddc
: https://github.com/brentp/slivar/wiki/data-driven-cutoffs - fix for slivar compound-hets reporting orphaned variants in some cases (thanks Steve Boyden for reporting)
- update to htslib >= 1.10
- change recommended setup for rare disease to use family-based approach, rather than trios as it is more
flexible and now well-tested. (https://github.com/brentp/slivar/wiki/rare-disease#full-analysis-for-trios-with-unaffected-parents) - add
VCF
object available in javascript which gets populated withVCF.CSQ == ["CONSEQUENCE", "CODONS","AMINO_ACIDS", "GENE", ...]
and VCF.ANN, VCF.BCSQ if any
or all of those annotations are available in the VCF (see #3)
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing.
users can also use slivar via docker at brentp/slivar:v0.1.9
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
bug-fix for subsetting vcfs
v0.1.8
- fix bug when subsetting vcfs (ped file contained subset of samples in vcf). Thanks Matt for reporting.
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing.
users can also use slivar via docker at brentp/slivar:v0.1.8
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
-
gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
-
gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
-
gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.8) here
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
bug-fixes and gnomad genomes v3
this release has small bugfixes and general usability experiments. see below for the list of changes and the links to the gnotate (zip) file for hg38 genomes from gnomad.
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing.
users can also use slivar via docker at brentp/slivar:v0.1.7
Full release notes:
v0.1.7
- slivar make-gnotate will error with warning if field is not a float or int
- compound-het provides a way to ignore some impacts (intron, non_coding, etc) and sets sane defaults. this removes
most candidates from WGS which were predominantly pairs of intronic or non-coding variants. - compound-het: remove -f flag, this is now discovered by default
- only samples with parents in the vcf (not just ped) are counted as trios
- expose sample.phenotype attribute which is a string taken directly from the pedigree file
- given an error on an INFO expression, don't bail, just report and continue
- fix for VCFs without index AND without contigs in header (#44)
- support for file sizes > 4.2GB and support for gnotate files with chr prefix annotating query files without
- fix gnotate bug where some long alleles were not annotated
gnotate annotation files
the gnotation files for fast annotation remain unchanged except for the addition of topmed.
Users can create their own gnotation files with slivar make-gnotate
, but we provide:
- gnomad for hg37 with AF popmax, numhomalts (total and controls only) here
- gnomad for hg38 with AF popmax, numhomalts (updated in release v0.1.2) here
+ gnomad genomes (71,702 samples) for hg38 with AF popmax, numhomalts (updated in release v0.1.7) here
This zip file is broken. Will be updated soon.
-
spliceai scores (maximum value of the 4 scores in spliceai) here
-
topmed allele frequencies (via dbsnp) these can be used with
INFO.topmed_af
. Useful when analyzing data in hg38 because some variants in hg38 are not visible in GRCh37
family-based filtering and `impactful`
until this release slivar
was heavily focused on trio expressions. but this left quads and more ragged families as possible but more difficult. this release allows things like:
--family-expr 'aff_only:fam.every(function(s) { s.het == s.affected && s.hom_ref == !s.affected && s.GQ > 15 })'
to extract variants where all affecteds in a family are heterozygotes and all unaffecteds are homozygous reference (and all samples have genotype quality > 15).
This release also introduces INFO.impactul
, which is extracted from the consequence field added by VEP, bcftools, or snpEff (or all of them if they are all present). This can be used to exclude, for example synonymous or non-coding or other variants as a way to quickly get candidate variants that are more likely functional. See the wiki for more details.
Installation
Just grab the binary
You can download it and use it without any other software. This is the recommended binary.
wget, chmod +x and start analyzing.
users can also use slivar via docker at brentp/slivar:v0.1.6
Full release notes:
v0.1.6
- add INFO.impactful (boolean) if CSQ/BCSQ/ANN is present. this value will be true of any of those annotation are high enough
- the
impactful
flag is automaticaly added to the output VCF for any variant passing through slivar and meeting the criteria (https://github.com/brentp/slivar/wiki/impactful) - better checks for length of AD field
- fix for empty groups (#38)
- add --family-expr for family expressions like:
fam.every(function(s) { s.het == s.affected && s.hom_ref == !s.affected && s.GQ > 5 })
- slivar tsv now outputs a sortable column for highest-impact. it makes use of a default
list of impact orderings from: https://uswest.ensembl.org/info/genome/variation/prediction/predicted_data.html and
supplemented with any type seen in bcftools or snpEff. - bug fix with gnotate (-g) when annotating a file with a empty chromsome followed by a non-empty (e.g. if
chr16_random
is mixed in with canonical chromosomes)