This is a test version of Biostars. For the public version, visit https://www.biostars.org.
multithread support for gatk

The below command is meant to use bqsr with multithreading:

gatk BaseRecalibrator --num_cpu_threads_per_data_thread=8 \
-I /home/cmccabe/Desktop/fastq/NA12878_marked_duplicates.bam \
-R /home/cmccabe/Desktop/fastq/resources/ucsc.hg19.fasta \
--known-sites /home/cmccabe/Desktop/fastq/fastq/resources/dbsnp_138.b37.vcf \
--known-sites /home/cmccabe/Desktop/fastq/resources/Mills_and_1000G_gold_standard.indels.hg19.vcf \
-O /home/cmccabe/Desktop/fastq/NA12878_recal_data.table

I have also tried nct 8 and both return an unrecognized exit code, so It does not appear to be supported or more likely I am not using it correctly? Thank you :).

gatk bqsr
• 3,891 views • link
written by bioguy24 bioguy24

What version of GATK?.

• 0 views • link
written by GenoMax GenoMax

no sorry forgot to mention, the version is 4.0.6.0.... thank you :).

• 0 views • link
written by bioguy24 bioguy24

1 answer

Update: Looks like you can't use nct with GATK 4.0. You have to use SPARK for parallelization.

• 0 views • link
written by GenoMax GenoMax

Thank you very much :).

• 0 views • link
written by bioguy24 bioguy24

Log in to answer this question.

More posts like this
  • Hey guys, I'm having a prob when using GATK4 BQSR . This dbsnp vcf file has chromosomes notated a…
    written by ANImaL • ANImaL
    anilkumar@ak-omen-laptop:~/NGStools/gatk-4.4.0.0$ gatk --java-options "-DGATK_STACKTRACE_ON_USER_EXCEPTION=true" BaseRecalibrator -I "/media/anilkumar/My Passport/CRC/fastq/C_4_mkdp.bam" -R "/media/anilkumar/My Passport/CRC/fastq/hg19.fa" --known-sites "/media/anilkumar/My Passport/CRC/fastq/dbsnp_138.b37.vcf" --known-sites "/media/anilkumar/My Passport/CRC/fastq/Mills_and_1000G_gold_standard.indels.b37.vcf" --known-sites "/media/anilkumar/My Passport/CRC/fastq/1000G_phase1.indels.b37.vcf" -O "/media/anilkumar/My Passport/CRC/fastq/C_4_bqsr.table" Using GATK …
  • finding compatible version of reference genome for use of GATK
    written by chrisclarkson100 chrisclarkson100
    I am mapping a series of fastq files, containing ChiP-Seq data, with bowtie2 (my mapping software of choice) to the GRCh37 genome: `https://www.gencodegenes.org/human/release_39lift37.html` . I …
  • So many variants detected.
    written by Mehmet Mehmet
    Dear All, I have done variant calling in Germline data that has single sample of each individual and two genes. I did following steps, but …
  • GATK BaseRecalibrator error.. how do I solve this?
    written by kwanghoon • kwanghoon
    I downloaded reference genome and known site vcf at GATK bundle (https://console.cloud.google.com/storage/browser/genomics-public-data/resources/broad/hg38/v0;tab=objects?pli=1&prefix=&forceOnObjectsSortingFiltering=false) I sorted and duduplicated bam file then tried to use baserecalibrator. But I …
  • gatk vqsr --recal-file error
    written by bioguy24 bioguy24
    In GATK 4.0.6.0 I am getting the below error using the command: gatk VariantRecalibrator \ -R /home/cmccabe/Desktop/fastq/resources/ucsc.hg19.fasta \ -V NA12878.vcf.gz \ --resource hapmap,known=false,training=true,truth=true,prior=15.0:/home/cmccabe/Desktop/fastq/resources/hapmap_3.3.hg19.sites.vcf.gz \ --resource …
  • error using bqsr and incompatible contigs
    written by bioguy24 bioguy24
    Using the below command in GATK 4.0.6.0, I am getting: A USER ERROR has occurred: Input files reference and features have incompatible contigs: No overlapping …
  • GATK MarkDuplicates output and bam
    written by bioguy24 bioguy24
    The below command executes and produces the dup_metrics below, I am trying to interpret them and don't really understand the metrics and in the bam …
  • illumina data analysis with bwa mem
    written by bioguy24 bioguy24
    I a new to illuminia data/analysis and just want to be sure I am processing correctly. I am using the below to pair each forward …
  • How about my parallel bwa+gatk pipeline?
    written by bluemonster0808 • bluemonster0808
    Normally, we run NGS analysis in the follwing steps: 1、 ~/bwa-0.7.16a/bwa mem -t 32 -M -R @RG\tID:1\tPL:Illumina\tSM:HCC1143 ~/reference/human_g1k_v37_decoy.fasta 1.fastq 2.fastq | samblaster --splitterFile my.split.sam --discordantFile …
  • Number of snps are very low
    written by waqasnayab waqasnayab
    <p>Hi,</p> <p>Dear All,</p> <p>I have WGS data of one healthy individual. I want to make a VCF file. I am running the following pipeline but …