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

About the INS sequence in vcf output #344

Open
ClarenceHsiang opened this issue Aug 12, 2023 · 4 comments
Open

About the INS sequence in vcf output #344

ClarenceHsiang opened this issue Aug 12, 2023 · 4 comments

Comments

@ClarenceHsiang
Copy link

image
Dear Delly team,

Good evening.
In the output vcf file,
Alt showed without the insertion sequence, is the INS sequence in the INFO column?

Thanks again

@tobiasrausch
Copy link
Member

This seems to be a VCF file after snpEff processing. Can you please check the original output VCF of delly if the ALT column contains symbolic values like or the inserted nucleotide sequence?

@ClarenceHsiang
Copy link
Author

ClarenceHsiang commented Aug 14, 2023

image

Hello, I checked the original vcf, BND showed this kind of information (i.e. ]Pmu3_1:7879959]T) in ALT column,
for other SVs (DEL, INS .. ) are the same to the previous file.

I did the process of variant calling as below,
By followering Germline SV calling, I converted the germline.bcf to vcf through bcftools view function.
Then, the output vcf showed as the attached figure.

@tobiasrausch
Copy link
Member

I see. Thanks for letting me know. I will try to fix this in the next version, for an insertion there should always be an ALT sequence. For the time being, you can then only use INFO/CONSENSUS which indeed includes the inserted sequence but also surrounding sequence.

@ClarenceHsiang
Copy link
Author

ClarenceHsiang commented Aug 14, 2023

Thank you so much. Is there any way that we can extract the insertion sequence?

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

2 participants