GATK Base Quality Score Recalibration (BQSR) for merged sequencing runs
1
0
Entering edit mode
7.1 years ago
ElUretsky ▴ 30

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 • 2.0k views
ADD COMMENT
2
Entering edit mode
7.1 years ago

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 COMMENT

Login before adding your answer.

Traffic: 2914 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6