Question: Adding read group information in a bam header
2
gravatar for gkuffel22
3.4 years ago by
gkuffel22170
United States
gkuffel22170 wrote:

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.

galaxy vcf bam • 1.5k views
ADD COMMENTlink modified 3.4 years ago by Jennifer Hillman Jackson25k • written 3.4 years ago by gkuffel22170
0
gravatar for Jennifer Hillman Jackson
3.4 years ago by
United States
Jennifer Hillman Jackson25k 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

ADD COMMENTlink written 3.4 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 132 users visited in the last hour