bowtie2 output all aligned reads from paired alignment (even unconcordant pairs or single reads)
0
0
Entering edit mode
5.9 years ago
JV ▴ 410

Hi,

I have assembled and binned a small metagenome.Now I want to extract all reads that map to the contigs of each bin and reassemble them and see if the assemblies improve significantly.

For this I have mapped the original reads back to the bins using bowtie2.

I know that bowtie can output all mapped reads in fastq-format which would be perfect for my needs.

I have a mixture of paired and unpaired reads after adapter-trimming and quality clipping and I want to keep the paired-end information when possible. The problem is, when running bowtie with paired-end settings, i only have the option to output the concordantly aligning reads (argument --al-conc). This does not include paired reads that align to seperate contigs. As it may be that such unconcordantly aligning pairs may connect different contigs to possible scaffolds, I would like to include these pairs as well.  Using the option "--un-conc", however,  produces all other pairs (those which align unconcordantly and those which don't align at all).

When running bowtie with unpaired reads it is quite straiightforward to output all aligned single reads (just add the "--al" argument), but then I would have to rather painfully reconstruct the read pairs and the orphaned reads from that output.

Is there any elegant way to solve this problem?

alignment Assembly bowtie2 BAM SAM • 2.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 1569 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