Jesper Eisfeldt
Jesper Eisfeldt
Hello! Sorry for the delay! Interesting error, usually I get this if the vcf is either corrupted or if there is issues with the header. Can you send the vcf...
Hello! Yes, I agree, refactoring is really needed, especially for collaboration! Thanks a lot! Feel free to tell us when the cleanup branch is ready, or if you have questions,...
Thanks, I'll work on that!
Hello! I have added SUPP_VEC similar to SURVIVOR: it is a list of 1 or 0 telling which callers detected the variant. The callers are listed in the header. Maybe...
Hello there! That's a tricky issue, it usually happens if the start position of the SV is the same, but the end position differs greatly. I have an example here...
Woopsiedaisy! And thanks for your suggestions! could you send me the variants you are trying to merge? Then I will look in to this! //JEsper
Hello! I have had a look, thanks for the test data! I notice you are correct, I have edited the readvcf function to this: posB=posA+abs(int(description["SVLEN"])) and now the variants merge...
Hello! Sorry for the delay! I have had a look, and implemented separate treatment of the INS variants. now PosB=PosA (i.e all ins will have lenght 0 in the reference);...
Hello! That is a good idea! Maybe we should add a "update" module? If we make it smart, we might even be able to use it for adding variants to...
Hello! Thanks for your feedback, and nice to hear that you found a solution! Have you checked the TIDDIT header? I believe the ##INFO= originates from an older version of...