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

Very low number of variants called by elprep in an heterozygous genome #67

Open
chklopp opened this issue Feb 14, 2023 · 0 comments
Open

Comments

@chklopp
Copy link

chklopp commented Feb 14, 2023

I'm trying to use elprep to call variants in a genome having one SNP every 50bp.
It runs but the g.vcf files have only hundred times less position where $5 != "<NON_REF>" as expected.
I use standard parameters in my command line without a reference vcf file because at this stage I would like to define the reference variation set by combining freebayes, bedtools and elprep results. freebayes and bedtools give high numbers of variants with most are concordant.
What should I modify in the standard elprep command line to increase the number of detected variants?

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