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

Similar Threads
Thread Thread Starter Forum Replies Last Post
How do aligners deal with ambiguous bases? luc Bioinformatics 0 10-05-2016 05:59 PM
Ambiguous bases should not be more than total 10% length or more than 14n's in a row. papori Bioinformatics 8 08-21-2014 12:59 AM
mapper which work with ambiguous bases dietmar13 Bioinformatics 2 10-30-2013 12:59 PM
Ambiguous bases vs mismatches in Tophat ameyer Bioinformatics 0 10-01-2012 10:59 AM
IUPAC ambiguous bases in vcf file? sdvie Bioinformatics 0 07-04-2012 02:47 AM

Thread Tools
Old 01-15-2018, 02:40 PM   #1
Location: Ithaca, NY

Join Date: Jun 2012
Posts: 37
Default Calling Ambiguous Bases with Read Mapping

Hi There,

I've used MeGAMerge to merge two assemblies of the same data (, but merging has produced many ambiguous bases. I have good read coverage, so I'd like to use read mapping to call the ambiguous bases (likely SNPs) based on majority rules. Can someone recommend a strategy for this? I've written a script to identify the location of all ambiguous bases and use 'samtools view' to grab all the reads that map to that location, and I could write a script to parse those reads at that specific position, but I hesitate because I'm sure someone has already solved this problem.

p.s. I've use pilon's fix '--amb' function, which should have done the trick, but the output is identical to my input. (

Thanks in advance,
roliwilhelm 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 08:40 PM.

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