Error when Phasing with Beagle 5.2
I’m having trouble phasing a multi-sample (9-samples) vcf file produced by gatk HaplotypeCaller with Beagle 5.2. I do not have a genetic map or reference panel. I am working with a very heterozygous group of organisms (sea urchins). When I run beagle with the following command,
java -Xmx100g -jar beagle_5.2.jar gt=filtered_calls.vcf.gz out=phased impute=false
I get this error:
Window 470 [NW_022145483.1:11416-24979]
Reference markers: 630
Study markers: 630
Burnin iteration 1: 0 seconds
Exception in thread "main" java.lang.IllegalArgumentException: 0
at main.RunStats.printEstimatedNe(RunStats.java:260)
at main.Main.phaseStage1Variants(Main.java:195)
at main.Main.phaseTarg(Main.java:181)
at main.Main.phaseAndImpute(Main.java:171)
at main.Main.main(Main.java:126
Beagle usually runs fine for about 15 min and outputs 1.4G of phased genotypes, then crashes. I’m not sure what this error code means. I have been playing around with memory usage and window size (anywhere between 5 and 40). Neither has seemed to help. When I change the window size, it crashes while processing different scaffolds. I’m fairly sure that there isn’t a problem with my input vcf file as other programs have run successfully using it as input.
Are there other programs that will phase without a genetic map and reference panel? Should I split my vcf files into separate files for each chromosome before phasing with Beagle? I have access to HPC so memory should not be an issue.
• 14 views