Read panaroos's output into a pagoo's R6 class object

This function handle conversion of panaroo's output files into a pagoo R6 class object. It takes the "gene_presence_absence.csv" file and (optionally but recommended) gff input file paths, and returns an object of class PgR6MS (or PgR6M if left empty the gffs argument). Panaroo identifies some genes with unusual lengths tagging them with 'stop', 'length', or 'refound' labels. In the current version, this function discards those genes.

Usage

panaroo_2_pagoo(gene_presence_absence_csv, gffs, sep = "__")

Arguments

gene_presence_absence_csv

character, path to the "gene_presence_absence.csv" file. (Do not confuse with the file with the same name but with .Rtab extension).

gffs

A character vector with paths to original gff files used as roary's input. Typically the return value of list.files() function. This parameter is optional but highly recommended if you want to manipulate sequences.

sep

character. Default: "__" (two underscores). See PgR6MS for a more detail argument description.

Value

A pagoo's R6 class object. Either PgR6M, if gffs

argument is left empty, or PgR6MS if path to gff files is provided.

References

Tonkin-Hill, G., MacAlasdair, N., Ruis, C. et al. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol 21, 180 (2020). https://doi.org/10.1186/s13059-020-02090-4

Examples

if (FALSE) {
gffs <- list.files(path = "path/to/gffs/",
                   pattern = "[.]gff$",
                   full.names = TRUE)
gpa_csv <- "path/to/gene_presence_absence.csv"

library(pagoo)
pg <- panaroo_2_pagoo(gene_presence_absence_csv = gpa_csv,
                    gffs = gffs)
}