SEQanswers

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 10:59 AM
Unique mapper on genome - multi mapper on transcriptome StephaniePi83 Bioinformatics 10 09-04-2012 01:49 AM
IUPAC ambiguous bases in vcf file? sdvie Bioinformatics 0 07-04-2012 02:47 AM
454 Mapper nitinkumar 454 Pyrosequencing 3 05-31-2011 12:57 PM
samtools rm ambiguous alignment EHC Bioinformatics 0 06-11-2010 09:50 AM

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

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

hello,

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


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?

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

Join Date: May 2011
Posts: 823
Default

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, 12:59 PM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

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
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 04:23 AM.


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