Please help, I have tried everything. So I have 200 unique samples, I have processed them accordingly [fastq groomer, trim, bwa]. I now have 100 bam files (the reads were paired end) and I would like to merge these bam files so that I can run FreeBayes (variant analysis). The only problem is I would like to avoid entering in unique read group data for each sample manually. Is there a tool or a script that can accomplish this? I have tried running samtools merge -r I was under the impression that this would attagen a RG tag to each alignment using a value inferred from the input file name, this did not work for me, still no RG info. In case anyone is wondering, the reason the RG group is important is ultimately I would like to have a VCF file that will list all variants found separated by unique sample. Any help is greatly appreciated.
Heads up! This is a static archive of our support site. Please go to help.galaxyproject.org if you want to reach the Galaxy community. If you want to search this archive visit the Galaxy Hub search
Question: Adding read group information in a bam header
2
gkuffel22 • 170 wrote:
ADD COMMENT
• link
•
modified 3.6 years ago
by
Jennifer Hillman Jackson ♦ 25k
•
written
3.6 years ago by
gkuffel22 • 170
0
Jennifer Hillman Jackson ♦ 25k wrote:
Hello,
If you do not wish to do this through the UI, then using the API for batch jobs is the solution. I believe that this info was shared with your for another tool, but I can let you know that the API works with all tools. More is in our wiki:
http://wiki.galaxyproject.org/Admin
http://wiki.galaxyproject.org/Develop/API
Best, Jen, Galaxy team
Please log in to add an answer.
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 16.09
Traffic: 174 users visited in the last hour