This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MOTHUR:16s analysis: maximum ambiguity

Hi friends!!! I am using mothur for analysing 16s rRNA gene sequence analysis. Can you please tell me how much ambiguity (maxambig=?) should I use in my protocol?

16s mothur ambiguity protocol

maxambig = x

x = maximum number of ambiguous bases in a given read?

Sorry, I can't understand. Can you please elaborate a little?

Bacterial genomes typically have multiple copies of the 16S rRNA gene, which are not necessarily identical. When you sequence DNA, the basecaller algorithm may e.g. see a position where A and G are as common and call the position as R, which is IUPAC code for A or G. I presume in mothur maxambig controls how many such bases you want to allow in a sequence before it's discarded

yes.... actually that is my question. How many ambiguous base should I keep in my sequence at most for doing taxonomic analysis in MOTHUR,

How many sequences are discarded if you allow 0, 1, 2, etc.?

0 answers

No answers yet.

Log in to answer this question.