Hello Everyone,
I have noticed an issue with my BAM file headers, where the @RG line is either mal-formed or is missing entirely. I think I can sed the files that are mal-formed, and add the sample names necessary to complete my further analyses in GATK. But is there a way to add an entire line to the headers?
Thanks,
2 answers
samtools addreplacerg -w -r ID:oldid -r LB:Lib -r SM:SAMPLE -O BAM -o output.bam bad.bam
You should be able to open and edit the header in plain text format, and you can just add it to the top of a headerless bam. Something like
samtools view -H bad.bam > bad_header.txt
<edit bad_header.txt>
samtools view bad.bam | cat good_header - | samtools view -hb - > good.bam
Hi,
I think that's close to the solution I want, but the BAM files I have are not headerless, their headers are just really poorly formatted. I simply want to add a line to the header, I think I can do what you said, but then use samtools reheader to edit the original BAM.
So:
samtools reheader edited_header.sam bad.bam > good.bam
Log in to answer this question.
-
Segmentation fault PLINK 1.90b
written by Alexandros.Frydas •
-
Problem assigning read group to a sample
written by taniamahmood38 •
-
edit headers of a fasta
written by Candela •
-
Corrupted FASTq files with missing "+" under some sequences.
written by akh22
-
Rename fasta headers using a variable already present in the header?
written by FreshBio •
-
How treat a Bam file with mutiple read groups (RG)
written by zamani_2012 •
-
Error parsing SAM header. @RG line missing SM tag
written by kwangchol.mun •
-
Read group header line
written by win
-
Add Information To @Rg Using Picard, Without Replacing Read Groups
written by DonJoe •
-
How To Set Rg Header For Tablerecalibration?
written by PeterPan •
