second pass base recalibration gatk4

SebBatSebBat BuffaloMember

hi there

I'm trying to run the second pass for the BaseRecalibrator with GATK4.
as seen in this post ( the previous -BQSR was substituted with -bqsr. also, it seems that there is no second pass available to do that?!

when I use it in:
GATK/gatk- BaseRecalibrator \
-R CanFam3.1.dna.toplevel.fa \
-I sorted_WES_.bam \
--known-sites GATK/final.Axelsson.SNPs_nochr.vcf \
-bqsr recal.table \
-O post_recal.table

I get indeed the error "b is not a recognized option", suggesting that the option is not there..

however, the documentation here ( still says that you can run it.

in the --help there is no mention of the second pass and I'm confused about what I can/should do.

also, in previous versions I used PrintReads with the recalibration table, while the first link above mentions, which one is the one to use?

Best Answers


Sign In or Register to comment.