Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Best way of Merging Multiple Lanes generated by HiSeq 2000, index BAM or VCF?
0
gravatar for venuraherath
2.1 years ago by
venuraherath • 20
venuraherath • 20 wrote:

Hi Folks, I just want to make sure I am on the right track. I have paired reads coming from three Lanes per each sample. Following is what I did so far; I mapped reads to reference genome using BWA in SAM tools.

Now I have indexed BAM files (3 per each sample).

I am thinking of merging bam files using Sam tools and then go for Variant calling.

I read that alternatively, variants can be called for each indexed bam file and then merge vcf files.

What is the more logical approach? Looking forward to get your suggestions. Thank you! 🙂
hiseq illumina mapping sam wgs • 811 views
ADD COMMENT • link •
written 2.1 years ago by venuraherath • 20
5
gravatar for ATpoint
2.1 years ago by
ATpoint ♦ 46k
ATpoint ♦ 46k wrote:

Given that these are lane replicates and variation between lanes is typically low, I would have merged already at the fastq level. Anyway, now you could merge them at the BAM level and then perform variant calling. Calling variants on the three technical replicates separately is probably not too informative as each of the three would have lower power to detect events than the merged file.
ADD COMMENT • link written 2.1 years ago by ATpoint ♦ 46k

Thank You!

I am going to use merge in SAMtools for merging index bam files. Do I need to keep bai files (this case 3) in the same folder as well? After merging do i need to index it again?

ADD REPLY • link written 2.1 years ago by venuraherath • 20
1

Merging does not require an index. Yes, you'll need to index the new file again, but you can do this in one command like:

samtools merge - file1.bam file2.bam file3.bam | tee merged.bam | samtools index - merged.bam.bai
ADD REPLY • link written 2.1 years ago by ATpoint ♦ 46k
1

Awesome! Thank you so much! Have a Great Day!

Venura

ADD REPLY • link modified 2.1 years ago • written 2.1 years ago by venuraherath • 20
Please log in to add an answer.

Similar posts • Search »

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