freebayes / freebayes

Bayesian haplotype-based genetic polymorphism discovery and genotyping.
http://arxiv.org/abs/1207.3907
MIT License
790 stars 265 forks source link

Freebayes stuck #774

Closed seynard closed 8 months ago

seynard commented 1 year ago

Hi,

I'm running freebayes for de novo variant calling on Oxford Nanopore Technology sequencing data (bam files produced by PromethION) for pool sequenced experiment. I'm testing freebayes on a chromosome for 15 pools and it stays stuck after 7Kb. It does not get killed but it also does not write anymore. I looked at my row data on IVG and I do not have a large coverage in this region (just a couple of X) as it as been mentioned in the issue #614 (https://github.com/freebayes/freebayes/issues/614). If I run freebayes with 3 pools the calling is performed but form 4 pools onwards (any of them) it gets stuck. Here is the command i'm running:
freebayes -f ${fasta} --min-alternate-count 1 --min-alternate-fraction 0.01 --use-best-n-alleles 2 --cnv-map ${bam_list} --pooled-discrete ${dirwork}/freebayes/*${c}.header.sorted.rename.filter.bam | vcffilter -f "QUAL > 20" -f "AF < 1" | vcffilter --or -f "TYPE = del" --or -f "TYPE = ins" -f "TYPE = snp" > ${dir_work}/freebayes/${c}.freebayes.vcf

I already tried to ask for more memory (up to 500G), without success. Would you have an idea where this issue is coming from?

friburgo-moc commented 1 year ago

I have published an issue that may be similar (https://github.com/freebayes/freebayes/issues/777), does it work better if you used a higher min-alternate-fraction? does it work better if you use a freebayes version lower than 1.3.1?

github-actions[bot] commented 8 months ago

This issue is marked stale because it has been open 120 days with no activity. Remove stale label or comment or this will be closed in 5 days

github-actions[bot] commented 8 months ago

This issue was closed for lack of activity. Feel free to re-open if someone feels like working on it.