SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Paired end alignments -- What to throw away and what to keep? Mike2188 Bioinformatics 0 07-30-2014 02:48 PM
--best --strata in paired end alignments using Bowtie carmeyeii Bioinformatics 0 12-06-2012 05:16 PM
bowtie- paired end - no alignments madsaan Bioinformatics 1 06-27-2011 01:24 PM
Paired end reads: line # does not match juan Illumina/Solexa 0 10-27-2010 05:56 PM
Bowtie got different alignments for fliped paired-end 狗熊趴趴 Bioinformatics 0 06-26-2010 07:39 AM

Reply
 
Thread Tools
Old 01-26-2017, 09:32 AM   #1
Azufre451
Junior Member
 
Location: Italy

Join Date: Jan 2017
Posts: 2
Default Bowtie2 paired end alignments do not match

Hi all,

Maybe someone already encoutered this issue... I am trying to use bowtie2 to map some paired-end reads against a set of assembled contigs (the contigs were built using those reads with SPAdes).

I am using bowtie2 in this way:
Code:
bowtie2 --threads 10 --end-to-end --no-unal -x contigs -1  R1.fastq -2 R2.fastq | samtools view -bS - > bams_bowtie2.3/my.bam
However, when inspecting the SAM/BAM (either with Tablet, IGV or manually) I see many cases in which the two mates do not align on the same contig.

I tried to look at the 3rd and 7th fields of the SAM files (the contig-names to which the two mates of the pair align), but the information are discordant. Let's say that I have an alignment that says:

R1/1 -> contig_1 | R1/2 -> contig_2

then I would expect to find another record where

R1/2 -> contig_2 | R1/1 -> contig_1

However, this isn't always true (and for some of the samples I am mapping, it rarely is!).

I also opened an issue/question on the official bowtie2 repository (here) with more detailed explaination of the problem.

I am pretty sure I am missing something important in how to use Bowtie2 on paired-end data... where is the problem?

Thank you very much for your help!
Moreno
Azufre451 is offline   Reply With Quote
Reply

Tags
bowtie2, illumina, mate alignment, paired end

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:27 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO