Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

reference letter case issue in VCF #433

Open
splaisan opened this issue Oct 12, 2023 · 1 comment
Open

reference letter case issue in VCF #433

splaisan opened this issue Oct 12, 2023 · 1 comment

Comments

@splaisan
Copy link

HI,

sniffles2 VCF results contain sequences in lower case while the reference genome is higher case in that location.
This leads to crashing when I phase the VCF using hiphase (PacificBiosciences/HiPhase#19).

Can I fix this issue at the source which seems to be here in future runs?

Thanks for your help

@splaisan
Copy link
Author

Hi,
Unfortunately, the incompatibility with hiphase does not stop here and also includes the ref and alt notation.

  • INS calls give N as ref base instead of the actual first reference base at the insertion point as instructed in the VCF format manifest.
  • DEL calls give N as alt sequence instead of the first nucleotide of the deletion as in the manifest.
    Being able to phase SV calls would be great to have but would require rewriting sniffles2 to comply with VCFV4.
    Please correct me if I missed something here or if there is a phasing tool that accepts sniffles2 data as input
    Thanks for the otherwise great tool
    Stephane

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant