Question: (Closed) Handling BAM files in technical replicate situations
gravatar for David_emir
3.3 years ago by
David_emir370 wrote:

Hello all,

I have a query and i hope you will be able to guide me to understand better, Please note: I am a non Tech Guy struggling hard to understand Life science(Out of my own intrest). I have received fastq files from my friend which are from Re-sequencing experiment. To minimize the likelihood of systematic bias in sampling, 2 paired-end libraries with insert size of 500bp are prepared for all samples and are then subjected to whole-genome sequencing four lanes of each library, resulting in at least 30-fold haploid coverage for each sample. Raw image files are processed by Illumina pipeline for base calling with default parameters and the sequences of each individual are generated as 100-bp-paired-end reads. So, I have around 3-4 BAM files for each Sample (Which i am assuming as Technical replicates) of the data that I am trying to analyse so far using GATK pipeline. I have aligned each file independent of each other with BWA-MEM now. I am at a step where I have a doubt. Should I combine (Merge) all Bam files from each sample before running the next step of Collect Alignment & Insert Size Metrics step or should i process each BAM file separately and call variants at last?

Process walk through: 1. Mapping to reference --> using BWA-MEM --> Output: .SAM file 2. Sorting the SAM file by coordinate and converting it to BAM --> Output: BAM file 3. Collect Alignment & Insert Size Metrics 4. Mark Duplicates 5. Build BAM Index 6. Create Realignment Targets 7. Realign Indels 8. Call Variants

If i need to merge these BAM files what would be the best practice and at what step it should be merged? it would be great if you can share the commands as well. Sorry for long post.

Thanks a lot and a very happy new year David Emir

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by David_emir370

Hello David_emir!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.


ADD REPLYlink written 3.3 years ago by WouterDeCoster43k
Please log in to add an answer.
The thread is closed. No new answers may be added.


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