GATK calling Merged Bam files

[jayce_ocean]

Dear All,

can share your idea? i have tried to merge my bam files (bam file1, bam file2...so on) using samtools and picard. After then, i used GATK to call the SNP and refer to the merged bam file. However, the GATK always throw out genotyp results with 2 columns or more pertaining to my number of bam files prior the merging instead of a single column result for the merged bam file.

Any idea?anyone can share the experience?many thanks in advance...

Regards
Jayce

[jayce_ocean]

i had done reheader and add the @RG tag referring to my merged bam file. But still, GATK always throw an error message stating the read error message or output each of every bam file results in each column prior the merging....any idea?i wanna results from GATK 1 column referring to the merged bam file results. thanks...

[jayce_ocean]

any trick?anyone have idea?or any settings in GATK i need to twist on?please help to share some idea or experience...thanks..

[Jose Blanca]

We usually merge bam files for GATK. Our code is exposed in a pipeline (http://bioinf.comav.upv.es/ngs_backbone/), but if you like python you could use the underlying library (https://github.com/JoseBlanca/franklin).

Regards.