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

VCF with both SNV and non SNV's #32

Open
Takadonet opened this issue Dec 15, 2017 · 3 comments
Open

VCF with both SNV and non SNV's #32

Takadonet opened this issue Dec 15, 2017 · 3 comments

Comments

@Takadonet
Copy link

I am looking to have both SNV called and non SNV's in the same vcf for downstream analyses.

My temporary workaround at the moment is to run the tool twice, once with the reference to generate the SNV's and other without to get non SNV. Then combine both of them using the SNV vcf's as first pick using vcf-vcf intersect tool

Is this possible?

@karel-brinda
Copy link
Owner

karel-brinda commented Mar 1, 2018

Hi Phil, I wanted to answer this question and put the answer to the FAQ section in the readme (will be created). However, I realized that I don't completely understand the trick with non-SNVs. Could you explain it in more details, please?

@Takadonet
Copy link
Author

Hey Karel,

The overall trick to have both SNV and non-SNV in the same vcf file is to run ococo 2 separate times. First time you run it, you pass the reference file, second time you don't.

After that I combine both vcf output files with vcf-vcf intersect and ensure that SNV only vcf is applied onto of the non-SNV version. Below is part of my Galaxy Workflow.

combineococo

My hope is to combine aligning, non-SNV calling , SNV calling and vcf vcf intersect into a SINGLE tool/job on Galaxy. Maybe even directly convert it to a bcf file format.

@Takadonet
Copy link
Author

@karel-brinda is there any plans to add this feature?

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