Comments (5)
Hello,
It took a while, but I finally found time to build this. The next release will have a collapse
flag for--median-info
. It adds three new info fields:
##INFO=<ID=CollapseStart,Number=1,Type=Integer,Description="Median start position of collapsed variants">
##INFO=<ID=CollapseEnd,Number=1,Type=Integer,Description="Median end position of collapsed variants">
##INFO=<ID=CollapseSize,Number=1,Type=Integer,Description="Median size of collapsed variants">
Note that I didn't test this extensively. So let me know how it goes.
from truvari.
Hello,
I'm happy to hear you find Truvari useful.
This is an interesting request. My first reaction is to say no, we won't do that exactly as described. Because truvari doesn't use any sequencing information, I don't feel like it really has the authority to change variant representations. And if we opened the door to changing variants, we'd maybe want to consider options like weighing the median breakpoints based on variant depth or AF. How that information is stored in a VCF can get really messy and I can imagine it becomes a huge development burden. Additionally, one of my main design goals is to have inputs/outputs be cleanly traceable (e.g. every variant representation given to collapse ends up either in the output VCF or the collapsed (i.e. removed) output VCF). So, if Truvari creates brand new variant representations, it gets a little tricky to figure out how to make that traceable - not to mention making it cleanly compatible across all the variant representations possible (e.g. what would the 'median' insertions's ALT sequence be?). Everyone's VCFs look a little different and I don't want to add more expectations on what their VCFs look like nor do I want to add more noise to that mix.
Having said that, I do think it may be reasonable to collect the median information and supplement the output VCF's INFO fields with it. Then, if users are interested in changing their variants, Truvari can help them work it out while taking on none of the variant representation responsibility. An initial design would be to add something like INFO/Median[Start/End/Len]
which could be calculated beside INFO/NumCollapsed
in the code.
Please let me know what you think of this idea of adding Median annotations.
Have a great day,
~/Adam English
from truvari.
Thank you so much for your detailed answer, that makes sense to me!
Indeed, adding INFO/Median[Start/End/Len]
could be really helpful to explore data from large cohorts.
May I request to also add the minimum start and maximum end in the cluster INFO/[MinStart/MaxEnd]
?
Thank you very much for Truvari,
Have a great week-end,
Véronique
from truvari.
Sounds good. I'm going to use this ticket to track the development so feel free to unsubscribe to the notifications.
I was close to cutting a v4.0.0 release with a bunch of changes so I'll use the weekend to figure out if I want to get this in for 4.0 or the inevitable 4.1.
Dev notes:
- add header definitions of new fields in
truvari.collapse.edit_header
. - populate the INFO fields inside
truvari.collapse.annotate_entry
- clean up the interface for
annotate_entry
since passing a parameter per-INFO isn't ideal - make a new method
collapse_boundaries
orcalculate_position
that takes the kept entry and the set of collapsed entries to return the Start/End/Len ranges which will need to be passed toannotate_entry
. - can call the new method inside
truvari.collapse.output_writer
- this change will break probably all of
repo_utils/sub_tests/collapse.sh
, so will need to update VCFs inrepo_utils/answer_key/collapse
from truvari.
Great!
Thank you
from truvari.
Related Issues (20)
- BED Region off-by-one error HOT 4
- Zero matches between base and comp HOT 4
- AttributeError: 'CollapsedCalls' object has no attribute 'consolidate' | version 4.2.1 HOT 4
- Calculate SNV HOT 7
- complex genotype problem HOT 3
- GT integrate HOT 1
- No TP or FP calls for CNV HOT 1
- merging different SV type? HOT 3
- No FP or TP calls HOT 2
- Unable to run MAFFT HOT 9
- md5sum FIPS issue HOT 1
- Support vector for intra-sample merge HOT 6
- some questions about the results in fp.vcf.gz
- some questions about the results in fp.vcf.gz HOT 1
- Getting same numbers of TP-base and TP-comp HOT 4
- Suggested minor documentation changes
- Truvari, STRs and Expansion Hunter - Query HOT 2
- Bug in benchmarking HOT 4
- Request: truvari collapse --keep option to mantain the ALT sequence HOT 1
- Inquiry on the Determination of Representative Structural Variants in Merged VCF Sets HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from truvari.