How do I specify the Recalibration_data.table for "PrintReads" in GATK 4.3, if the tool does not admit -BQSR as an option anymore? · Biostar

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

How do I specify the Recalibration_data.table for "PrintReads" in GATK 4.3, if the tool does not admit -BQSR as an option anymore?

I executed BaseRecalibrator on my bam and obtained a Recalibration_data.table file. Now, When I apply PrintReads, the option -BQSR is not acknowledged by PrintReads. Even online or in the list of arguments there is nothing I can relate to this, the how does PrintReads apply the recalibration if I don't provide this file?

Keep in mind I'm very new in the bioinformatics field. Thanks a lot!

baserecalibrator bqsr printreads

• 1,060 views • link

updated 3.3 years ago by Ram 4.5K • | written 3.3 years ago by javiflaja 5

1 answer

do you want ApplyBQSR "Apply base quality score recalibration " ? https://gatk.broadinstitute.org/hc/en-us/articles/360036856671-ApplyBQSR

• 0 views • link

written 3.3 years ago by Pierre Lindenbaum 17K

Log in to answer this question.

More posts like this

BSQR in GATK without known variants
written 1.8 years ago by fuwamozu 0

Hi, to preface: - I am new to bioinformatics (not my field, just using bioinformatic tools to process some data) - I found some similar …
Exome alignment and preprocessing: When I perform an exome alignment, should I use a Ref Genome or …
written 3.3 years ago by javiflaja 5

I'm learning the basis of preprocessing, and I can't find anywhere a source that would tell me what's the difference between preprocessing a genome for …
Base recalibration in normal vs. tumor somatic variant calling in WXS data?
written 4.9 years ago by rebeliscu 6

Hi there, I have a tumor and a normal BAM file and am preparing to run base recalibration. I was planning on calling variants on …
Recalibrating variants in GATK4.1 as compared to GATK3.7
written 5.7 years ago by vinayjrao 26

Hello everyone, I had been using GATK3.7 until last Septemeber to identify variants from human whole exome sequencing samples. However, recently, our institute has upgraded …
Scatter / Gather for BaseRecalibrator on a single human WES dataset?
written 6.6 years ago by asg 0

Hello everyone! Please excuse me if this question is a bit naïve: I'm new to bioinformatics in general and GATK in particular. I am using …
GATK quality recalibration to get recalibrated BAM files before SNP calling
written 7.3 years ago by alicia.poppy 5

I'm having problems with **quality score recalibration** using **GATK** as most tutorials/examples use an old version of GATK which has different syntax and arguments to …
GATK PrintRead Producing Empty File
written 7.5 years ago by BAGeno 19

Hi, I am running gatk command on human chr21 which I downloaded from UCSC. I have performed alignment with `bwa`. Then sorted, marked duplicates and …
DNAseq pipeline steps
written 8.2 years ago by shuksi1984 6

Following are my DNAseq pipeline steps keeping GATKBestPrac in mind: Step1-Quality check of raw data (fastqc) Step2-Alignment with reference genome(bwa) Step3-SAM to BAM conversion and …
[ERROR] Malformed VCF: empty alleles are not permitted in VCF records
written 10.6 years ago by umn_bist 39

I am running BaseRecalibrator for my RNA-seq: java -jar -Xmx120g ${GATK} -T BaseRecalibrator \ -R "${reference}" \ -I "${file4}" \ -knownSites "${gerVar}" \ -knownSites "${somVar}" …
Base Recalibration Of Stampy Alignments In Gatk
written 13.4 years ago by jake 2

<p>I am struggling to produce a good base recalibration in GATK (version 2.4-9-g532efad) on multi-sample, single-end illumina data sets aligned with Stampy. After mapping reads …