Question: inconsistency in aligning metagenomic data in Bowtie2
gravatar for oatirosh
3.6 years ago by
oatirosh0 wrote:

I have reads from Illumina HiSeq, paired end, 126 bp each reads (after adaptor trimming) from metagenomic samples. reads are post filtration against human genome (hg19). I want to filter these reads against all known bacterial/fungal/viral genomes. I doing this step by step, each time outputting the un-aligned reads (using the --un-conc flag). I use --very-sensitive Bowtie2. at the end of the process I get reads that still can be mapped to viral/bacterial databases. It seems like the filtration is inconsistent and only portion of the reads are being filtered out, according to there abundance in the library.

Does anyone have an idea why this is happening and how to solve it?

reads some threads here, I though that maybe is because one mate-pair aligns to a database and the other does not, so the paired end reads is passing as unaligned.

I would be happy to hear any thoughts and solutions about this.

metagenome alignment bowtie2 • 1.0k views
ADD COMMENTlink written 3.6 years ago by oatirosh0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 829 users visited in the last hour