SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Can small RNA library be sequenced paired end? rnaeye Illumina/Solexa 3 07-06-2017 05:13 AM
Need help in mapping Single end and paired-end reads Lizex RNA Sequencing 0 02-17-2013 10:48 PM
Can paired-end mapping produce more reads than single-end ? warrenemmett Bioinformatics 13 03-20-2012 11:10 PM
Paired End Mapping using Bowtie adarshjose Bioinformatics 8 01-26-2012 07:41 AM
Can we extract f3 reads while f5 reads are being sequenced in paired end Raa Bioinformatics 2 12-25-2011 08:46 PM

Reply
 
Thread Tools
Old 01-17-2015, 04:53 PM   #1
jeni
Member
 
Location: UK

Join Date: Oct 2014
Posts: 17
Default Mapping of Paired end sequenced Amlicon?

Hi guys,

I am puzzled how to do mapping of my amplicons PE seq reads. Library was prepared in such a way that PE seq has been done from both sides of amplicons. Now I am slightly confused, How to map these reads over the reference amplicons using bowtie?
jeni is offline   Reply With Quote
Old 01-17-2015, 06:06 PM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Is there a particular reason you need to use bowtie?
Brian Bushnell is offline   Reply With Quote
Old 01-18-2015, 11:31 AM   #3
jeni
Member
 
Location: UK

Join Date: Oct 2014
Posts: 17
Default

No there is no perticular reason. I could use other tool as well.

My question is related with the mapping. If it is single end, I can map it easily over the reference amplicon but how to map the Paired reads over the amplicons?
jeni is offline   Reply With Quote
Old 01-18-2015, 12:17 PM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,976
Default

NGS aligners are able to account for the paired-end reads. You would have no problems as long as you properly indicate input files with read pairs being mapped. What size amplicons do you have?
GenoMax is offline   Reply With Quote
Old 01-19-2015, 05:19 AM   #5
jeni
Member
 
Location: UK

Join Date: Oct 2014
Posts: 17
Default

I have reads which are around 200 bp long.
jeni is offline   Reply With Quote
Old 01-19-2015, 05:22 AM   #6
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,976
Default

See BBMap and the example for how to use paired end reads for mapping.
GenoMax is offline   Reply With Quote
Old 01-19-2015, 05:29 AM   #7
jeni
Member
 
Location: UK

Join Date: Oct 2014
Posts: 17
Default Merge overlapping PE reads

In continuation of the question, I have another idea of merging the overlapping reads of the PE sequenced fastq data.

I found a very interesting page related to this: http://thegenomefactory.blogspot.com...aired-end.html

They have suggested few tools. eg. fastq-join, FLASH etc. Have anybody tried any of such tools? So that I can have a review.


I would appreciate if any one could suggest me for the PE sequenced data,
1) Should I directly go for Mapping using PE option (eg. bowtie1) OR
2) merge the PE reads R1 and R2 together to make a super-long SE read, with extra confidence of the middle bases from consensus of the overlapping sections of R1 and R2?


Thanks all in advance.
jeni is offline   Reply With Quote
Old 01-19-2015, 05:32 AM   #8
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,976
Default

If you know that your reads are overlapping then you can merge them into a longer representation before mapping. If you find that the reads are not merging well (because insert size was bigger than you expected) then you should proceed with mapping them as paired-end.

BBMap suite has a tool that does the merging called BBMerge. FLASH is another well regarded tool for merging reads.

Last edited by GenoMax; 01-19-2015 at 05:34 AM.
GenoMax is offline   Reply With Quote
Old 01-19-2015, 05:54 AM   #9
jeni
Member
 
Location: UK

Join Date: Oct 2014
Posts: 17
Default

Thanks for your response.

I have ~350bp long DNAs fragments and after PE sequencing the length of reads are 200 (FastQC). I would expect an overlap of ~50 bases. And that is the reason, I am thinking of merging not directly mapping with PE parameter.

Last edited by jeni; 01-19-2015 at 07:03 AM.
jeni is offline   Reply With Quote
Old 01-19-2015, 10:16 AM   #10
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Both approaches are valid, just bear in mind that not all PE reads will merge even if they overlap, particularly if are low-complexity in the overlapping region or if the reads have high error rates in their tails. Merged reads are typically easier to work with, but unmerged reads may have lower bias if a large percentage do not merge successfully.

It's also possible to map paired reads, merge them based on mapping coordinates, then remap the merged reads.

Is your goal variant calling, quantification, or something else?
Brian Bushnell is offline   Reply With Quote
Old 01-20-2015, 02:50 AM   #11
jeni
Member
 
Location: UK

Join Date: Oct 2014
Posts: 17
Default

Basically quantification.
jeni is offline   Reply With Quote
Old 01-20-2015, 10:16 AM   #12
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

For quantification, I recommend mapping the unmerged reads as pairs, to prevent bias against low-complexity sequence.
Brian Bushnell is offline   Reply With Quote
Reply

Tags
amplicons, illumina, sequencing

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 01:00 PM.


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