SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
mapping on to custom reference sequence PratikC Bioinformatics 14 02-06-2019 04:41 PM
BSMAP: whole genome Bisulfite Sequence MAPping program wei Epigenetics 4 03-20-2014 01:13 PM
Mapping Bisulphite converted sequence simonandrews Epigenetics 35 04-20-2012 08:27 PM
Samtools (sequence alignment/mapping) wisosonic Bioinformatics 2 03-28-2011 04:31 AM
whole genome Bisulfite Sequence MAPping program wei Bioinformatics 0 08-07-2009 02:46 PM

Reply
 
Thread Tools
Old 08-23-2012, 04:50 PM   #1
Poni
Junior Member
 
Location: Australia

Join Date: May 2012
Posts: 7
Smile Sequence mapping

Hi,
I want to map sequence reads into a reference gene. However, it's possible to have a higher number of SNPs and few insertions and/or deletions for this gene in the genome that I'm analyzing. May I please know which mapping software is suitable for this purpose and whether its useful to allow many mismatches?
Your advise is highly appreciated.

Poni
Poni is offline   Reply With Quote
Old 08-24-2012, 04:03 AM   #2
TiborNagy
Senior Member
 
Location: Budapest

Join Date: Mar 2010
Posts: 329
Default

I think there is no specific software for this job. Try BWA and decrease the mismatch penalty (-M).
TiborNagy is offline   Reply With Quote
Old 08-24-2012, 04:25 AM   #3
qiongyi
Member
 
Location: Brisbane

Join Date: Nov 2010
Posts: 10
Default

Try BWA or bowtie2.
qiongyi is offline   Reply With Quote
Old 08-24-2012, 03:14 PM   #4
Poni
Junior Member
 
Location: Australia

Join Date: May 2012
Posts: 7
Default

Thanks for the suggestions.
Poni
Poni is offline   Reply With Quote
Old 08-25-2012, 10:40 PM   #5
Torst
Senior Member
 
Location: The University of Melbourne, AUSTRALIA

Join Date: Apr 2008
Posts: 275
Default

Quote:
Originally Posted by Poni View Post
Hi,
I want to map sequence reads into a reference gene. However, it's possible to have a higher number of SNPs and few insertions and/or deletions for this gene in the genome that I'm analyzing. May I please know which mapping software is suitable for this purpose and whether its useful to allow many mismatches?
Is your reference a whole genome, or just a single gene?

Nearly all short read aligners trade sensitivity for speed. What you need is increase the sensitivity OR use a true dynamic programming aligner which implements Smith-Waterman or Needleman-Wunsch alignment. Tools which can do this are the FASTA 36 package (http://faculty.virginia.edu/wrpearson/fasta/), Exonerate with the affine model (http://www.ebi.ac.uk/~guy/exonerate/) and Needle/Water (http://emboss.sourceforge.net/).

Personally, I would use Shrimp 2.x (http://compbio.cs.toronto.edu/shrimp/) which does well with indels.
Torst is offline   Reply With Quote
Old 08-27-2012, 07:58 PM   #6
Poni
Junior Member
 
Location: Australia

Join Date: May 2012
Posts: 7
Default

Dear Torst,
Thanks for this discriptive reply. my reference is a single gene.
I will try Shrimp 2.x and see.

Thanks
Poni
Poni is offline   Reply With Quote
Reply

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:24 AM.


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