Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Ambiguous bases vs mismatches in Tophat ameyer Bioinformatics 0 10-01-2012 09:59 AM
Unique mapper on genome - multi mapper on transcriptome StephaniePi83 Bioinformatics 10 09-04-2012 12:49 AM
IUPAC ambiguous bases in vcf file? sdvie Bioinformatics 0 07-04-2012 01:47 AM
454 Mapper nitinkumar 454 Pyrosequencing 3 05-31-2011 11:57 AM
samtools rm ambiguous alignment EHC Bioinformatics 0 06-11-2010 08:50 AM

Thread Tools
Old 10-29-2013, 10:09 PM   #1
Senior Member
Location: Vienna

Join Date: Mar 2010
Posts: 107
Default mapper which work with ambiguous bases


which (RNA-seq) mappers can use ambiguous bases in the reference genome?

ambigous bases are e.g.
W for A,T
S for G,C
Y for C,T
B for non-A

and is there a script which can use a SNP dataset and modify the genome accordingly: replace all SNP positions with the correct ambiguous bases?

dietmar13 is offline   Reply With Quote
Old 10-30-2013, 11:49 AM   #2
David Eccles (gringer)
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838

Tophat/Bowtie works with 'N' (i.e. any base), and I think any non-ACGTU are converted to Ns. To accept all ambiguous bases in a DNA database, you need to double the minimum number of bits used to 4 (i.e. presence/absence of each base). The search algorithms also need to change quite a bit, which will probably at least double the time taken. I'm not yet convinced the added complexity is worth it.
gringer is offline   Reply With Quote
Old 10-30-2013, 11:59 AM   #3
Devon Ryan
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480

Gsnap can map to a genome with SNPs though I don't think it takes the exact input type you have in mind.
dpryan is offline   Reply With Quote

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 07:43 PM.

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