Question: GATK Base Quality Score Recalibration (BQSR) for merged sequencing runs
0
gravatar for ElUretsky
3.1 years ago by
ElUretsky20
Germany
ElUretsky20 wrote:

I have a BAM file with 10 different read groups all from the same sample:

  • 2 of those were test libraries sequenced in a MiSeq.
  • The remaining 8 are from a HiSeq and are much more abundant.

Can GATK recognize that the input quality scores will have a different profile for each read group? Should I split my BAM file with respect to sequencing platform or even on a per read group basis?

sequencing rna-seq • 1.2k views
ADD COMMENTlink modified 3.1 years ago by Kevin Blighe67k • written 3.1 years ago by ElUretsky20
2
gravatar for Kevin Blighe
3.1 years ago by
Kevin Blighe67k
Republic of Ireland
Kevin Blighe67k wrote:

Can GATK recognize that the input quality scores will have a different profile for each read group?

The BQSR calibration is indeed read-group aware, so, the answer to your question is 'yes'.

Should I split my BAM file with respect to sequencing platform or even on a per read group basis?

Keep them in the same BAM. The GATK will identify them by read-group.

[Source: Geraldine at GATK Forum]

ADD COMMENTlink modified 3.1 years ago • written 3.1 years ago by Kevin Blighe67k
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: 1550 users visited in the last hour