SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools merge bair Bioinformatics 5 10-10-2012 01:51 PM
samtools merge frymor Bioinformatics 4 10-26-2011 05:12 AM
Samtools Merge and rmdup AvinashP Bioinformatics 2 06-28-2010 02:16 AM
Samtools merge wangzkai Bioinformatics 1 05-01-2010 01:35 PM
samtools merge bair Bioinformatics 4 03-05-2010 01:23 AM

Reply
 
Thread Tools
Old 07-20-2011, 05:52 AM   #1
nans_bn
Member
 
Location: United Kingdom

Join Date: Mar 2011
Posts: 12
Default samtools merge

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.

Last edited by nans_bn; 07-20-2011 at 05:55 AM.
nans_bn is offline   Reply With Quote
Old 07-21-2011, 08:26 AM   #2
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

Try using the merge facility within Picard. It should handle merging headers properly.
nilshomer is offline   Reply With Quote
Old 07-22-2011, 04:45 AM   #3
nans_bn
Member
 
Location: United Kingdom

Join Date: Mar 2011
Posts: 12
Default

Great, thanks a lot. I'll try that out.
nans_bn is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:38 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO