Hey to all , we have a problem with the align.seqs command , we would to analyze function genes, and we try many different parameters , but mothur is inserting to many unnecessary gaps despite a correct reference alignment without gaps , How we can fix this problem, and wich parameters can we used for it ?


Can you explain how you have a reference alignment that doesn’t have gaps?

The reference file just have gaps , but mothur insert new gaps in false postions instead of the alignment corectly of the reference alignment wich paramters do we have to change to avoid these problems.

Best Regards


If you could send your reference file and an example sequence to mothur.bugs at with the url for this thread we can take a look.