SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
pairs out of sync!!! emolinari Bioinformatics 1 06-29-2015 12:40 PM
Lots of broken pairs matth431 Illumina/Solexa 5 08-18-2014 11:26 AM
FPKM and discordant pairs westeros Bioinformatics 0 09-26-2013 05:14 AM
filling in the pairs Sarah_Tulin Bioinformatics 1 08-08-2013 11:40 AM
Orientation of mate-pairs JackieBadger Bioinformatics 6 09-18-2012 05:07 PM

Reply
 
Thread Tools
Old 01-11-2016, 04:13 PM   #1
wisense
Member
 
Location: Shanghai,China

Join Date: Sep 2011
Posts: 30
Default Consider dovetail pairs as concordant pairs?

Hi, all

Recently I have a paired-end data from Illumina/Solexa sequencing machine. My sequencing DNA fragments are about 100-300 bp, and my read length is 125 bp for one end of the pairs. My question is, if the dovetail pairs are general in my data? Do I need to consider these pairs as concordant pairs in the Bowtie2 alignment?

I have tried to specify "--dovetail" parameter in Bowtie2, and it has 10% more concordant pairs in the alignment compared to the results without "--dovetail".

Without "--dovetail":
Code:
30047395 reads; of these:
  30047395 (100.00%) were paired; of these:
    7364630 (24.51%) aligned concordantly 0 times
    18847070 (62.72%) aligned concordantly exactly 1 time
    3835695 (12.77%) aligned concordantly >1 times
75.49% overall alignment rate
With "--dovetail":
Code:
30047395 reads; of these:
  30047395 (100.00%) were paired; of these:
    4190935 (13.95%) aligned concordantly 0 times
    21610905 (71.92%) aligned concordantly exactly 1 time
    4245555 (14.13%) aligned concordantly >1 times
86.05% overall alignment rate
Hope to receive your suggestions, thanks.
wisense
wisense is offline   Reply With Quote
Old 02-08-2016, 08:40 PM   #2
dcameron
Member
 
Location: Australia

Join Date: Mar 2013
Posts: 25
Default

Quote:
Originally Posted by wisense View Post
My question is, if the dovetail pairs are general in my data? Do I need to consider these pairs as concordant pairs in the Bowtie2 alignment?
Have you done any adapter trimming? Alignment of full-length reads is likely to result in soft clipped dovetailing reads. Since a significant portion of your library contains fragments shorter than the read length, dovetailing fragments are not unexpected. What you should do with them depends on your pipeline and analysis intent, but yes it is likely that the dovetailing reads originated from fragments shorter than the read length.

Be aware that the untrimmed reads can overalign due to microhomology between the adapter and the reference at the alignment position thus resulting in reads that appear to have the incorrect orientation (which for structural variant calling, could be a big issue).
dcameron is offline   Reply With Quote
Reply

Tags
bowtie2, dovetail, hts, illumina, 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 10:22 AM.


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