Biostar Beta. Not for public use.
Replace @RG group by the one of another BAM
0
Entering edit mode
15 months ago
da.de • 0

Hi, I want to replace the @RG info from one bam file (output.bam) by the @RG info from another (input.bam) creating a new.bam I tried it as follows:

samtools view -H input.bam | grep '@RG' > RG_header

samtools addreplacerg -r RG_header -o new.bam output.bam

I get the following message:

[parse_args] The supplied RG line lacks an ID tag.

RG_header looks as follows:

@RG ID:CAUPUANXXnnn PU:170913_7001253F  PL:ILLUMINA LB:55426    SM:55426

Thanks for your help.

Best, Dagmar

ADD COMMENTlink
1
Entering edit mode
13 months ago
Belgium

-r expects a string, not a file.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1