Two Fastq Files And Alignment With Bwa.
I have two fastq files. I want to use bwa to produce two sam files as this is more efficent for me than merging the fastq files and running bwa on the bigger file. If I run bwa separately on the fastq files and merge the sam files later, will the result be different? (I don't know if this is an important detail, but I know a priori that the fastq files have reads from separate chromosomes.)
1 answer
Individual reads are mapped independently of each other. So mapping the reads in two chunks should produce the same mappings as a single file. However, I wouldn't expect mapping runs to always be identical as there is a random element to mapping repetitive hits.
• 0 views
•
link
written
by
Damian Kao
Log in to answer this question.
More posts like this
-
BBmerge output question
written by ocortes •
-
Bcftools merge taking too much time and producing large file
written by BAGeno
-
peak calling of ChIP-seq
written by Ben
-
workflow for variant calling using GATK in multiple paired end reads
written by reza
-
I have same samples in multiple lanes. What are all the steps to be taken before downstream analysi…
written by nalandaatmi •
-
Is combining fastq files before running tophat the same as combining bam file output from tophat?
written by colin.kern
-
BWA MEM alignment output of splited fastq files differ from the original(unsplit) fastq file
written by manojkumar_bhosale
-
BWA alignment using paired end reads from different lanes
written by deepue
-
Merge 4 Fastq Files
written by bashwin.u •
-
The Role Of Alignment In Processing Sequencing Data
written by KCC
