.uc
) to OTU table-cluster_otu
option flag.
For details about installing and running usearch, please refer to the
.uc
extension
in usearch documentation.import_usearch_uc(ucfile, colRead = 9, colOTU = 10, readDelimiter = "_",
verbose = TRUE)
connection
corresponding to the file that contains the usearch output table.
This is passed directly to read.table
.
Please see its file
argument documentation for further
links and details.9
.10
.7.0.109
.
Hopefully later versions of usearch maintain this function and format,
but the phyloseq team has no way to guarantee this,
and so any feedback about this will help maintain future functionality.
For instance, it is currently
assumed that the 9th and 10th columns of the .uc
table
hold the read-label and OTU ID, respectively;
and it is also assumed that the delimiter between sample-name and read
in the read-name entries is a single "_"
.
If this is not true, you may have to update these parameters,
or even modify the current implementation of this function.Also note that there is now a UPARSE-specific output file format,
import_uparse()
function.
import
usearchfile <- system.file("extdata", "usearch.uc", package="phyloseq")
import_usearch_uc(usearchfile)
Run the code above in your browser using DataLab