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

No variants found during rerun #211

Open
dmjask opened this issue Jan 23, 2025 · 0 comments
Open

No variants found during rerun #211

dmjask opened this issue Jan 23, 2025 · 0 comments

Comments

@dmjask
Copy link

dmjask commented Jan 23, 2025

Hi,
I have previously successfully used PEPPER-Margin-DeepVariant to obtain haplotagged reads (--only-haplotag flag) for a chromosome of interest. Recently I needed to repeat this for different chromosome, which yielded empty vcf; no candidates or variants found according to PEPPER.
To troubleshoot I repeated the haplotagging for the original chromosome, which also resulted in empty vcf and no candidates/variants.

[01-23-2025 09:23:48] INFO: STARTING CANDIDATE FINDING. [01-23-2025 09:23:49] INFO: FINISHED PROCESSING, TOTAL CANDIDATES FOUND: 0 [01-23-2025 09:23:49] INFO: FINISHED PROCESSING, TOTAL VARIANTS IN PEPPER: 0 [01-23-2025 09:23:49] INFO: FINISHED PROCESSING, TOTAL VARIANTS SELECTED FOR RE-GENOTYPING: 0 [01-23-2025 09:23:49] INFO: FINISHED PROCESSING, TOTAL SNP VARIANTS SELECTED FOR RE-GENOTYPING: 0 [01-23-2025 09:23:49] INFO: FINISHED PROCESSING, TOTAL INDEL VARIANTS SELECTED FOR RE-GENOTYPING: 0

`Running OpenMP with 64 threads.

Parsing model parameters from file: /opt/margin_dir/params/phase/allParams.haplotag.ont-r94g507.snp.json
Parsed 0 total VCF entries from input.vcf.gz; kept 0 HETs, skipped 0 for region, 0 for not being PASS, 0 for being homozygous, 0 for being INDEL
No valid VCF entries found!`

I am running the software within a conda env.

singularity exec $pepper_margin_deepvariant run_pepper_margin_deepvariant call_variant \ -b "$indir/${ID}_${CHR}_${START_POS}_${END_POS}.bam" \ -f "${REF}" \ -o "${outdir}" \ -p "${POUT}" \ -t 64 \ --ont_r9_guppy5_sup \ --only_haplotag \ -v "$ref":"$ref"

Any suggestion what may be causing this issue?

Thanks!

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