View Single Post
Old 10-19-2011, 12:11 AM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

In addition to the @RG line in the SAM header, every single read belonging to that read group has to declare this in its tags - and that part means adding this requires a lot of IO and will be comparatively slow.

If all you want to do is add some information to an existing read group, you just need to edit the SAM header. Using samtools reheader let's you modify this efficiently in a BAM file.

Last edited by maubp; 10-19-2011 at 12:12 AM. Reason: Typo
maubp is offline   Reply With Quote