To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

Memory problem with HaplotypeCaller -L

Morning,

I´m using HaplotypeCaller program spliting by chromosomes in sheep genome.

On 22 chromosomes the process has been completed correctly, but on the others there have been several errors. The next paragraph is the command used for all 158 samples. We use nodes of 16 cores (-ntc 16) and 32 Gb of memory RAM.

java -Xmx28g -Djava.io.tmpdir=/scratch/dpaule_1/dpaule_1_3/tmp -jar /home/dpaule_1/COMUNES/WGS_software/gatk/3.8.0/GenomeAnalysisTK.jar -T HaplotypeCaller -R /home/dpaule_1/
COMUNES/reference_genome/reference_genome_Oar-v3.1-r90_toVariantCalling/Ovis_aries_v3.1-r90.fa -L 9 -I /home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AFS32_SRR501853/AFS32_SRR50
1853_recal_reads.bam -I /home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AFS33_SRR501871/AFS33_SRR501871_recal_reads.bam -I /home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AU
MEHM000000000001/AUMEHM000000000001_recal_reads.bam -I
/home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AUMEHM000000000002/AUMEHM000000000002_recal_reads.bam -I /home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AUMEHM000000000003/AUMEHM000000000003_recal_reads.bam -I /home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AUMEHM000000000004/AUMEHM000000000004_recal_reads.bam -I
home/dpaule_1/dpaule_1_3/HECTOR/resultados_WGS_20171128_1704/AUMEPM000000000005/AUMEPM000000000005_recal_reads.bam....

ERROR MESSAGE: An error occurred because you did not provide enough memory to run this program. You can use the -Xmx argument (before the -jar argument) to adjust the maximum heap size provided to Java. Note that this is a JVM argument, not a GATK argument.

Best Answers

Answers

  • Hi SkyWarrior,

    Thank you for the answer. I can not run GATK without -ntc, because we are trying to reduce the analysis time as much as possible. But we will prove it as a last resort.

    Thank you for everything.

    Héctor

Sign In or Register to comment.