HaplotypeCaller using one sample split in two bam files
AnsweredI have a sample that has been split and ran in two different machines. For that sample I have reads with two different read groups and in 2 pairs of fastq files. Would be okay to align them separately so I will obtain 2 different bam files, each of the bam files will have the appropriate read group. Would be then possible to pass those 2 bam files to the base quality recalibration tool and then to the haplotypecaller? Otherwise, could you please suggest how to deal with this?
-
Thank you for your post, elena vigorito ! I want to let you know we have received your question. We'll get back to you if we have any updates or follow up questions.
Please see our Support Policy for more details about how we prioritize responding to questions.
-
Hi elena vigorito,
Thank you for writing to the GATK forum! I hope that we can help you sort this out.
In this case, you can use the MergeSamFiles (Picard) tool. This tool will allow you to merge your BAM files from different read groups into a single file. You’ll want to combine these before inputting them into HaplotypeCaller. Please find an example command below.
java -jar picard.jar MergeSamFiles \
I=input_1.bam \
I=input_2.bam \
O=output_merged_files.bamI hope this helps! Please let me know if you find success with this. If you have any other questions, please do not hesitate to reach out.
Best,
Anthony -
Hi elena vigorito,
We haven’t heard from you in a while so we will be closing out your ticket in our system. If you still require assistance, you need only respond to this thread, and we’ll create a follow-up ticket to pick up where we left off.
Thank you again for contributing to our GATK forum!
Best,
Anthony
Please sign in to leave a comment.
3 comments