Question: Error in GATK BaseRecalibrator arguments
0
gravatar for manubiomed20
11 months ago by
manubiomed2010
manubiomed2010 wrote:

I tried to recalibrate the bases of my complete genome assembly, but it gave the following error:

 invalid argument `I`

I used this command:

 -gatk BaseRecalibrator \
 -I my_reads.bam \
 -R reference.fasta \
 --known-sites sites_of_variation.vcf \
 --known-sites outro / opcional / setOfSitesToMask.vcf \
 -O recal_data.table

I wrote:

gatk BaseRecalibrator \I- /home/ims.santos06/bam/SRR062634.bam -- home/ims.santos06/bam/Hg38knownindels.vcf -R \home/ims.santos06/reference/hg38.fa   -O  /home/ims.santos06/bam/SRR062634recal_data.table
gatk assembly • 705 views
ADD COMMENTlink modified 11 months ago • written 11 months ago by manubiomed2010

I see two command , what did you wrote in the end ?

what is that \I- ?

what was the error message ?

ADD REPLYlink written 11 months ago by Pierre Lindenbaum129k

This message appeared:

First pass of the Base Quality Score Recalibration (BQSR) -- Generates recalibration table based on various
    user-specified covariates (such as read group, reported quality score, machine cycle, and nucleotide context).
    Version:4.1.2.0
    Required Arguments:

    --input,-I:String             BAM/SAM/CRAM file containing reads  This argument must be specified at least once.
                                  Required. 

    --known-sites:FeatureInput    One or more databases of known polymorphic sites used to exclude regions around known
                                  polymorphisms from analysis.  This argument must be specified at least once. Required. 

    --output,-O:File              The output recalibration table file to create  Required. 

    --reference,-R:String         Reference sequence file  Required. 

    A USER ERROR has occurred: Invalid argument 'I-'
ADD REPLYlink modified 11 months ago by RamRS27k • written 11 months ago by manubiomed2010

From what I understand I should add the covariates to the input, but I don't know how.

ADD REPLYlink modified 11 months ago by RamRS27k • written 11 months ago by manubiomed2010

Could you post the exact command you used. FYI you should use -I and not I-

ADD REPLYlink written 11 months ago by Nicolas Rosewick8.8k

I used this command:

gatk BaseRecalibrator -I /home/ims.santos06/bam/SRR062634.bam -R /home/ims.santos06/reference/hg38.fa  --known-sites /home/ims.santos06/bam/Hg38knownindels.vcf -O  /home/ims.santos06/bam/SRR062634recal_data.table
ADD REPLYlink modified 11 months ago by RamRS27k • written 11 months ago by manubiomed2010

This message Appeared:

A USER ERROR has occurred: Input /home/ims.santos06/bam/Hg38knownindels.vcf must support random access to enable queries by interval. If it's a file, please index it using the bundled tool IndexFeatureFile
ADD REPLYlink modified 11 months ago by RamRS27k • written 11 months ago by manubiomed2010

manubiomed20 : Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized. SUBMIT ANSWER is only for new answers to the original question.

ADD REPLYlink modified 11 months ago • written 11 months ago by genomax85k
0
gravatar for manubiomed20
11 months ago by
manubiomed2010
manubiomed2010 wrote:
thanks. You must index vcf and then use the command below:
gatk IndexFeatureFile.

gatk is running recalibration. I believe now it worked
ADD COMMENTlink written 11 months ago by manubiomed2010
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1926 users visited in the last hour