Adding read groups to mark duplicated bam file
1
0
Entering edit mode
7 weeks ago
reza ▴ 270

hi

I have a mark duplicated bam file and I want to use it to variant calling using GATK. GATK needs read groups and I want to add read groups to mark duplicated bam file. Now my question is: Can I add read groups via Picard to the mentioned file or should I do it immediately after mapping reads to reference?

Picard GATK • 132 views
ADD COMMENT
1
Entering edit mode
7 weeks ago
vdauwera ★ 1.1k

You have several options depending on what is most convenient to you: during alignment (you can provide read group information to BWA mem), or at any point afterward (before GATK processing) as a separate step with Picard AddOrReplaceReadGroups. The GATK development team recommends doing it as early as possible for data management reasons, but there is no penalty in terms of quality if you do it after marking duplicates, as long as you are doing alignment and duplicate marking separately per read group (assuming you have multiple read groups per sample -- if not, then you can ignore this point).

ADD COMMENT

Login before adding your answer.

Traffic: 2514 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