This is a test version of Biostars. For the public version, visit https://www.biostars.org.
why bwa-mem only return one read with pair-read end

Hi all,

I have a very basic question. We introduce read1 and read2 into the bwa-mem; why does in the bam file output, they only give 1 sequence in return?

sequencing bwa-mem

Have you tried name sorting the alignment file and/or viewing the alignment in IGV? Paired-end reads are not guaranteed to be next to each other in the file. There may be a discordant alignment for the read pair (because the read captures a chromosomal breakpoint). Or the second read has not mapped.

0 answers

No answers yet.

Log in to answer this question.