adding RG tags to sorted bam file
Hello,
Could anyone help me with adding read groups to the sorted bam files, I have tired with samtools merge but failed and I have found a one line script
perl -e 'print "\@RG\tID:Disc1\tSM:'$x'\tPL:Illumina\n"' >rg.txt
samtools merge -rh rg.txt - index10.aln.sort.bam
adds the read group but later calling for SNP with freebayes causes an error in the vcf file.
Thank you all
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
READ GROUP in GATK
written by Payal
-
Using shell command in Nextflow DSL=2 for BWA mem aligner
written by rohitsatyam102
-
freebayes error for SNP calling
written by evafinegan •
-
Pileup or merge?
written by storm1907 •
-
samtools fixmate: Couldn't write to output file: Broken pipe
written by Ric
-
Variants calling by bcftools
written by Mbillah
-
Dealing with multiple read groups, from BWA to sorted .bam
written by jamiedm •
-
bwa mem: Passing a variable to read group
written by firestar
-
Running GATK DepthOfCoverage on BAM files with multiple RG's
written by stefano.iantorno •
-
Sam/Bam File Samfilereader
written by mad.cichlids
samtools merge adds the RG from the filename, not from the header you're adding, so I would assume that you're adding one thing in the header and something else on the end of the individual alignment lines. Also make sure you have an intact bam header, not just that @RG line.
Is x defined?
X is given as sample name
Have you tried Picard tools? It has a convenient function for this.
I can recommend Picard:
I'm trying Picard, hope it works ..,
Thank you