You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hello Adrian
Apparently the decompose_blocksub tool does not produce proper VCF's in respect to REF types. You can see this demonstrated here. The filtered.vcf has been produced by processing multiple samples from a control DNA korell with freebayes, breaking them with breakmulti and filtering them. However decompose_blocksub will attempt to break REF fields further, which they should not (a R type should always contain 2 elements).
Further it will throw an exception on MISSING samples, which might or might not be desired behaviour depending on your implementation.
Thanks a lot.
Fohlen
The text was updated successfully, but these errors were encountered:
Hello Adrian
Apparently the
decompose_blocksub
tool does not produce proper VCF's in respect toREF
types. You can see this demonstrated here. Thefiltered.vcf
has been produced by processing multiple samples from a control DNA korell withfreebayes
, breaking them withbreakmulti
and filtering them. Howeverdecompose_blocksub
will attempt to breakREF
fields further, which they should not (aR
type should always contain 2 elements).Further it will throw an exception on
MISSING
samples, which might or might not be desired behaviour depending on your implementation.Thanks a lot.
Fohlen
The text was updated successfully, but these errors were encountered: