Get.sharedseqs, unique otus

Commands in mothur
1

Hello,

I’m running get.sharedseqs to get a list of otus that are unique to certain group in a shared file. My shared file has only two groups, “a” and “an”. I am getting otus that are definitely present in both groups listed as unique. For instance, otu0001 is present in both groups, but is reported as unique for both groups.

get.sharedseqs(shared=CC.0.binary.merge_o2.shared, uniquegroups=a, output=accnos)

get.sharedseqs(shared=CC.0.binary.merge_o2.shared, uniquegroups=an, output=accnos)

I’ve used this with other data sets with more than two groups and it seemed to work fine, is there something else I’m doing wrong here?

Thanks,
Audrey

2

Could you send your shared and log file to mothur.bugs@gmail.com?

3

This topic was automatically closed 10 days after the last reply. New replies are no longer allowed.