Question: GATK Base Quality Score Recalibration (BQSR) for merged sequencing runs
0
gravatar for ElUretsky
23 months 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 • 856 views
ADD COMMENTlink modified 23 months ago by Kevin Blighe48k • written 23 months ago by ElUretsky20
2
gravatar for Kevin Blighe
23 months ago by
Kevin Blighe48k
Kevin Blighe48k 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 23 months ago • written 23 months ago by Kevin Blighe48k
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: 1211 users visited in the last hour