This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Running an alignment with merged paired end reads, can I use bowtie2 or bwa mem?

I want to align some paired end reads that I merged and trimmed with BBmerge and BBduk to a referene genome. Bwa mem and bowtie2 seem to be popular options but I can't work out if merged paired end reads can be used as input for these programs. Has anyone done this? If not are there any other programs that can take merged paired end reads as input?

alignment next-gen sequencing sequence

1 answer

Sure, a merged PE read is simply a single-end read and both tools support that.

For bowtie2 would I then flag the reads with -U (Comma-separated list of files containing unpaired reads to be aligned)? The only other options are -1,-2 for mate pairs and -interleaved

Yes, -U would be the choice.

Thanks for your help!

Log in to answer this question.