Hwo to identify that BQSR is performed on CRAM file

Question

[Sd]

Hi, I have a bunch of CRAM files of WGS that I want to check if Base Quality Score Recalibration (BQSR) has been done or not. Does anyone can help me how can I check it?

Answer 1

[raphael.B]

If you run samtools view -H on your alignment file, you should have an entry starting with @PG ID:GATK ApplyBQSR .... in your CRAM/BAM header.

[Sd]

Do you have an idea that how long on average does it take to run gatk ApplyBQSR on a whole genome sequencing (WGS) data? Does it necessary to perform this step on CRAMs?

[raphael.B]

I don't know what are your computing resources but compared to the alignement itself it's quite short.