Read the output of merops and extract abundance profile

read_merops calculates the abundance of each protease within the bins based on the merops output from run_merops.

read_merops(merops_path, write=FALSE, profile=TRUE)

Arguments

merops_path

a path where Merops output data are. They should have the extension .txt and all files in the path are the ones that need to be read. Output data should have 6 columns with the bin names followed by the Genes obtained in every algorithm (HMMER,Hotpep,DIAMOND), column 'Signalp' indcating if a Peptide signal is found and a column '#ofTools" indicating the number of algorithms that found this Gene.

write

a logical value indicating to save the data imported as a formatted table with .tsv extension with a time stamp and it will be located in your current workin directory

profile

a logical value indicating if you want to print a profile or not.

Examples

if (FALSE) { # \dontrun{
read_merops("C:/Users/bins/")
} # }