Tweet
Hello,
I know there are many threads on samtools merge command but unfortunately, I am yet to find my answer. I hope someone can help me out in this!
I have just one sample (SOLiD data) with two different libraries - A and B.
A has 3 mate pair reads and 2 fragments. B has 1 paired end read.
I wish to merge all the above bam files into one for snp calling.
I tried it with
samtools merge -r out.bam in1.bam in2.bam in3.bam ...etc
I obtained a out.bam file. Now I tried to call for the snps but when I run the command it gives me an error stating
[group_smpl] Read group <B_sample.hg19> used in file <out.bam> but absent from the header or an alignment missing read group.
I would appreciate any suggestions..
Thank you.
I know there are many threads on samtools merge command but unfortunately, I am yet to find my answer. I hope someone can help me out in this!
I have just one sample (SOLiD data) with two different libraries - A and B.
A has 3 mate pair reads and 2 fragments. B has 1 paired end read.
I wish to merge all the above bam files into one for snp calling.
I tried it with
samtools merge -r out.bam in1.bam in2.bam in3.bam ...etc
I obtained a out.bam file. Now I tried to call for the snps but when I run the command it gives me an error stating
[group_smpl] Read group <B_sample.hg19> used in file <out.bam> but absent from the header or an alignment missing read group.
I would appreciate any suggestions..
Thank you.
Comment