SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
input BAM files for GATK Jane M Bioinformatics 26 07-30-2015 09:58 PM
Convert merged BAM back to per lane BAM or FASTQ file danielsbrewer Bioinformatics 6 10-03-2013 07:29 AM
RG tags in merged BAM files, muliple lanes one sample cow_girl Bioinformatics 17 08-27-2013 11:28 PM
GATK pileup with merged BAM files yxl Bioinformatics 0 04-22-2011 07:07 PM
Calling multiple BAM files for SNPs and vcf newbietonextgen Bioinformatics 3 04-19-2011 11:29 AM

Reply
 
Thread Tools
Old 03-03-2011, 11:40 PM   #1
jayce_ocean
Junior Member
 
Location: location

Join Date: Sep 2010
Posts: 5
Default GATK calling Merged Bam files

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 is offline   Reply With Quote
Old 03-03-2011, 11:44 PM   #2
jayce_ocean
Junior Member
 
Location: location

Join Date: Sep 2010
Posts: 5
Default

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 is offline   Reply With Quote
Old 03-04-2011, 12:07 AM   #3
jayce_ocean
Junior Member
 
Location: location

Join Date: Sep 2010
Posts: 5
Default

any trick?anyone have idea?or any settings in GATK i need to twist on?please help to share some idea or experience...thanks..
jayce_ocean is offline   Reply With Quote
Old 03-16-2011, 12:15 AM   #4
Jose Blanca
Member
 
Location: Valencia, Spain

Join Date: Aug 2009
Posts: 70
Default

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.
Jose Blanca is offline   Reply With Quote
Reply

Tags
gatk, merged bam

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 06:48 AM.


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