Running an alignment with merged paired end reads, can I use bowtie2 or bwa mem?
1
0
Entering edit mode
4.2 years ago
m.radz ▴ 10

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.9k views
ADD COMMENT
2
Entering edit mode
4.2 years ago
ATpoint 81k

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

ADD COMMENT
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode

Yes, -U would be the choice.

ADD REPLY
0
Entering edit mode

Thanks for your help!

ADD REPLY

Login before adding your answer.

Traffic: 2628 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6